纳维–斯托克斯方程

1. 为什么重要（千禧年难题）

绕过船身的水和绕过机翼的空气遵循同一种局部平衡：质量守恒，力改变动量。纳维—斯托克斯方程把这一思想写成连续介质模型。写出方程，远比证明每一个解可能具有怎样的行为容易。

研究状态核查日期：2026年10月13日。 9月8日，OpenAI宣布构造了具有光滑外力的三维不可压缩纳维—斯托克斯方程的有限时间爆破解，并发布论文和Lean代码。该结果针对克雷千禧年难题的选项 C和D。这是一项声称取得重大进展的成果，但它并没有证明无外力纳维—斯托克斯方程会发生爆破。另行公布的无外力欧拉方程结果对应零黏性，不应与之混淆。OpenAI公告（9月10日更新）。

这里必须区分结果的发布、形式化命题的检查、可供人理解的证明，以及数学界的认可。第14节按日期介绍评估状态，包括10月出现的审视与讨论。下面先从物理上说明，为什么能量有界仍可能伴随越来越剧烈的局部运动。

---

2. 连续介质场与守恒定律

设 $\rho (\mathbf{x},t)$ 为质量密度，$\mathbf{v}(\mathbf{x},t)$ 为速度，$p(\mathbf{x},t)$ 为压强，$\mathbf{\sigma }$ 为柯西应力。单位质量所受的体力记为 $\mathbf{f}$，例如重力。

场为每个位置和时刻赋予一个数值，代表对大量分子的平均。这个模型假定，所关注的长度尺度足够大，使这种平均和局部应力定律具有意义。数学奇点讨论的是连续介质方程，并不预言现实中的某个分子会达到无限速度。

2.1 质量守恒

$${\partial }_t\rho +\nabla \cdot (\rho \mathbf{v})=0.$$

第一项表示局部质量积累的速率，散度项表示向外的净质量通量。密度恒定时，质量守恒化为速度散度为零。流体微团可以拉伸、改变形状，同时保持体积不变。

2.2 动量守恒（柯西方程）

$$\rho \left({\partial }_t\mathbf{v}+(\mathbf{v}\cdot \nabla )\mathbf{v}\right)=\nabla \cdot \mathbf{\sigma }+\rho \mathbf{f}.$$

将应力分解为各向同性的压强部分与黏性部分：

$$\mathbf{\sigma }=-p\mathbf{I}+\mathbf{\tau }.$$

物质导数 $D_t={\partial }_t+\mathbf{v}\cdot \nabla$ 跟随运动中的流体微团。即使速度场不随时间变化，微团也可能加速：进入狭窄通道的水，会移动到速度不同的区域。这就是非线性输运项 $(\mathbf{v}\cdot \nabla )\mathbf{v}$ 的含义。

---

3. 本构关系（牛顿流体）

对于各向同性的牛顿流体，

$$\mathbf{\tau }=2\mu \mathbf{S}+\lambda (\nabla \cdot \mathbf{v})\mathbf{I},\mathbf{S}\equiv \frac{1}{2}\left(\nabla \mathbf{v}+\nabla {\mathbf{v}}^{\top }\right),$$

其中 $\mu$ 为动力黏度，$\lambda$ 为第二黏性系数；体积黏度为 $\zeta =\lambda +2\mu /3$。

当黏性系数在空间中恒定时，代入动量平衡可得可压缩纳维—斯托克斯方程：

$$\rho \left({\partial }_t\mathbf{v}+\mathbf{v}\cdot \nabla \mathbf{v}\right)=-\nabla p+\mu {\nabla }^2\mathbf{v}+(\mu +\lambda )\nabla (\nabla \cdot \mathbf{v})+\rho \mathbf{f}.$$

当 $\rho$ 为常数时，定义运动黏度 $\nu \equiv \mu /\rho$。

黏性响应的是形变，而不是匀速平移。其量纲为 $[\mu ]=\mathrm{Pa}\mathrm{s}$ 和 $[\nu ]={\mathrm{m}}^2/s$。在平板流动示意图中，若间距为 $h$、上板速度为 $U$，则 $v_x(y)=Uy/h$，剪切应力为 ${\tau }_{xy}=\mu U/h$。这个简单的精确解是光滑的。奇点存在的主张针对某些特定数据，并不是说每一种流动都会出现奇点。

