热流体耦合 — CAE术语解释

主要是「流体的流动」和「固体的热传导」相互影响的现象。例如，在发动机气缸盖上，高温燃烧气体（流体）加热壁面，热量在金属（固体）内部传导，最后被冷却水（另一种流体）带走。求解这整个过程就是热流体耦合分析。

边界条件的处理方式决定性地不同。分别求解时，需要假设流体-固体界面的温度或热流量（比如「壁面温度恒定300℃」）。但实际上，界面温度分布由流体的对流热传递与固体内部热传导的平衡决定。耦合分析直接求解这个平衡。误差可能会超过几十摄氏度。

正是这样。流体区域使用Navier-Stokes方程和能量守恒方程，固体区域使用热传导方程（傅里叶定律）。耦合的关键是通过界面处的两个条件连接这些方程：
1. 温度连续条件：$$ T_f = T_s $$（界面处流体侧温度等于固体侧温度）
2. 热流量连续条件：$$ -k_f \frac{\partial T_f}{\partial n} = -k_s \frac{\partial T_s}{\partial n} $$（进入流体的热流量等于流出固体的热流量）

基本相同，但上下文会影响含义。「热流体耦合」是更宽泛的概念，可能包括流体-固体以外的现象，也可能包含辐射传热或相变。而「共轭热传递（Conjugate Heat Transfer, CHT）」更狭义，特别强调流体和固体热传递的耦合。Ansys Fluent的菜单里就标注为「Conjugate Heat Transfer」。

主要有两种方式：「单一求解」和「分割求解」。单一求解是将流体和固体的全部控制方程合并成一个大矩阵一次性求解，内存消耗大但收敛性好。而实际应用中常用的是分割求解，特别是「部分重叠型交替迭代法」很普遍。

例如，在某个时间步骤（或迭代步骤）中，步骤如下：
1. 以当前的界面温度作为边界条件，运行流体分析，获得界面处的热流量分布。
2. 以步骤1得到的热流量作为边界条件，运行固体热传导分析，获得新的界面温度分布。
3. 重复1和2，直到界面温度变化量小于允许值（例如0.01 K）。
为了防止发散，通常会使用「松弛系数」（例如0.7）来调整界面条件的更新。

通常分别建立。界面有两种处理方式：「匹配网格」和「非匹配网格」。匹配网格要求界面处节点位置和单元分割完全一致，数据转移很精确，但限制了网格划分的自由度。非匹配网格则允许不同的网格，通过补间函数传递热流量和温度。由于形状复杂，Siemens Star-CCM+和COMSOL Multiphysics等软件主要使用「非匹配网格」。

影响很大。最简单的方法是「逐次耦合」，某时间步先完全求解流体，再求解固体，最后进入下一步。但当流体和固体的时间尺度差异很大时（例如：空气响应快，金属热响应慢），这个方法效率低且容易不稳定。因此应该采用「强耦合」算法，在一个时间步内进行多次界面数据交换，确保时间积分的精度和稳定性。Ansys CFX默认就用强耦合算法。

首先要正确定义「耦合界面」。流体和固体接触的所有表面都要登记为接口，不能遗漏。容易忽略的是螺栓孔和微小缝隙。其次，各区域的物性值。固体的热导率是否考虑温度依赖性？例如不锈钢SUS304可以用线性近似 $$ k = 14.6 + 0.0127 \times T $$（T单位为℃），或从JIS和厂商数据直接输入表格。

仅看残差下降不够。对耦合分析来说，「界面热流量和温度的平衡」是最重要的收敛指标。具体来说，要监测跨界面的正净热流量不平衡（Net Heat Flux Imbalance）。在Ansys Fluent中，用「Report → Fluxes」可以检查，通常以输入热量的0.1~1%以下为目标。同时还要确认界面及固体内代表点的温度已经稳定到定常值。

不需要完全一致，但关键是能否分辨「热边界层」。流体侧壁面附近温度梯度很陡峭，需要极细的边界层网格，第一个单元到壁面的距离y+要小于1。而固体内部温度梯度通常较平缓，可以用相对粗的网格。不过，如果界面两侧网格大小差异太大（例如流体0.1mm，固体5mm），非匹配网格的数据补间误差会增大，需注意。

