不可压缩与可压缩Navier-Stokes：密度、压力与正则性

"不可压缩与可压缩"归结为密度问题。密度保持恒定,还是会变化?

试试看。用水填充一支注射器并推动活塞。水会移动,但在日常条件下它不会明显压缩。水对压缩的抵抗如此强烈,以至于将其视为不可压缩是一个极好的近似。这是不可压缩的。现在用空气填充注射器并密封末端。推动活塞,你会感觉到空气让步,相同质量的空气现在被压缩到更小的体积中,在你的拇指下压缩。这就是可压缩流。

在不可压缩流中,密度 $\rho$ 在整个流体中是恒定的,每个微小的流体团在空间中移动时都保持其体积。可压缩流不同。密度成为一个变量,可以在不同位置和不同时刻自由变化。喷气发动机周围的空气、爆炸中的气体、大尺度的大气:都是可压缩的,都由密度变化驱动。

为什么要关心?因为这种区分重塑了纳维-斯托克斯方程的形式、它们的预测内容,以及分析和求解它们的难度。

不可压缩假设将密度场$\rho$设定为整个流动区域内的正常数。从物理上讲,这意味着每个物质体积元在流动映射下保持其体积不变。可压缩情形将$\rho(x,t)$提升为完整的未知量,由其自身的演化方程控制。

这不是同一系统的两种记号。它们在未知量的数目、约束方程的结构以及压力的性质上都有所不同。不可压缩系统有$d+1$个标量未知量($u$和$p$在$\mathbb{R}^d$中);可压缩系统通常有$d+2$个主要场(速度、密度和内能或温度),压力通过状态方程确定。

这一区别也改变了压力的偏微分方程类型:在不可压缩模型中为椭圆型,其中压力由瞬时空间求解确定;而在有粘性的可压缩情形中则是混合双曲-抛物结构。声学模态在无粘或双曲部分以有限声速传播,而粘性和热传导又会在数学上引入抛物扩散。

不可压缩不可压缩Navier-Stokes方程描述密度恒定的流体。它们是出现在Clay千禧年大奖难题中的版本,也是本网站关注的版本。

该系统有两个部分。动量方程:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f$$

以及不可压缩性约束:

$$\nabla \cdot u = 0$$

该约束$\nabla \cdot u = 0$表示速度场是无散的:流体既不会在任何地方堆积,也不会变稀薄。流入微小区域的流体必须以相同速率流出。这一条件取代了整个密度方程。密度不变,因此不需要方程来追踪它。

压力在这里起着特殊作用。它不是由热力学定律(如理想气体定律)确定的。相反,它在各处瞬时调整以保持流动无散。数学上,$p$求解从约束导出的Poisson方程。压力变化以无限快的速度传播。不可压缩流动中没有"声速"。

不可压缩Navier-Stokes系统有两个未知场:速度$u$和压力$p$。这种简单性是有欺骗性的。非线性项$(u \cdot \nabla)u$仍然使该系统在三维中极其困难。

不可压缩Navier-Stokes系统在$\mathbb{R}^3$上,运动粘性系数为$\nu > 0$:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f, \qquad x \in \mathbb{R}^3,\; t > 0,$$

$$\nabla \cdot u = 0, \qquad u(x,0) = u_0(x).$$

无散条件$\nabla \cdot u = 0$是一个逐点约束,而非演化方程。它编码了局部体积守恒:流映射$\Phi_t$满足$\det(D\Phi_t) = 1$对所有$t$。

对动量方程应用散度算子并利用不可压缩性,可得压力Poisson方程:

$$-\Delta p = \partial_i \partial_j (u_i u_j) - \nabla \cdot f.$$

这是一个关于椭圆型方程,在每个固定时刻求解$p$。压力并非独立的热力学变量;在标准的偏微分方程解释中,它充当强制执行无散约束的Lagrange乘子,由速度场全局且瞬时地确定。信息通过压力以无限速度传播,这是与可压缩系统的结构性差异,无法被掩盖。

未知量是$u : \mathbb{R}^3 \times [0,T) \to \mathbb{R}^3$和$p : \mathbb{R}^3 \times [0,T) \to \mathbb{R}$。对于Clay问题的表述(Fefferman, 2000),$u_0 \in C^\infty(\mathbb{R}^3)$是无散的,问题是$u$是否保持在$C^\infty(\mathbb{R}^3 \times [0,\infty))$中且能量有界。