---

4. 不可压缩形式

对于恒定密度 $\rho ={\rho }_0$，且满足

$$\nabla \cdot \mathbf{v}=0,$$

时，可得

$${\partial }_t\mathbf{v}+(\mathbf{v}\cdot \nabla )\mathbf{v}=-\frac{1}{{\rho }_0}\nabla p+\nu {\nabla }^2\mathbf{v}+\mathbf{f},\nabla \cdot \mathbf{v}=0.$$

可以把这条方程读成：加速度等于压强引起的加速度、黏性扩散和外加加速度之和。

项	含义
${\partial }_t\mathbf{v}$	在固定位置测得的时间变化
$(\mathbf{v}\cdot \nabla )\mathbf{v}$	流体微团穿过速度场而产生的变化
$-\nabla p/{\rho }_0$	压力使动量重新分配
$\nu {\nabla }^2\mathbf{v}$	黏性使空间变化趋于平滑
$\mathbf{f}$	外加的单位质量所受力

不可压缩并不要求速度或加速度为零。压强作为约束，使速度场保持无散；在这个恒密度模型中，它不是由独立的状态方程确定的。

对动量方程取散度，得到压强泊松方程，并需配合适当的边界条件：

$${\nabla }^{2}p=-{\rho }_0\nabla \cdot [(\mathbf{v}\cdot \nabla )\mathbf{v}]+{\rho }_0\nabla \cdot \mathbf{f}.$$

---

5. 涡量与螺度

定义涡量 $\mathbf{\omega }\equiv \nabla \times \mathbf{v}$。 对于不可压缩流动（$\nabla \cdot \mathbf{v}=0$）：

$${\partial }_t\mathbf{\omega }=\nabla \times (\mathbf{v}\times \mathbf{\omega })+\nu {\nabla }^2\mathbf{\omega }+\nabla \times \mathbf{f}.$$

分量形式为：

$${\partial }_t{\omega }_i+v_j{\partial }_j{\omega }_i={\omega }_j{\partial }_j{v}_i+\nu {\partial }_{jj}{\omega }_i+(\nabla \times \mathbf{f}{)}_i.$$

其中 ${\omega }_j{\partial }_j{v}_i$ 是涡旋拉伸项，仅在三维中出现。 螺度 $H\equiv \int \mathbf{v}\cdot \mathbf{\omega }dV$ 可用来考察涡线的拓扑结构。

沿轴向被拉伸的涡管可能变细，旋转也可能增强。在涡量方程中，拉伸可以放大旋转，而黏性使它扩散。仅凭其中任何一种效应，都无法判断完整的耦合动力学是否始终光滑。

对于真正的平面流动 $\mathbf{v}=(v_x(x,y),v_y(x,y),0)$，涡量垂直于平面，且各场在该方向上没有变化。因此 $(\mathbf{\omega }\cdot \nabla )\mathbf{v}=0$。这使二维流动少了一种主要的放大机制，却不意味着每个三维涡旋都必然爆破。

---

6. 能量平衡（不可压缩情形）

对于光滑解，将动量方程与 $\mathbf{v}$ 点乘，并在具有静止无滑移边界的区域 $\Omega$ 上积分。在周期盒中，或具有充分衰减的全空间中，也有相同的恒等式：

$$\frac{1}{2}\frac{d}{dt}{\int }_{\Omega }|\mathbf{v}{|}^{2}dV+\nu {\int }_{\Omega }|\nabla \mathbf{v}{|}^{2}dV={\int }_{\Omega }\mathbf{f}\cdot \mathbf{v}dV.$$

当 $\mathbf{f}=0$ 时，动能单调下降；黏性耗散率为 $\nu \Vert \nabla \mathbf{v}{\Vert }_2^2$。

这里的能量已除以恒定密度；物理动能是速度平方积分的 ${\rho }_0/2$ 倍。输运项和压强项只贡献边界通量，而在这些假设下，边界通量为零。分部积分给出非负的黏性项。运动的壁面或开放边界则会带来额外的功或通量项。

