Euler方程与Navier-Stokes方程的区别：黏性、耗散与正则性

欧拉方程描述零内摩擦的流体。完全没有粘性。纳维-斯托克斯方程描述包含粘性的相同流体。

在数学上,全部差异在于一项:$\nu \Delta u$,即粘性扩散项。去掉它,纳维-斯托克斯就变成欧拉。保留它,方程就获得了一种平滑机制,这种机制以你从单个额外项无法预料的方式改变了物理和分析。

正是这一项解释了为什么烟雾会消散,为什么边界层沿表面形成,以及为什么纳维-斯托克斯千禧年问题与相应的欧拉问题具有完全不同的性质。

在$\mathbb{R}^3$上的不可压缩Euler方程为

$$\partial_t u + (u \cdot \nabla)u = -\nabla p, \qquad \nabla \cdot u = 0.$$

不可压缩Navier-Stokes方程为

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0,$$

其中$\nu > 0$为运动粘度。形式上,在Navier-Stokes中令$\nu = 0$可恢复Euler方程。但这种形式代换掩盖了一个事实:$\nu \to 0$极限是奇异的:粘性项$\nu \Delta u$含有系统中最高阶的空间导数,去掉它会改变偏微分方程的类型以及解所在的函数空间。

两个方程的标准不可压缩形式,写法便于对比:

Euler:

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

Navier-Stokes:

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

相同的左边:速度的变化率加上非线性自输运项$(u \cdot \nabla)u$。两者都通过$\nabla \cdot u = 0$强制不可压缩性。唯一的结构差异是Navier-Stokes右边的粘性项$\nu \Delta u$。

参数$\nu$是运动粘度,流体的物理常数。蜂蜜的$\nu$较大。空气的较小。Euler方程对应于$\nu = 0$:一个完全无摩擦的理想化,可以在远离边界处近似某些高Reynolds数流动,但在任何真实流体中都不存在。

两个系统都共享双线性形式$(u \cdot \nabla)u$以及由无散度约束隐式确定的压力。对动量方程取散度并使用$\nabla \cdot u = 0$可得压力Poisson方程

$$-\Delta p = \nabla \cdot ((u \cdot \nabla)u) = \partial_i \partial_j (u_i u_j),$$

对两个系统都相同。在任一情况下,压力都是$u$的非局部泛函。

粘性项$\nu \Delta u$是作用于每个速度分量的二阶线性椭圆算子。它是抛物型正则化:它的存在使Navier-Stokes成为半线性抛物型系统,而不可压缩Euler方程是具有非局部压力耦合的一阶非线性输运方程。偏微分方程类型上的这一差异驱动了正则性理论中几乎所有后续的差异。

粘性是流体相邻层之间的摩擦。快层紧邻慢层?粘性在它们之间传递动量,平滑速度差异。概念很简单。后果是巨大的,它们以三种方式将纳维-斯托克斯与欧拉分开。

• 耗散。动能转化为热量。搅拌咖啡,然后停止。它最终会静止,因为粘性将运动作为热能消耗掉。欧拉完全无法预测这一点,因为方程中没有将动能耗散为热量的机制。
• 边界层。真实流体粘附在表面上(无滑移条件),在壁面附近产生速度快速变化的薄层。这些产生飞机机翼上的阻力、管道中的摩擦损失以及高速下的湍流发生。欧拉流满足滑移条件,因此它们完全忽略了粘性壁面阻力。
• 小尺度平滑。粘性消除最尖锐的速度梯度。没有它?没有什么能阻止流动发展出无限精细的结构,永远越来越尖锐。这种平滑正是使纳维-斯托克斯的正则性问题与欧拉成为不同野兽的原因。

Navier-Stokes在$\mathbb{R}^3$(或周期域)上的能量恒等式为

$$\frac{d}{dt}\frac{1}{2}\|u\|_{L^2}^2 = -\nu \|\nabla u\|_{L^2}^2,$$

因此动能单调耗散。对于Euler($\nu = 0$),$L^2$的$u$范数形式上守恒。

在边界层面,Navier-Stokes使用无滑移条件($u|_{\partial \Omega} = 0$),而Euler只要求不可穿透性($u \cdot n = 0$)。当$\nu \to 0$时,这些条件之间的不匹配产生Prandtl边界层。这是一个奇异摄动现象,人们自Prandtl 1904年的论文以来一直在与之搏斗。

物理上,粘度充当高频滤波器:它以$\nu |k|^2$的速率阻尼Fourier模式,优先抑制小尺度。这种谱阻尼是湍流中Kolmogorov耗散尺度$\eta \sim (\nu^3 / \varepsilon)^{1/4}$背后的机制。完整的标度图景见Reynolds数与湍流。

形式上?是的。设 $\nu = 0$ 就得到欧拉。但停止思考这一点是很糟糕的。

