结果族 376：Universal computation in forced Navier–Stokes flows · 学科：Partial differential equations · 验证状态：暂无形式化证明，请以社区核验为准
想象一个全自动的水下仓库：水流是搬运工，把一个个密封箱从 A 区挪到 B 区，箱子既压不扁也撞不散。这篇论文证明，只要外力设计得当，这样一股三维水流就能”跑程序”——计算机的每一步计算被翻译成箱子的搬运，程序算完的那一刻，一个预先放好的示踪粒子会被水流恰好冲进指定的”答案方块”。
关键词卡片
看个具体例子
机器的三个栈被编码成三个实数（进制 @@M@@B@@）。比如一条规则在左栈压入一位（@@M@@x\mapsto 3x+2@@）、在右栈弹出一位（@@M@@y\mapsto (y-d)/3@@）、历史栈不动——一个方向放大 3 倍、另一方向缩小 3 倍，三个栈合起来的行列式恰为 1，正是流体搬得动的”盒子”。最终判据：粒子 @@M@@x_*=(4,0,0)@@ 进入开立方体 @@M@@(-1,2)^3@@ 当且仅当机器停机；不停机时它全程待在 @@M@@x_1>2@@ 一侧。
为什么值得关心
它把”通用计算”完整地藏进了最经典的流体方程：理论上，追踪一个粒子的行踪就能读出任意机器的停机答案，而所需外力只是查有限表格、从不真的执行那台机器。
暂无形式化证明（AI 结果待核验）
第一级台阶
导读里的水下仓库可以改述成半数学的语言：对象是三维空间中成批的实心盒子，操作是把盒子整体平移并沿三个方向各自伸缩，且任何一步都不改变体积——水池挤不出也吸不进水，这一条由流体本身担保。”跑程序”即按指令表把每个盒子从源区搬到指定的目标区，示踪粒子只负责报告终局。
第二级台阶
用本科术语精确化：每条机器指令诱导一个保体积仿射映射——线性部分是沿坐标轴的缩放 diag(s₁,s₂,s₃) 且 s₁s₂s₃=1，再加一个平移。流体一侧的主角是无散（solenoidal）光滑速度场：散度处处为零，等价于其流的时间一映射保体积，恰好搬得动这类盒子。图灵机的三个栈（左带、右带、历史）用奇数基 B 的进位展开写成三个实数，控制状态与停机态安置在固定的格点位置，”机器执行一步”就成了”某个盒子被搬走”。
第三级台阶
主定理在受迫 Navier–Stokes 方程的语境下陈述：粘性 ν>0 事先给定，流体从静止出发，外力 f=f₀+νf₁ 光滑、只有一个紧空间支集、在加载时段 [0,1] 之后以 1 为周期。定理断言：固定粒子 x*=(4,0,0) 进入开立方体 (−1,2)³ 当且仅当被编码的机器停机，且解在声明的比较类中唯一。盒子的有理参数与正间隙挡住了构造中求逆和取极限的麻烦；”有效性”——外力求值只需查有限规则表、从不执行机器——则挡住了循环定义的陷阱。
证明有限族”正对角、行列式为 1”的有理实体盒仿射映射——源族与像族各自不交、但两族可互相重叠——都能被有效的光滑不可压流实现为时间一映射；据此用平衡三栈记录器构造周期外力下从静止出发的三维 Navier–Stokes 流，固定粒子进入固定开立方体当且仅当图灵机停机。
一条机器指令对应盒子上的保体积仿射映射，但流体实现必须搬动整个盒子的邻域而不只是编码点、在盒外保持不可压，并且在源盒与目标盒互相重叠时仍可执行；报告答案的粒子还必须在一切非停机过渡中避开观测区域。Cardona、Miranda、Peralta-Salas 与 Presas（2021）用收缩–搬运–扩张加 Moser 面积校正实现了平面圆盘上的版本；Moore 的广义移位（generalized shift）提供了从符号规则到位置展开的框架，Landauer 与 Bennett 的历史保留原理解决可逆性。本文把问题推进到三维平直空间中的实体盒（solid box），并把有效性（可计算的系数与求值程序）和全时段观测证明一并纳入定理陈述：几何固定、粘性 @@M@@\nu>0@@ 给定、初速为零，机器只进入外力。
几何主定理（盒子路由）：设有理实体盒 @@M@@Q_i=a_i+\prod_j[-h_{ij},h_{ij}]@@ 与 @@M@@Q_i'=b_i+\prod_j[-k_{ij},k_{ij}]@@（半宽度均为正，两族各自两两不交），尺度比 @@M@@s_{ij}=k_{ij}/h_{ij}@@ 满足 @@M@@\prod_js_{ij}=1@@。则存在有效的紧支集光滑无散（solenoidal）速度场，在时段 @@M@@[1/4,3/4]@@ 外为零，其时间一映射在每个 @@M@@Q_i@@ 的整邻域上恰为 @@M@@a_i+y\mapsto b_i+\operatorname{diag}(s_{i1},s_{i2},s_{i3})\,y@@。
计算主定理（平衡三栈实现）：对每台机器与有限输入，存在有效系数 @@M@@f_0,f_1@@ 使 @@M@@f=f_0+\nu f_1@@ 在 @@M@@\mathbb R^3@@ 上有显式的零初速解 @@M@@(U,0)@@；力与速度具有单一紧空间支集、各阶混合导数有界，并在 @@M@@t>1@@ 后以 1 为周期；固定粒子 @@M@@x_*=(4,0,0)@@ 进入开立方体 @@M@@(-1,2)^3@@ 当且仅当机器停机；解在比较类 @@M@@u\in CH^2\cap C^1L^2@@（速度、时间导数与二阶空间导数连续，压力模常数属 @@M@@CH^1@@）中唯一。对固定机器，周期场与输入无关，输入只进入加载段。文中还给出变体：固定环面图表版本、从零时刻起即周期的初始化、持久与守卫式历史域、有界盒与平板（slab）观测者，以及由满慢钟（onto slow clock）导出的导数时间平方可积的衰减外力。
几何证明分三层。第一层是局部旋度（localized curl）：设运动盒子 @@M@@K_i(t)=c_i(t)+D_i(t)K_i^0@@，其中 @@M@@D_i@@ 正对角、行列式为一，对势 @@M@@\chi_i\{\tfrac12\dot c_i\times y+\tfrac13(L_iy)\times y\}@@（@@M@@L_i=\dot D_iD_i^{-1}@@ 无迹）取旋度，在平台上恰得仿射速度 @@M@@\dot c_i+L_iy@@；不交的平台支撑使各盒子的场可直接叠加。第二层是调度：先以公因子 @@M@@\Lambda@@ 把全部源中心扩张到彼此分离，在分离中心处完成形状的几何插值（半宽按 @@M@@h_{ij}^{1-\theta}k_{ij}^{\theta}@@ 单调变化），再把盒子逐个撤往远离两族的停车中心，全部停妥后逐个交付到扩张后的目标中心，最后公因子收缩还原。第三层是避障：单个盒子的直线路径若穿过静盒中心周围的 @@M@@4R@@ 立方体障碍，就沿其表面绕行——进入点、顶点、棱、出口点连成有理折线——得到与一切静中心距离至少 @@M@@4R@@ 的路径。源族与像族可以重叠的根本原因在时间调度：交付开始前所有源已全部撤离。
计算部分的关键是”平衡”。以基数 @@M@@B=2|\Sigma|+1@@ 的奇数字把三个栈 @@M@@L@@（头左侧带）、@@M@@R@@（头右侧带与边界 @@M@@\star@@）、@@M@@J@@（历史记录）编码为实数三元组，控制状态置于 @@M@@o_s=(4i_s,0,0)@@、停机态放在原点；一条规则替换三个有限前缀，前缀长度总和不变，故仿射映射的行列式为 @@M@@B^{\sum|\alpha_j|-\sum|\alpha'_j|}=1@@，恰好落入几何定理的适用范围。可逆历史这样取得：先标记并右扫，从工作带边界之外取一个新空白推上 @@M@@J@@，回扫还原工作栈，最后才用记录符号 @@M@@\kappa@@ 替换该空白并执行工作指令；一步普通计算耗 @@M@@2m+2@@ 条辅助规则，记录还能复原磁头位移。停机检测的全时段正确性来自坐标估计：一切源码与非停机目标的第一个坐标至少为 4，而扩张、变形与停车路径的界表明非停机分支全程 @@M@@x_1>2@@，中途不可能误入 @@M@@(-1,2)^3@@。最后在 @@M@@[0,1]@@ 用平移场把 @@M@@x_*@@ 加载到初始有理码，从 @@M@@t=1@@ 起逐周期重复盒场；残差恒等式 @@M@@f=U_t+(U\cdot\nabla)U-\nu\Delta U@@ 给出外力（并按 @@M@@\nu@@ 线性分解为 @@M@@f_0+\nu f_1@@），能量比较引理给出比较类中的唯一性。力的求值只查有限规则表，从不执行机器。
本文主结果未形式化，请以社区核验为准（族元数据提及配套 Lean 文档页，但本篇标记为未形式化）。本文与族内姊妹篇共享编译器、进制编码、比较类唯一性与慢钟等框架，并在引言中与速度场积分检测计算的工作互相对照观测事件，文末还显式导入姊妹篇的平面处理器与空间钟完成输入偏移应用。按 OpenAI 官方声明，未经形式化的结果可能有问题。
问 1：为什么每条指令必须”平衡”到行列式恰为 1？松一点不行吗？
答：不可压速度场的时间一映射保体积，流体只能实现体积不变的搬运。三栈编码下，一条规则替换的三个前缀总长度不变，仿射映射的行列式才为 B 的零次幂即 1，恰好落进几何定理的适用范围（引言明确陈述这一点）。去掉平衡，整套编译就断了。
问 2：源盒族与目标盒族允许互相重叠，为什么不出乱子？
答：靠时间调度：先把全部源中心扩张分离，完成形状插值，再把盒子逐个撤往远离两族的停车中心，全部停妥后才逐个交付（引言概括的 separates centers—changes shapes—evacuates—delivers 流程）。重叠在空间上被允许、在时间上被错开，这是全文最脆弱的一环。
问 3：全时段观测证明防的是什么？最容易在哪里出错？
答：防的是粒子在两条规则之间的搬运途中误入观测立方体。论文靠坐标估计排除：源码与非停机目标的首坐标至少为 4，扩张、变形与停车路径的界保证非停机分支全程 x₁>2。这类跨文件拼接的显式常数最易出 AI 式失配：各引理各自成立，接缝处差一个因子就前功尽弃。
问 4：外力”有效”究竟何意？会不会偷偷执行了机器？
答：有效指一个有限程序能把外力的任意阶混合导数算到任意有理误差，并给出证明所用的界；系数由机器与输入决定，但求值过程不运行机器（引言 effectiveness 小节）。论文同时明确不主张对无界计算的一致精度容忍度。
问 5：与 Cardona–Miranda–Peralta-Salas–Presas 等前人相比，真正新意在哪？
答：前人的普适计算或用自适应几何与平稳 Euler 流，或用无限能量的 Beltrami 流；本文在固定平直度量、给定粘性、零初速下把机器整个放进外力，并把有效性、全时段观测与比较类唯一性一并写进定理。三维实体盒上的直接正对角实现与有界盒、平板、环面等观测变体也是新的。
问 6：唯一性为何只在比较类 CH²∩C¹L² 中陈述？
答：论文自己声明不主张任意流体数据的正则性定理；能量比较引理针对声明的类，配合外力的光滑性与紧支集给出该类中的唯一性。大外力下三维方程的整体强唯一性本非已知，限缩比较类是诚实之举。
问 7：这个结果最可能以什么方式出错，谁会用它？
答：最可能在调度、绕行与编码交界处的常数与边界情形。下游：它把”固定粒子是否进入固定立方体”变成停机问题的物理化身，与姊妹篇速度场积分检测的工作互为对照；受迫流体不可判定性方向的工作会直接引用。物理可实现性与数值方面论文未展开。
总评
把通用计算完整装进固定平直空间、给定粘性、从静止出发的受迫 Navier–Stokes 流，且有效性、观测与唯一性一并交付，分量十足；最大风险在多文件构造接缝处的常数与完全未形式化的证明，须经逐引理的人工核验。