6.1 为什么有限能量不能保证光滑性

积分的界控制的是运动总量，而不是局部最大值。作为尺度变换的例子，取一个固定的、光滑紧支撑无散场 $\mathbf{\varphi }$，定义

$${\mathbf{v}}_{\ell }(\mathbf{x})={\ell }^{-1}\mathbf{\varphi }(\mathbf{x}/\ell ).$$

则当 $\ell \to 0$ 时，$\Vert {\mathbf{v}}_{\ell }{\Vert }_{\infty }={\ell }^{-1}\Vert \mathbf{\varphi }{\Vert }_{\infty }$ 增大，而 $\Vert {\mathbf{v}}_{\ell }{\Vert }_2^2=\ell \Vert \mathbf{\varphi }{\Vert }_2^2$ 减小。这些只是缩放后的场的快照，并非方程的解，也不是爆破证明。它们说明，速度在越来越小的体积中变得越来越大，并不一定违反能量上界。正则性问题还需要解决导数与集中的控制。

---

7. 无量纲化与控制参数

选取特征长度 $L$ 和特征速度 $U$，定义 ${\mathbf{x}}^{\prime }=\mathbf{x}/L$，$t^{\prime }=tU/L$，${\mathbf{v}}^{\prime }=\mathbf{v}/U$，$p^{\prime }=p/(\rho U^2)$。

不可压缩方程变为

$${\partial }_{t^{\prime }}{\mathbf{v}}^{\prime }+({\mathbf{v}}^{\prime }\cdot {\nabla }^{\prime }){\mathbf{v}}^{\prime }=-{\nabla }^{\prime }{p}^{\prime }+\frac{1}{\mathrm{Re}}{\nabla }^{\prime 2}{\mathbf{v}}^{\prime }+{\mathbf{f}}^{\prime },{\nabla }^{\prime }\cdot {\mathbf{v}}^{\prime }=0,$$

其中雷诺数为

$$\mathrm{Re}=\frac{UL}{\nu }.$$

这里 ${\mathbf{f}}^{\prime }=L\mathbf{f}/U^2$。扩散时间为 $L^2/\nu$，平流时间为 $L/U$，两者之比就是 $\mathrm{Re}$。雷诺数大，表示在所选尺度上，平流比黏性平滑作用快。它既不是普适的湍流阈值，也不是奇点判据：几何形状、扰动和测量尺度同样重要。例如，$U=0.1m/s$、$L=0.01\mathrm{m}$、$\nu =10^{-6}{\mathrm{m}}^2/s$ 时，$\mathrm{Re}=1000$。

其他按需使用的无量纲数包括：

• 马赫数 $\mathrm{Ma}=U/c_s$（可压缩性）。
• 弗劳德数 $\mathrm{Fr}=U/\sqrt{gL}$（重力）。
• 普朗特数 $\mathrm{Pr}=\nu /\kappa$（动量扩散与热扩散之比）。

---

8. 可压缩纳维—斯托克斯—傅里叶方程概述

加入内能 $e$ 或温度 $T$，以及傅里叶热通量 $\mathbf{q}=-k_{\mathrm{th}}\nabla T$。其中 $k_{\mathrm{th}}$ 是热导率，不同于上文的热扩散率 $\kappa$：

$${\partial }_t(\rho e)+\nabla \cdot (\rho e\mathbf{v})=-p\nabla \cdot \mathbf{v}+\mathbf{\tau }:\nabla \mathbf{v}-\nabla \cdot \mathbf{q}.$$

再用状态方程使系统闭合，例如理想气体的 $p=(\gamma -1)\rho e$。

此时，压缩功以及动能向内能的转化都显式地发挥作用。不可压缩千禧年难题的结果，并不能解决这个不同系统中的激波形成、热输运或正则性理论。

---

9. 数学理论（存在性与唯一性）

局部存在指解在初始时刻之后的某个时间区间存在；全局存在指解在以后所有时刻都存在。经典解具有足够的导数，可以逐点满足方程；这里的光滑指任意阶可微。弱解满足积分形式的方程，把导数转移到光滑测试函数上。这些是不同层次的存在性概念。