有几个办法。如果是定常问题，不要直接开启耦合，而是先用假设的壁面条件单独计算流体，待流场接近收敛后再打开耦合。再就是尽量简化固体区域网格。除了热流量集中的地方（如散热翅片根部）外，其他部分可以粗化。在求解器设置上，Ansys Fluent可以减少每次耦合迭代内部迭代次数（例如从10次减到5次），然后增加外层耦合迭代次数。并行计算是必须的，考虑到界面数据通信效率来规划区域分割。

哲学不同。Ansys（Fluent/Mechanical）和Siemens（Star-CCM+）属于「专用求解器耦合型」。分别为流体和结构优化的不同求解引擎协同进行耦合计算。特别是Star-CCM+，在一个GUI内同时处理两套网格和物理场，「多物理场耦合」是其卖点。而COMSOL Multiphysics是「单一求解型」，所有物理场方程都整合在一个框架内，用单一求解器求解。对小规模复杂耦合有优势，但大型乱流计算不太适合。

有的。Fluent与Mechanical（稳态热分析）的耦合通过名为「System Coupling」的专用模块进行分割求解，数据映射功能丰富。而CFX原本就内置「流动和热传导一体求解」功能。用TurboGrid或Meshing生成网格后，只需在其中定义流体和固体区域，就能相对简便地建立耦合分析。但它不能直接耦合复杂的固体应力分析。

用OpenFOAM等开源软件可以做。例如，标准提供「chtMultiRegionFoam」求解器，支持多区域热流体耦合。但没有GUI，需要直接编辑配置文件（如`constant/regionProperties`），学习曲线陡峭。商业软件中，Autodesk CFD也有热流体耦合功能，但在复杂乱流模型和高级耦合控制方面不如高端软件全面。

这两者是「基于部件的热流体分析」工具。FloTHERM（Mentor Graphics，现为Siemens）和Ansys Icepak内置印刷电路板、风扇、散热器等专用部件库，特别优化了电子设备冷却的工作流程。内部也是求解热流体耦合，但与通用CFD不同，不是任意形状网格，而是简化形状（如积木）的快速设置。在筐体整体温度预测的「系统热设计」阶段效率很高。

最常见原因是耦合迭代的「松弛系数」设得太大。初始值如果设为1.0（完全更新）很容易发散。从0.3～0.5这样小的值开始，待收敛后再逐步提高到0.7～0.9。另外要检查初始条件。如果固体初始温度与流体入口温度差距很大，最初的热冲击可能导致发散。最好先用单独流体分析得到的壁面温度作为固体初始温度。

怀疑界面网格分辨率不足，特别是流体侧边界层网格。热流量与 $$ q = h (T_{fluid} - T_{wall}) $$有关，热传递系数h强烈依赖于壁面附近的速度梯度。y+太大会低估h，平衡难以取得。同时检查固体侧热流量集中处（如转角）的网格。考虑「边界层适配网格」和「自适应网格细分」。

这是数据映射（补间）产生的数值误差。对策有三种。1. 改变映射方法：从默认的「投影法」改为「搜索球法」等，增加补间点数（5点改9点）。2. 调整网格尺度：消除极端的尺寸差异，特别是较粗一侧网格适度细化。3. 平滑处理结果：COMSOL中有「数据转移平滑选项」，Ansys System Coupling可以切换「Conservative」和「Non-Conservative」映射模式。要确保选择不损害热流量守恒的设置。

首先怀疑「热损失遗漏」。实机上，固体外表面一直存在「自然对流」和「辐射」放热。分析中如果这些边界条件设为「绝热」，热无处逃逸，内部温度应该升高，反之说明——外部流体（周围空气）区域没有建模，或固体表面没有设置适当的热传系数（自然对流目标值：3～10 W/m²K）和辐射率（铝约0.1，涂装面约0.8）。还有接触热阻（螺栓连接部等）的遗漏也是常见原因。