We study the global stability of large solutions to the compressible isentropic magnetohydrodynamic equations in a three-dimensional (3D) bounded domain with Navier-slip boundary conditions. It is shown that the solutions converge to an equilibrium state exponentially in the $L^2$-norm provided the density is essentially uniform-in-time bounded from above. Moreover, we also obtain that the density and magnetic field converge to their equilibrium states exponentially in the $L^\infty$-norm if additionally the initial density is bounded away from zero. These greatly improve the previous work in (J. Differential Equations 288 (2021), 1-39), where the authors considered the torus case and required the $L^6$-norm of the magnetic field to be uniformly bounded as well as zero initial total momentum and an additional restriction $2μ>λ$ for the viscous coefficients. This paper provides the first global stability result for large strong solutions of compressible magnetohydrodynamic equations in 3D general bounded domains.
Global stability for the compressible isentropic magnetohydrodynamic equations in 3D bounded domains with Navier-slip boundary conditions 论文ID : 2505.19749标题 : Global stability for the compressible isentropic magnetohydrodynamic equations in 3D bounded domains with Navier-slip boundary conditions作者 : Yang Liu (长春师范大学), Guochun Wu (厦门理工学院), Xin Zhong (西南大学)分类 : math.AP (数学分析-偏微分方程)发表时间 : 2025年1月 (arXiv预印本)论文链接 : https://arxiv.org/abs/2505.19749 本文研究三维有界区域内具有Navier-slip边界条件的可压缩等熵磁流体动力学方程大解的全局稳定性。研究表明，在密度本质上一致有界的条件下，解在L 2 L^2 L 2 L ∞ L^∞ L ∞ L 6 L^6 L 6 2 μ > λ 2\mu>\lambda 2 μ > λ
磁流体动力学(MHD)方程描述了导电流体在磁场作用下的运动，广泛应用于天体物理、等离子体物理和工程领域。可压缩等熵MHD方程组合了可压缩Navier-Stokes方程和Maxwell电磁方程，具有重要的理论价值和实际意义。
本文要解决的核心问题是：在三维有界区域内，具有Navier-slip边界条件的可压缩等熵MHD方程大解的长时间行为和全局稳定性。
区域限制 ：之前的工作主要集中在环面T 3 T^3 T 3 R 3 \mathbb{R}^3 R 3 技术条件过强 ：需要磁场L 6 L^6 L 6 2 μ > λ 2\mu>\lambda 2 μ > λ 边界条件 ：缺乏对Navier-slip边界条件下全局稳定性的系统研究建立更一般区域内MHD方程的全局稳定性理论 去除不必要的技术限制条件 为实际应用提供更符合物理实际的边界条件处理方法 首次建立 了三维一般有界区域内可压缩MHD方程大强解在Navier-slip边界条件下的全局稳定性结果显著改进 了现有理论，去除了多个技术限制：不需要磁场L 6 L^6 L 6 不需要零初始总动量假设 不需要粘性系数限制2 μ > λ 2\mu>\lambda 2 μ > λ 建立了指数收敛率 ：在L 2 L^2 L 2 在额外条件下，密度和磁场在L ∞ L^{\infty} L ∞ 提供了完整的数学框架 ，包括强解的定义、能量估计方法和关键不等式考虑三维有界区域Ω ⊂ R 3 \Omega \subset \mathbb{R}^3 Ω ⊂ R 3
{ ρ t + div ( ρ u ) = 0 ( ρ u ) t + div ( ρ u ⊗ u ) + ∇ P − μ Δ u − ( μ + λ ) ∇ div u = ( ∇ × H ) × H H t − ∇ × ( u × H ) = − ∇ × ( ν ∇ × H ) div H = 0 \begin{cases}
\rho_t + \text{div}(\rho u) = 0 \\
(\rho u)_t + \text{div}(\rho u \otimes u) + \nabla P - \mu\Delta u - (\mu+\lambda)\nabla \text{div} u = (\nabla \times H) \times H \\
H_t - \nabla \times (u \times H) = -\nabla \times (\nu \nabla \times H) \\
\text{div} H = 0
\end{cases} ⎩ ⎨ ⎧ ρ t + div ( ρ u ) = 0 ( ρ u ) t + div ( ρ u ⊗ u ) + ∇ P − μ Δ u − ( μ + λ ) ∇ div u = ( ∇ × H ) × H H t − ∇ × ( u × H ) = − ∇ × ( ν ∇ × H ) div H = 0 配以初始条件和Navier-slip边界条件：
u ⋅ n = 0 , curl u × n = 0 , H ⋅ n = 0 , curl H × n = 0 u \cdot n = 0, \quad \text{curl} u \times n = 0, \quad H \cdot n = 0, \quad \text{curl} H \times n = 0 u ⋅ n = 0 , curl u × n = 0 , H ⋅ n = 0 , curl H × n = 0
定义强解( ρ , u , H ) (ρ, u, H) ( ρ , u , H )
0 ≤ ρ ∈ C ( [ 0 , T ] ; W 1 , q ) 0 \leq \rho \in C([0,T]; W^{1,q}) 0 ≤ ρ ∈ C ([ 0 , T ] ; W 1 , q ) ρ t ∈ C ( [ 0 , T ] ; L q ) \rho_t \in C([0,T]; L^q) ρ t ∈ C ([ 0 , T ] ; L q ) ( u , H ) ∈ C ( [ 0 , T ] ; H 2 ) ∩ L 2 ( 0 , T ; W 2 , q ) (u,H) \in C([0,T]; H^2) \cap L^2(0,T; W^{2,q}) ( u , H ) ∈ C ([ 0 , T ] ; H 2 ) ∩ L 2 ( 0 , T ; W 2 , q ) ( ρ u t , H t ) ∈ L ∞ ( 0 , T ; L 2 ) (\sqrt{\rho}u_t, H_t) \in L^{\infty}(0,T; L^2) ( ρ u t , H t ) ∈ L ∞ ( 0 , T ; L 2 ) 密度上界条件 ：sup t ≥ 0 ∥ ρ ( ⋅ , t ) ∥ L ∞ ≤ ρ ^ \sup_{t \geq 0} \|\rho(\cdot,t)\|_{L^{\infty}} \leq \hat{\rho} sup t ≥ 0 ∥ ρ ( ⋅ , t ) ∥ L ∞ ≤ ρ ^
初始密度下界条件 （用于L ∞ L^{\infty} L ∞ inf x ∈ Ω ρ 0 ( x ) ≥ ρ ∗ > 0 \inf_{x \in \Omega} \rho_0(x) \geq \rho_* > 0 inf x ∈ Ω ρ 0 ( x ) ≥ ρ ∗ > 0
(1) 有效粘性通量方法 定义有效粘性通量：
F = ( 2 μ + λ ) div u − ( P − P ˉ ) − 1 2 ∣ H ∣ 2 F = (2\mu + \lambda)\text{div} u - (P - \bar{P}) - \frac{1}{2}|H|^2 F = ( 2 μ + λ ) div u − ( P − P ˉ ) − 2 1 ∣ H ∣ 2
利用F F F
(2) Bogovskii算子应用 使用Bogovskii算子处理压力项，建立能量-耗散不等式：
d d t E ~ ( t ) + D ~ ( t ) ≤ 0 \frac{d}{dt}\tilde{E}(t) + \tilde{D}(t) \leq 0 d t d E ~ ( t ) + D ~ ( t ) ≤ 0
(3) 分层证明策略 第一步 ：证明∥ ( ρ − ρ ˉ 0 , ρ u , H ) ∥ L 2 \|(ρ-\bar{ρ}_0, \sqrt{ρ}u, H)\|_{L^2} ∥ ( ρ − ρ ˉ 0 , ρ u , H ) ∥ L 2 第二步 ：证明∥ ( ∇ u , ∇ H ) ∥ L 2 \|(\nabla u, \nabla H)\|_{L^2} ∥ ( ∇ u , ∇ H ) ∥ L 2 第三步 ：证明∥ ( ρ u ˙ , H t ) ∥ L 2 \|(\sqrt{ρ}\dot{u}, H_t)\|_{L^2} ∥ ( ρ u ˙ , H t ) ∥ L 2 第四步 ：证明∥ ρ − ρ ˉ 0 ∥ L ∞ \|ρ-\bar{ρ}_0\|_{L^{\infty}} ∥ ρ − ρ ˉ 0 ∥ L ∞ 关键创新在于处理边界项∫ ∂ Ω F t ( u ˙ ⋅ n ) d S \int_{\partial\Omega} F_t(\dot{u} \cdot n)dS ∫ ∂ Ω F t ( u ˙ ⋅ n ) d S ( u ⋅ n ) ∣ ∂ Ω = 0 (u \cdot n)|_{\partial\Omega} = 0 ( u ⋅ n ) ∣ ∂ Ω = 0 u ⋅ ∇ u ⋅ n = − u ⋅ ∇ n ⋅ u on ∂ Ω u \cdot \nabla u \cdot n = -u \cdot \nabla n \cdot u \quad \text{on } \partial\Omega u ⋅ ∇ u ⋅ n = − u ⋅ ∇ n ⋅ u on ∂ Ω
构造时间能量泛函：
M i ( t ) = D i ∫ [ 基本能量项 ] + 高阶项 − 耦合项 M_i(t) = D_i \int[\text{基本能量项}] + \text{高阶项} - \text{耦合项} M i ( t ) = D i ∫ [ 基本能量项 ] + 高阶项 − 耦合项
建立Lyapunov型不等式实现指数衰减。
通过粒子轨迹方法和Grönwall不等式，证明密度的正下界在时间演化中得以保持。
本文为纯理论研究，不涉及数值实验，主要通过严格的数学分析建立结果。
Sobolev空间理论 Gagliardo-Nirenberg不等式 Poincaré不等式 div-curl系统的L p L^p L p 能量方法和Grönwall不等式 引理2.1 ：Gagliardo-Nirenberg不等式的特殊情况引理2.2 ：Bogovskii算子的性质引理2.3 ：广义Poincaré不等式引理2.4 ：有效粘性通量和涡度的先验估计设初始数据满足K = ∥ ρ 0 − ρ ˉ 0 ∥ L 2 + ∥ ρ 0 u 0 ∥ L 2 + ∥ ∇ u 0 ∥ H 1 + ∥ ∇ H 0 ∥ H 1 < ∞ K = \|\rho_0 - \bar{\rho}_0\|_{L^2} + \|\sqrt{\rho_0}u_0\|_{L^2} + \|\nabla u_0\|_{H^1} + \|\nabla H_0\|_{H^1} < \infty K = ∥ ρ 0 − ρ ˉ 0 ∥ L 2 + ∥ ρ 0 u 0 ∥ L 2 + ∥∇ u 0 ∥ H 1 + ∥∇ H 0 ∥ H 1 < ∞ C 1 , η 1 C_1, η_1 C 1 , η 1
∥ ( ρ − ρ ˉ 0 , ρ u , H , ∇ u , ∇ H , ρ u ˙ , H t ) ( ⋅ , t ) ∥ L 2 ≤ C 1 e − η 1 t \|(\rho-\bar{ρ}_0, \sqrt{ρ}u, H, \nabla u, \nabla H, \sqrt{ρ}\dot{u}, H_t)(\cdot,t)\|_{L^2} \leq C_1 e^{-η_1 t} ∥ ( ρ − ρ ˉ 0 , ρ u , H , ∇ u , ∇ H , ρ u ˙ , H t ) ( ⋅ , t ) ∥ L 2 ≤ C 1 e − η 1 t
若额外假设inf x ∈ Ω ρ 0 ( x ) ≥ ρ ∗ > 0 \inf_{x \in Ω} ρ_0(x) \geq ρ_* > 0 inf x ∈ Ω ρ 0 ( x ) ≥ ρ ∗ > 0 C 2 , η 2 C_2, η_2 C 2 , η 2 ∥ ( ρ − ρ ˉ 0 , H ) ∥ L ∞ ≤ C 2 e − η 2 t \|(ρ-\bar{ρ}_0, H)\|_{L^{\infty}} \leq C_2 e^{-η_2 t} ∥ ( ρ − ρ ˉ 0 , H ) ∥ L ∞ ≤ C 2 e − η 2 t
如果初始密度包含真空，即inf x ∈ Ω ρ 0 ( x ) = 0 \inf_{x \in Ω} ρ_0(x) = 0 inf x ∈ Ω ρ 0 ( x ) = 0 t ≥ 0 t \geq 0 t ≥ 0 inf x ∈ Ω ρ ( x , t ) = 0 \inf_{x \in Ω} ρ(x,t) = 0 inf x ∈ Ω ρ ( x , t ) = 0
论文建立了一系列关键的能量估计：
基本能量不等式 （方程3.11）：
d d t M 1 ( t ) + M 1 ( t ) D 1 + ∥ div u ∥ L 2 2 + ∥ curl u ∥ L 2 2 + ∥ ∇ H ∥ L 2 2 D 1 ≤ 0 \frac{d}{dt}M_1(t) + \frac{M_1(t)}{D_1} + \frac{\|\text{div} u\|^2_{L^2} + \|\text{curl} u\|^2_{L^2} + \|\nabla H\|^2_{L^2}}{D_1} \leq 0 d t d M 1 ( t ) + D 1 M 1 ( t ) + D 1 ∥ div u ∥ L 2 2 + ∥ curl u ∥ L 2 2 + ∥∇ H ∥ L 2 2 ≤ 0 高阶能量估计 （方程3.28）：
d d t M 2 ( t ) + M 2 ( t ) D 3 + ∥ ρ u ˙ ∥ L 2 2 + ∥ curl 2 H ∥ L 2 2 D 3 ≤ 0 \frac{d}{dt}M_2(t) + \frac{M_2(t)}{D_3} + \frac{\|\sqrt{ρ}\dot{u}\|^2_{L^2} + \|\text{curl}^2 H\|^2_{L^2}}{D_3} \leq 0 d t d M 2 ( t ) + D 3 M 2 ( t ) + D 3 ∥ ρ u ˙ ∥ L 2 2 + ∥ curl 2 H ∥ L 2 2 ≤ 0 Kawashima 13 ：建立了接近非真空平衡态小扰动的全局光滑解Hu-Wang 10 ：基于Lions-Feireisl框架证明了有限能量弱解的全局存在性Li-Xu-Zhang 15 ：得到了小初始能量但可能有大振荡和真空的经典解Wu-Zhang-Zou 23 ：获得了间断初始数据弱解的最优时间衰减率Zhang-Zhao 27 ：建立了稳定参考态小扰动的时间衰减估计Chen-Huang-Peng-Shi 3 ：在有界区域内得到了指数衰减结果，但需要小初始能量相比于**Zhu-Zi 28 **在环面T 3 T^3 T 3
处理一般有界区域而非环面 去除技术限制2 μ > λ 2\mu > \lambda 2 μ > λ 不需要磁场L 6 L^6 L 6 不要求零初始总动量 首次建立 了三维一般有界区域内可压缩MHD方程大强解的全局稳定性证明了指数收敛 ：解在L 2 L^2 L 2 给出了L ∞ L^{\infty} L ∞ ：在额外条件下密度和磁场也指数收敛显著放松了技术条件 ，使结果更具一般性发展了处理Navier-slip边界条件的新技术 建立了有效的能量估计方法 提供了处理边界积分的创新技巧 密度上界假设 ：需要密度一致有界，这在某些物理情况下可能不现实初始数据限制 ：需要初始数据满足一定的正则性条件边界几何 ：要求区域是单连通光滑有界区域研究更一般的边界条件（如部分滑移条件） 考虑非等熵情况 研究临界指标情况下的长时间行为 发展数值方法验证理论结果 理论突破 ：首次在一般有界区域建立了MHD方程大解的全局稳定性技术创新 ：发展了处理复杂边界条件的新方法结果完整 ：既有L 2 L^2 L 2 L ∞ L^{\infty} L ∞ 条件放松 ：去除了多个不必要的技术限制证明严谨 ：数学推理严密，估计精确物理假设 ：密度上界条件在某些物理情况下可能过强构造性 ：证明是存在性的，没有给出构造方法数值验证 ：缺乏数值实验支持理论结果应用指导 ：对实际应用的指导意义有限学术价值 ：为MHD方程理论研究提供了重要工具方法价值 ：技术方法可推广到其他流体方程理论完善 ：填补了有界区域MHD稳定性理论的空白等离子体约束装置的理论分析 磁流体发电机的稳定性研究 天体物理中磁场-流体相互作用的建模 核聚变装置中等离子体行为的理论预测 论文引用了28篇相关文献，涵盖了MHD方程的全局存在性、长时间行为、数值方法等各个方面，为研究提供了坚实的理论基础。关键参考文献包括Kawashima的开创性工作、Lions-Feireisl的弱解理论框架，以及近年来在稳定性分析方面的重要进展。