极限 $\nu \to 0$ 是奇异的。粘性携带方程中的最高阶导数,因此去除它并不是做一个小调整。它完全改变了你正在处理的偏微分方程类型。边界层不会优雅地变薄。它们可能爆发成湍流。在纳维-斯托克斯下完全光滑的解可以在欧拉下发展出完全不同的行为。

是的,两个方程共享它们的数学DNA。但零粘度极限是整个流体动力学中最深刻的开放问题之一,不是餐巾纸计算。

无粘极限$\nu \to 0$是奇异摄动:$\nu \Delta u$含有系统中最高阶的空间导数,因此令$\nu = 0$会降低偏微分方程的阶数。在有边界的区域上,该极限与Prandtl边界层展开的有效性相关,后者可能惨遭失败(Grenier 2000, Gérard-Varet & Dormy 2010)。

在$\mathbb{R}^3$或$\mathbb{T}^3$(无边界)上,情况变得更清晰。如果Euler解$u^E$在$[0,T]$上保持光滑,那么Navier-Stokes解$u^\nu$在$u^E$范数下收敛到$L^2$,当$\nu \to 0$时(Kato 1972)。更清晰,但并不干净:Euler解在三维情况下是否甚至能全局保持光滑本身就是未解决的,而这正是整个问题所在。

该极限也与湍流理论正面碰撞。Kolmogorov的图景需要$\nu > 0$来定义耗散尺度。然而反常耗散,即当$\nu \to 0$时能量耗散的持续性,几十年来一直是一个未解决的猜想。Onsager猜想(现在是定理:Isett 2018, 由Buckmaster, De Lellis, Székelyhidi和Vicol在2019年改进)精确刻画了Euler解何时可以在完全没有粘度的情况下耗散能量。

每当粘性与其他作用力相比可以忽略时。这比你想象的更常见:

• 远离表面的高速空气动力学。远离机翼,气流几乎是无粘的。工程师通常为主体流使用欧拉求解器,并在壁面附近打补丁加入边界层修正。
• 天体物理流动。星际气体云、恒星内部、黑洞周围的吸积盘。在这些尺度下,分子粘度完全无关紧要(尽管湍流有效粘度可能不是)。
• 可压缩气体动力学。激波。爆轰。超音速飞行。主导的物理是压力和惯性,而不是摩擦。
• 纯理论。欧拉本身值得研究,不仅仅是通向纳维-斯托克斯的垫脚石。它连接到黎曼几何、涡动力学和关于湍流结构本身的深刻问题。

但对于摩擦、阻力或边界行为重要的任何情况(管道流动、表面附近的车辆空气动力学、血液循环、人类尺度的天气),你需要纳维-斯托克斯。句号。

欧拉方程在高雷诺数 $\mathrm{Re} = UL/\nu \gg 1$ 时支配领先阶行为。在这些雷诺数下,粘性效应被限制在薄边界层和内部剪切层中,而远离壁面的主体流由欧拉很好地近似。

可压缩欧拉方程,一个具有有限传播速度的双曲系统,是气体动力学的标准模型,包括激波形成和黎曼问题。这些不同于上面讨论的不可压缩欧拉方程:可压缩欧拉是真正的双曲型,而不可压缩欧拉具有非局部压力耦合和无限传播速度。

在数学分析中,欧拉既作为消失粘度问题的极限对象,又作为一个丰富的偏微分方程系统,具有自己的正则性理论、守恒量(通过微分同胚群上的欧拉-阿诺德框架的螺旋度、卡西米尔量)以及与几何力学的联系。

这是差距最重要的地方,也是事情真正变得有趣的地方。

纳维-斯托克斯千禧年问题提出了一个听起来几乎太简单的问题:如果你从三维中光滑、表现良好的流动开始,解是否永远保持光滑,还是会爆炸?地球上没有人知道答案。

欧拉的同样问题在3D中也是开放的。但这两个问题感觉完全不同:

• 纳维-斯托��斯站在粘性这一边。总是平滑,总是耗散能量,总是阻尼最尖锐的梯度。真正的问题是这种平滑是否足够强以在非线性项产生奇点之前压制它。
• 欧拉什么都没有。零平滑。零耗散。非线性项可以在绝对没有对抗力的情况下放大速度梯度,从光滑的3D初始数据是否真的产生有限时间奇点是偏微分方程理论中最大的开放问题之一。

在2D中,对于光滑初始数据,两个方程都是全局适定的。解决了。完成了。谜团完全存在于三维中,对于两个方程都是如此,但原因根本不同。

正则性图景:

2D:两个系统的光滑解的全局存在性和唯一性都是已知的。对于具有光滑数据的2D欧拉,Wolibner(1933)在赫尔德空间中证明了全局存在性;Yudovich(1963)为有界涡度数据建立了唯一性。对于2D纳维-斯托克斯,全局正则性源于拉迪任斯卡娅不等式和涡度最大值原理。

