软件显示“收敛”,为什么结果还是可能错?

计算运行了几个小时，残差曲线终于降下来，屏幕上出现“Converged”。李工松了一口气，把结果截图放进报告：“计算已收敛，方案满足要求。”

这句话中间，其实跳过了不少证据。

追问一步：即使图上的数值已经算出来，软件也显示收敛，我们就能相信它吗？

收敛通常说明某一层求解满足了设定的数值停止条件。它不自动说明模型符合实际，更不自动说明结果足以支持当前决策。

一、软件说“收敛”，究竟在说什么

把连续的物理问题转成计算机能够处理的方程后，求解器往往需要反复迭代，使方程的不平衡量或相邻迭代的变化量逐渐减小。达到设定条件，软件便可能报告收敛。

这里有三个限定：当前模型、当前离散方式、当前判据。

以残差为例，可以把它理解为当前近似解代入离散方程后留下的不平衡量。残差减小，是重要的数值证据；但残差不是结果相对真实世界的误差百分比。残差达到10⁻⁵，不能直接解释成“结果误差只有十万分之一”。

不同软件、不同方程的残差可能采用不同的归一化方式，数值也不能简单横向比较。工程上没有一个脱离问题、对所有软件都适用的“万能收敛线”。

还要看清提示属于哪一层：线性方程求解完成，不代表外层非线性迭代完成；某个时间步收敛，不代表整个瞬态过程已经计算完毕。达到最大迭代次数而停止，更不能直接写成收敛。

二、最容易混淆的，是四个不同的问题

这些证据各有分工，不能互相替代。

同一张粗网格上的方程可以解得很充分，但它对真实物理过程的近似仍然可能不够好。

NASA的CFD验证教程也明确区分迭代误差与离散误差，并通过网格收敛研究检查后者。

三、模型设错了，也可能收敛得很漂亮

设想一个设备散热计算。实际设备通过有限接触面积把热量传到支架，接触界面存在热阻。建模时，却把整个底面设置为恒定低温。

求解器可能很快收敛，最高温度也很低。但它回答的是“底面始终被维持在这个温度时，设备有多热”，与实际散热条件未必一致。

类似问题并不罕见：结构计算把可变形连接设为完全固定；流动计算漏掉重要泄漏通道；材料参数使用了不适合当前温度的数值；载荷单位或作用方向录入错误。它们都可能得到稳定的数值结果。

所以，看残差之前，先把问题说清楚：算的对象是什么，边界在哪里，载荷从哪里来，哪些物理过程被忽略了。

软件能够求解输入的模型，却不会仅凭“收敛”替我们确认这个模型就是实际问题。

四、残差降了，真正关心的指标未必稳了

工程师关心的通常不是残差本身，而是压降、流量、最高温度、位移、反力或其他决策指标。

例如，风道计算中残差持续降低，但进出口压差仍在缓慢漂移；某个散热模型平均温度变化很小，局部热点却还在变化。只盯一条残差曲线，容易漏掉这些现象。

因此，应同步监控与决策有关的量，并记录其在一定迭代区间内的变化。Ansys的CFD说明同样强调，除残差下降外，还需要监测预先定义的目标量。

质量、能量或力的平衡检查也很重要。不过要按具体问题列全收支项：瞬态问题可能存在积累，结构动力学需要计入惯性，热问题可能包含体热源。不能只把进出口或外力简单相减，就认定所有模型都应等于零。

对于本来就存在涡脱落、周期激励等非定常现象的问题，结果随时间波动可能是物理行为。此时要检查时间分辨率、周期或统计量，而不是一味调整参数，把曲线“压平”。

五、一个算例：三次都收敛，为什么结论还会变

假设要评估某设备的稳态最高温度，教学设定的限值为85℃。保持物理模型和边界条件一致，采用逐步细化的三套网格，并充分完成各自迭代，得到下表。

只看第一次，似乎离限值还有7℃。但随着网格细化，热点温度升高，最后甚至超过限值。

三次计算都“收敛”，只说明各自通过了相应迭代检查。它们并没有因此消除网格造成的差异。

也不能看到最后两次相差2℃，就直接宣布“误差是2℃”，更不能把细网格结果当成真值。还需要检查细化方式、关键区域分辨率与收敛趋势；在满足相应条件时，可用外推或网格收敛指数等方法估计离散不确定性。[3]

这个算例最重要的启示是：当结果接近限值时，看起来不大的数值差异，也可能改变结论。可接受误差要结合用途和裕度确定，而不能仅凭“变化百分比很小”。

六、网格检查通过，也还需要现实证据

即使网格细化后结果变化很小，迭代、守恒检查也没有明显问题，仍然可能存在模型形式或输入参数方面的偏差。

例如，简化的湍流模型能否描述当前流动？热接触参数是否适合实际装配状态？边界条件的测量误差会不会明显改变结论？这些问题不会因为增加网格数量就自动消失。

与解析解或可信基准对比，可以帮助检查数值实现和求解过程；与独立试验、实测数据对比，则为模型在指定条件下的现实预测能力提供证据。两者用途应区分。用于调参的数据，也不应直接当作独立验证数据。[4]

验证不能只写“趋势一致”。还应说明比较了哪些量、在哪些位置和工况比较、偏差多少，以及测量本身有什么不确定性。一个工况下吻合，并不意味着可以无条件外推到更高载荷或另一种流动状态。

七、看到“收敛”之后，再问六句话

• 停在哪里？是正常收敛、达到迭代上限，还是仅完成了某个子步骤？是否完成全部目标时间或载荷过程？
• 算的是什么？几何、材料、单位、载荷、约束、初始条件与实际问题是否对应？
• 关注量稳了吗？除了残差，压降、温度、位移等决策指标有没有持续漂移？
• 离散影响查了吗？网格与时间步改变后，关键结果和方案排序是否稳定？
• 物理上说得通吗？收支平衡、数量级、变化方向和极端工况是否合理？
• 凭什么用于决策？有哪些独立验证证据？误差与不确定性相对允许裕度是否可接受？

这些检查应围绕用途安排。用于方案初筛与用于接近限值的最终判断，对证据充分程度的要求并不相同。

报告中，与其只写“模型已收敛，结果可靠”，不如写清楚：“本工况已满足所述迭代判据，关键指标在末段迭代中的变化已记录；网格和时间步敏感性检查结果见附件。当前模型的验证范围为……，在该范围之外使用仍需补充证据。”尚未完成的检查，应直接注明。

计算完成，是软件运行状态；结果可信，是针对具体用途、有证据支持的判断。

“收敛”值得重视，但真正让结果能够用于工程决策的，是从输入、求解到验证的一整条证据链。

参考资料

以下资料用于核对数值收敛与模型验证概念。具体收敛定义和默认阈值以所用软件、求解器及版本为准。

[1]NASA / NPARC：CFD仿真中的不确定性与误差

[2]Ansys CFX：CFD仿真中的残差与目标量监控

[3]NASA / NPARC：空间网格收敛检查

[4]NASA / NPARC：模型验证评估