共轭传热（CHT）分析方法概述

这个认识对了一半也错了一半。关键在于“在求解时，直接满足共享边界界面上的热通量和温度的连续性”。传统的顺序分析中，是将流体侧求得的对流传热系数（h）作为固体侧的边界条件给出，而CHT则无需此假设。例如，在汽车刹车盘中，当盘表面温度超过400°C时会发生热膨胀，从而改变流体（空气）的流道形状。捕捉这种强耦合现象正是CHT的本质。

在数学上，它被视为一个联立方程组。将区域分离为流体（f）和固体（s）并进行离散化后，可以用如下的分块矩阵形式表示。

最基本的条件有两个。第一是温度的连续性：

软件内部使用“多物理场耦合”框架。例如，COMSOL Multiphysics会构建一个整合了各物理场弱形式的单一雅可比矩阵。另一方面，像Ansys Fluent这样的分离式求解器，则采用“分区耦合”方法，即分别求解流体和固体区域，同时在每个迭代步骤中更新界面条件。后一种情况下，固体侧的热传导方程通常在流体求解器内部用FVM（有限体积法）进行离散化。

不，会重复直到收敛。典型的工作流程如下：1. 迭代求解流体方程数次（例如5次）。2. 将得到的界面温度/热通量传递给固体区域。3. 完全求解固体的热传导方程。4. 将固体侧的界面条件传递给流体侧。重复步骤1〜4，直到界面热通量的残差达到例如1e-06 W/m²以下。为了提高这个“外循环”的收敛性，可以设置类似于SIMPLE算法的松弛系数（0.5〜0.8）。

这取决于软件。在Abaqus/Standard与CFD（Fluent/Star-CCM+）进行协同仿真时，即使网格不匹配，也可以通过插值来传输数据。然而，对于Ansys Fluent内部的“固体区域”功能或Siemens Star-CCM+的“共轭传热”模型，建议在界面处匹配网格（共形网格）。如果使用非匹配网格，则需要特殊的插值方法（如GGI：通用网格界面）来保证热通量守恒，计算成本会增加约5〜10%。

首先是判断“是否需要CHT”。计算无量纲数毕渥数

界面附近的分辨率至关重要。在固体侧，从产生陡峭温度梯度的界面开始，至少铺设3〜5层边界层网格（膨胀层）。其厚度的参考标准是热渗透深度

仅监视残差是不够的。必须监视的是“界面处的净热通量”。监视从流体流入固体的热通量与从固体流出到流体的热通量（计算上符号相反）的差值，并确认其处于整个计算域的能量收支误差范围内（例如0.1%以下）。此外，应绘制代表点（固体内部的热点、界面上的一点）的温度历程图，确认其是否已达到稳态。在Ansys中可使用“Surface Report”中的“Heat Transfer Rate”，在Star-CCM+中可使用“Interface Flux Balance”报告。

存在架构上的根本差异。Ansys的主流方式是在“Fluent”内部统一求解基于FVM的流体和固体。另一方面，Abaqus（FEM）与Fluent（FVM）协作的“协同仿真”则使用MPI通信并行运行各求解器，并定期交换界面数据。这在发动机支架等复杂的结构热应力分析中表现出色。Siemens Star-CCM+专注于FVM，标榜在同一FVM框架内求解流体、固体甚至电磁场的“多物理场耦合”。COMSOL Multiphysics基于FEM，采用将所有物理场作为一个联立方程组进行公式化和求解的“完全耦合”方法。

差异很大。通常，基于单一代码库在内存内进行耦合的Fluent或Star-CCM+，数据传递开销少，在大规模并行（1000核以上）下的扩展效率高。特别是Star-CCM+在其并行化技术方面有优势。协同仿真（Abaqus+Fluent）由于两个求解器间的数据交换（通常每个时间步一次）成为瓶颈，并行效率往往下降。COMSOL的完全耦合精度高，但雅可比矩阵容易变得巨大且稠密，有时需要使用直接法求解器，内存消耗和计算时间成为挑战。有基准测试显示，在1000万单元的问题中，Fluent的CHT比Star-CCM+快约15%（取决于硬件）。

开源软件是可能的。OpenFOAM的`chtMultiRegionFoam`求解器支持基于FVM的CHT分析。但是，前处理（网格创建、边界条件设置）和后处理比商业软件更费功夫。此外，也可以组合使用CalculiX（FEM结构求解器）和Elmer（多物理场FEM求解器）。商业软件中，Autodesk CFD以相对较低的价格提供CHT功能，但在湍流模型和复杂物理模型的丰富性上不如高端软件。首先用OpenFOAM尝试学术基准测试（例如，T型管中的CHT）是性价比较高的选择。

最常见的原因是流体与固体之间的“热容不匹配”和“时间尺度的差异”。固体（例如钢）的热扩散率

该警告意味着数值误差超过了允许值。主要原因有三个。1. 网格不匹配：界面处网格不匹配，插值误差大。请应用GGI等守恒性插值，或使网格匹配。2. 收敛不足：在外循环（流体←→固体迭代）充分收敛之前停止了计算。请检查界面热通量的残差历程。3. 遗漏边界条件：例如，在固体侧设置了发热（焦耳热等），但流体侧的能量方程中没有对应的项（例如辐射）。首先，应该在没有发热和外部热通量的简单稳态热传导问题中，与解析解比较，验证界面平衡是否成立。

那是“数值振荡”或“条件稳定”的问题。特别是在显式解法，或者在协同仿真中将数据交换间隔（耦合频率）设置得比时间步长更粗的情况下会发生。例如，如果流体的时间步长为1e-4秒，数据交换间隔为1e-3秒，那么流体推进10步期间，固体一直使用旧的边界条件。这会产生不稳定性。对策：1. 使数据交换间隔与流体时间步长一致（“每步交换”）。2. 使用隐式解法。3. 尝试关闭对固体和流体使用不同时间步长的“子循环”功能。在Star-CCM+中，启用“Implicit Coupling”选项可以缓解此问题。

这是并行化中的一个重大陷阱。在CHT中，构成界面的单元/面不能被分割到不同的处理器（CPU核心）上。因为计算界面条件需要两侧（流体和固体）的信息位于同一内存空间内。对策：在前处理器中设置区域分割时，确保“不跨越界面”。在Ansys Fluent中，在“Mesh/Partition/Domain”中，先将界面指定为`boundary`再进行分割。在Star-CCM+中使用“Partition Surface”来保护界面。如果使用自动分割，请务必进行此项检查。否则会导致界面热通量局部出现2倍或0倍的误差。