分数阶 Burgers 方程时间周期弱解的唯一性与渐近稳定性
Uniqueness and Asymptotic Stability of Time-Periodic Solutions for the Fractional Burgers Equation
Received: 2022-10-8 Revised: 2023-04-10
Fund supported: |
|
作者简介 About authors
徐斐,E-mail:
该文主要讨论实轴上分数阶 Burgers 方程时间周期弱解的性质. 首先利用伽辽金逼近和傅里叶展开的方法, 得到了分数阶 Burgers 方程对应的线性化问题时间周期弱解的存在性. 然后通过构造合适的压缩映射, 得到了非线性方程时间周期弱解的存在唯一性, 并进一步证明了该时间周期弱解是渐近稳定的. 此外, 该文中的方法也可以用来处理有界区间上的问题, 并进一步提升了文献 [
关键词:
The paper is concerned with the time-periodic (T-periodic) problem of fractional Burgers equation on the real line. Based on the Galerkin approximates and Fourier expansion, we first prove the existence of T-periodic solution to a linearized version. Then, the existence and uniqueness of T-periodic solution for the nonlinear equation are established by constructing a suitable contraction mapping. Furthermore, we show that the unique T-periodic solution is asymptotically stable. In addition, our method can be extended to the classical forced Burgers equation in a bounded region, which improves the known result.
Keywords:
本文引用格式
徐斐, 张勇.
Xu Fei, Zhang Yong.
1 引言和主要结果
其中 κ>0 表示黏性系数, 0<α<2 表示耗散指数. 分数阶耗散项 Λαu=(−∂xx)α2u 由下面傅里叶变换所定义
与经典的整数阶模型相比, 分数阶扩散具有在动态发展过程中保持记忆特性和遗传特性的能力, 这使得在一些物理现象中用分数阶模型更加准确. 此外, 分数阶拉普拉斯还可以理解为稳态 Lévy 扩散过程的无穷小生成子, 它们也会出现在等离子体物理[17]、拟地转流[4]和湍流的异常扩散中[3]. 从数学层面来考虑, 上述分数阶 Burgers 方程的适定性严格依赖于耗散指数 α 的取值范围 (参见文献 [8,15,18]). 在超临界耗散情形 (0<α<1), 该分数阶 Burgers 方程是局部适定的, 且它的解会在有限时刻发展为梯度爆破. 在临界耗散情形 (α=1) 和次界耗散情形 (1<α<2), 方程的解整体适定, 有限时间内不会出现奇性. 此外, 在文献 [9] 中考虑了次临界情形下 Burgers 方程弱解的整体正则化效应, 在文献 [1] 中考虑了超临界情形下 Burgers 方程弱解的非唯一性. 在文献 [13] 中还得到了临界情形下方程解的解析性和大时间行为.
其中 f(t,x) 时间周期外力, 即 f(t+T,x)=f(t,x). 在文献 [5] 中, 作者在假设条件 f∈C(0,T;L2(S)∩C1(0,T;L2(S)) 成立的前提下, 可以在空间 C(0,T;C1(S)∩C1(0,T;L2(S)) 中找到唯一的时间周期弱解. 这里我们的结果不仅可以提升该结论, 而且还可以处理如下全空间分数阶外力驱动下的 Burgers 方程
其中 1<α≤2, f(t,x) 表示时间周期外力. 在物理上, 该方程可用于描述一维空间中时间周期外力作用下粘性非均匀扩散液体的流动. 当对流体施加一个时间周期的外力, 可以预期之后的流动会经历一个相同时间周期的运动. 这是研究时间周期分数阶 Burgers 方程的一个主要动机. 因此, 我们很自然地会提出疑问, 方程 (1.1) 在数学上是否存在时间周期解. 事实上, 时间周期解的存在性受到了许多数学家的关注, 例如文献 [10,12,16,20]. 但是据我们所知, 全空间中分数阶 Burgers 方程 (1.1) 的时间周期解的性质很少有人研究. 在这篇文章中, 我们将给出该问题的部分研究结果. 当 f≤cg 其中 c 为常数, 简记 f≲g, 令 Hs(R) 和 ˙Hs(R) 分别表示非齐次和齐次索伯列夫空间. 这篇文章的主要结果是有关于方程 (1.1) 时间周期弱解的存在唯一性和渐近稳定性, 定理陈述如下.
定理 1.1 假设 1<α<32, f∈X=H1(0,T;˙H−α2(R))∩L2(0,T;L2(R))∩L2(0,T;˙H−α(R)) 且外力 f(t,x) 满足
其中 C(\kappa,\alpha)=\min\{ \frac{2\kappa^{2}\sqrt{2(3-2\alpha)}}{9(\kappa+1)^{2}(\sqrt{2(3-2\alpha)}+1)}, \frac{\kappa^{2}}{6(\kappa+1)^{2}} \}, 那么问题 (1.1) 存在唯一时间周期弱解
并在 H^{\frac{\alpha}{2}}(\mathbb{R}) 范数意义下满足 u(0)=u(T), 而且不等式
成立. 此外, 上述得到的时间周期弱解 u(t,x) 是渐近稳定的. 也就是说, 如果 U(t,x) 是初值为 U_{0}\in L^{2}(\mathbb{R}) 柯西问题 (1.1) 的弱解, 那么
注 1.1 基于文献 [20] 中的结果, 我们知道当 U_{0}\in L^{2}(\mathbb{R}) 且 f\in X, 那么柯西问题 (1.1) 一定存在一个全局弱解 U\in L^{\infty}([0,T),L^{2}(\mathbb{R}))\cap L^{2}([0,T),\dot{H}^{\frac{\alpha}{2}}(\mathbb{R})), 并且当 T>0 时该弱解满足
另外值得注意的是, 证明上述定理 1.1 的方法也可以用来提升有界区间上经典 Burgers 方程的结果, 具体陈述如下.
推论 1.1(当 \alpha=2) 假设 f\in X=H^{1}(0,T;\dot{H}^{-1}(\mathbb{S}))\cap L^{2}(0,T;L^{2}(\mathbb{S})) 且外力 f(t,x) 满足
那么带 Dirichlet 边界条件的方程 (1.1) 存在唯一的时间周期解
且该解 u 满足 \|u(0)\|_{H^{1}(\mathbb{S})}=\|u(T)\|_{H^{1}(\mathbb{S})} 以及
此外该时间周期解 u(t,x) 也是渐近稳定的. 也就是说,
本文写作思路如下, 在第 2 节中, 我们利用对时间变量的傅里叶展开和对空间变量的伽辽金逼近的方法证明了方程 (1.1) 对应的线性化问题解的存在唯一性. 在第 3 节中, 主要通过构造合适的压缩映射来证明定理 1.1. 最后一节中, 我们主要在有界区间处理 \alpha=2 的情况, 即证明推论 1.1.
2 线性化问题
在这节中, 我们主要考虑方程 (1.1) 对应的线性化问题
定理 2.1 假设 f\in X=H^{1}(0,T;\dot{H}^{-\frac{\alpha}{2}}(\mathbb{R}))\cap L^{2}(0,T;L^{2}(\mathbb{R}))\cap L^{2}(0,T;\dot{H}^{-\alpha}(\mathbb{R})) 且在 \dot{H}^{-\frac{\alpha}{2}}(\mathbb{\mathbb{R}}) 范数意义下 f(0)=f(T), 那么问题 (2.1) 存在唯一的时间周期弱解
且满足在 H^{\frac{\alpha}{2}}(\mathbb{R}) 意义下 u(0)=u(T) 以及
证 为了更清晰的展示证明过程, 我们首先给出定理 2.1 的证明思路. 对于任给 \varphi\in C^{\infty}_{0,T}((0,T)\times \mathbb{R})=:\{\varphi\in C^{\infty}_{0}| \varphi(0)=\varphi(T)\}, 我们首先尝试找到函数 u\in H^{1}(0,T;\dot{H}^{\frac{\alpha}{2}}(\mathbb{R})) 以及 u_{t}\in L^{2}(0,T;L^{2}(\mathbb{R})) 满足
然后证明
最后, 我们用傅里叶分析的方法提升解的正则性, 也就是说建立估计
和
下面我们分三步来完成定理 2.1 的证明
第一步 证明周期逼近序列 \{u_{nm}\} 的存在性和一致有界性.
令 \{\psi_{j}\}\subset C^{\infty}_{0}(\mathbb{R}) 是空间 H^{\frac{\alpha}{2}}(\mathbb{R}) 的一组基且满足
因为想证明时间周期解的存在性, 我们不妨假设方程 (2.1) 的解可以写为
其中
对于每个固定的 n\in N, 我们首先构造伽辽金逼近序列 \{u_{nm}\}_{m\in N}
在分布意义 (2.3) 下满足线性方程 (2.1).
现在我们证明逼近解序列 u_{nm}(t,x) (即系数 a(n,m)_{lk} 和 b(n,m)_{lk}) 是存在且唯一的. 为了实现这一目标, 我们选择测试函数 \varphi\in C^{\infty}_{0,T}((0,T)\times \mathbb{R}) 且有如下形式
把 (2.12) 和 (2.13) 式带入到 (2.3) 式, 我们可得
和
考虑到三角函数如下的正交性质
则 (2.14) 和 (5.15) 式可约化为
和
那么对于固定的 i\in \{1,2,\cdots,n\}, 我们可以把 (2.16) 和 (2.17) 式写成如下的代数方程组
其中 A 是如下 (2m+2)\times(2m+2) 的矩阵
y 是如下 (2m+2)\times 1 的向量
F 也是如下 (2m+2)\times 1 的向量
易证矩阵 A 是可逆的, 那么由 (2.18) 式可得
是可以被唯一确定的.
现在我们建立伽辽金逼近解 u_{nm}(t,x) 的一致估计. 用 a(n,m)_{ij} 作用 (2.14) 式, 用 b(n,m)_{ij} 作用 (2.15) 式, 并把他们相加可得
由于 u_{nm}(t) 是时间周期的, 这表明
类似地, 用 -b(n,m)_{ij}\frac{2i\pi}{T} 作用 (2.14) 式, 用 a(n,m)_{ij}\frac{2i\pi}{T} 作用 (2.15) 式, 并把他们相加可得
由于 \Lambda^{\frac{\alpha}{2}}u_{nm}(t) 是时间周期的, 那么
最后用 -a(n,m)_{ij}(\frac{2i\pi}{T})^{2} 作用 (2.14) 式, 用 -b(n,m)_{ij}(\frac{2i\pi}{T})^{2} 作用 (2.15) 式, 并把他们相加可得
又由 \Lambda^{\frac{\alpha}{2}}u_{nm}(t), \partial_{t}u_{nm}(t) 和 f(t) 是时间周期的, 那么
因此, 当 n 固定, 上述得到的一致界 (2.19)-(2.21) 式表明存在一个函数列 u_{n}\in H^{1}(0,T;\dot{H}^{-\frac{\alpha}{2}}(\mathbb{R})) 且 \partial_{t}u_{n}\in L^{2}(0,T;L^{2}(\mathbb{R})) 使得
由 (2.22) 可得, 对于所有 \varphi\in C^{\infty}_{0,T}((0,T)\times \mathbb{R}),
成立.
第二步 证明序列 \{u_{n}\} 的周期性和一致有界性.
为了证明序列 u_{n} 在空间 H^{\frac{\alpha}{2}}(\mathbb{R}) 是时间周期的, 定义下面子空间
由于 S_{m} 是有限维的, 易得下面嵌入
是紧的. 那么估计 (2.21) 式表明我们可以选择合适的子序列 \{u_{nm}\}_{m\in N} 使得 \Lambda^{\frac{\alpha}{2}}u_{nm}(t)\rightarrow \Lambda^{\frac{\alpha}{2}}u_{n}(t) 在 [T] 是一致的. 又由于 u_{nm}(t) 是时间周期的, 那么在空间 \dot{H}^{\frac{\alpha}{2}}(\mathbb{R})中 u_{n}(0)=u_{n}(T). 类似的, 由 (2.20) 式可得在空间 L^{2}(\mathbb{R}) 中 u_{n}(0)=u_{n}(T). 因此 u_{n}(t,x) 在空间 H^{\frac{\alpha}{2}}(\mathbb{R}) 中是时间周期的, 那么可写为 (2.21) 式这样的形式.
对序列 \{u_{n}\} 执行类似第一步中的过程, 可得下面一致的界
那么存在一个函数 u\in H^{1}(0,T;\dot{H}^{\frac{\alpha}{2}}(\mathbb{R})) 且 u_{t}\in L^{2}(0,T;L^{2}(\mathbb{R})) 使得
因此, 由 (2.23) 和 (2.25) 式可推得 (2.3) 式, 由 (2.11) 和 (2.24) 式可推得 (2.5)-(2.7) 式.
第三步 证明 u(t,x) 的时间周期性和缺失的正则性
值得注意的是, 这里证明 (2.4) 式中 u 的时间周期性和建立 (2.8) 和 (2.9) 式中 u 的正则性不仅可以完成定理 2.1 的证明, 而且正则性 (2.8) 和 (2.9) 式在下一节处理非线性问题时也发挥着重要的作用. 现在我们首先证明 u(t,x) 是时间周期的, 对于 \varphi\in C^{\infty}_{0,T}((0,T)\times \mathbb{R}), 由于 u(t,x) 是方程 (2.3) 的弱解, 可得
此外, 又由 (2.23) 式可得
考虑到 u_{n}(t) 是时间周期的, 则 (2.25) 和上式表明
因此由 \varphi 的任意性, (2.26) 和 (2.27) 式可得 u(t,x) 是时间周期的.
最后我们用傅里叶分析的技巧来证明 (2.8) 和 (2.9) 式. 对方程 (2.1) 两边取傅里叶变换, 可得
由 (2.5)-(2.7) 式可知 \||\xi|^{\frac{\alpha}{2}}\hat{u}\|_{L^{2}(0,T;L^{2}(\mathbb{R}))} 是有定义的. 对任意的 r>0, 我们有
和
那么 \|\hat{u}\|_{L^{2}(0,T;L^{2}([-r,r]^{c}))} 和 \||\xi|^{\alpha}\hat{u}\|_{L^{2}(0,T;L^{2}([-r,r]))} 也是有定义的. 对 (2.28) 式两边同时平方, 我们可得
又由 (2.31) 式可得
考虑到 (2.29) 式和 u(t) 的时间周期性, 对上式两端在 [T]\times [-r,r]^{c} 上积分可得
那么在 (2.32) 式中令 r\rightarrow 0 可得 (2.8) 式. 类似地, 我们可以从 (2.31) 式推出
考虑到 (2.30) 式和 u(t) 的时间周期性, 对上式两端在 [T]\times [-r,r] 上积分可得
那么在 (2.33) 式中令 r\rightarrow \infty 可得 (2.9) 式. 此外, 容易验证线性问题 (2.1) 的解是唯一的, 证毕.
3 定理 1.1 的证明
在这一节中, 我们主要通过构造压缩映射的方法来证明非线性方程 (1.1) 时间周期解的存在唯一性.
证 为了方便起见, 我们首先引入两个空间
和
根据 Gagliardo-Nirenberg 不等式和文献 [14] 中的 Kato-Ponce 型交换子估计, 接下来我们将证明对于任意 v\in Y_{T} 都有 vv_{x}\in X_{T}. 容易证明当 v\in Y_{T}, 那么(vv_{x})(0)=(vv_{x})(T). 此外, 由插值定理可得
和
这里我们用到 H^{\frac{\alpha}{2}}(\mathbb{R})\hookrightarrow L^{\infty}(\mathbb{R}) 和 Hölder 不等式.
最后, 我们只需估计 \|vv_{x}\|_{L^{2}(0,T;\dot{H}^{-\alpha}(\mathbb{R}))} 的界, 这时限制 1<\alpha<\frac{3}{2} 是本质的. 事实上, 我们希望做到 1<\alpha<2, 但是目前的技术只能限制在 1<\alpha<\frac{3}{2}. 由于 vv_{x}=\frac{1}{2}\partial_{x}(v^{2}), 取傅里叶变换可得
且
那么
因此, 我们有
又由 (3.1)-(3.4) 式, 可得
基于第二节的结果, 我们发现方程 (1.1) 的解可以看作映射 \Gamma: Y_{T}\rightarrow Y_{T} 的不动点, 也就是说对于任给的 v\in Y_{T}, 则 u=\Gamma(v) 是下面问题
的时间周期解. 但是, 通过计算发现很难直接证明映射 \Gamma 是从空间 Y_{T} 到 Y_{T} 的压缩映射. 为了解决这一问题, 我们在空间 Y_{T} 中取一些限制, 那么定义Y_{T}的子空间
对于给定的 v\in Y_{T}^{*}, 由定理 2.1 和 (3.5) 式可得, 方程 (3.6) 存在唯一解 u\in Y_{T}, 且满足
下面我们将证明映射 \Gamma 是从空间 Y^{*}_{T} 到 Y^{*}_{T} 的压缩映射. 对于任意给定的 v\in Y^{*}_{T}, 那么 (3.7) 式和 \|f\|_{X} 的小性表明
也就是说
此外, 如果假设 v_{1}, v_{2}\in Y^{*}_{T} 且 u_{1}=\Gamma(v_{1}), u_{2}=\Gamma(v_{2}), 那么定理 2.1和 \|f\|_{X} 的小性表明
那么 \Gamma 是一个压缩映射, 则时间周期解 u(t,x) 的存在唯一性证得.
最后, 我们来证明该时间周期弱解 u(t,x) 是渐近稳定的. 令 w(t,x)=U(t,x)-u(t,x), 其中 U(t,x) 是注 1.1 中得到的方程 (1.1) 更弱的解, 那么 w 满足
用 w 作用在方程 (3.9) 两端并在 \mathbb{R} 上积分可得
此外, 易得
那么由 (3.10) 和 (3.11) 式可得
对不等式 (3.12) 两端关于时间变量在 (0,T) 上积分, 可得
那么又由 (3.13), (1.3) 式和 \|f\|_{X}\leq\frac{\kappa^{2}}{6(\kappa+1)^{2}}< \frac{\kappa^{2}}{6(\kappa+1)}, 则可推得 (1.4) 式, 证毕.
注 3.1 根据定理 1.1 的证明过程, 估计 (1.3) 可以做如下改进
其中 \beta>1. 此外, 我们期望该结果可以推广到 1< \alpha\leq2, 但是用来构造压缩映射的估计 (3.4) 式限制了该指标, 即目前只能做到 1< \alpha< \frac{3}{2}.
4 推论 1.1 的证明
在本小节中, 我们来证明推论 1.1. 基于上面的观察和分析, 我们发现在有界区间 \mathbb{S} 上且 \alpha=2 时, 当外力 f\in H^{1}(0,T;\dot{H}^{-1}(\mathbb{S}))\cap L^{2}(0,T;L^{2}(\mathbb{S})), 也可以在空间 H^{1}(0,T;H^{1}(\mathbb{S})) 找到时间周期弱解. 此时, 在构造压缩映射时无需建立类似 (3.4) 型的估计. 值得注意的是, 这节中我们考虑的是带零迪利克雷边界条件的方程, 我们在使用分部积分时可以忽略边界项, 但傅里叶变换的技巧就不可用了.
证 根据定理 1.1 的证明, 当 \alpha=2 时, 我们容易建立类似 (2.5)-(2.8) 式的如下估计
因此可以得到类似的线性化问题的时间周期弱解. 此外, 只需证明对于任意 v\in H^{1}(0,T;H^{1}(\mathbb{S})), 有 vv_{x}\in H^{1}(0,T;\dot{H}^{-1}(\mathbb{S}))\cap L^{2}(0,T;L^{2}(\mathbb{S})), 便可构造从 Z_{T} 到 Z_{T} 的压缩映射, 其中 Z_{T}:=\{\phi\in H^{1}(0,T;H^{1}(\mathbb{S}))| \phi(0)=\phi(T) in H^{1}(\mathbb{S})\}, 也就是说, 建立如下类似 (3.1)-(3.3) 式的估计
事实上, 估计 (4.5) 成立是由于
估计 (4.6) 成立是由于
估计 (4.7) 成立是由于
参考文献
Non-uniqueness of weak solutions for the fractal Burgers equation
DOI:10.4171/aihpc URL [本文引用: 1]
Anomalous diffusion in disordered media, Statistical mechanics, models and physical applications
DOI:10.1016/0370-1573(90)90099-N URL [本文引用: 1]
Correlation problem in a one-dimensional model of turbulence
Drift diffusion equations with fractional diffusion and the quasi-geostrophic equation
DOI:10.4007/annals URL [本文引用: 1]
Asymptotic stability and bifurcation of time-periodic solutions for the viscous Burger's equation
DOI:10.1016/j.jmaa.2016.08.018 URL [本文引用: 6]
On a quasi-linear parabolic equation occurring in aerodynamics
DOI:10.1090/qam/1951-09-03 URL [本文引用: 1]
Fractional diffusion on bounded domains
DOI:10.1515/fca-2015-0023 URL [本文引用: 1]
Finite time singularities and global well-posedness for fractal Burgers equations
DOI:10.1512/iumj.2009.58.3505 URL [本文引用: 1]
Global solution and smoothing effect for a non-local regularization of a hyperbolic equation
DOI:10.1007/s00028-003-0503-1 URL [本文引用: 1]
Existence and uniqueness of time-periodic physically reasonable Navier-Stokes flow past a body
DOI:10.1007/s00205-004-0306-9 URL [本文引用: 1]
The partial differential equation u_{t}+uu_{x}=\mu u_{xx}
DOI:10.1002/cpa.v3:3 URL [本文引用: 1]
On the asymptotic stability analysis and the existence of time-periodic solutions of the primitive equations
DOI:10.1512/iumj.2013.62.4902 URL [本文引用: 2]
Analyticity and large time behavior for the Burgers equation and the quasi-geostrophic equation, the both with the critical dissipation
DOI:10.4171/aihpc URL [本文引用: 1]
Commutator estimates and the Euler and Navier-Stokes equations
DOI:10.1002/cpa.v41:7 URL [本文引用: 1]
Blow up and regularity for fractal Burgers euation
DOI:10.4310/DPDE.2008.v5.n3.a2 URL [本文引用: 1]
Time periodic solutions of the Navier-Stokes equations with the time periodic Poiseuille flow under (GOC) for a symmetric perturbed channel in R^{2}
Fractional diffusion limit for collisional kinetic equations
DOI:10.1007/s00205-010-0354-2 URL [本文引用: 1]
Global well-posedness of the critical Burgers equation in critiacal Besov spaces
DOI:10.1016/j.jde.2009.03.028 URL [本文引用: 1]
Forced oscillation of viscous Burgers' equation with a time-periodic force
Uniqueness and stability of steady-state solution with finite energy to the fractal Burgers equation
/
〈 |
|
〉 |
