跳转至

验证

本站出现的每个标志性数字都来自一道真实运行过的验收门,而不是设计意图。三道头条门锚定整个项目——其中第二道本身捆着四道内部门:每道门都把桥与一个独立基线相比,而且只有当双方运行同一个钉死的 Julia 1.12.7 时才算数——跨 Julia 版本的比较测的是版本漂移,不是桥。

门一:CJ 爆速

反应流收官示例(reacting_flow)模拟 Chapman–Jouguet 爆轰:Fortran Euler 驱动器 + 内嵌 Clawpack 核,经 Strang 分裂与 Catalyst 化学耦合。

量 值
仿真 CJ 爆速 7.12391
解析 CJ 爆速 7.12470
偏差 0.011%

仿真波速与解析 CJ 值的偏差仅 0.011%,意味着算子分裂、分裂内部的桥穿越、以及 Julia 化学积分器,全部在离散误差范围内保持物理不损。

门二:桥接对原生——四道内部门零偏差

GRI 点火算例(gri_ignition)运行两遍:一遍过桥,一遍作为原生 Julia 基线。四道内部验收门比较输出。

设定 结果
双方 Julia 版本 1.12.7(钉死)
比较的门数 4
发现的偏差 0
残差 < 1e-12

低于 1e-12 的残差已到浮点非结合性消失、精确一致开始的量级:桥没有改变任何算术。

为什么版本钉死是指标的一部分

早先一次对拍失败,原因与桥毫无关系:一个版本垫片把 Julia 静默解析到了与桥内嵌版本不同的 minor 版本上,基线随之漂移。自那以后,每次比较的双方都钉死同一个 1.12.7 二进制。

门三:GRI 点火延迟

同一套机构计算 GRI-Mech 3.0 均相混合气(53 物种、325 反应)的点火延迟:

  • t_ign = 1.1075 ms,取 argmax |dT/dt|——最大放热率时刻,这是不依赖阈值约定的标准锐利定义。

从 325 反应的网络算出物理上合理的点火延迟,是一次端到端检查:机理解析(ChemMechSim.jl,已开源;本次验证使用 v1.2.0)、网络构建(Catalyst)、刚性积分(OrdinaryDiffEq)、以及数据回流 Fortran 数组的通路。

汇总

指标 值 验收门
CJ 爆速 7.12391 对 7.12470,偏差 0.011% reacting_flow,Clawpack1d + Strang 分裂
桥接对原生 4 门 0 偏差,残差 < 1e-12 gri_ignition,双方同钉 1.12.7
GRI 点火延迟 t_ign = 1.1075 ms(argmax abs(dT/dt)) gri_ignition,GRI-Mech 3.0,ChemMechSim.jl v1.2.0

复现步骤、计时方法与交付形态阶梯(JIT 17.3 s → sysimage 1.90 s → AOT 首调约 12 μs)见用户指南与开发者附件。