固体-流体界面处理（共轭传热分析）
理论与物理
概述
老师，我听说在CHT分析中“界面的处理”很重要，但为什么界面会成为问题呢？
问得好。在CHT（共轭传热）分析中，需要同时求解固体的热传导和流体的对流，但在这两个区域相接的“界面”处，必须正确地传递温度和热通量。这里是决定CHT分析精度的最关键点。
原来如此。具体在什么情况下会成为问题呢？
例如，在燃气轮机叶片的冷却分析中，叶片内部（固体）用几毫米的网格就足够了，但外侧的高温气体（流体）为了解析边界层，需要几十微米的网格。这样一来，固体和流体的网格在界面处就无法匹配。这时，选择“匹配网格”还是“非匹配网格”，会极大地影响工作量和精度。
界面连续条件的数学公式化
在界面处需要遵守的数学条件具体是怎样的？
在界面 $\Gamma$ 处需要满足两个条件。首先是温度连续条件（Dirichlet条件）：
其次是热通量连续条件（Neumann条件）：
这里 $k_s$, $k_f$ 分别是固体和流体的热导率，$n$ 是界面的法线方向。同时满足这两个条件，就能保证界面处的能量守恒。
原来光是温度相等还不够啊。热量的进出也必须匹配。
没错。如果只匹配温度而热通量不匹配，就会在界面处产生能量凭空出现或消失的“人工热源”。在实际工作中，有时温度分布看起来合理，但整体热平衡却出现2~3%的偏差，原因往往就是这个。
匹配网格与非匹配网格
请再详细讲解一下匹配网格和非匹配网格的区别。
总结起来大概是这样。
| 项目 | 匹配（Matching）网格 | 非匹配（Non-matching）网格 |
|---|---|---|
| 界面的节点 | 完全一致 | 不一致 |
| 插值误差 | 零（直接传递） | 取决于插值方法 |
| 网格生成 | 需同时考虑固体和流体 | 各区域可独立生成 |
| 网格效率 | 一方会产生不必要的细化 | 各区域可优化 |
| 实际应用中的采用 | 简单形状、高精度要求时 | 复杂形状时为主流 |
实际上，像汽车发动机缸体冷却分析这样固体形状复杂的情况，分别制作固体和流体的网格，再用非匹配方式连接，是业界的常识。
原来实际工作中非匹配网格用得很多啊。但是插值的精度有点让人担心…
这个担心很合理。所以，“使用哪种插值方法”就成了非匹配网格CHT的生命线。
数值解法与实现
插值方法的比较
在非匹配网格中传递温度或热通量时，有哪些插值方法呢？
我来介绍四种代表性的方法。
1. 最近邻点插值（Nearest Neighbor）
直接使用最近节点的值。实现简单，但当网格尺寸差异大时会出现阶梯状伪影。精度约为 $O(h)$ 量级。
2. 线性插值（Linear Interpolation）
从周围节点进行线性内插。精度提高到 $O(h^2)$，但能量守恒不能严格保证。
3. RBF插值（Radial Basis Function）
使用径向基函数构建平滑的插值面。精度高，但计算成本为 $O(N^2)$，在大规模模型中可能难以扩展。
4. 守恒型插值（Conservative Interpolation）
也称为GGI（General Grid Interface）或AMI（Arbitrary Mesh Interface）。通过面积加权分配热通量，严格满足能量守恒。在CHT分析中，这是首选。
守恒型插值（GGI/AMI）的原理
守恒型插值是什么机制？我想直观地理解一下。
思路很简单。几何计算界面固体侧网格面 $A_s^i$ 与流体侧网格面 $A_f^j$ 的重叠面积 $A_{ij}$。然后创建权重系数：
使用这个权重，流体侧单元 $j$ 的界面温度计算如下：
因为是基于面积的权重，所有的热量都能毫无遗漏地保存下来。可以把它想象成“把披萨切开，根据拿到几片来决定味道的分配比例”。
因为是按面积分摊，所以没有遗漏啊！不过各求解器里的名称不同，有点容易混淆…
是啊。在OpenFOAM中是AMI，在STAR-CCM+中是In-place Interface，在Ansys Fluent中是Coupled Wall，在COMSOL中是Identity Pair。名称虽然不同，但根本思路都是“基于面积的守恒型映射”。
实践指南
界面网格设计的要点
在实际工作中创建界面网格时，请告诉我一些教科书上没有的“现场智慧”！
我列举5条现场积累的经验。
- 网格尺寸比控制在1:5以内 — 固体侧和流体侧的界面单元尺寸比过大时，一个大的固体单元会对应很多流体单元，导致插值精度下降。经验上1:5是边界线。
- 流体侧务必加入边界层网格 — 界面附近的温度梯度由流体侧主导。结合y+管理，至少设置10层以上的膨胀层。
- 界面形状尽可能简化 — 复杂曲面的重叠面积计算精度会下降。在CAD阶段去除圆角等简化操作是有效的。
- 统一界面的法线方向 — 务必确认固体侧和流体侧的法线方向是否相反。因符号错误导致热通量反向的事故相当多。
- 接触热阻要显式建模 — 实际的界面存在油脂或空气层。忽略这些会导致温差被低估，与实测结果不符。
主要求解器中的设置
请告诉我各求解器中具体的设置要点。
| 求解器 | 界面设置 | 注意事项 |
|---|---|---|
| Ansys Fluent | Coupled Wall边界条件 | 确认Wall thickness=0。务必检查Shadow face是否存在 |
| STAR-CCM+ | In-place Interface → CHT | 可设置Contact Resistance。推荐Mapped连接 |
| COMSOL | Heat Transfer模块 → Pair条件 | 使用Non-conforming mesh自动插值。推荐Identity Pair |
| OpenFOAM | chtMultiRegionFoam + AMI | 在createBaffles前确认网格质量。regionProperties设置是必须的 |
OpenFOAM的chtMultiRegionFoam配置文件好像很多，挺麻烦的…
门槛确实高，但tutorials里的heatTransfer/chtMultiRegionFoam总结得很好，建议先从那里开始。区域间的边界条件定义（compressible::turbulentTemperatureCoupledBaffleMixed）是最关键的部分，理解了那里，后面就会顺利了。
软件比较
到底哪个求解器在处理界面上最强呢？
各有利弊，根据用途选择才是最好的。
| 角度 | Ansys Fluent | STAR-CCM+ | COMSOL | OpenFOAM |
|---|---|---|---|---|
| 非匹配网格支持 | 良好 | 非常强 | 良好 | 通过AMI支持 |
| GUI操作便捷性 | 直观 | 直观 | 非常简单 | 无GUI |
| 接触热阻 | 手动设置 | 可在GUI内设置 | 通过Pair条件设置 | 通过字典文件设置 |
| 大规模并行 | 擅长 | 非常擅长 | 稍弱 | 擅长 |
无论哪个求解器，“务必确认界面处的能量平衡”这个原则是共通的。后处理时输出Flux Report，检查固体输入热量与流体接收热量的差值是否在0.1%以内。
先进技术
在界面处理的研究中，最近有什么值得关注的趋势吗？
介绍三个值得关注的趋势。
1. 重叠网格（Chimera）的CHT应用
在旋转机械中固体运动的情况下，使用重叠网格处理界面的方法正在走向实用化。比传统的滑移网格更能应对大位移、大变形。
2. 通过机器学习加速插值
利用神经网络预测界面温度场的“数据驱动型插值”研究正在推进。在形状优化循环中，有可能省去每次网格变更后重新进行GGI计算的麻烦。
3. 浸入边界法（IBM）
不显式地划分界面网格，而是在背景网格上作为惩罚项处理的方法。网格生成量急剧减少，但界面附近的高精度化尚处于研究阶段，实际应用还在发展中。
故障排除
CHT分析中界面周围的故障，哪些比较常见？
症状1: 界面处温度发生不连续跳跃
原因：界面法线方向反转，或Shadow face未定义。
对策：在网格查看器中可视化法线向量，确认两面的对应关系。
症状2: 整体热平衡偏差超过数%
原因：使用了非守恒型插值，或网格尺寸比超过1:10。
对策：切换到守恒型插值（GGI/AMI），将尺寸比调整到1:5以内。
症状3: 界面附近出现振荡性温度分布
原因：在分割型耦合的显式解法中，界面数据更新滞后一个步长。
对策：将界面的松弛系数设为0.5~0.7，或改为隐式耦合。
症状4: 计算发散
原因：固体与流体的热导率差异极端（例：铜400 W/(m·K) vs 空气0.026 W/(m·K)，相差15,000倍）。
对策：从解析解设置初始值，并逐步提高耦合的松弛系数。
界面问题从外观上很难看出来，真可怕。首先要养成确认热平衡的习惯！
这是正确的做法。即使温度云图看起来很漂亮，如果热平衡误差超过0.5%，也一定要追究原因。特别是界面的通量报告，每次都必须确认。仅此一项就能预防80%的故障。
相关主题
なった
詳しく
報告