结果族 376：Universal computation in forced Navier–Stokes flows · 学科：Partial differential equations · 验证状态：暂无形式化证明，请以社区核验为准
能不能把一段计算机程序”编译”进一杯水里？这篇论文造出了这样的编译器：输入任意图灵机，输出一套搅水的外力。水里一个固定的小颗粒就是指示灯——程序停机，它迟早漂进指定的探测带；程序不停机，它就永远不进去。于是”这杯水会不会亮灯”和停机问题一样无法预判。
关键词卡片
看个具体例子
把机器状态铺在一条横轴上：左端 x<3/16 是”还在算”区，右端 3/4≤x<13/16 是”停机”方块，中间 1/2<x<1 是探测带。不停机的颗粒全程 x≤9/32，只在左区打转；停机的颗粒则在有限时间后被精确搬进右端方块。每条指令由四个剪切实现 (r,s)→(λr, s/λ)：先把矩形抬到专属高度，改形后平移、落下，各指令互不打架。
为什么值得关心
它把纯逻辑的”不可判定”装进了最经典的物理方程：没有任何算法能预判这杯水里的颗粒是否会进入某片区域。构造还格外干净——压强恒为零、外力无散度且自初始化后周期运行，机器与输入全部只藏在外力里。
暂无形式化证明（AI 结果待核验）
第一级台阶
导读说这篇论文造了一台”把程序编译进水里”的编译器。翻译成半数学语言：对象是一杯不可压缩的黏性流体，装在一个固定不动的三维环面”杯子”里，从完全静止出发；操作是一套预先写好的光滑外力，像多只看不见的手按时刻表分层错动液体；被观测的对象是其中一个指定颗粒走的路线。程序的全部信息只藏进这份外力时刻表里。
第二级台阶
用本科术语精确化：速度场 u 满足强迫 Navier–Stokes 方程 u_t+(u·∇)u=−∇p+νΔu+f，在平坦环面上无散度、初速为零，粘度 ν 取任一给定正数。”剪切”指速度只依赖一个坐标、只推动另一个坐标的场，它把每层平移一段与层高成正比的距离；四个剪切按时刻表复合，恰好实现”横坐标缩 λ 倍、纵坐标缩 1/λ 倍”的保面积仿射变形。颗粒轨迹是常微分方程 Ẋ=u(t,X) 的解。
第三级台阶
主定理在此陈述：存在一个终止的算法，把任一带头移动限于左、停、右的单带图灵机及有限输入，编译成光滑、无散度、零均值的外力 f；方程有全局光滑解，压强恒为零，且在古典解类中唯一；一切混合导数全局有界，初始化之后外力以一为周期。假设各自的作用：可逆记录器补回被机器抹掉的历史，因为光滑流的有限时间映射必可逆；每条指令专属的”高度层”挡住源矩形与目标矩形重叠的冲突；偏移估计保证中间轨线始终不越出 2h 邻域，挡住”路过探测带造成误报”。结论：固定颗粒进入固定开带，当且仅当机器停机。
在固定平坦环面、固定正粘度、零初速的强迫 Navier–Stokes 流中，本文给出一个终止的编译器，把任意图灵机与输入变成光滑、无散度、零均值的外力，使固定粒子在有限时刻进入固定开带当且仅当机器停机，从而该”固定粒子事件”不可判定。
停机问题（halting problem）的不可判定性是计算的边界；能否让一个具体物理系统”演示”这一边界？Cardona–Miranda–Peralta-Salas–Presas 在适配的黎曼三球面上构造了模拟图灵机的定常欧拉流，Dyhr 等人又用 Hodge 粘性得到无外力定常 Navier–Stokes 的图灵完备性，但两者都要改动几何、且用非零初速。悬而未决的是：在固定平坦环面（flat torus）@@M@@\T^3@@、固定正粘度 @@M@@\nu@@、流体从静止出发的强迫体系里能否实现通用计算。障碍有二：机器指令可能抹除信息，而光滑流的有限时间映射必可逆；一条指令的源区域也可能与另一条指令的目标区域重叠。本文以”保留历史＋分层高度”逐一化解，给出一条从普通机器到流体的完整构造链。
主定理（定理 1.1，initialized computation by a solenoidal force）：固定正的可计算粘度 @@M@@\nu@@，存在算法把带头移动限于 @@M@@\{-1,0,1\}@@ 的单带确定图灵机 @@M@@M@@ 与有限输入 @@M@@w@@ 编译为有效光滑外力 @@M@@f_{M,w}@@，满足：(i) @@M@@f_{M,w}@@ 无散度（solenoidal）且零空间均值，方程 @@M@@u_t+(u\cdot\nabla)u=-\nabla p+\nu\Delta u+f@@、@@M@@u(0)=0@@ 有全局光滑解 @@M@@(u,p)=(U,0)@@——压强恒为零——且在古典类中唯一；(ii) @@M@@U@@ 与 @@M@@f_{M,w}@@ 的一切混合导数全局有界，动能一致有界，两场自 @@M@@t\ge1@@ 起周期为一，周期机制只依赖 @@M@@M@@；(iii) 固定粒子标签 @@M@@a_*=(1/8,3/8,1/2)@@ 在某有限时刻进入固定开带 @@M@@O=\{(x,y,z):1/2<x<1\}@@，当且仅当 @@M@@M@@ 停机于 @@M@@w@@（含初始即停机的情形）。推论：对该编译器产出的一切力，无算法能判定固定粒子事件。
构造分四步。第一步可逆化：仿照 Bennett 的可逆模拟（reversible simulation），设计七行转移表的单头记录器（recorder），其带符号 @@M@@(a,\ell,m)@@ 同时携带工作符号、历史与返回标记；每执行一条原始指令，就把其编号抄写到不断右移的历史前沿上，第 @@M@@n@@ 个检查点（checkpoint）精确代表原机第 @@M@@n@@ 步后的构形，每步耗时 @@M@@2(r_0+n-h')+4@@ 个微步，且该表在其整个定义域上可由”目标状态＋写下的符号”唯一回溯前驱。第二步位置编码：取 Moore 广义移位（generalized shifts）视角，把带头两侧的带子写成 @@M@@B@@ 进制小数，状态方块按 @@M@@x@@ 分置（非终态在 @@M@@1/8\le x<3/16@@，终态在 @@M@@3/4\le x<13/16@@），每条记录器指令化为互逆仿射映射（reciprocal affine map），线性部分为 @@M@@\diag(\lambda_i,\lambda_i^{-1})@@、@@M@@\lambda_i=B^{-d_i}@@；引理证明源矩形族与目标矩形族各自两两分离，但源与目标允许重叠。第三步剪切实现：七个时间上依次分离的脉冲完成”抬升—四剪切—平移—下降”：先把源矩形抬到该分支专属高度 @@M@@z_i@@，再用四个剪切 @@M@@\mathsf Y(-\lambda),\ \mathsf X(\lambda^{-1}-1),\ \mathsf Y(1),\ \mathsf X(\lambda-1)@@ 精确实现 @@M@@(r,s)\mapsto(\lambda r,s/\lambda)@@，平移至目标中心后降回 @@M@@z_0@@。关键的偏移估计（excursion estimate）断言中间各坐标偏离不超过 @@M@@2h@@ 且与 @@M@@\lambda@@ 无关；每段速度都横向于自身变化的方向，故对流项（convective term）@@M@@(V\cdot\nabla)V\equiv0@@，压强恒为零，@@M@@f=\partial_tV-\nu\Delta V@@ 自动无散度、零均值；源与目标的重叠因使用时刻不同、中间隔有专属高度而互不干扰。第四步装载与归纳：一次性装载流把 @@M@@a_*@@ 沿两段直线送到有理的输入编码，此后周期机制重复；归纳证明每个微步粒子恰落在一个源矩形内并被送到正确的后继编码。不停机的运行始终处于非终态，全程 @@M@@x\le 9/32<1/2@@，绝不触发探测带；停机则在有限个周期后进入 @@M@@x\ge3/4@@ 的终态方块。解的唯一性由 Leray 式能量比较（energy comparison）加 Gronwall 不等式得到。
本文暂无形式化证明；按 OpenAI 官方声明，未经形式化的结果可能有问题，请以社区核验为准。同族两篇姊妹篇分别以”普通机器＋第三坐标历史”和”可逆记录表＋空间时钟”独立实现同类事件等价，三篇共享”编译器→预设速度→残差力→能量法唯一性”的骨架，互为印证；本文独有的亮点是压强恒为零、外力无散度且初始化后周期。结论针对精确物质轨线（material trajectory）的可达性，编码使用精确实坐标，不含有限精度鲁棒性。
问 1：证明中最关键也最脆弱的一步是什么？
（答：定理 thm:shears 里”精确轨迹”一段：作者写出候选路径，验证它满足局部化常微分方程，再用 ODE 唯一性把它认同为真实轨迹；其前提是偏移引理保证四剪切全程不越出 2h、从而水平截止函数的平台条件成立。若某个 λ 极端值下偏移界失效或平台判断有漏，”精确实现仿射映射”便崩塌，停机等价随之失效。）
问 2：为什么必须先做可逆化，直接模拟原始机器不行吗？
（答：不行。光滑流的有限时间映射是单射，而原始图灵机的转移可以合并不同输入（不可逆）。论文用七行转移表的单头记录器把每条已执行指令的编号抄到右移的历史前沿上（Bennett 技巧），并在整个定义域上显式证明”目标状态＋写下的符号”可唯一回溯前驱——这是后文目标矩形彼此分离（间隔 ≥κ/B²）的组合基础。）
问 3：与已有流体计算工作相比，真正的新意在哪里？
（答：Cardona–Miranda–Peralta-Salas–Presas 的图灵完备欧拉流要改造成适配的黎曼三球面且用非零初速；Dyhr 等人的无外力黏性计算要对度量做形变并引入 Hodge 粘性。本文在固定的平坦环面、固定正粘度、零初速下，把机器全部编进外力，且压强恒为零；论文明言未引用前人任何几何或解析结论，链条自足。）
问 4：压强恒为零是怎么白拿的，会不会太巧？
（答：每个速度场都”横向于自身变化的方向”，故对流项 (V·∇)V 恒等于零；各脉冲时序分离，交叉对流项也不出现，于是 f=∂ₜV−νΔV 自动无散度、零均值，(V,0) 直接解方程。这是构造的漂亮之处，也意味着它只覆盖这一族预设速度，论文同样明示：这是强迫可达性构造，不对任意 Navier–Stokes 数据的正则性作任何断言。）
问 5：这个结果最可能以什么方式出错？
（答：AI 生成证明的典型风险是组合细节算错：七行表在”任意合法带子”上的逐情形反演、B 进制区间间隔的 κ/B 与 κ/B² 估计、每原始步 2(r₀+n−h′)+4 个微步的计数、以及 δ 邻域论证里五个 margin 的极小值。任何一处缺一个情形或差一个因子，等价性就可能只在”好带子”上成立。且本文未形式化。）
问 6：下游推论是什么，谁会用到？
（答：直接推论是该编译器产出力类上”固定粒子事件”的不可判定性（把判定算法与终止编译器复合即得停机判定，论文并援引图灵的符号打印问题）。用家是可计算性分析与流体可计算性方向；但需注意结论针对精确实数坐标上的物质轨线，论文明确不做任何有限精度鲁棒性声明，物理可实现性仍隔一层。）
总评
这是”机器到流体”完整编译链的一次干净实现，亮点是平坦几何、零初速、零压强与初始化后周期性全部同时成立；最大风险在于长组合构造中未经形式化核验的区间间隔与微步计数细节。