🔢 mathematics

Nonuniqueness of Leray-Hopf solutions to the unforced incompressible 3D Navier-Stokes Equation

本文通过开发创新的数值方法并结合严格的计算机辅助证明，在自相似框架下构造了无强迫不可压缩三维纳维 - 斯托克斯方程的勒雷 - 霍普夫解，并通过验证线性化算子存在不稳定特征对，首次严格证明了该方程解的非唯一性。

原作者： Thomas Hou, Yixuan Wang, Changhe Yang

发布于 2026-03-20

📖 1 分钟阅读🧠 深度阅读

原作者： Thomas Hou, Yixuan Wang, Changhe Yang

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

这篇论文讲述了一个流体力学领域困扰了人类近一个世纪的“大谜题”的突破。为了让你轻松理解，我们可以把这篇论文的内容想象成一场关于**“水流是否唯一”**的侦探游戏。

1. 核心谜题：水流只有一种走法吗？

想象你在浴缸里放了一团墨水，或者在河里扔了一块石头。水会怎么流动？

传统观点（Leray-Hopf 解）： 数学家们认为，只要初始条件（比如你扔石头的位置和力度）确定了，水流的路径就应该是唯一的。就像你按下一个开关，灯只会亮一种模式。
新发现（非唯一性）： 这篇论文的作者（Thomas Hou 等人）发现，在某些极端情况下，同一个初始状态，竟然可以演化出无数种完全不同的水流路径！ 就像你按下一个开关，灯可能亮红灯，也可能亮绿灯，甚至可能同时亮两种颜色，而且这两种状态都是完全符合物理定律的。

这在数学上被称为**“非唯一性”**，是 3D 纳维 - 斯托克斯方程（描述流体运动的核心方程）中最大的未解之谜之一。

2. 他们是怎么做到的？（三个关键步骤）

作者没有用纯笔算（因为太难了），而是开发了一套**“超级计算机辅助证明”的方法。我们可以把这个过程比作“寻找并验证一个不稳定的平衡点”**。

第一步：寻找“走钢丝”的平衡点（自相似解）

想象你在走钢丝。通常，如果你稍微歪一点，你会掉下来（不稳定）。但在某些特殊位置，你可能刚好能保持平衡。

作者首先用超级计算机在数学空间里“寻找”一个特殊的流体状态（称为自相似解）。这个状态就像是一个完美的、静止的“走钢丝”姿势。
难点： 这个姿势必须非常精确，稍微有点误差，整个证明就崩塌了。

第二步：测试“推一把”会发生什么（线性化算子与特征值）

找到了这个平衡点后，作者问：“如果我轻轻推它一下，它会怎样？”

如果它只是晃一下然后回到原位，那就是稳定的（唯一解）。
如果它像多米诺骨牌一样，一旦推倒就彻底改变方向，那就是不稳定的。
作者通过复杂的数学计算（特征值分析），发现这个特殊的流体状态确实极其不稳定。只要有一点点扰动，它就会分裂成不同的路径。

第三步：用计算机做“数学公证”（计算机辅助证明）

这是这篇论文最牛的地方。以前大家只能算出“大概是这样”，但数学要求绝对精确。

传统方法： 计算机算出结果，但可能有舍入误差（就像用尺子量东西，永远有 0.01 毫米的误差）。
作者的方法： 他们发明了一种**“区间算术”和“严格误差控制”**的方法。
- 想象一下，他们不是告诉计算机“答案是 5"，而是告诉计算机“答案在 4.999 到 5.001 之间”。
- 然后，他们把整个数学证明过程拆解成无数个小块，用计算机严格验证每一个小块都在误差允许范围内。
- 这就好比盖房子，他们不仅画了图纸，还拿着游标卡尺，精确到微米级地检查了每一块砖、每一根钢筋，确保房子绝对不会塌。

3. 一个生动的比喻：分岔路口

想象你开车在一个分岔路口：

以前的理论： 只要你的车（初始条件）一样，你只能走左边那条路，或者只能走右边那条路，不可能同时走两条。
这篇论文的发现： 作者发现了一个特殊的路口（那个不稳定的自相似解）。在这个路口，如果你以某种特定的速度开过去，你可以选择向左转，也可以向右转，甚至可以在路口原地打转。
更神奇的是，作者不仅找到了这个路口，还严格证明了在这个路口，向左和向右都是合法的驾驶行为（都符合物理定律），而且这两种行为可以无限细分，产生无数种可能的驾驶轨迹。

