制造解法 (MMS) 代码验证模拟器

参数设置

计算域长度 L

一维计算域 x ∈ [0, L] 的长度

粗网格划分数 N

最粗网格的单元数。后续分四阶细分

细化比 r

相邻网格间h缩小1/r的比率

期望收敛阶数 p_design

格式的理论阶数（中心差分为2）

扩散系数 D

−D·u'' = f(x) 的物性系数

制造解 u_mms(x)

本工具采用 u_mms(x) = sin(πx/L)，反演得 f(x) = D·(π/L)²·sin(πx/L) 作为代码输入。

计算结果

—

粗网格误差 E₁

—

细网格误差 E₄

—

观测收敛阶数 p_obs（平均）

—

期望阶数 p_design

—

MMS验证结果

—

误差比 E₁/E₄

—

网格细分与数值解对比制造解

4阶网格（粗→细）求解 u_mms(x) = sin(πx/L)，将数值解叠加绘制。误差范围随网格细化而缩小的过程可视化。

收敛图 — log E vs log h

网格间观测收敛阶数 p_obs

理论与关键公式

$$f_{mms}(x) = L[u_{mms}(x)],\quad E_k = \|u_h^{(k)} - u_{mms}\|_2,\quad p_{obs} = \frac{\log(E_1/E_2)}{\log r}$$

L为控制方程的差分算子，u_mms为任意光滑函数，p_obs ≈ p_design则验证成功。

$$-D\,u''(x) = f(x),\qquad u_{mms}(x) = \sin\!\left(\tfrac{\pi x}{L}\right),\qquad f(x) = D\!\left(\tfrac{\pi}{L}\right)^{2}\sin\!\left(\tfrac{\pi x}{L}\right)$$

本工具求解的一维定常扩散方程，以及从制造解反演得到的源项f(x)。边界条件为u(0)=u(L)=0。

$$E_k \;\propto\; h_k^{\,p_{design}},\qquad \frac{E_k}{E_{k+1}} = r^{\,p_{obs}}$$

二阶中心差分下E ∝ h²，r=2时相邻网格误差比的理论值为2² = 4。本工具从该比值反演p_obs。

通过制造解法 (Method of Manufactured Solutions, MMS) 进行代码验证

🙋

这个"制造解法"名字很有趣。通常是去"求解"问题得到解，为什么要"制造"解呢？

🎓

问得很好。思路是反过来的。常规验证采用"与基准问题解析解对比"的方式，但具有已知解析解的问题很有限。MMS则反其道行之，先随意选定一个光滑函数u_mms(x)，比如u_mms = sin(πx/L)。将其代入控制方程 -D·u''(x) = f(x)，左边可以计算，所以我们能"逆推"为了满足这个u需要多大的源项f。这样得到f(x) = D·(π/L)²·sin(πx/L)。

🙋

原来如此，先定答案再出题。这样怎么就能发现代码的问题呢？

🎓

把这个f输入你的代码，让它求解数值解u_h。理论上u_h应该完全等于u_mms。但实际上会有离散化截断误差，数值解有偏差。我们把误差E = ‖u_h − u_mms‖₂逐步测量，同时不断细化网格。二阶中心差分理论上E ∝ h²，所以h缩小一半，误差应该缩小到1/4。本工具右上角显示的E₁/E₄是64（=2⁶），正是这个缘故——h被缩小了3次，每次缩小2，所以2²×2²×2²=64。

🙋

所以说"观测收敛阶数"与设计阶数一致，就能证明代码是对的？

🎓

有很强的说服力。观测阶数p_obs = log(E_k/E_{k+1})/log(r)与p_design完美吻合，意味着"离散化按设计精度实现""边界条件正确施加""源项符号无误"这些条件全部满足。在实际工程中，常见的问题是"代码改了以后能运行，但精度下降了"。如果验证只用现有的基准算例，可能看不出来；但用MMS的话，p_obs会从2明确掉到1，你就知道引入了bug。Roy和Oberkampf的论文以及ASME V&V 10/40标准为什么推荐MMS作为金标准，就是因为它这种高灵敏度。

🙋

反过来，如果p_obs与设计阶数对不上，应该怎样排查问题呢？

🎓

试试把"期望收敛阶数"改成1，看看收敛曲线的斜率怎样。观测阶数低于设计阶数的典型原因有4个。第一是边界条件实现有误——边界一点精度低会拖累全域，导致p_obs降到1阶。第二是源项符号错了——f的正负反掉了，误差就与网格无关、卡在某个常数。第三是迭代求解器没收敛够——线性方程组的残差大于截断误差，p_obs就接近0。第四是还没进入渐近域——粗网格那边误差出现平台。如果p_obs偏离设计阶数太远，首先核查边界条件和f的符号。

常见问题

常规验证采用"对比基准问题的解析解"方式，只能处理有限的问题类型。而MMS反其道行之，先自由"制造"任意光滑解析解u_mms(x,t)，将其代入控制方程反演源项f。将该源项输入代码，即可为任意复杂的方程和边界条件创建包含精确解的测试用例。这样可以比较理论收敛阶数（设计阶数p）与观测阶数p_obs，定量证明代码实现的正确性，是V&V的黄金标准。

