2. ROM误差来源分类

降阶模型(ROM)的误差很少来自单一环节,它更像一条“误差链”:从高保真解流形被采样、被投影、被插值、被积分,每一环都会引入偏差,并在后续环节被放大或掩盖。若只盯着最终相对误差曲线,往往无法判断问题出在“基不够”还是“采样点没覆盖关键区域”。因此,本节按误差产生的物理与数值位置,把ROM误差来源分为六类:投影截断、基模态缺失、采样不足、参数插值、非线性处理、数值求解。理解这六类的区别与耦合关系,是后续误差评估与定位的基础。

核心判断原则:先定位误差“发生在哪一步”,再判断它“是否可被后续步骤补偿”。投影截断与基模态缺失属于“表示能力”问题,采样不足与参数插值属于“信息获取”问题,非线性处理与数值求解属于“计算实现”问题。三类问题的排查手段和修正代价完全不同。
ROM 误差来源分类与传播链 高保真解 → 采样 → 降维投影 → 参数插值 → 非线性处理 → 数值积分 → ROM 输出 高保真解 u(μ) 全阶求解 / 实验数据 采样 / 快照集 参数点选择 降维投影 POD / 贪心 / 模态基 参数插值 基/系数随 μ 变化 非线性处理 超缩减 / DEIM 数值求解 时间积分 ① 投影截断误差 保留 r 个模态,丢弃尾部 ε ∝ Σ σᵢ² (i>r) ② 基模态缺失 快照未覆盖的物理模态 如激波、局部高频 ③ 采样不足 参数空间欠采样 快照数 < 自由度需求 ④ 参数插值误差 基/系数插值偏离真实 插值阶次、节点密度 ⑤ 非线性处理误差 超缩减采样点不足 DEIM 近似非线性项 ⑥ 数值求解误差 时间步长、刚性、舍入 降阶系统条件数 误差耦合关系:表示能力 × 信息获取 × 计算实现 • 采样不足 → 基模态缺失 → 投影截断误差被低估 • 参数插值误差在参数区间边界处最显著,常与采样密度共同作用 • 非线性处理误差随非线性强度增大而放大 • 数值求解误差在长时间积分或刚性系统中累积

一、投影截断误差:表示能力的“硬上限”

投影截断误差来自降维投影本身。以POD为例,快照矩阵的奇异值 σ₁ ≥ σ₂ ≥ … 衰减越快,说明低维子空间能捕捉的能量越多。若保留前 r 个模态,截断误差的能量上界为 ε² = Σi>r σi²。这个误差是“表示能力”的硬上限:无论后续采样多密、积分多准,都无法恢复被丢弃模态所携带的信息。

实际判断时,不应只看累计能量占比(如99.9%),还要看被丢弃模态的空间形态。若尾部模态集中在某个局部区域(如壁面附近、激波下游),即使能量占比很小,也可能对关注量(如壁面热流、局部应力)产生显著影响。此时应针对关注量单独评估截断误差,而不是用全局能量指标代替。

二、基模态缺失:快照没“看见”的物理

基模态缺失与投影截断不同:截断是“有模态但没保留”,缺失是“快照里根本没有”。常见情形包括:参数采样点未覆盖分岔区间,导致某个失稳模态从未出现在快照中;瞬态快照只取了稳态段,遗漏了启动阶段的低频模态;或者高保真求解本身分辨率不足,无法解析小尺度结构。

识别基模态缺失的实用方法是留出验证:用未参与训练的“新参数点”高保真解投影到已有基上,若投影残差远大于训练集内的截断误差,说明基的覆盖范围不足。此时增加模态数无效,必须补充采样点或更换采样策略(如贪心采样、自适应采样)。

三、采样不足:参数空间的信息盲区

采样不足有两层含义:一是快照数量不足,导致基向量估计不准;二是参数点分布不合理,导致某些区域“过采样”、某些区域“欠采样”。对于参数化ROM,后者往往更致命。例如参数 μ 在 [0,1] 内变化,若只在 μ=0.1,0.2,…,0.9 均匀采样,而系统在 μ≈0.95 附近发生模态切换,则ROM在该区域必然失效。

工程上可用参数敏感度指导采样:先以少量高保真解估计关注量对 μ 的梯度,在梯度大的区域加密采样。采样数 N 的经验下界可参考 N ≥ r·ln(r) 量级,但更可靠的判据是:当继续增加采样点后,基的投影残差不再显著下降,说明采样已趋于充分。

四、参数插值误差:基与系数的“中间态”

许多ROM在离线阶段只计算若干参数点上的基和系数,在线阶段通过插值获得任意 μ 的近似。插值误差来自两方面:一是插值阶次与节点密度不匹配,例如用线性插值处理强非线性依赖;二是插值对象选择不当,直接插值基向量可能破坏正交性,而插值系数或插值矩阵更稳定。

一个常见误区是:训练点处误差很小,就认为插值误差可控。实际上插值误差在节点之间最大,尤其在参数区间边界外推时急剧增大。评估时应专门在节点中点和区间边界设置验证点,观察误差是否出现峰值。

五、非线性处理误差:超缩减的代价

非线性ROM为降低计算量,通常采用超缩减或DEIM类方法,只在一小部分采样点上计算非线性项。这引入两类误差:一是采样点数量不足,无法准确重构非线性项的全场分布;二是采样点位置由训练快照确定,当参数变化导致非线性活跃区域迁移时,原采样点可能失效。

判断非线性处理误差是否主导,可做消融对比:分别用全阶非线性项和缩减非线性项计算同一降阶系统,若两者输出差异显著,则误差主要来自非线性处理,而非投影截断。此时应增加采样点或采用自适应采样点更新策略。

六、数值求解误差:被忽视的“最后一公里”

降阶系统虽然规模小,但可能继承原系统的刚性,或因投影导致条件数变差。时间步长过大、隐式求解迭代不收敛、浮点舍入累积,都会在长时间积分中显现。这类误差的典型特征是:短时间仿真结果正常,长时间后出现漂移或能量不守恒。

排查时建议固定降阶模型,仅改变时间步长或求解器容差,观察输出变化。若误差随步长减小而显著下降,说明数值求解误差主导;若几乎不变,则应回到投影截断或采样环节查找原因。

六类误差并非孤立:采样不足会同时导致基模态缺失和投影截断误差被低估;参数插值误差在边界处常与采样密度共同作用;非线性处理误差随非线性强度增大而放大。实际评估时,建议按“表示能力 → 信息获取 → 计算实现”的顺序逐层排查,避免在错误的环节反复调参。

误差来源速查表

误差来源 典型表现 快速判断方法 优先修正手段
投影截断 训练点内误差随 r 增大而下降后趋平 绘制误差–模态数曲线 增加保留模态数或加权模态
基模态缺失 新参数点投影残差远大于训练点 留出验证 + 投影残差对比 补充采样点或改用贪心采样
采样不足 误差在参数空间分布不均 参数扫描 + 误差热图 按敏感度加密采样
参数插值 节点中点或边界处误差峰值 节点间验证点测试 提高插值阶次或改插系数
非线性处理 非线性项重构偏差大 全阶/缩减非线性项消融对比 增加采样点或自适应更新
数值求解 长时间积分漂移、能量不守恒 步长/容差敏感性测试 减小步长、改善条件数