4. 为什么这很重要？

打破认知： 这推翻了我们对流体运动“确定性”的直觉。它告诉我们，在极端条件下，流体可能具有内在的“随机性”或“多义性”。
数学界的里程碑： 这是人类第一次用计算机辅助的方式，严格证明了 3D 纳维 - 斯托克斯方程存在非唯一解。这就像在数学界投下了一颗重磅炸弹。
方法创新： 他们开发的这套“严格验证”的数学工具，未来可以用来解决其他超级复杂的科学问题（比如预测天气、设计飞机等），因为这套方法能保证计算结果绝对可信。

总结

简单来说，Thomas Hou 和他的团队利用超级计算机和严密的数学逻辑，证明了在 3D 空间中，同样的水流起点，可以演化出无数种不同的未来。他们不仅找到了这个“奇迹”，还像最严格的法官一样，用计算机给这个结论盖上了“绝对真实”的印章。

这就好比他们发现了一个魔法开关，按下去后，世界不再只有一种走向，而是充满了无限的可能性。

这篇论文《无强迫不可压缩三维纳维 - 斯托克斯方程 Leray-Hopf 解的非唯一性》（Nonuniqueness of Leray–Hopf solutions to the unforced incompressible 3D Navier–Stokes Equation）由 Thomas Hou, Yixuan Wang 和 Changhe Yang 撰写，发表于 2026 年 3 月。该论文在数学流体力学领域取得了突破性进展，首次通过严格的计算机辅助证明（Rigorous Computer-Assisted Proof），证明了无强迫三维纳维 - 斯托克斯（Navier-Stokes, NS）方程的 Leray-Hopf 弱解具有非唯一性。

以下是该论文的详细技术总结：

1. 研究问题 (Problem)

三维不可压缩纳维 - 斯托克斯方程解的存在性与唯一性是千禧年大奖难题之一。虽然 Leray 和 Hopf 证明了全局弱解（Leray-Hopf 解）的存在性，但这些解是否唯一长期以来是一个悬而未决的开放问题。

核心挑战：在无强迫（unforced）且初始数据为光滑紧支集的情况下，证明存在多个不同的 Leray-Hopf 解。
背景：之前的非唯一性结果多依赖于凸积分（Convex Integration）方法（构造低正则性解）或引入外部强迫项（Forced NS）。对于无强迫且满足能量不等式的“合适”解（Suitable solutions），其非唯一性尚未被严格证明。

2. 方法论 (Methodology)

作者采用了一种结合自相似变量（Self-similar variables）、谱分析和严格数值验证的创新策略。

2.1 自相似框架与谱假设

自相似解构造：研究形式为 $u(t, x) = \frac{1}{\sqrt{t}} \bar{U}(\frac{x}{\sqrt{t}})$ 的自相似解。将 NS 方程转化为关于 $\bar{U}$ 的椭圆型方程组。
非唯一性机制：根据 Jia-Sverak 和 Albritton-Brue-Colombo 的理论，如果线性化算子 $L_{\bar{U}}$ 在自相似变量下存在一个实部为正的离散特征值（即不稳定的特征模态），则可以通过在该不稳定流形上扰动，构造出第二个解，从而证明非唯一性。
目标：找到一个精确的自相似剖面 $\bar{U}$ ，使得其线性化算子 $L_{\bar{U}}$ 具有不稳定的特征值。

2.2 严格的计算机辅助证明框架

为了克服纯数值模拟无法保证数学严谨性的问题，作者开发了一套完整的计算机辅助证明框架：

高精度数值近似：
- 使用有限元方法（FEM）在截断域上计算候选剖面 $\bar{U}$ 和特征向量 $v$ 。
- 利用坐标变换 $r = \tan \beta$ 将无穷远域映射到有限域，并使用三角张量积基底进行插值，确保解的光滑性和散度为零（Divergence-free）的精确性。
- 通过梯度下降法优化参数，寻找具有负特征值（对应线性化算子的不稳定模态）的构型。
后验误差分析与存在性证明：
- 残差控制：计算数值解的残差 $E$ ，并证明其足够小。
- 线性算子分解：将线性化算子 $L$ $L$ 分解为强制部分（Coercive part, $L_1$ $L_{1}$ ）和紧扰动部分（Compact perturbation, $L_2$ $L_{2}$ ）。
  - 利用粘性项提供的强耗散性， $L_1$ 在 $L^2/H^1$ 范数下具有显著的强制系数（Coercivity coefficients）。
