三维可压缩磁流体力学方程的重要估计 The Important Estimates for Compressible MHD Equations

doi:10.12677/PM.2022.128143

Pure Mathematics
Vol. 12 No. 08 ( 2022 ), Article ID: 54822 , 7 pages
10.12677/PM.2022.128143

三维可压缩磁流体力学方程的重要估计

王帅，陈菲^*，王传宝

●How to Cite this Article

青岛大学，数学与统计学院，山东青岛

收稿日期：2022年7月11日；录用日期：2022年8月10日；发布日期：2022年8月17日

摘要

本文主要研究了三维可压缩磁流体力学方程其线性方程解的重要估计。当频率分别满足 $| ξ | \leq k_{0}$ 和 $k_{1} \leq | ξ | \leq k_{2}$ 时，通过傅里叶变换的方法，我们建立了两个重要的估计。从频率分解角度来看， $ξ$ 的两个取值范围对应着低频率和中频率，这对于研究三维可压缩磁流体力学非线性方程解的时间衰减率有重要作用。

关键词

频率分解，傅里叶变换，可压缩磁流体力学方程

The Important Estimates for Compressible MHD Equations

Shuai Wang, Fei Chen^*, Chuanbao Wang

School of Mathematics and Statistics, Qingdao University, Qingdao Shandong

Received: Jul. 11^th, 2022; accepted: Aug. 10^th, 2022; published: Aug. 17^th, 2022

ABSTRACT

In this paper, the important estimates of solution about linear equations to three-dimensional compressible magnetohydrodynamic equations are mainly studied. By Fourier splitting method, we establish two important estimates in case of the frequency satisfying $| ξ | \leq k_{0}$ and $k_{1} \leq | ξ | \leq k_{2}$ respectively. From the perspective of frequency decomposition, the two value ranges of $ξ$ correspond to low frequency and medium frequency, which plays an important role in the study of the time decay rates of solution about nonlinear equations to three-dimensional compressible magnetohydrodynamic equations.

Keywords:Frequency Decomposition, Fourier Transform, Compressible Magnetohydrodynamic Equations

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

1. 引言

对于 $(x, t) \in R^{3} \times R^{+}$ ，可压缩磁流体力学(magnetohydrodynamic，简记为MHD)方程为：

