2.3 离散格式与精度控制

2.3 离散格式与精度控制

数值引擎的最后一根柱子落在每个控制体的小小界面上：面值到底怎么插值。这个看似细节的选择直接决定两件事——你的结果有多准（精度阶数）、你的计算有多稳（发散风险）。2.1 节手算案例里"剖面被算肥"的元凶，本节正式审判；瞬态计算的时间离散与库朗数约束也在这里一并交代。

对流项离散：信息从上游来

对流输运的方向性是它的灵魂：信息沿流动方向传播，因此面上的对流值应该反映上游状态。问题是"反映多少"：

一阶迎风。面值完全取上游单元值。绝对稳定、永不振荡，代价是只有一阶精度，数值扩散随网格佩克莱数增长——上节的管道案例已经见过它把剖面抹肥的手艺。它的现代定位是"攻坚工具"：初始收敛阶段用一阶把流场拉进合理区间，随后切高阶精修。

二阶迎风。用上游方向两个单元做线性外推，精度二阶。Fluent 的默认主力，绝大多数稳态问题的正确选择。强间断附近仍可能轻微振荡，需要时用限制器兜底。

QUICK 与中心类格式。QUICK 用上风加权的二次插值，对结构化网格上的对流占优问题表现好；纯中心差分精度高但对流主导时振荡风险大，主要留给湍流直接数值模拟这类特殊场景。工程里最常见的组合还是二阶迎风加限制器。

幂律与 BODY FORCE WEIGHTED 等老格式如今出场率很低，认识即可，新项目不必启用。

图 2-3：迎风与中心类格式的剖面行为对比

时间离散：显式、隐式与库朗数

瞬态模拟把时间轴也切成步长，槽位里的瞬态项用时间导数近似。显式格式用当前时刻已知量直接外推新时刻，每步便宜，但稳定性受库朗数约束——库朗数等于速度乘时间步除以网格尺度，必须小于一附近，网格密的地方时间步被迫很小。隐式格式把新时刻未知量联立求解，时间步可以放大数十倍仍稳定，代价是每步要解联立方程。Fluent 瞬态默认隐式；显式主要见于密度基高速可压缩求解器。

库朗数在压力基瞬态里的用法相反：Fluent 让你设上限（常见取一到五，耦合瞬态可放宽到数十），求解器自适应调节实际步长。它同时是收敛稳定器——瞬态残差振荡时把库朗数上限降半，往往立竿见影。

精度控制的全程视角

离散格式只控制了一处精度，整条链路上还有四道闸门。梯度重构（Least Squares 与 Green Gauss 之别）影响二阶精度的成色；压力插值（Standard、Linear、PRESTO、Body Force Weighted）在浮升力、旋流与多孔介质场景各有讲究，PRESTO 是旋流与强浮力的常客；网格自身质量是所有格式的前提，歪网格让名义二阶退化成一阶甚至零阶；时间步长与每步内迭代数（瞬态每步二十次迭代左右保线性化收敛）决定瞬态解的时间精度。验收单语言总结：格式、梯度、压力插值、网格、时间步，五项全部在案，精度声明才有出处。

事故复盘：被一阶格式伪造的掺混

背景。某燃烧室改造评估，两股不同温度的射流在混合段掺混，工程师赶进度全程一阶迎风，温度云图显示掺混迅速、均匀、漂亮，报告初稿据此得出"掺混性能达标"。

操作。复核时按验收单流程加密网格一倍并切换二阶迎风重算，同时提取混合段中心线的温度标准差作为掺混度的定量指标。

结果。二阶解显示掺混明显更慢，混合段出口温度标准差比一阶解高出近四成——一阶格式的数值扩散凭空制造了掺混假象。用一阶解做的设计结论全部作废，所幸复核发生在评审之前。

解读。数值扩散的物理后果是"假装黏性更强"，凡是结论依赖掺混、边界层、火焰前锋这类陡峭结构的场景，一阶格式的解都不具证据效力。经典防御动作是"一阶起跑、二阶冲刺"：前期用一阶把流场拉稳，中后期切二阶继续收敛，最终报告只认二阶解。

变式。若几何复杂到网格质量只能勉强及格，二阶格式可能反复振荡，此时优先改善网格而非降格式；实在要保收敛，至少在报告中声明"关键区域采用一阶，结论为保守下界"。

格式选择速查与高频疑问

一张速查表覆盖九成日常决策：

三个高频疑问。限制器开了就一定不振荡吗：限制器压制的是格式过冲，物理上不连续（激波、界面）附近仍可能轻微波动，判读时区分数值波动与物理波动。二阶瞬态时间格式（二阶隐式）何时需要：时间精度敏感的统计量（脱落频率、瞬态峰值）建议启用；准周期演化的粗看可以一阶隐式起步。库朗数上到多少算安全：压力基瞬态耦合算法数十都常见，判据不是数字而是每步内迭代收敛——线性化没收敛的大步长等于没走。

格式与网格的联动提醒：再高阶的格式也建立在面几何信息可信的前提上，歪斜单元的高阶重构可能比低阶更糟。次序永远是网格关在前、格式关在后——这也是验收单把网格放第一栏的原因。

补一个容易被忽略的时间维度精度问题：稳态计算没有时间精度可言，但瞬态计算的时间离散误差与空间离散误差独立存在、可以互相掩盖。判读方法是分别做收敛检查——空间用网格加密、时间用步长减半，两个误差带在报告里分列。只做网格无关性就宣称"误差带百分之二"的瞬态报告，隐含了"时间误差可忽略"的未证明假设；步长减半的对照在关键统计量上做一次，成本不高，堵的是评审会上最刁钻的那个追问。瞬态验收单因此比稳态多一栏：空间 GCI 加时间步收敛性，两栏都签才完整。

本节验收点

• 最终解采用二阶格式（或已书面声明一阶的适用理由）。
• 瞬态计算记录了时间步、库朗数上限与每步迭代数。
• 压力插值方案与物理场景匹配：旋流浮升力场景用了 PRESTO。
• 掺混、分离、边界层类结论没有仅凭一阶格式解作出。