收敛判定 — CAE用语解说

是的。残差是「当前解满足支配方程的程度」的指标，当它降低到设定的阈值（tolerance）以下时，就判定为「已收敛」并停止计算。在OpenFOAM中，通过fvSolution的solvers块中的设置如 p: tolerance 1e-6 来配置。压力残差降到1e-6以下时停止迭代，这就是意思。但这只是「局部收敛的指标」，并不能保证整体解的质量——「残差低≠解正确」这是个陷阱。

我来介绍两个典型陷阱。①过度松弛（Over-relaxation）导致振荡，残差无法下降的解——松弛系数太大时残差会在1e-3附近一直震荡。②全局收敛和局部收敛的不一致——入口到出口的积分量（总热流量或总升力系数）已经稳定，但局部涡结构仍在变化。在实际应用中，①不能只看残差日志，②还要加上目标量（阻力系数、最大温度等）的监测，等这些也稳定了才认为收敛，这是标准的双重检查方法。

概念相同，但指标略有不同。非线性FEM的Newton-Raphson法需要同时检查「力的残差（力的不平衡）」和「位移增量（位移的更新量）」。Abaqus默认设置是节点力的残差Rmax < 0.005 * F_max（荷载的0.5%以下）且位移增量cmax < 0.01 * u_max（位移的1%以下）时判定为收敛。有接触或材料非线性时收敛会变慢，会出现「DIVERGENCE」提示——此时标准处理是减小time increment（荷载增量）后重新尝试。

计算时间会无谓地增加，或者根本收敛不了。比如在OpenFOAM中设tolerance为1e-10，会碰到数值误差（舍入误差）的墙，残差无法再降低——计算就会永远进行下去。在实际应用中应根据问题的目标精度来设置：设计筛选阶段1e-3就足够了，最终精密评估才用1e-5～1e-6这样的设置。「每次都用最高精度」往往只是浪费时间。