Perelman熵的构造 · Böchner公式的障碍 · 人工几何化 · 熵单调性
要构造 NS 的 Perelman-style 熵，必须先理解 Perelman 的 W-熵在 Ricci 流中做了什么。
Ricci 流：∂gij/∂t = -2Rij
Perelman 定义了 W-熵泛函（对度量 g 和标量函数 f）：
其中 τ 是后向时间参数，R 是标量曲率，n 是维数。
Perelman 的关键发现：如果 f 满足后向热方程 ∂f/∂τ = Δf - |∇f|2 + R，则：
W 沿 Ricci 流单调不减。
Perelman 的 W-熵有效，依赖三个深层结构：
这三个结构都是几何的——它们依赖 Riemann 流形的结构。
NS 方程没有 Riemann 结构。速度场 u 是 R3 上的向量场，不是某个流形的度量。没有 Riemann 曲率，没有 Böchner 公式，没有 Bianchi 恒等式。
飞轮说："也许可以构造一个人工几何结构——将速度场解释为某个黎曼流形上的向量场。"让我试试。
在 Riemann 几何中，Christoffel 符号 Γkij 编码了度量的导数。也许可以把 u 的某个导数解释为 Γ。
具体地，设 g 是 R3 上的某个度量，使得：
如果这个对应成立，则 NS 方程 ∂u/∂t + (u·∇)u = -∇p + νΔu 可能对应某个"度量流"方程。
障碍：Γkij 是对称的（Γkij = Γkji），但 ∂iuk 不对称。∂iuk 可以分解为对称部分 Sik 和反对称部分 Aik：
换一个思路。不是把 u 解释为度量的导出量，而是从 u 构造一个度量。
定义：gij = δij + ε·∂iuj + ε·∂jui = δij + 2ε·Sij
其中 ε 是小参数。这给出了一个度量，其偏离欧氏度量的部分由应变率张量决定。
这个度量的 Ricci 曲率（到 ε 的一阶）：
而 NS 方程的耗散项是 νΔu。如果我们选择 ε 使得 -εΔSij 与 νΔu 相关联...
障碍：Ricci 曲率涉及二阶导数，而 NS 的耗散项也是二阶导数。但 NS 的非线性项 (u·∇)u 是一阶导数的二次项，而 Ricci 曲率的非线性部分也是二阶导数的二次项。两者结构不匹配。
从 u 构造度量 g 的方法遇到了根本困难：
既然几何化方法失败了，也许可以直接在 NS 方程上构造一个类似 W-熵的泛函——不需要度量，只需要 u 本身。
受 Perelman W-熵启发，定义：
其中 f 是辅助函数（类似 Perelman 的 f），τ 是后向时间参数。
这里用 |∇u|2 替代了 Perelman 的 |∇f|2 + R。|∇u|2 在 NS 中的作用类似于"能量密度"。
计算 dWNS/dτ：
从 NS 方程：∂u/∂t = -(u·∇)u - ∇p + νΔu。
代入后：
展开各项：
Perelman 的 W-熵有效，是因为所有项组合后恰好是一个完全平方（|Ric + ∇²f - g/(2τ)|² ≥ 0）。
NS 版本中：
更根本的问题：Perelman 的熵是单调递增的（dW/dτ ≥ 0），这给出了熵的上界，从而控制了曲率的增长。如果我们构造的 NS 熵是单调递减的，它给出的是熵的下界——不能阻止能量（或梯度）增长。
要得到单调递增的熵，需要熵的被积函数中耗散项的符号与对流项"合作"。但在 NS 中，耗散项总是阻尼（减少能量），而对流项总是放大（增加能量）。两者的效应相反。
Perelman 的 W-熵灵感来自 Boltzmann 的 H-定理。也许直接迁移 H-定理更可行。
Boltzmann 方程：∂tf + v · ∇xf = Q(f,f)
H-定理：H(t) = ∫ f log f dv dx 满足 dH/dt ≤ 0，等号仅当 f 是 Maxwellian 分布时成立。
H-定理有效的原因：碰撞算子 Q(f,f) 满足微观可逆性（detailed balance），导致熵产生项非负。
NS 与 Boltzmann 的关系：NS 方程是 Boltzmann 方程在 Knudsen 数 → 0 时的宏观极限（Chapman-Enskog 展开）。在这个极限下，Boltzmann 的 H-函数退化为 NS 的能量泛函。
在 Chapman-Enskog 展开中，NS 的能量方程是：
这正是 NS 的能量不等式——能量单调递减。
但这不是我们需要的。能量递减是已知的（Leray 1934）。我们需要的是某种"更高级"的量的单调性——比如梯度的某种函数。
考虑：
这是 H1 Sobolev 范数的平方。如果 H1[u] 有界，则 u 在 H1 中有界。但 H1 在 3D 中不是临界空间——需要控制更高阶的范数。
计算 dH1/dt：
对流项 -∫ ∇u · ∇[(u·∇)u] dx 展开后包含 ∫ ∂iu · ∂i(uj∂ju) dx 类型的项。分部积分后得到 ∫ ∂iuk · ∂iuj · ∂juk dx——这是三个一阶导数的乘积，符号不确定。
dH1/dt 中的对流项 ∫ ∂u · ∂u · ∂u dx 用 Hölder 不等式估计：
在 3D 中，用 Gagliardo-Nirenberg 不等式：
所以对流项 ≤ C · ‖u‖L23/4 · ‖∇²u‖L29/4 = C · E3/4 · ‖∇²u‖L29/4
而耗散项 = -ν‖∇²u‖L22
要使总和非正，需要：ν‖∇²u‖2 ≥ C · E3/4 · ‖∇²u‖9/4
即：‖∇²u‖-1/4 ≥ (C/ν) · E3/4
这等价于：‖∇²u‖2 ≤ (ν/C)4 · E-3
但这正好是反向的不等式！我们需要 ‖∇²u‖ 足够大来控制对流项，但实际得到的是 ‖∇²u‖ 必须足够小——这与我们的目的矛盾。
这又回到了那个老问题：在临界情况下，非线性项和耗散项的尺度相同，无法用耗散压制非线性。
五次尝试，五次都失败了。
| 尝试 | 思路 | 失败原因 |
|---|---|---|
| 26 | 理解 Perelman W-熵 | 依赖 Böchner 公式——NS 没有 |
| 27 | 从 u 构造度量 g | 涡量无法编码；维数不匹配；非线性结构不同 |
| 28 | 直接构造 NS 的 W-熵 | 耗散项符号反了（熵递减而非递增）；对流项符号不确定 |
| 29 | Boltzmann H-定理迁移 | 能量递减已知；H1 估计遇到循环论证 |
| 30 | 综合评估 | — |
Perelman 的方法之所以成功，是因为 Ricci 流是纯几何的——它只依赖 Riemann 流形的内在结构。Böchner 公式、Bianchi 恒等式、曲率演化方程——这些都是几何恒等式，不依赖具体的方程形式。
NS 方程是物理的——它描述流体运动，不描述空间几何。速度场 u 不是度量，不是联络，不是曲率。它是 R3 上的向量场。
要把 NS 方程"几何化"，需要找到一个几何对象（度量、联络、曲率、或其他），其演化方程恰好是 NS 方程。但 NS 方程的物理结构（对流 + 压力 + 耗散）与已知的几何流方程（Ricci 流、平均曲率流、Yamabe 流等）的结构完全不同。
除非：存在一个全新的几何框架，在其中 NS 方程是一个"自然"的几何流方程。这需要发明新的几何语言——而这不是我们在这个探索系列中能做的事。
Ricci 流类比和 Perelman-style 熵泛函方法不可行——结构性障碍太深。
但有一个收获：H-定理的迁移尝试（第二十九步）揭示了能量估计在 3D 中遇到循环论证的具体机制——对流项的三次积分与耗散项的二次积分在临界尺度上恰好平衡。这为理解 NS 困难提供了一个清晰的定量框架。
飞轮应转向第三方向：Elgindi C1,α 爆破裂值的 C∞ 磨光策略。