9.1 弱解（勒雷—霍普夫）

在标准的全空间或周期情形中，对于无散、有限能量的初始数据及允许的外力，勒雷—霍普夫理论给出满足能量不等式的全局弱解。光滑解在保持光滑的时间内，在相应的能量类中是唯一的，这称为弱—强唯一性。因此，即使可能失去光滑性，也不意味着所有广义解都不再存在。费弗曼撰写的官方问题说明。

唯一性依赖解的类别、初始数据和外力。Albritton、Brué和Colombo 证明了某些外力下的非唯一性。Hou、Wang和Yang 报告了针对奇异初始数据、在无外力情形下的计算机辅助非唯一性结果。这类结果不能改称为光滑初始数据在无外力下产生爆破的证明。

9.2 正则性判据（三维）

• 涡量延拓判据（Beale—Kato—Majda型）： 在通常的全空间或周期情形中，假设具有充分正则性，并且外力直到时刻 $T$ 及跨越该时刻都足够正则，那么经典解在 $T$ 破裂必须满足

$${\int }_0^T\Vert \mathbf{\omega }(\cdot ,t){\Vert }_{L^{\infty }}dt=\infty .$$

Luo关于纳维—斯托克斯延拓判据的讨论陈述了无外力情形的判据；上面对外力的限定排除了由奇异外部输入强行导致的破裂。

• Prodi—Serrin判据： 对于具有适当外力的能量类弱解，若 $\mathbf{v}\in L^p(0,T;L^q)$，且

$$\frac{2}{p}+\frac{3}{q}\le 1,q>3,$$

则解在 $(0,T]$ 上光滑。

• Caffarelli—Kohn—Nirenberg定理： 在定理对外力的假设下，适当弱解的奇点集合具有零一维抛物型豪斯多夫测度。这是部分正则性，并不证明奇点集合为空。官方问题说明中的背景与参考文献。

这些都是条件性命题。它们说明什么样的控制能够阻止破裂，并没有为所有三维流动提供这种控制。上面的指数反映了自然尺度变换 ${\mathbf{v}}_a(\mathbf{x},t)=a\mathbf{v}(a\mathbf{x},a^{2}t)$；在三维中，单靠能量范数并不具有临界尺度性质。

9.3 二维全局正则性

在标准的全空间或周期情形中，对于光滑且满足相容条件的数据及足够正则的外力，二维不可压缩流动在每个有限时刻都保持光滑。无外力时，涡量平方积分 $\int |\mathbf{\omega }{|}^{2}dV$ 不增；外力会增加源项。由于没有涡旋拉伸，其控制比三维情形强得多。

9.4 可压缩流动

全局有限能量弱解在一定条件下存在，例如满足 $\gamma$ 律的气体；但正则性与唯一性仍大多未解决。

9.5 克雷的四个选项究竟问什么

官方表述要求在正黏性下，证明以下四个命题中的任意一个：

选项	区域	外力与问题
A	${\mathbb{R}}^3$	零外力；所有允许的光滑初始数据是否都具有全局光滑有限能量解
B	周期区域 ${\mathbb{R}}^3/{\mathbb{Z}}^3$	零外力；所有允许的光滑周期初始数据是否都具有全局光滑解
C	${\mathbb{R}}^3$	是否存在某些允许的光滑数据与光滑快速衰减外力，使能量一致有界的全局光滑解不存在
D	周期区域 ${\mathbb{R}}^3/{\mathbb{Z}}^3$	是否存在某些允许的光滑周期数据与外力，使全局光滑解不存在

在全空间中，初始数据与外力须满足指定的衰减条件；周期外力也须满足时间衰减条件。C和D允许非零外力，因此证明其中之一，并不从逻辑上否定无外力命题A或B。“解决克雷的一个选项”和“解决所有版本的正则性问题”是不同的主张。

---

10. 湍流与能谱

在高 $\mathrm{Re}$ 的均匀各向同性湍流理论（K41）中，惯性区间的能谱为

$$E(k)=C_K{\epsilon }^{2/3}{k}^{-5/3},$$