可压缩Navier-Stokes方程控制密度变化的流动。更大的系统。更多未知量。更多方程。

你仍然有一个动量方程,但现在密度$\rho$显式出现:

$$\partial_t (\rho u) + \nabla \cdot (\rho u \otimes u) = -\nabla p + \nabla \cdot \tau + \rho f$$

约束$\nabla \cdot u = 0$不再存在。取而代之的是一个连续性方程,它跟踪密度如何演化:

$$\partial_t \rho + \nabla \cdot (\rho u) = 0$$

这表示质量守恒:密度改变是因为流动压缩或膨胀了流体微团。

该系统还需要一个能量方程和一个状态方程,即热力学关系式,如$p = \rho R T$(理想气体定律),它将压力与密度和温度联系起来。压力不再是约束条件的被动执行者。它的声学部分以声速传播,而完整的有粘系统还包含扩散效应。

可压缩系统对于高速空气动力学、天体物理气体动力学、燃烧以及任何密度变化重要的流动都是必不可少的。但它是一个与不可压缩方程完全不同的数学对象。更多未知量,更多方程,完全不同的偏微分方程结构。

可压缩Navier-Stokes系统耦合了速度$u(x,t)$、密度$\rho(x,t)$、压力$p(x,t)$和比内能$e(x,t)$(或温度$\theta$)。守恒形式为:

连续性方程: $$\partial_t \rho + \nabla \cdot (\rho u) = 0$$

动量方程: $$\partial_t (\rho u) + \nabla \cdot (\rho u \otimes u) + \nabla p = \nabla \cdot \tau + \rho f$$

能量方程: $$\partial_t (\rho E) + \nabla \cdot ((\rho E + p)u) = \nabla \cdot (\tau \cdot u) + \nabla \cdot (\kappa \nabla \theta) + \rho f \cdot u$$

其中$E = e + \tfrac{1}{2}|u|^2$是总比能,$\tau$是粘性应力张量(对于牛顿流体,$\tau = \mu(\nabla u + \nabla u^T) + \lambda (\nabla \cdot u)I$其中体积粘度$\lambda$),$\kappa$是热传导系数。

封闭方程组需要状态方程,例如$p = (\gamma - 1)\rho e$,其中$\gamma$是理想气体的绝热指数。

在可压缩系统的无粘或双曲部分,声扰动以有限速度传播,特征声速为$c = \sqrt{\partial p / \partial \rho |_s}$。完整的有粘Navier-Stokes-Fourier系统则是混合双曲-抛物型,所以粘性和热扩散在数学上会引入无限传播。可压缩方程支持激波、稀疏波和接触间断,这些在不可压缩流动中没有对应物。

何时可压缩性重要?一个数字决定:Mach数。

$$\text{Ma} = \frac{|u|}{c}$$

$|u|$是流动速度。$c$是声速。它们的比值告诉你流动相对于压力扰动在介质中传播的速度有多快,这种比较决定了你是否可以安全地忽略密度变化,还是密度变化将主导物理现象。

当$\text{Ma} < 0.3$时,密度变化小于约5%。不可压缩方程适用。室内空气、管道中的水、建筑物周围的风:所有这些低Mach数流动中,压力扰动的传播速度远快于流动本身,因此密度几乎不变。

超过$\text{Ma} \approx 0.3$时,可压缩性开始显现,在$\text{Ma} \approx 1$附近会进入跨声速区域,局部超声速区域出现并形成激波。战斗机。火箭喷管。再入航天器。

这不是二元开关。大多数日常流体流动以及Clay千禧年难题,都稳固地处于低Mach数区域,不可压缩方程适用。

马赫数$\text{Ma} = |u|/c$ 参数化了可压缩性的重要程度,其中$c = \sqrt{\partial p / \partial \rho |_s}$ 是等熵声速。形式上,不可压缩方程作为可压缩系统的低马赫数极限出现。

关于$\text{Ma}^2$ 幂次的渐近展开(见 Klainerman & Majda, 1981, 1982; Schochet, 1986)表明,当$\text{Ma} \to 0$ 且初始数据适当时,可压缩解收敛到不可压缩解。在具有良好准备数据的标准低马赫数无量纲化中,通常将压力写成几乎空间均匀的热力学背景加上一个较小的动力学修正,该修正在极限情况下强制执行不可压缩性约束。

