涡量-应变对齐 · 信息熵方法 · 文献地图
第六步结束时我们到了一个有趣的观察：数值模拟显示涡量 ω 倾向于与应变率张量 S 的中间特征向量对齐，而非最大特征向量。这个现象如果严格成立，可能会抑制涡量拉伸。
这是对齐研究方向最早的严格结果之一。
如果涡量方向 ξ = ω/|ω| 在空间上是 Lipschitz 连续的，且 Lipschitz 常数足够小，则解保持光滑。
直观含义：如果涡量的方向变化不太剧烈（即涡线不太"弯曲"），则涡量拉伸被抑制，解保持光滑。
为什么有效：涡量拉伸项 ω · S · ω 不仅依赖 |ω| 和 ‖S‖，还依赖 ω 相对于 S 的"方向"。如果方向变化平缓，涡量拉伸的"相干性"被破坏，能量无法在一点集中。
Constantin-Fefferman 的结果要求 先验地 假设涡量方向是 Lipschitz 的。但我们不知道这是否成立——这正是我们要证明的事情。
所以这不是一个完整的正则性判据，而是一个"如果...则..."的条件结果。
后续改进：
让我更仔细地看看为什么 ω 会与 S 的中间特征向量对齐。
S 是 3×3 对称矩阵，有三个实特征值 λ1 ≤ λ2 ≤ λ3，满足 λ1 + λ2 + λ3 = 0（因为 ∇·u = 0 → tr(S) = 0）。
因此必有：λ1 ≤ 0 ≤ λ3。λ2 的符号不确定。
涡量拉伸速率：
如果 ω 完全沿 e3（最大拉伸方向），则 ω · S · ω = λ3|ω|2，最大增长。
如果 ω 完全沿 e1（最大压缩方向），则 ω · S · ω = λ1|ω|2 ≤ 0，涡量被压缩。
如果 ω 完全沿 e2（中间方向），则 ω · S · ω = λ2|ω|2，增长速率取决于 λ2 的符号。
数值模拟表明：在湍流中，ω 倾向于与 e2 对齐。但 λ2 的符号分布是：正的概率略大于负的概率。
如果 λ2 > 0，则 ω · S · ω > 0，涡量仍然增长，但速率比沿 e3 时小得多。
问题：即使对齐到 e2，只要 λ2 > 0，涡量仍在增长。我们需要的不是"增长变慢"，而是"增长的积分有界"。
所以对齐现象可能 延缓 爆破，但不一定 阻止 爆破。
前一步提出的信息论框架：把 NS 方程理解为能量谱在尺度间的信息传递。现在让我更具体地尝试。
定义归一化的能量谱：
其中 Etot(t) = ½∫ |u|2 dx 是总动能。在无外力无耗散情况下，Etot 守恒。有耗散时 Etot 递减。
定义能量谱的 Shannon 熵：
H(t) 度量能量在不同尺度上的"分散程度"。H 越大，能量越均匀分布在各个尺度上；H 越小，能量越集中在少数尺度上。
对 H(t) 求导：
从 Fourier 空间的 NS 方程可以得到 ∂E/∂t 的方程：
其中 T(k) 是能量转移项（来自非线性卷积），满足 ∫ T(k) dk = 0（非线性项不改变总能量，只是重新分配）。
因此：
由于 ∫ T(k) dk = 0 且 ∫ log p · T(k) dk 的符号不确定，这个方程无法简化为一个确定的符号判断。
问题在于能量转移项 T(k)：它既可以把能量从低频推到高频（增熵），也可以从高频拉回低频（减熵）。净效应取决于 u 的具体结构，而不是一个通用不等式。
在 Kolmogorov 湍流的统计稳态中，能量从大尺度注入、小尺度耗散，中间有一个能量级串（energy cascade）区域。在这个区域中，T(k) 的净效应是向高频传递，熵增加。
但 NS 存在性问题不是在统计意义下问的，而是在逐点意义下问的。即使平均而言熵增加，也不排除在某个特定时刻能量突然向某个小尺度极端集中。
结论：信息熵方法在统计/平均意义下有用，但对于逐点正则性问题，不够精确。
Shannon 熵太"粗"了。试试 Fisher 信息，它对分布的"尖锐程度"更敏感：
Fisher 信息在 p(k) 高度集中时很大（尖峰），在 p(k) 均匀分布时很小。如果 I(t) 在有限时间内发散，意味着能量在某个波数 k* 处形成 delta 函数——这对应物理空间中的奇点。
De Bruijn 恒等式连接了 Shannon 熵和 Fisher 信息：
这意味着：如果耗散足够强，I(t) 有界，则 H(t) 单调递减（能量趋于均匀分布）。如果耗散不够强，I(t) 可能增长，能量可能集中。
Fisher 信息方法比 Shannon 熵更有希望，因为它直接探测"尖峰形成"。
但问题依然存在：要证明 I(t) 有界，需要对 T(k)（非线性能量转移）做精细估计——而这又回到了老问题。
可能的突破口：如果能证明 T(k) 有一个内在的"展宽效应"（即非线性卷积天然地使能量谱展宽而非尖锐化），那么即使初始条件很尖锐，演化也会自动展宽，阻止 delta 函数形成。
但这需要证明一个关于 NS 非线性结构的具体性质，而不是通用估计——这正是 Tao 说必须的。
在上述探索中，我们重新发现了很多前人已知的结果。现在系统地画一张地图。
| 路线 | 代表人物/年份 | 核心思想 | 进展 | 障碍 |
|---|---|---|---|---|
| 能量方法 | Leray (1934), Hopf (1951) | L2 能量不等式，弱解存在性 | 弱解存在（能量有限） | 弱解是否唯一/光滑？未知 |
| 正则性判据 | Prodi (1959), Serrin (1962), Ladyzhenskaya (1967) | 如果 u ∈ Lp(Lq) 满足 Serrin 条件，则光滑 | 条件正则性：满足某些范数有界则光滑 | 无法先验地证明这些范数有界 |
| 涡量方向 | Constantin-Fefferman (1993), Cordoba-Fefferman (1998) | 如果涡量方向变化不太剧烈，则光滑 | 条件正则性：方向 Lipschitz 则光滑 | 无法先验地证明方向正则性 |
| BKM判据 | Beale-Kato-Majda (1984) | 如果 ∫‖ω‖L∞ dt < ∞ 则光滑 | 最简洁的正则性判据 | 无法先验地控制 ‖ω‖L∞ |
| 小初值 | Fujita-Kato (1964), Cannone (1995) | 小初值时全局光滑解存在 | 小初值已解决 | 大初值？未知 |
| 轴对称 | Ukhovskii-Yudovich (1968), Hou-Lei (2009) | 轴对称无 swirl 时全局光滑 | 轴对称无 swirl 已解决 | 轴对称有 swirl？未知 |
| 调和分析 | Constantin (1990), Lemarie-Rieusset (2002) | Littlewood-Paley分解，Besov空间 | 在临界 Besov 空间中部分进展 | 临界空间中的估计不够紧 |
| 平均化模型 | Tao (2016) | 构造共享 NS 对称性但会爆破的模型 | 证明了通用方法不可能成功 | 不直接解决真实 NS |
| 有限时间爆破候选 | Hou-Chen (2014), Elgindi (2019) | 数值模拟/分析构造爆破候选 | Elgindi 证明了 C1,α 速度的轴对称 NS 在特定边界条件下爆破 | C1,α 不是光滑初值 |
| 深度学习 | 近期（2020s） | PINNs、神经算子、数据驱动正则性检测 | 在数值模拟中有帮助 | 不能替代数学证明 |
Tarek Elgindi (2019, Duke Math. J.) 证明了：对于具有 C1,α 速度场的轴对称无 swirl Navier-Stokes 方程，在特定边界条件下，存在有限时间爆破的解。
这是一个真正的数学定理——不是数值模拟，不是条件结果。
但是：C1,α 初值不是光滑初值（光滑意味着 C∞）。C1,α 是 Hölder 连续的导数，但不是无限可微的。
所以 Elgindi 的结果告诉我们：在"接近光滑但不够光滑"的初值条件下，NS 方程确实可以爆破。但光滑初值是否爆破？仍然未知。
意义：这至少说明 NS 方程的"可爆破性"不是遥不可及的——它在 C1,α 边界上就已经发生了。光滑初值的情况可能只差"一步"。
Thomas Hou 和 Guo Luo (2014) 在三维欧拉方程（ν = 0 的 NS，即无黏性情况）的数值模拟中发现了可能的有限时间爆破候选。随后：
从第一步到第十步，让我画一张完整的"探索地图"。
| 步骤 | 内容 | 状态 |
|---|---|---|
| 1 | 涡量方程：找到三维灾难源 (ω·∇)u | ✅ 已知 |
| 2 | Lp 能量估计：遇到循环论证 | 🚨 撞墙 |
| 3 | 信息论框架：能量放大器 vs 衰减器 | ⚠️ 框架建立，精确估计困难 |
| 4 | 涡量-应变几何分析：S 的分解 | ✅ 标准推导 |
| 5 | Tao 平均化模型：通用方法不可能 | ✅ 已排除一条死路 |
| 6 | 诚实总结：墙的形状 | ✅ 认知收获 |
| 7 | 涡量-应变对齐：Constantin-Fefferman | ⚠️ 条件结果，无法先验验证 |
| 8 | 信息熵方法：Fisher 信息 | ⚠️ 统计有用，逐点不够 |
| 9 | 文献地图：10条技术路线 | ✅ 地图画出来了 |
| 10 | 综合评估 | — |
我们没有证明任何东西。但我们比出发时更清楚地知道：