其中 $\epsilon$ 是平均耗散率，$C_K\sim \mathcal{O}(1)$ 是柯尔莫哥洛夫常数。

柯尔莫哥洛夫微尺度：

$$\eta \sim {\left(\frac{{\nu }^3}{\epsilon }\right)}^{1/4}.$$

使用K41估计 $\epsilon \sim U^3/L$ 可得 $L/\eta \sim {\mathrm{Re}}^{3/4}$。在三维中解析这一尺度，意味着仅空间网格点数就约为 ${\mathrm{Re}}^{9/4}$ 阶，还未计入时间推进的成本。这种估算促使人们采用不同分辨率与建模选择，包括直接数值模拟（DNS）、**大涡模拟（LES）和雷诺平均纳维—斯托克斯（RANS）**方法。LES和RANS需要闭合假设来处理未解析的相关量。

$-5/3$ 能谱描述的是惯性区间的统计模型，并不适用于每一种流动和每个尺度。间歇性与边界会使这一图景更复杂。专门构造的奇异解不会直接提供普适的湍流模型，也不会让数值误差估计变得多余。

---

11. 边界条件与几何

• 无滑移壁面： $\mathbf{v}{|}_{\partial \Omega }=0$。
• 自由滑移： $\mathbf{n}\cdot \mathbf{v}=0$，$\mathbf{n}\cdot \mathbf{\tau }\cdot \mathbf{t}=0$。
• 用于理论研究与DNS的周期盒（环面）。
• 流入／流出边界，配以适当的应力或压强条件。

这些选择影响适定性、能量收支和数值稳定性。

周期盒将相对的面视为同一面，并没有物理壁面。在全空间中，无穷远处的衰减代替了壁面条件。这些几何下的结果不能自动移用于飞机表面或管道：边界条件本身就是数学问题的一部分。

---

12. 数值方法（一页指南）

• 投影法： 推进 ${\mathbf{v}}^{\ast }$，求解压强 $p$ 的泊松方程，再投影到 $\nabla \cdot {\mathbf{v}}^{n+1}=0$。
• 能量稳定格式： 采用斜对称形式的对流项来保持离散能量。
• 显式平流与扩散格式的CFL限制： $\Delta t\lesssim \min \{\Delta x/U,\Delta x^2/\nu \}$。
• 散度控制： 使用无散谱基，或使用满足inf–sup稳定性的速度—压强空间对的混合有限元。严格无散需要合适的空间，单有inf–sup稳定性并不能保证这一点。

模拟可以显示峰值增长，但有限网格本身无法确立无限极限。要区分真正的集中和数值不稳定，需要网格加密研究、离散化误差控制及解析论证。计算机辅助证明提供严格的界；形式化证明助手检查从编码的定义和假设出发的推导。它们都不只是高分辨率动画。

---

13. 联系与推广

• 欧拉方程： 令 $\nu =0$ 就去掉了黏性扩散。取极限 $\nu \to 0$ 本身也可能很微妙，尤其是在壁面附近。OpenAI另行宣布的无外力欧拉方程结果讨论的是这个无黏方程，而不是无外力、正黏性的问题。公告与另行发布的欧拉方程论文。
• 相对论流体力学： 朴素的纳维—斯托克斯理论不满足因果性；二阶Israel–Stewart 理论在对输运系数和流体状态施加适当条件时，可以恢复因果性（含体积黏性的严格例子）。
• 量子流体： 采用双流体模型和量子化涡旋的超流体纳维—斯托克斯类似系统。

---

14. 研究现状

截至2026年10月13日： 光滑外力驱动的三维构造已被宣布，正接受检验；无外力三维纳维—斯托克斯方程的正则性问题仍未解决。二维正则性和光滑流动解在日常应用中的价值没有被推翻。

14.1 精确的主张与构造机制

OpenAI论文的定理1.1陈述：对于每个 $\nu >0$，存在一个在空间和时间上均紧支撑的光滑外力，以及从静止出发、定义在 ${\mathbb{R}}^3\times [0,1)$ 上的光滑不可压缩流动，使得

