We investigate the global stability of large solutions to the compressible isentropic Navier-Stokes equations in a three-dimensional (3D) bounded domain with Navier-slip boundary conditions. It is shown that the strong 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 obtain that the density converges to its equilibrium state exponentially in the $L^\infty$-norm if additionally the initial density is bounded away from zero. Furthermore, we derive that the vacuum states will not vanish for any time provided vacuum appears (even at a point) initially. This is the first result concerning the global stability for large strong solutions of compressible Navier-Stokes equations with vacuum in 3D general bounded domains.
Global stability for compressible isentropic Navier-Stokes equations in 3D bounded domains with Navier-slip boundary conditions 论文ID : 2504.17136标题 : Global stability for compressible isentropic Navier-Stokes equations in 3D bounded domains with Navier-slip boundary conditions作者 : Yang Liu¹, Guochun Wu², Xin Zhong³
¹长春师范大学数学学院 ²厦门理工学院数学与统计学院 ³西南大学数学与统计学院 分类 : math.AP (数学分析-偏微分方程)提交时间 : 2025年4月论文链接 : https://arxiv.org/abs/2504.17136v4 本文研究三维有界区域上带有Navier滑移边界条件的可压缩等熵Navier-Stokes方程大解的全局稳定性。证明了在密度本质一致有界的条件下，强解在L 2 L^2 L 2 L ∞ L^∞ L ∞
可压缩Navier-Stokes方程描述粘性可压缩流体的运动，是流体力学中的基本方程组。研究其解的长时间行为和稳定性对理解流体运动规律具有重要意义。
边界条件限制 : 以往研究主要集中在全空间R 3 \mathbb{R}^3 R 3 T 3 T^3 T 3 边界条件类型 : 大多数研究考虑无滑移边界条件，而Navier滑移边界条件的研究相对缺乏解的大小限制 : 许多结果要求初始数据较小，对大解的稳定性研究不够充分真空状态 : 带真空的情形增加了数学分析的复杂性Navier滑移边界条件在物理上更符合某些实际情况，如粗糙边界、多孔边界等 三维有界区域的几何复杂性需要新的数学技巧 大解的全局稳定性分析具有重要的理论价值 首次结果 : 这是首个关于三维一般有界区域中带真空的可压缩Navier-Stokes方程大强解全局稳定性的结果指数收敛性 : 证明了在适当条件下解的指数收敛性，给出了具体的衰减率真空持续性 : 证明了真空状态的持续性——一旦出现真空就不会消失技术创新 : 发展了处理Navier滑移边界条件的新技术方法考虑三维有界区域Ω ⊂ R 3 \Omega \subset \mathbb{R}^3 Ω ⊂ R 3
{ ρ t + div ( ρ u ) = 0 ( ρ u ) t + div ( ρ u ⊗ u ) + ∇ P = μ Δ u + ( μ + λ ) ∇ div u \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
\end{cases} { ρ t + div ( ρ u ) = 0 ( ρ u ) t + div ( ρ u ⊗ u ) + ∇ P = μ Δ u + ( μ + λ ) ∇ div u
其中P = ρ γ P = \rho^\gamma P = ρ γ γ > 1 \gamma > 1 γ > 1 u ⋅ n = 0 , curl u × n = 0 在 ∂ Ω 上 u \cdot n = 0, \quad \text{curl} u \times n = 0 \quad \text{在} \partial\Omega \text{上} u ⋅ n = 0 , curl u × n = 0 在 ∂ Ω 上
对T > 0 T > 0 T > 0 q ∈ ( 3 , 6 ] q \in (3,6] q ∈ ( 3 , 6 ] ( ρ , u ) (ρ,u) ( ρ , u )
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 ∈ C ( [ 0 , T ] ; H 2 ) ∩ L 2 ( 0 , T ; W 2 , q ) u \in C([0,T]; H^2) \cap L^2(0,T; W^{2,q}) u ∈ C ([ 0 , T ] ; H 2 ) ∩ L 2 ( 0 , T ; W 2 , q ) ρ u t ∈ L ∞ ( 0 , T ; L 2 ) \sqrt{\rho}u_t \in L^∞(0,T; L^2) ρ u t ∈ L ∞ ( 0 , T ; L 2 ) u t ∈ L 2 ( 0 , T ; H 1 ) u_t \in L^2(0,T; H^1) u t ∈ L 2 ( 0 , T ; H 1 ) 构造Lyapunov泛函：
M 1 ( t ) = D 1 ∫ ( 1 2 ρ ∣ u ∣ 2 + G ( ρ ) ) d x − ∫ ρ u ⋅ B [ ρ − ρ ˉ 0 ] d x M_1(t) = D_1 \int \left(\frac{1}{2}\rho|u|^2 + G(\rho)\right)dx - \int \rho u \cdot B[\rho - \bar{\rho}_0]dx M 1 ( t ) = D 1 ∫ ( 2 1 ρ ∣ u ∣ 2 + G ( ρ ) ) d x − ∫ ρ u ⋅ B [ ρ − ρ ˉ 0 ] d x
其中G ( ρ ) G(\rho) G ( ρ ) B B B
定义有效粘性通量：
F = ( 2 μ + λ ) div u − ( P − P ˉ ) F = (2\mu + \lambda)\text{div} u - (P - \bar{P}) F = ( 2 μ + λ ) div u − ( P − P ˉ )
利用其满足的椭圆方程：
{ Δ F = div ( ρ u ˙ ) 在 Ω 中 ∂ F ∂ n = ρ u ˙ ⋅ n 在 ∂ Ω 上 \begin{cases}
\Delta F = \text{div}(\rho \dot{u}) & \text{在} \Omega \text{中} \\
\frac{\partial F}{\partial n} = \rho \dot{u} \cdot n & \text{在} \partial\Omega \text{上}
\end{cases} { Δ F = div ( ρ u ˙ ) ∂ n ∂ F = ρ u ˙ ⋅ n 在 Ω 中 在 ∂ Ω 上
对涡度ω = ∇ × u \omega = \nabla \times u ω = ∇ × u ( ω × n ) ∣ ∂ Ω = 0 (\omega \times n)|_{\partial\Omega} = 0 ( ω × n ) ∣ ∂ Ω = 0 ∥ ∇ u ∥ L q ≤ C ( ∥ div u ∥ L q + ∥ ω ∥ L q ) \|\nabla u\|_{L^q} \leq C(\|\text{div} u\|_{L^q} + \|\omega\|_{L^q}) ∥∇ u ∥ L q ≤ C ( ∥ div u ∥ L q + ∥ ω ∥ L q )
边界积分处理 : 通过恒等式u ⋅ ∇ u ⋅ n = − u ⋅ ∇ n ⋅ u u \cdot \nabla u \cdot n = -u \cdot \nabla n \cdot u u ⋅ ∇ u ⋅ n = − u ⋅ ∇ n ⋅ u 分层能量估计 : 采用三层能量估计策略：
第一层：∥ ( ρ − ρ ˉ 0 , ρ u ) ∥ L 2 \|(\rho - \bar{\rho}_0, \sqrt{\rho}u)\|_{L^2} ∥ ( ρ − ρ ˉ 0 , ρ u ) ∥ L 2 第二层：∥ ∇ u ∥ L 2 \|\nabla u\|_{L^2} ∥∇ u ∥ L 2 第三层：∥ ρ u ˙ ∥ L 2 \|\sqrt{\rho}\dot{u}\|_{L^2} ∥ ρ u ˙ ∥ L 2 粒子轨道方法 : 利用特征线方法分析密度的下界假设初始数据满足K = ∥ ρ 0 − ρ ˉ 0 ∥ L 2 + ∥ ρ 0 u 0 ∥ L 2 + ∥ ∇ u 0 ∥ H 1 < ∞ K = \|\rho_0 - \bar{\rho}_0\|_{L^2} + \|\sqrt{\rho_0}u_0\|_{L^2} + \|\nabla u_0\|_{H^1} < \infty K = ∥ ρ 0 − ρ ˉ 0 ∥ L 2 + ∥ ρ 0 u 0 ∥ L 2 + ∥∇ u 0 ∥ H 1 < ∞ sup t ≥ 0 ∥ ρ ( ⋅ , t ) ∥ L ∞ ≤ ρ ^ \sup_{t \geq 0}\|\rho(\cdot,t)\|_{L^∞} \leq \hat{\rho} sup t ≥ 0 ∥ ρ ( ⋅ , t ) ∥ L ∞ ≤ ρ ^ C 1 , η 1 C_1, \eta_1 C 1 , η 1
∥ ( ρ − ρ ˉ 0 , ρ u , ∇ u , ρ u ˙ ) ( ⋅ , t ) ∥ L 2 ≤ C 1 e − η 1 t \|(\rho - \bar{\rho}_0, \sqrt{\rho}u, \nabla u, \sqrt{\rho}\dot{u})(\cdot,t)\|_{L^2} \leq C_1 e^{-\eta_1 t} ∥ ( ρ − ρ ˉ 0 , ρ u , ∇ u , ρ u ˙ ) ( ⋅ , t ) ∥ L 2 ≤ C 1 e − η 1 t
若另外inf x ∈ Ω ρ 0 ( x ) ≥ ρ ∗ > 0 \inf_{x \in \Omega}\rho_0(x) \geq \rho_* > 0 inf x ∈ Ω ρ 0 ( x ) ≥ ρ ∗ > 0 C 2 , η 2 C_2, \eta_2 C 2 , η 2 ∥ ρ − ρ ˉ 0 ∥ L ∞ ≤ C 2 e − η 2 t \|\rho - \bar{\rho}_0\|_{L^∞} \leq C_2 e^{-\eta_2 t} ∥ ρ − ρ ˉ 0 ∥ L ∞ ≤ C 2 e − η 2 t
在定理1.1的条件下，若inf x ∈ Ω ρ 0 ( x ) = 0 \inf_{x \in \Omega}\rho_0(x) = 0 inf x ∈ Ω ρ 0 ( x ) = 0 t ≥ 0 t \geq 0 t ≥ 0 inf x ∈ Ω ρ ( x , t ) = 0 \inf_{x \in \Omega}\rho(x,t) = 0 inf x ∈ Ω ρ ( x , t ) = 0
通过乘子法和Bogovskii算子，建立能量-耗散不等式：
d d t M 1 ( t ) + M 1 ( t ) D 1 + ∥ div u ∥ L 2 2 + ∥ curl u ∥ L 2 2 D 1 ≤ 0 \frac{d}{dt}M_1(t) + \frac{M_1(t)}{D_1} + \frac{\|\text{div} u\|_{L^2}^2 + \|\text{curl} u\|_{L^2}^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 ≤ 0
利用动量方程和有效粘性通量的性质，得到：
d d t Ψ ( t ) + 1 2 ∥ ρ u ˙ ∥ L 2 2 ≤ C ( ∥ ∇ u ∥ L 2 2 + ∥ ∇ u ∥ L 2 6 + ∥ ρ − ρ ˉ 0 ∥ L 2 2 ) \frac{d}{dt}\Psi(t) + \frac{1}{2}\|\sqrt{\rho}\dot{u}\|_{L^2}^2 \leq C(\|\nabla u\|_{L^2}^2 + \|\nabla u\|_{L^2}^6 + \|\rho - \bar{\rho}_0\|_{L^2}^2) d t d Ψ ( t ) + 2 1 ∥ ρ u ˙ ∥ L 2 2 ≤ C ( ∥∇ u ∥ L 2 2 + ∥∇ u ∥ L 2 6 + ∥ ρ − ρ ˉ 0 ∥ L 2 2 )
通过改进的能量估计，控制∥ ρ u ˙ ∥ L 2 \|\sqrt{\rho}\dot{u}\|_{L^2} ∥ ρ u ˙ ∥ L 2
利用特征线方法和密度的下界估计，证明L ∞ L^∞ L ∞
边界项处理 : Navier滑移边界条件导致速度在边界上不为零，需要特殊处理边界积分非线性估计 : 处理∫ ∂ Ω F t ( u ⋅ ∇ n ⋅ u ) d S \int_{\partial\Omega} F_t(u \cdot \nabla n \cdot u)dS ∫ ∂ Ω F t ( u ⋅ ∇ n ⋅ u ) d S 真空区域 : 密度可能为零增加了分析的复杂性恒等式应用 : 巧妙利用( u ⋅ n ) ∣ ∂ Ω = 0 (u \cdot n)|_{\partial\Omega} = 0 ( u ⋅ n ) ∣ ∂ Ω = 0 分层Lyapunov方法 : 构造多层能量泛函实现递进式估计粒子轨道分析 : 结合Lagrange观点分析密度演化Matsumura-Nishida (1980s) : 小扰动情形的经典结果Lions (1998) : 大初值弱解的存在性Huang-Li-Xin (2012) : Cauchy问题的大解存在性He-Huang-Wang (2019) : 全空间大解的稳定性Cai-Li (2023) : 有界区域小能量解的存在性相比现有工作，本文首次在三维有界区域的Navier滑移边界条件下建立了大解的全局稳定性理论。
证明了三维有界区域Navier滑移边界条件下可压缩Navier-Stokes方程大强解的指数稳定性 建立了密度在L 2 L^2 L 2 L ∞ L^∞ L ∞ 揭示了真空状态的持续性机制 数学理论 : 丰富了可压缩流体方程的稳定性理论方法论 : 发展了处理滑移边界条件的新技术几何效应 : 揭示了有界区域几何对解行为的影响密度上界假设 : 需要假设密度一致有界，这在实际中可能难以验证强解假设 : 结果基于强解的存在性，而强解存在性本身是困难问题特殊边界条件 : 仅考虑了特定形式的Navier滑移条件弱化假设 : 研究如何减弱密度上界等技术假设一般边界 : 考虑更一般的滑移边界条件数值验证 : 发展相应的数值方法验证理论结果原创性强 : 首次解决该设定下的全局稳定性问题技术深度 : 发展了处理边界条件的新技术方法结果完整 : 不仅证明收敛性，还给出了具体衰减率物理意义 : 真空持续性结果具有重要物理含义方法创新 : Lyapunov泛函构造和边界项处理技巧具有创新性估计精细 : 多层能量估计策略设计合理论证严谨 : 数学论证逻辑清晰，技术细节处理得当假设较强 : 密度一致上界的假设在物理上较为限制性适用范围 : 结果局限于特定的边界条件和区域类型数值验证缺失 : 缺乏数值实验支持理论结果理论贡献 : 为可压缩流体稳定性理论做出重要贡献方法价值 : 所发展的技术方法可能适用于相关问题后续研究 : 为进一步研究奠定了理论基础理论研究 : 偏微分方程理论和流体力学数学理论数值分析 : 为数值方法的长时间行为分析提供理论支撑工程应用 : 为涉及滑移边界的流体问题提供理论指导论文引用了44篇相关文献，涵盖了可压缩Navier-Stokes方程理论的主要发展历程，从经典的Matsumura-Nishida工作到近期的相关研究，文献调研较为全面。
总体评价 : 这是一篇高质量的理论数学论文，在可压缩流体方程稳定性理论方面取得了重要进展。虽然技术假设较强，但所解决的问题具有重要理论意义，所发展的方法具有创新性和推广价值。
