4. 高保真参考解生成

在 ROM 精度验证中，参考解（Reference Solution）的可信度直接决定了后续误差评估是否有意义。如果参考解本身存在显著数值误差，那么 ROM 与参考解之间的偏差就无法区分是降阶模型引入的，还是参考解不准确造成的。因此，高保真参考解的生成不是“跑一次 FOM 就行”，而是需要系统控制求解精度、网格分辨率、时间步长和求解器设置四个关键因素。

一、FOM 求解精度：参考解的“基准线”

FOM（Full Order Model）是参考解的唯一来源。它的求解精度由离散格式、收敛容差和数值耗散共同决定。以有限体积法为例，空间离散通常采用二阶迎风格式，压力-速度耦合使用 SIMPLE 或 PISO 算法。若将动量方程的残差收敛标准从 1e-4 收紧到 1e-6，参考解在回流区或激波附近的局部误差可能下降一个数量级。

实际操作中，建议先做一次“精度敏感性扫描”：固定网格和时间步长，仅改变收敛容差，观察关键监测量（如升力系数、出口平均温度、最大速度）的变化。当容差从 1e-5 收紧到 1e-6 时，若监测量变化小于 0.1%，可认为该容差已足够。示例参数：某二维圆柱绕流，Re=100，收敛容差 1e-6，监测量为升力系数幅值，两次容差下差异为 0.03%，则参考解精度可接受。

二、网格收敛：从“能算”到“算准”

网格收敛性验证是参考解可信度的核心证据。推荐采用网格收敛指数（GCI）方法，至少使用三套系统加密的网格。以特征尺寸 h 为基准，网格加密比 r = h_coarse / h_fine 通常取 1.3～2.0。计算三套网格上的监测量 f1、f2、f3，用 Richardson 外推估计精确值 f_exact，并计算 GCI。

示例：某换热器通道，三套网格单元数分别为 12 万、38 万、120 万，对应 h 比为 1.5。监测量为压降 Δp，结果分别为 245.3 Pa、251.7 Pa、253.9 Pa。计算得观测精度阶 p ≈ 1.8，外推值 254.6 Pa，最细网格 GCI ≈ 0.28%。若 GCI 小于 1%，可认为网格收敛。若 GCI 大于 5%，必须继续加密或检查网格质量。

三、时间步长：瞬态参考解的“时间分辨率”

对于非定常问题，时间步长 Δt 必须足够小以解析最高关注频率。经验规则：每个周期至少 20～40 个时间步。若关注频率为 f_max，则 Δt ≤ 1/(20 f_max)。同时，CFL 数应控制在 1 以下（显式格式）或 5 以下（隐式格式）。

时间步长收敛性验证方法与网格类似：将 Δt 减半，观察监测量（如 Strouhal 数、涡脱落频率、峰值压力）的变化。若变化小于 1%，可认为时间步长收敛。示例：某圆柱绕流，Re=200，预期涡脱落频率 f ≈ 0.2 Hz。取 Δt = 0.05 s（每周期 100 步），再取 Δt = 0.025 s，Strouhal 数从 0.198 变为 0.199，变化 0.5%，可接受。

四、求解器设置：隐式与显式的选择

求解器设置包括时间推进格式、线性求解器容差、松弛因子和并行策略。对于刚性系统，隐式格式（如 BDF2）允许更大时间步长，但需保证线性求解收敛。线性求解器相对容差建议设为 1e-8～1e-10，避免因线性求解不充分引入“伪收敛”。

松弛因子对稳态问题影响显著。压力松弛因子通常取 0.3，动量取 0.7。若松弛因子过大，残差振荡不收敛；过小则收敛缓慢。建议先以默认值运行 100 步，观察残差曲线，再微调。并行计算时，分区数不宜过多，否则交界面的插值误差会累积。示例：16 核并行与 4 核并行的监测量差异应小于 0.5%。

五、综合判断与交付物

四个因素并非独立，而是相互耦合。例如，网格加密后，时间步长可适当放宽；收敛容差收紧后，线性求解器容差也需同步收紧。建议按以下顺序执行：先做网格收敛，确定空间分辨率；再做时间步长收敛，确定时间分辨率；然后收紧求解器容差，最后做一次全局敏感性扫描。

为便于实际操作，下表给出参考解可信度检查表，可直接用于项目评审或自检。

最后，参考解的可信度应形成书面记录，包括：网格量、时间步长、收敛容差、线性求解器设置、监测量收敛曲线、GCI 计算结果。这些记录不仅是 ROM 误差评估的依据，也是后续复现和审计的基础。只有参考解本身经过严格验证，ROM 的精度结论才具有工程意义。