基于修正方程的解差异特征可检测对流换热求解器中隐藏的稳定化变化

《International Communications in Heat and Mass Transfer》：Modified-equation-informed solution-discrepancy signatures detect hidden stabilization changes in convective heat-transfer solvers

　　

**结论**

• **解-差异特征签名仅从求解器输出即可识别隐藏的 SUPG 变化。**
• **在 1% 白噪声下，仅特征签名在匹配的换热读数中达到 95% 准确率。**
• **在通道审计中，故障以单位速率触发；物理对照返回无故障。**
• **在 OpenFOAM v2506 的三维瞬态和 RANS 运行中识别出散度格式。**
• **低通和乘法衰减极限达到 0.8；其余保持审查未决。**

**引言**

在工程 Peclet 数下的对流换热，是有限元求解器无法使用其纯 Galerkin 形式的设置之一。在无法分辨热边界层的网格上，对流-扩散能量方程的非稳定 Galerkin 离散化会产生虚假振荡，包括加热壁面旁的非物理负温度，而标准补救措施是迎风方向流稳定化（SUPG）[1]、[2]、[3]、[4]、[5]，即加入受控的、沿流线方向对齐的数值扩散，并按参数 τ 缩放。图 1 展示了加热通道上的效果：稳定化后，壁面法向温度分布保持单调且物理合理，而在同一粗网格上的 Galerkin 分布下冲至 T<0。

该稳定化的强度是一个数值选择，而非物理选择，在实际生产实践中，它很少是用户会设置或查看的内容。商业和库形式的换热求解器公开了边界条件、材料属性和网格控制，但稳定化公式及其参数是内部的，它们可能在求解器版本之间静默变更，或作为配置变更的副作用而改变 [3]、[6]。当稳定化改变时，求解器实际积分的离散能量算子随之改变，工程师从解中读取的换热参数也随之改变。使这种情况危险而不仅仅是令人不满的，在于这种改变是静默的：图 2 表明，将 SUPG 参数减半或加倍在测试范围内仅使平均 Nusselt 数偏移最多约 1.5%，峰值壁温偏移最多约 0.6%，偏差小到足以隐藏在现实观测噪声之下，同时仍引入系统误差，且这些误差保持在常规输出监控容限之下。

确认求解器按其预期数值格式计算是代码验证的范畴 [7]、[8]、[9]、[10]、[11]、[12]、[13]，其工具基于不同的访问假设，值得精确陈述。制造解方法需要将控制方程表达为可制造的形式。Richardson 和网格收敛指标类型的常规网格细化估计可直接应用于求解器输出，因此不需要离散化；它们返回的是给定输出的观测阶数和不确定度范围，而非产生它的格式标识。通过 τ-估计和误差输运方法进行截断误差估计则需要更多条件，即离散算子本身或一系列网格上的解 [14]、[15]。通过单一输出场观测的封闭换热求解器无法满足其中任何一种对稳定化的访问要求，而这恰恰是用户无法检查的格式部分。我们提出的问题是，仅从换热输出可以进行的验证还有哪些，具体而言，隐藏的稳定化变更是否可以被检测并区分于换热模拟在工况、网格和测量噪声中表现出的常规变化。

我们的起点是标准的。一致的离散化并不精确求解能量方程；它求解一个相邻的修正方程，其前导截断项由产生它们的格式逐项确定 [16]、[17]、[18]。在这种语言下，稳定化是一个有意添加的截断项，改变其强度会以特征性的方式改变修正方程。这种改变在整体换热输出中不可见，但存在于求解器场与可信参考之间的差异中，且可以在其中读取。我们使用的材料没有新的。修正方程分析是经典的 [16]、[17]、[18]；从计算解估计截断误差是已确立的验证实践，通过 τ-估计 [14] 和本地截断误差引导的代码验证 [15]；通过回归与导数库对比恢复微分算子在数据驱动发现中是标准的 [19]、[20]、[21]，并已专门应用于截断误差 [22]；在无需代码模型的情况下比较模拟输出本身是一个已被研究的问题 [23]；机器学习已被用于调整 SUPG 稳定化而非审计它 [24]。据我们所知，此前未被组装和测量的是此处测试的特定组合：以全局嫌疑-参考差异取代本地残差，归一化到系数方向使总体系数幅度消失，将该方向作为隐藏 SUPG 变化的检测器进行有监督分类并提供误差界，在操作条件变化时保持格式固定的物理负对照，以及当网格不可控时的明确弃权。我们声称的是该组装及其测量的极限，而非其中的材料。