3D纳维-斯托克斯:Leray(1934)证明了 $L^2$ 中弱解的全局存在性,但唯一性和正则性仍然开放。Caffarelli、Kohn和Nirenberg定理(1982)表明奇异集具有零的一维抛物豪斯多夫测度,因此任何爆破,如果发生,都是极其稀疏的。粘性项提供了关键的先验估计 $\int_0^T \|\nabla u\|_{L^2}^2 \, dt \leq C(u_0)$,但这种能量层面的控制低于三维的临界标度,不足以闭合自举论证。参见为什么纳维-斯托克斯很难了解超临界间隙。

3D欧拉:对于光滑数据不存在全局理论。在索伯列夫空间 $H^s$,$s > 5/2$ 中的局部适定性是经典的(Kato 1972,Kato和Ponce 1988)。Beale、Kato和Majda准则(1984)将爆破检测归结为 $\int_0^T \|\omega\|_{L^\infty} \, dt$:当且仅当这个积分有限时,解在 $[0,T]$ 上保持光滑。爆破要求涡度增长得足够快,以至于在时间上不可积。Elgindi(2021,数学年刊)证明了 $C^{1,\alpha}$ 数据的有限时间奇点形成。真正的突破,但低于光滑($C^\infty$)阈值。光滑欧拉解是否在3D中爆破仍然是开放的。

湍流。这是欧拉与纳维-斯托克斯比较在物理上变得生动、几乎有形的地方。

在湍流中,能量在大尺度(管道的大小、机翼、风暴)进入,并级联到越来越小的涡旋。这就是能量级联,它是所有物理学中最引人注目的现象之一。在级联的最底部,粘性最终将动能转化为热量。终点线。

欧拉捕捉了惯性范围动力学:由非线性驱动的跨尺度能量传递。但它没有粘性截止。级联没有底部。没有在任何确定尺度上将动能转化为热量的机制。能量是否仍然可以在无粘极限中耗散,即所谓的异常耗散,仍然是一个深刻的开放问题。

这就是为什么湍流建模几乎总是使用纳维-斯托克斯。雷诺数 $\mathrm{Re} = UL/\nu$ 告诉你级联有多宽:高 $\mathrm{Re}$ 意味着从能量输入到粘性消耗之间有许多数量级的尺度。真正的湍流存在于向下倾泻能量的无粘级联和在最小尺度上破坏能量的粘性截止之间的张力中。

Kolmogorov 1941年的理论预测惯性区间中的能谱为$E(k) \sim \varepsilon^{2/3} k^{-5/3}$,该区域既不受大尺度强迫也不受粘性耗散支配。耗散尺度$\eta = (\nu^3/\varepsilon)^{1/4}$设定了该区间的下界。在$\eta$以下,粘度占优。

Euler方程是形式上的$\nu = 0$模型。但将湍流状态下的Navier-Stokes解解释为收敛到Euler解确实很微妙,能量耗散是否在此极限中持续存在(反常耗散)的问题正是Onsager猜想的内容。刚性方面,即对于$u \in C^{0,\alpha}$且$\alpha > 1/3$没有耗散,由Constantin, E和Titi在1994年证明。柔性方面,即在$C^{1/3}$以下构造耗散的Euler解,由Isett在2018年完成,建立在De Lellis和Székelyhidi的凸积分程序之上。

对于Navier-Stokes,级联图景直接连接到能量平衡:$\varepsilon = \nu \langle |\nabla u|^2 \rangle$。未解决的问题很明确。Navier-Stokes解是否保持足够长时间的光滑性,使得Kolmogorov的统计理论在数学上得到证明?存在性与光滑性问题在某种程度上正是在问这个问题:粘性耗散是否足够强大,能在每个尺度上、在所有时间里驯服级联。

欧拉和纳维-斯托克斯之间的差异是一项:$\nu \Delta u$。那一项改变了一切。

欧拉不是简化的纳维-斯托克斯。它是一个根本不同的系统,恰好共享其大部分结构。选择在实践中很重要:选错模型(在粘性重要时用欧拉,在粘性不重要时用纳维-斯托克斯)可能完全破坏模拟。完整方程参见什么是纳维-斯托克斯方程?。障碍参见为什么它很难。奖项参见千禧年问题。不可压缩与可压缩流动参见不可压缩与可压缩纳维-斯托克斯。

粘性项 $\nu \Delta u$ 将一阶非线性系统(欧拉)转换为半线性抛物系统(纳维-斯托克斯)。你得到的是:能量耗散、高阶先验估计,以及Leray在1934年建立并由Fujita和Kato在1964年扩展的纳维-斯托克斯理论的解析半群结构。

然而,这种正则化在3D中不足以闭合全局正则性论证。差得远。克雷千禧年问题精确地问纳维-斯托克斯中的抛物平滑是否可以控制所有尺度上的非线性能量传递,在所有时间。对于欧拉,平行问题同样鲜明:完全缺乏平滑是否保证从光滑数据的有限时间爆炸?没有人知道。

两个问题都位于流体数学理论的中心,比较它们可以准确揭示粘性给你买到了什么和没有买到什么。对于任一系统的3D正则性问题,仍然是整个分析中最难的开放问题之一。