${\begin{cases} \partial_{t} σ + div (σ u) = 0, \\ \partial_{t} (σ u) + div (σ u \otimes u) - μ Δ u - (μ + λ) \nabla div u + \nabla P = curl M \times M, \\ \partial_{t} M - ν Δ M = curl (u \times M), div M = 0, \\ (σ, u, M) (x, 0) = (σ_{0}, u_{0}, M_{0}) (x) \to (1, 0, 0), | x | \to + \infty, \end{cases}$ (1)

其中， $σ \in R^{+}$ 表示流体密度， $u \in R^{3}$ 表示流体速度场， $M \in R^{3}$ 表示磁场和 $P = P (σ) = σ^{γ}$ ( $γ \geq 1$ 为绝热指数)表示压力。粘性系数 $μ$ 和 $λ$ 满足 $μ > 0$ ， $2 μ + 3 λ > 0$ ，以及正常数 $ν$ 表示磁扩散系数。

三维可压缩MHD方程不但具有重要的物理意义和广泛的实际应用，研究MHD方程的数学理论也有重要的意义。研究其数学理论的方法有很多，例如能量方法、调和分析方法(傅里叶变换、频率分解等)和RSAV方法等(参考文献 [1] - [7] 等)。本文主要借助于调和分析方法来建立三维可压缩MHD方程线性方程的中低频部分的估计，此类估计在后续建立三维可压缩MHD非线性方程某种大初值解的最优时间衰减率时将发挥重要作用。

下面，我们来回顾一些关于三维可压缩MHD方程解的时间衰减率的结果。Shi-Zhang [4] 研究了初值低频部分在 ${\dot{B}}_{2, \infty}^{- b}$ 范数意义下小时，解在 $L^{r}$ 范数下的时间衰减率。当初值满足 $L^{l} (1 \leq l < \frac{6}{5}) \cap H^{3}$ 时，Chen-Tan [8] 建立了小初值解的整体存在性和时间衰减率，类似的时间衰减率结论可以参考 [9] [10]，进一步， [11] [12] 得到了解的高阶导数的最优时间衰减率。除此之外，当初值分别在 $H^{S} (S \geq 3) \cap {\dot{B}}_{2, \infty}^{- s} (0 \leq s \leq \frac{5}{2})$ 和 $H^{L} (L \geq 3) \cap {\dot{H}}^{- n} (0 \leq n < \frac{3}{2})$ 时，Huang-Lin-Wang [13] 和Tan-Wang [14] 得到了对应的小初值解的时间衰减率。近来，Chen-Huang-Xu [15] 研究了 $L^{1} \cap H^{2}$ 中的大初值解的时间衰减率： ${‖ (σ - 1, u, M) ‖}_{H^{1}} \leq C {(1 + t)}^{- \frac{3}{4}}$ ，Gao-Wei-Yao [16] 进一步优化了磁场的时间衰减率： ${‖ M_{t} ‖}_{L^{2}} + {‖ \nabla M ‖}_{H^{1}} \leq C {(1 + t)}^{- \frac{5}{4}}$ 。

基于文章 [15] [16] 中解的二阶导数的时间衰减率并不是最优的，因此，我们试图借助Wang-Wen [6] 的方法，通过能量估计、傅里叶变换和频率分解等方法建立大初值解的最优时间衰减率。而在频率分解之后，最基本的一步就需要建立对MHD方程其线性方程解的低频率和中频率的估计，所以，我们得到了本文的主要结论：定理1.1和定理1.2，分别对应着对于方程(1)傅里叶变换的线性方程(11)解的低频率和中频率的估计。基于两个定理，可以得到方程(1)线性方程解的中低频的能量估计，进而建立非线性方程解的中低频的能量估计，最终就可以得到非线性方程解的时间衰减率。

本文结构主要为：首先是主要定理的描述，其次是定理的证明过程，包括改写方程、定理1.1证明和定理1.2证明。

2. 主要结论

定理1.1.

对于方程(1)傅里叶变换的线性方程(11)，对任意的 $| ξ | \leq k_{0}$ ， $k_{0}$ 见，存在常数 $C_{l}$ 使得

$l (ξ, t) \leq e^{- c_{l} {| ξ |}^{2} t} l (ξ, 0),$ (2)

其中， $l (ξ, t) ~ {| \hat{q} |}^{2} + {| \hat{ϒ} |}^{2} + {| \hat{M} |}^{2}$ ， $q = σ - 1$ ， $ϒ = {(- Δ)}^{- \frac{1}{2}} div u$ 。

定理1.2.

对于方程(1)傅里叶变换的线性方程(11)，任意的常数 $0 < k_{1} < k_{2}$ ，存在常数 $ς > 0$ (依赖于 $η, ν, k_{1}, k_{2}$ )，使得对任意的 $k_{1} \leq | ξ | \leq k_{2}$ 和 $t \in R^{+}$ ， $| e^{- t Q} | \leq C e^{- ς t}$ 成立，Q见，进一步可得到

$| (\hat{q}, \hat{ϒ}, \hat{M}) (ξ, t) | \leq C e^{- ς t} | (\hat{q}, \hat{ϒ}, \hat{M}) (ξ, 0) |,$ (3)

其中，常数 $C > 0$ 与时间无关。

备注

定理1.1和定理1.2的具体应用将会在后续对解的时间衰减率的研究中得以体现。

3. 定理证明过程

3.1. 改写方程

定义 $q = σ - 1$ ，改写(1)

${\begin{cases} \partial_{t} q + div u = f_{1}, \\ \partial_{t} u + P^{'} (1) \nabla q - μ Δ u - η \nabla div u = f_{2}, \\ \partial_{t} M - ν Δ M = f_{3}, div M = 0, \\ (q, u, M) (x, 0) = (q_{0}, u_{0}, M_{0}) (x) \to (0, 0, 0), | x | \to + \infty, \end{cases}$ (4)

令 $η = μ + λ$ ，此处的非线性项为

${\begin{cases} f_{1} = - q div u - u \cdot \nabla q, \\ f_{2} = - u \cdot \nabla u - (\frac{P^{'} (1 + q)}{1 + q} - P^{'} (1)) \nabla q - \frac{q}{q + 1} (μ Δ u + η \nabla div u) \\ + \frac{1}{q + 1} (M \cdot \nabla M + M \cdot \nabla^{t} M), \\ f_{3} = M \cdot \nabla u - (div u) M - u \cdot \nabla M . \end{cases}$ (5)

定义 $Λ = {(- Δ)}^{\frac{1}{2}}$ ， $ϒ = Λ^{- 1} div u$ ，进一步， $u = - Λ^{- 1} \nabla ϒ + Λ^{- 1} curl (Λ^{- 1} curl u)$ ，由此根据(4)，可以验证

${\begin{cases} \partial_{t} q + Λ ϒ = f_{1}, \\ \partial_{t} ϒ - P^{'} (1) Λ q - (μ + η) Δ ϒ = Λ^{- 1} div f_{2}, \\ \partial_{t} M - ν Δ M = f_{3}, \end{cases}$ (6)

且 $Γ u = Λ^{- 1} curl u$ 满足

${\begin{cases} \partial_{t} Γ u - μ Δ Γ u = Γ f_{2}, \\ Γ u (x, 0) = Γ u_{0} (x), \end{cases}$ (7)

其中， $div Γ u = 0$ 。

基于上面的分析，对u的估计便转化为对 $ϒ$ 和 $Γ u$ 的估计。

首先，对(7)作傅里叶变换，其线性部分为

$\partial_{t} \hat{Γ u} + μ {| ξ |}^{2} \hat{Γ u} = 0,$ (8)

因而，对任意的 $| ξ | \geq 0$ ，有

${| \hat{Γ u} (ξ, t) |}^{2} \leq C e^{- μ {| ξ |}^{2} t} {| \hat{Γ u} (ξ, 0) |}^{2} .$ (9)

其次，对(6)作傅里叶变换，得到

${\begin{cases} \partial_{t} \hat{q} + | ξ | \hat{ϒ} = {\hat{f}}_{1}, \\ \partial_{t} \hat{ϒ} - P^{'} (1) | ξ | \hat{q} + (μ + η) {| ξ |}^{2} \hat{ϒ} = \hat{Λ^{- 1} div f_{2}}, \\ \partial_{t} \hat{M} + ν {| ξ |}^{2} \hat{M} = {\hat{f}}_{3}, \end{cases}$ (10)

显然，其线性方程为

${\begin{cases} \partial_{t} \hat{q} + | ξ | \hat{ϒ} = 0, \\ \partial_{t} \hat{ϒ} - P^{'} (1) | ξ | \hat{q} + (μ + η) {| ξ |}^{2} \hat{ϒ} = 0, \\ \partial_{t} \hat{M} + ν {| ξ |}^{2} \hat{M} = 0, \end{cases}$ (11)

事实上，(11)可写成

$\frac{d}{d t} (\begin{matrix} \hat{q} \\ \hat{ϒ} \\ \hat{M} \end{matrix}) + Q (\begin{matrix} \hat{q} \\ \hat{ϒ} \\ \hat{M} \end{matrix}) = 0,$ (12)

此处

$Q = (\begin{matrix} 0 & | ξ | & 0 \\ - P^{'} (1) | ξ | & (μ + η) {| ξ |}^{2} & 0 \\ 0 & 0 & ν {| ξ |}^{2} \end{matrix}) .$ (13)

3.2. 定理1.1证明

根据(11)，计算得到

$\frac{d}{d t} (\frac{P^{'} (1) {| \hat{q} |}^{2} + {| \hat{ϒ} |}^{2} + {| \hat{M} |}^{2}}{2}) + (μ + η) {| ξ |}^{2} {| \hat{ϒ} |}^{2} + ν {| ξ |}^{2} {| \hat{M} |}^{2} = 0.$ (14)

对(11)的第一个式子乘 $\bar{\hat{ϒ}}$ ，第二个式子先求共轭再乘 $\hat{q}$ ，然后，把二者所得的结果相加并两边取实部得到

$2 \frac{d}{d t} Re (\hat{q} \bar{\hat{ϒ}}) + | ξ | {| \hat{ϒ} |}^{2} - P^{'} (1) | ξ | {| \hat{q} |}^{2} = - (μ + η) {| ξ |}^{2} Re (\hat{q} \bar{\hat{ϒ}}) .$ (15)

选取一个小的固定常数 $δ_{*} > 0$ ，把 $- δ_{*} | ξ | \times$ (15)和(14)相加，并根据杨不等式得

$\begin{array}{l} \frac{d}{d t} (\frac{P^{'} (1) {| \hat{q} |}^{2} + {| \hat{ϒ} |}^{2} + {| \hat{M} |}^{2} - 4 δ_{*} | ξ | Re (\hat{q} \bar{\hat{ϒ}})}{2}) \\ + (μ + η) {| ξ |}^{2} {| \hat{ϒ} |}^{2} + ν {| ξ |}^{2} {| \hat{M} |}^{2} - δ_{*} {| ξ |}^{2} {| \hat{ϒ} |}^{2} + P^{'} (1) δ_{*} {| ξ |}^{2} {| \hat{q} |}^{2} \\ = δ_{*} (μ + η) {| ξ |}^{3} Re (\hat{q} \bar{\hat{ϒ}}) \\ \leq \frac{δ_{*} {(μ + η)}^{2}}{2 P^{'} (1)} {| ξ |}^{4} {| \hat{ϒ} |}^{2} + \frac{P^{'} (1) δ_{*}}{2} {| ξ |}^{2} {| \hat{q} |}^{2} . \end{array}$ (16)

选取

$0 < δ_{*} \leq \min {\frac{1}{2}, \frac{μ + η}{2}},$ (17)

可以验证

$\begin{array}{l} \frac{d}{d t} (\frac{P^{'} (1) {| \hat{q} |}^{2} + {| \hat{ϒ} |}^{2} + {| \hat{M} |}^{2} - 4 δ_{*} | ξ | Re (\hat{q} \bar{\hat{ϒ}})}{2}) \\ + \frac{μ + η}{2} {| ξ |}^{2} {| \hat{ϒ} |}^{2} + ν {| ξ |}^{2} {| \hat{M} |}^{2} + \frac{P^{'} (1) δ_{*}}{2} {| ξ |}^{2} {| \hat{q} |}^{2} \leq \frac{{(μ + η)}^{2}}{4 P^{'} (1)} {| ξ |}^{4} {| \hat{ϒ} |}^{2} . \end{array}$ (18)

选取常数 $k_{0}$ 满足

$| ξ | \leq k_{0} \leq \min {\sqrt{\frac{P^{'} (1)}{μ + η}}, \frac{1}{2}, \frac{P^{'} (1)}{2}},$ (19)

使得

$\frac{d}{d t} l (ξ, t) + \frac{μ + η}{4} {| ξ |}^{2} {| \hat{ϒ} |}^{2} + ν {| ξ |}^{2} {| \hat{M} |}^{2} + \frac{P^{'} (1) δ_{*}}{2} {| ξ |}^{2} {| \hat{q} |}^{2} \leq 0,$ (20)

此处

$l (ξ, t) = \frac{P^{'} (1) {| \hat{q} |}^{2} + {| \hat{ϒ} |}^{2} + {| \hat{M} |}^{2} - 4 δ_{*} | ξ | Re (\hat{q} \bar{\hat{ϒ}})}{2},$ (21)

由于 $δ_{*} k_{0} \leq \min {\frac{1}{4}, \frac{P^{'} (1)}{4}}$ ，根据杨不等式可以得到

$l (ξ, t) ~ {| \hat{q} |}^{2} + {| \hat{ϒ} |}^{2} + {| \hat{M} |}^{2} .$ (22)

因此，对任意 $| ξ | \leq k_{0}$ ，存在常数 $C_{l}$ 使得

$C_{l} {| ξ |}^{2} l (ξ, t) \leq \frac{μ + η}{4} {| ξ |}^{2} {| \hat{ϒ} |}^{2} + ν {| ξ |}^{2} {| \hat{M} |}^{2} + \frac{P^{'} (1) δ_{*}}{2} {| ξ |}^{2} {| \hat{q} |}^{2} .$ (23)

结合(20)和(23)，得到

$l (ξ, t) \leq e^{- c_{l} {| ξ |}^{2} t} l (ξ, 0) .$

定理1.1的证明完毕。

3.3. 定理1.2证明

根据矩阵Q，计算其特征多项式

$\begin{matrix} Q_{λ_{0}} = (λ_{0}^{2} - (μ + η) {| ξ |}^{2} λ_{0} + P^{'} (1) {| ξ |}^{2}) (λ_{0} - ν {| ξ |}^{2}) \\ = a_{0} λ_{0}^{3} - a_{1} λ_{0}^{2} + a_{2} λ_{0} - a_{3}, \end{matrix}$ (24)

此处， $a_{0} = 1$ ， $a_{1} = (μ + η) {| ξ |}^{2} + ν {| ξ |}^{2}$ ， $a_{2} = ν (μ + η) {| ξ |}^{4} + P^{'} (1) {| ξ |}^{2}$ ， $a_{3} = ν P^{'} (1) {| ξ |}^{4}$ 。

由于 $a_{0}$ 至 $a_{3}$ 均大于0，根据Routh-Hurwitz定理，所有的根有正实部当且仅当

$Q_{1} = a_{1} = (μ + η) {| ξ |}^{2} + ν {| ξ |}^{2} > 0$ ，

$Q_{2} = | \begin{matrix} a_{1} & a_{0} \\ a_{3} & a_{2} \end{matrix} | = ν {(μ + η)}^{2} {| ξ |}^{6} + P^{'} (1) (μ + η) {| ξ |}^{4} + ν^{2} (μ + η) {| ξ |}^{6} > 0 ，$

在文章 [1] 中，通过根据Routh-Hurwitz定理，分析得到了相应特征多项式的所有根有正实部，进而建立了关于中频率的估计( [1] 中(55)式)，据此便有，对于任意的常数 $0 < k_{1} < k_{2}$ ，存在常数 $ς > 0$ ， $ς$ 仅依赖于 $η, ν, k_{1}, k_{2}$ ，使得对任意的 $k_{1} \leq | ξ | \leq k_{2}$ 和 $t \in R^{+}$

$| e^{- t Q} | \leq C e^{- ς t}$ (25)

成立，进一步，结合和，即可得到

$| (\hat{q}, \hat{ϒ}, \hat{M}) (ξ, t) | \leq | e^{- t Q} (\hat{q}, \hat{ϒ}, \hat{M}) (ξ, 0) | \leq C e^{- ς t} | (\hat{q}, \hat{ϒ}, \hat{M}) (ξ, 0) | .$

定理1.2的证明完毕。

基金项目

国家自然科学基金资助项目(编号：12101345)，山东省自然科学基金资助项目(编号：ZR2021QA017)。

文章引用

王帅,陈菲,王传宝. 三维可压缩磁流体力学方程的重要估计
The Important Estimates for Compressible MHD Equations[J]. 理论数学, 2022, 12(08): 1305-1311. https://doi.org/10.12677/PM.2022.128143

参考文献

1. Danchin, R. and Ducomet, B. (2014) On a Simplified Model for Radiating Flows. Journal of Evolution Equations, 14, 155-195. https://doi.org/10.1007/s00028-013-0211-5

2. Hu, X.P. and Wang, D.H. (2010) Global Existence and Large-Time Behavior of Solutions to the Three-Dimensional Equations of Compressible Magnetohydrodynamic Flows. Archive for Rational Mechanics and Analysis, 197, 203-238. https://doi.org/10.1007/s00205-010-0295-9

3. Jiang, M.S., Zhang, Z.Y. and Zhao, J. (2022) Improving the Ac-curacy and Consistency of the Scalar Auxiliary Variable (SAV) Method with Relaxation. Journal of Computational Physics, 456, Article ID: 110954. https://doi.org/10.1016/j.jcp.2022.110954

4. Shi, W.X. and Zhang, J.Z. (2021) A Remark on the Time-Decay Estimates for the Compressible Magnetohydrodynamic System. Applicable Analysis, 100, 2478-2498. https://doi.org/10.1080/00036811.2020.1745779

5. Suen, A. and Hoff, D. (2012) Global Low-Energy Weak Solutions of the Equations of Three-Dimensional Compressible Magnetohydrodynamics. Archive for Rational Mechanics and Analysis, 205, 27-58. https://doi.org/10.1007/s00205-012-0498-3

6. Wang, W.J. and Wen, H.Y. (2022) Global Well-Posedness and Time-Decay Estimates for Compressible Navier-Stokes Equations with Reaction Diffusion. Science China Mathematics, 65, 1199-1228. https://doi.org/10.1007/s11425-020-1779-7

7. Zhu, L.M. and Zi, R.Z. (2021) Decay Estimates of the Smooth Solution to the Compressible Magnetohydrodynamic Equations on T3. Journal of Differential Equations, 288, 1-39. https://doi.org/10.1016/j.jde.2021.04.010

8. Chen, Q. and Tan, Z. (2010) Global Existence and Convergence Rates of Smooth Solutions for the Compressible Magnetohydrodynamic Equations. Nonlinear Analysis: Theory Methods and Applications, 72, 4438-4451. https://doi.org/10.1016/j.na.2010.02.019

9. Li, F.C. and Yu, H.J. (2011) Optimal Decay Rate of Classical Solu-tions to the Compressible Magnetohydrodynamic Equations. Proceedings of the Royal Society of Edinburgh. Section A: Mathematics, 141, 109-126. https://doi.org/10.1017/S0308210509001632

10. Pu, X.K. and Guo, B.L. (2013) Global Existence and Conver-gence Rates of Smooth Solutions for the Full Compressible MHD Equations. Zeitschrift für angewandteMathematik und Physik, 64, 519-538. https://doi.org/10.1007/s00033-012-0245-5

11. Gao, J.C., Chen, Y.H. and Yao, Z.A. (2015) Long-Time Behavior of Solution to the Compressible Magnetohydrodynamic Equations. Nonlinear Analysis, 128, 122-135. https://doi.org/10.1016/j.na.2015.07.028

12. Gao, J.C., Tao, Q. and Yao, Z.A. (2016) Optimal Decay Rates of Classical Solutions for the Full Compressible MHD Equations. Zeitschrift für Angewandte Mathematik und Physik, 67, Article No. 23. https://doi.org/10.1007/s00033-016-0616-4

13. Huang, W.T., Lin, X.Y. and Wang, W.W. (2021) Decay-in-Time of the Highest-Order Derivatives of Solutions for the Compressible Isentropic MHD Equations. Journal of Mathematical Analysis and Applications, 502, Article ID: 125273. https://doi.org/10.1016/j.jmaa.2021.125273

14. Tan, Z. and Wang, H.Q. (2013) Optimal Decay Rates of the Compressible Magnetohydrodynamic Equations. Nonlinear Analysis: Real World Applications, 14, 188-201. https://doi.org/10.1016/j.nonrwa.2012.05.012

15. Chen, Y.H., Huang, J.C. and Xu, H.Y. (2019) Global Stability of Large Solutions of the 3-D Compressible Magnetohydrodynamic Equations. Nonlinear Analysis: Real World Appli-cations, 47, 272-290. https://doi.org/10.1016/j.nonrwa.2018.11.001

16. Gao, J.C., Wei, Z.Z. and Yao, Z.A. (2020) Decay Rate of Strong Solution for the Compressible Magnetohydrodynamic Equations with Large Initial Data. Applied Mathematics Letters, 102, Article ID: 106100. https://doi.org/10.1016/j.aml.2019.106100

期刊菜单