$$\mathbf{v}(\cdot ,0)=0,\underset{0\le t<1}{\sup }\Vert \mathbf{v}(t){\Vert }_{L^2}<\infty ,\underset{t\uparrow 1}{\mathrm{limsup}}\Vert \mathbf{v}(t){\Vert }_{L^{\infty }}=\infty .$$

在时刻1之前，速度与压强都具有固定的紧空间支撑；外力则光滑延伸并跨越该时刻。论文声称证明选项C，并通过推论10.6的周期构造证明D。论文的压强约定吸收了恒定密度因子。

提出的机制使一个向内螺旋运动、沿轴向拉伸的涡旋发生集中。振荡修正及其动量输运，抵消本来会变得奇异的残余力。随意选择一个爆破的速度，再由方程定义外力，并不足够：外力及它的所有阶导数必须在跨越奇异时刻时仍保持光滑。这只是对所声称构造的概念说明，并非独立证明。

“对每种黏性，都存在某个外力”允许外力依赖黏性。它不意味着同一个固定外力能克服任何黏性，不意味着任意小的外力就足够，也不意味着最初静止的流体能在没有能量输入时爆破。

14.2 Lean、数学审视与奖金

已发布的Lean仓库包含全空间与周期情形的目标命题，以及使用Comparator独立检查的说明。形式化检查器验证的是针对精确编码的命题与假设的推导。定义是否表达了原本要研究的问题，以及论证如何解释文字证明，是额外需要考察的问题。本笔记并不声称已经重新构建或独立验证了这个Lean项目。

10月6日，Bastounis、Circelli和Hansen发表文章，批评自然语言阐述与Lean之间的对应关系。其中的例子包括：引理8.6对导数的要求不同，引理10.5中的压强估计采用了不同的论证。这质疑了把形式化视为文字证明每个步骤都已自动获证的做法。但仅凭这一点，不能反驳最终的有外力爆破定理，也不能说明最终形式化命题是错误的。

9月11日，克雷数学研究所的回应欢迎这一看似重大的突破，同时强调其刻意不急于作出判断的评估过程。奖金规则要求符合条件的出版、在该出版之后至少经过两年，以及获得数学界的普遍认可，之后才会考虑授奖。不能仅凭公告日期就认定奖金评估的计时已开始。这里并不声称已经授奖；OpenAI表示无意申领奖金。

14.3 后续研究说明了什么、尚未说明什么

9月的后续研究同样需要注意适用范围。Petrillo和Glimm明确区分有外力的结果与仍然存在的无外力问题。Lei和Ren的阐述研究剖面构造，把进一步的振荡抵消留给后续工作，因此不是完整的独立验证。Constantin、Ignatova和Vicol在特定的几何、增长与外力假设下得到正则性结果，从而阻止某些奇点形成。这与尝试从这一机制中移除外力有关，并不是所有解析外力流动都全局正则的定理。

现在有价值的问题变得具体：所声称的假设和结论是否得到核验，能否透明地解释这一机制，以及移除外力后还有什么能够保留？湍流预测、所有解类别中的弱解唯一性，以及可压缩流动，仍是另外的问题。

15. 检查你的理解

1. 匀速平移： 速度和压强恒定、外力为零时，不可压缩方程中的哪些项消失？所有空间与时间导数都消失，因此任意正黏性都与这种运动相容。
2. 能量与最大速度： 在第6.1节中作变量替换 $\mathbf{y}=\mathbf{x}/\ell$。为什么能量按 $\ell$ 缩放？因为振幅的平方贡献 ${\ell }^{-2}$，三维体积贡献 ${\ell }^3$。
3. 改变尺度： 固定 $U$ 和 $\nu$，若 $L$ 减半，雷诺数怎样变化？它也减半；在更小的尺度上，黏性的相对重要性增加。
4. 阅读定理： 为什么选项C成立时，A仍可能悬而未决？C允许选取一个非零光滑外力，而A只讨论零外力。
5. 证明检查器： Lean接受证明后，还须核查什么？须核查编码的定义、假设和最终结论是否符合原本的数学问题；与另行书写的阐述之间的对应关系也需要审阅。

---