这是一个奇异极限:声速$c \to \infty$ 且可压缩系统的双曲特性退化为不可压缩流动的椭圆型压力方程。声学模态变得无限快并与涡旋动力学解耦。

区域$\text{Ma} < 0.3$ 是一个工程启发式准则。它反映了相对密度变化$\delta\rho / \rho \sim \text{Ma}^2 / 2$ 在此范围内保持低于$\sim$5% 的经验观察。数学上的合理性依据是低马赫数极限的收敛定理,该定理需要良好准备的初始数据和相容的边界条件。

克雷千禧年大奖难题提出了一个精确的问题:给定 $\mathbb{R}^3$ 上的光滑无散度初始速度,不可压缩纳维-斯托克斯系统是否总能产生一个对所有时间都存在的光滑解?

为什么特别是不可压缩?三个原因。

首先,它已经足够困难。自1934年勒雷的基础性工作以来,不可压缩三维方程的全局正则性证明一直未能解决。添加可变密度、热力学和激波会使��题变得更加困难,而不是更易处理。

其次,困难是纯流体力学的。不可压缩系统分离出核心数学挑战,即非线性平流 $(u \cdot \nabla)u$ 与粘性耗散 $\nu \Delta u$ 之间的竞争,没有热力学或声学复杂性。这是提出正则性问题的最干净场所。

第三,物理是干净的。不可压缩方程模拟最常见的日常流动。它们是否能从光滑数据产生奇异性是关于经典流体力学数学一致性的基本问题。

可压缩系统有其自己深刻的开放问题(大数据的全局解存在性、激波的形成和相互作用),但那些是具有不同结构的不同问题。克雷奖针对不可压缩情形,因为那是Fefferman为三维纳维-斯托克斯方程表述的特定正则性问题。

官方的Clay问题表述(Fefferman, 2000)指定了不可压系统在$\mathbb{R}^3$上:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0, \qquad u|_{t=0} = u_0,$$

其中$u_0 \in C^\infty(\mathbb{R}^3)$无散,问题是是否$u \in C^\infty(\mathbb{R}^3 \times [0,\infty))$且$\int_{\mathbb{R}^3} |u(x,t)|^2\,dx$对所有$t \geq 0$有界。

选择不可压系统有数学上的动机。关键的开放问题,即Leray-Hopf弱解(全局存在但可能不唯一或不光滑)与经典光滑解(局部存在但可能爆破)之间的鸿沟,是三维不可压方程所特有的。在二维情况下,光滑不可压Navier-Stokes解的全局正则性已知;三维情况仍未解决。

可压缩系统引入了性质上不同的困难:激波形成(即使对于Euler方程和光滑数据也会发生)、真空状态($\rho \to 0$),以及涡量与声学模态之间的耦合。这些都是重要的开放问题,但它们在结构上与不可压正则性问题截然不同。

不可压问题隔离了能量超临界非线性与粘性耗散之间的竞争。自然的能量估计给出$u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$,在三维情况下比标度临界空间$L^\infty_t \dot{H}^{1/2}_x$差半个导数。弥合这一鸿沟,或证明它无法弥合,正是千禧年问题的核心。

从这里开始。想要拆解不可压缩系统中的每一项,并从头解释每一部分的物理意义和数学作用?什么是纳维-斯托克斯方程?

这个系统从何而来?纳维-斯托克斯方程的推导。

去掉粘性,你会得到欧拉方程,它早了一个世纪,在纸面上看起来更简单,但在某些方面在数学上甚至更难理解,因为你失去了扩散项的平滑效应。欧拉方程与纳维-斯托克斯方程。

大奖。纳维-斯托克斯存在性与光滑性问题。

从这里建议的路径:

• 什么是纳维-斯托克斯方程?,完整的不可压缩系统及逐项分析
• 纳维-斯托克斯方程的推导,从连续介质力学公理到偏微分方程系统
• 欧拉方程与纳维-斯托克斯方程,无粘极限 $\nu \to 0$ 以及粘性在正则性中的作用
• 纳维-斯托克斯存在性与光滑性问题,精确的克雷表述、勒雷-霍普夫理论,以及弱解与强解之间的差距