- 有限秩近似：将紧算子 $L_2$ 近似为一个有限秩算子 $L_{ap}$ ，并严格估计近似误差。
- Schauder 不动点定理：在强制部分和紧扰动部分的基础上，构建不动点映射，证明在数值解附近存在一个精确解。
区间算术（Interval Arithmetic）：
- 在计算积分、特征值、范数上界等关键步骤时，使用区间算术处理舍入误差，确保所有数值界限在数学上是严格成立的。

3. 主要贡献 (Key Contributions)

首个无强迫情形的严格证明：这是首次针对无强迫三维 NS 方程，在光滑紧支集初始数据下，严格证明 Leray-Hopf 解的非唯一性。之前的类似工作（如 Albritton-Brue-Colombo）需要引入奇异强迫项，而本文完全去除了这一限制。
创新的数值验证框架：
- 提出了一种将数值计算与严格分析紧密结合的新范式。
- 利用粘性的强制性质（Coercivity），使得线性算子的谱验证比欧拉方程（Euler equation）的奇点验证（如 Chen-Hou 的工作）更加简单和稳健。
- 成功处理了无穷远域边界条件和散度约束的严格满足问题。
构造了具体的不稳定剖面：通过数值优化，找到了一个具体的自相似剖面 $\bar{U}$ ，并严格证明了其线性化算子存在一个实部为负的特征值（在方程符号约定下对应不稳定模态），从而打破了唯一性。

4. 主要结果 (Key Results)

定理 1（非唯一性）：存在无穷多个不同的、合适的 Leray-Hopf 解，它们对应于同一个具有紧支集的光滑无散初始条件 $u_{in}$ 。这些解在 $t>0$ 时是光滑的，但违反了保证唯一性的 Prodi-Serrin 条件。
定理 2（存在性）：存在一个弱解 $(\bar{U}, \bar{v}, \bar{\lambda})$ 满足自相似方程组，其中 $\bar{\lambda}$ 是线性化算子的不稳定特征值。
数值常数验证：
- 验证了残差 $\|E_U\|_{L^2} \le 6.5 \times 10^{-7}$ 和 $\|E_v\|_{L^2} \le 1.9 \times 10^{-5}$ 。
- 验证了线性算子的强制系数（Coercivity coefficients）： $U$ 方程约为 0.051， $v$ 方程约为 0.164。
- 验证了特征值 $\bar{\lambda} \approx -0.113$ （对应不稳定模态）。
代码公开：所有验证代码已公开在 GitHub 上，确保了结果的可复现性。

5. 意义与影响 (Significance)

理论突破：该结果从根本上挑战了关于纳维 - 斯托克斯方程解唯一性的传统认知，表明即使在无强迫、初始数据光滑且满足能量不等式的情况下，解的唯一性也可能失效。
方法论示范：展示了计算机辅助证明（Computer-Assisted Proof, CAP）在解决非线性偏微分方程（PDE）核心难题中的巨大潜力。它证明了通过严格的数值验证，可以处理传统解析方法难以触及的复杂非线性问题。
对奇点理论的启示：虽然本文关注的是非唯一性，但其构造的自相似剖面与潜在的有限时间奇点（Blow-up）机制密切相关。理解这种不稳定性对于最终解决 NS 方程的奇点问题（即解是否会在有限时间内爆破）提供了新的视角和工具。
对比欧拉方程：相比于无粘欧拉方程的奇点证明，粘性 NS 方程的强制性质使得数值验证更加可行，这为未来研究粘性流体中的其他开放问题提供了新的思路。

总结

这篇论文通过结合高精度的数值计算、严格的误差分析和不动点理论，首次为无强迫三维纳维 - 斯托克斯方程 Leray-Hopf 解的非唯一性提供了确凿的数学证据。这不仅解决了流体力学中的一个核心开放问题，也标志着计算机辅助证明在数学物理领域迈上了一个新的台阶。