具体贡献如下：

1. 一个由修正方程引导的差异特征签名，从嫌疑换热场和可信基准中读取隐藏 SUPG 变化，无需访问求解器实现，并附带对其拟合系数是什么及不是什么明确说明（第 2.3 节）。

2. 在已验证加热通道上留出检测极限，在测量误报水平下，针对同一噪声观测和同一分类器驱动的联合和三个单常规换热读数（第 3.3 节）。

3. 对观测网格、导数库和噪声结构的敏感性分析，定位方法最锐利边界——空间结构噪声，而非观测网格、导数库或嫌疑求解器运行网格（第 3.4 节）。

4. 对标注工况变化的选择性，网格混淆的弃权规则，以及电池冷却液筛查模型中在 1% 噪声下无法观测的微妙失调，且该漏检对覆盖十六倍范围的电导率无系统依赖（第 3.5 节选择性：物理对照审计，第 3.9 节电池冷却液筛查模型中的稳定化变更漏检）。

5. 通过 OpenFOAM v2506 和 scalarTransportFoam 的三维瞬态审计，仅从输出场以 0.95 分组准确率识别散度格式，并附带随机水平负对照，以及通过 simpleFoam 的湍流 k-ω SST 审计，在所有案例中识别格式（第 3.11 节外部有限体积包的三维瞬态案例，第 3.12 节湍流雷诺平均计算中的格式识别）。

6. 一种主动多分辨率方法，移除网格混淆以在规范问题上进行家族级区分，以及对公开求解器包中两个已记录时间积分器的审计，我们未编写该包，其中在每一个案例中识别字段背后的积分器，并从输出单独恢复其收敛速率（第 3.8 节）。

此类热模拟支持广泛的当前工程研究，包括光伏热和相变冷却、流-固耦合、太阳蒸发器和热脱盐、电子和发光二极管热管理、纳米流体对流以及磁流体动力学边界层 [25]、[26]、[27]、[28]、[29]、[30]、[31]、[32]、[33]、[34]、[35]、[36]。我们引用这些作为求解器输出在实际中被依赖的示例。它们不是本文验证主张的证据，后者仅基于下文报告的控制数值实验。

证据的范围应从开始就清楚。加热通道和电池模型的所有换热结果均为稳态、二维、层流且具规定速度场。我们未测试商业包，因为无法获得许可。第 3.11 节和第 3.12 节报告了外部有限体积包的三维瞬态标量输运审计和湍流雷诺平均审计；该湍流为模型化而非解析，且湍流案例审计动量而非温度。规范瞬态、间断 Galerkin、耦合和二维非结构案例（第 3.10 节和补充材料）支持基础机制的广度，而非换热主张的广度。

**方法部分片段**

本节规定两个换热模型、特征签名构造及其含义、基准要求、网格混淆及其主动修复、决策规则及其误差控制、观测模型和评估。

**结果**

贯穿全文，受控对比是数值格式而非物理：工况、网格和观测噪声作为对照固定，数值方面仅被赋予其通过它们后所余的内容。

**讨论**

差异特征签名比匹配的换热读数更早解析隐藏稳定化变化，在采样的失调 |Δα|=0.4 减弱和 0.5 增强处达到 95% 目标，而联合对、单标量和全场距离保持审查未决。机制直截了当：标量输出将离散化误差的方向结构坍塌为单一幅度，而归一化导数组成保留相对权重。

**结论**

求解器解-差异场的修正方程引导方向提供了对流换热中隐藏稳定化变化的基于输出的指标。在已验证加热通道上，1% 白观测噪声下，它在采样的 |Δα|=0.4 减弱和 0.5 增强处解析失调，而同一观测和同一分类器驱动的联合和单常规换热读数及全场距离保持右删失在最大声明值之上。

**生成性 AI 和 AI 辅助技术在手稿制备过程中的声明**

在本文准备期间，作者使用 Claude（Anthropic）以改善手稿的语言和清晰度。使用该工具后，作者根据需要审查和编辑了内容，并对发表文章的内容承担全部责任。

**资金**

本研究未获得公共、商业或非营利部门的任何特定拨款。

**利益冲突声明**

作者声明不存在已知的可能影响本文工作的竞争经济利益或个人关系。

Arnav Garg