对于格子k的代表网格尺度h_k，L2误差定义为E_k = sqrt(Σ(u_h − u_mms)² · h)。对于相邻两格子，当网格细化比为r时，从误差比得到p_obs = log(E_k / E_{k+1}) / log(r)。本工具在4阶网格（粗→中→细→最细）中，以h缩小1/2（r=2），计算3组p_obs并取平均。若p_obs与设计阶数p（如二阶中心差分为2.0）在±0.15内相符，则代码实现的离散化阶数无误。

常见原因有4种：(1)边界条件实现有误——边界处精度低会主导全域误差，使p_obs降至1阶。(2)源项符号错误——通常f的正负反了，导致误差与网格无关而保持常数。(3)迭代求解器收敛不足——残差超过截断误差时，p_obs看似接近0。(4)网格还未进入渐近区——粗网格侧的误差出现平台。若p_obs与设计阶数偏离较大，首先核查边界条件和f的符号。

不是。MMS的最大优势在于适用于非线性、非定常、多维和耦合问题。Navier-Stokes、对流扩散、结构非线性、流固耦合、辐射换热等都可应用。选定u_mms，通过符号计算（如SymPy）展开其算子L[u_mms]即得f。制造解无需具有物理意义，关键是选择能激活代码全部项的光滑函数（多项式、三角函数、指数函数组合）。本工具以线性一维扩散为例展示MMS的基本逻辑，相同步骤可直接应用于大规模CFD代码验证。

实际应用领域

商用CFD/FEM求解器验证：ANSYS Fluent、OpenFOAM、CFD++、Code_Saturne、ABAQUS、LS-DYNA等主流求解器的开发团队在发版前必做MMS验证。每次新增乱流模型、输运方程或边界条件类型，都要对应生成u_mms和f，通过自动化回归测试检查观测阶数。美国国家航空航天局的CFL3D和FUN3D，美国桑迪亚国家实验室的SIERRA等用于安全关键领域的代码，其MMS测试套件是质量保证的核心。

学术论文审查与公开评审：近年来，顶级学术期刊（JCP、CMAME、IJNMF等）在刊登新型数值方法论文时，越来越多地要求用MMS提供收敛阶数验证。若论文声称"我们提出了二阶精度格式"，就应该附上观测阶数2.00±0.05左右的曲线。对于审稿人而言，相比单纯的解析解对比，MMS的缺陷检测能力更强，被认为是更可信的验证手段。

核电、航空等安全关键领域的评估：原子炉热工水力（CFD/CFD-NS）、火箭燃烧室、返回舱气动加热等安全规制直接相关的分析，都遵循ASME V&V 10/20/40或ANS-2.27的三层框架：代码验证→计算验证→确认。代码验证阶段的核心技术就是MMS，用制造解得到的观测阶数PASS是后续Validation结果能否被信任的前提。

企业内部开发脚本的质量把关：大型企业工程部门越来越多采用MMS作为内部Python/MATLAB分析脚本投入实际工程前的质量闸门。因为MMS计算成本小（仅需数个网格），易于集成到CI/CD流水线。本工具的交互模式，正是这种"最小化MMS测试运行器"的微缩版。

常见误解与注意事项

最大的误解是，"MMS通过就意味着代码正确=结果可信"。MMS只能证明"代码验证(Code Verification)"这一件事——程序确实在正确求解所声称的数学模型。至于这个模型与实际物理现象的符合度（Validation），需要另外用实验数据验证。比如RANS湍流解析的MMS通过，不代表模型本身与实流相符——模型如果选错了，结果还是不对。MMS保证"方程求解无误"，但不保证"方程本身对"，这个区分必须始终铭记。

其次，把制造解"优化"得过于物理合理而失败的情况。制造解无需有任何物理意义，只要光滑易算就行。反而，若u_mms过于简单（比如纯多项式），可能会漏掉代码的某些分支（对流项、非线性项、特殊边界条件处理）。实际中的做法是"刻意选择让方程全项都有非零贡献的人工函数"——多项式、三角、指数组合。本工具采用sin(πx/L)，既满足边界条件u(0)=u(L)=0，又足够简洁，是最小化示例。

第三，网格未进入渐近域就输出观测阶数的陷阱。p_obs只有在"误差由主项截断误差支配"的区域才有意义。网格太粗会混入高阶项和非渐近效应，p_obs偏离设计阶数。判别方法很简单：看4阶网格的3个p_obs值是否稳定。若稳定则在渐近域，若波动则非渐近。本工具若把N（粗网格数）调到4，有时会看p_obs轻微波动。实际代码中也要验证——不仅看粗网格侧，必须至少3对以上网格确认p_obs的收敛性。

通过制造解法 (Method of Manufactured Solutions, MMS) 进行代码验证

常见问题

实际应用领域

常见误解与注意事项

使用指南

具体计算示例

工程实践注意点