整体型（单一系统）CHT求解

根本区别在于支配方程的离散化和求解阶段。分割法分别用不同求解器交替求解流体和固体方程，在边界交换数据。而整体法从一开始就将其定义为一个连立方程组。例如，将非定常热传导和流体能量方程，界面处温度和热流连续条件作为内部边界条件直接组织在一起。这样可以从根本上消除分割法中界面数据交换滞后和小时间步长时的不稳定性。

很好的问题。变量需要统一。很多实现中，整个计算区域（流体+固体）都采用「温度」作为主要变量。流体区域求解Navier-Stokes方程和连续方程，加上温度形式的能量方程。固体区域求解热传导方程。把它们组织成一个矩阵，就会形成我说的块结构：

是的，在界面单元（跨越流体和固体的单元）的离散化中，将弱形式的边界积分项适当处理就可实现。具体来说，在界面Γ上，

这正是实现的关键。整体矩阵会变得很大且非对称。直接法不现实，所以必须采用Krylov子空间法（如GMRES）配合预条件。具体的预条件方式是，对流体块

这正是整体法的一大优势。作为一个连立方程组处理，就能对整体应用隐式积分法（如二阶向后差分公式BDF2）。这样就不会受流体CFL条件或固体Fourier数的限制，可以根据物理现象选择时间步长。例如，汽车发动机部件的热响应分析需要同时处理燃烧循环（数毫秒）和块的热渗透（数秒），用整体法就可以用1ms到10ms的时间步长稳定求解。Abaqus/Standard的共轭热传导分析也采用了类似的隐式时间积分。

不一定非要一致网格。处理非匹配网格需要在界面积分时，利用「映射」技术，将一个网格侧（如流体侧）节点的温度或热流，用另一网格侧（固体侧）的形状函数进行补间。这样

首先要检查「物性值的数量级」。特别是热导率和体积热容量。空气热导率约0.026 W/mK，钢约50 W/mK，差3个数量级。整体法中这些值直接进矩阵，导致条件数恶化，求解器可能发散。对策是进行无量纲化，或启用软件的缩放选项。用ANSYS Mechanical APDL做CHT时，`MP`命令输入物性时必须严格统一单位系统。

完全同意。至少要进行以下两点的定量检查。

必须包括界面一起做。因为整体法的精度强烈依赖于界面附近的网格分辨率。具体步骤是，先适当生成边界层网格（流体侧）和考虑固体热渗透深度的网格。然后系统地细化整体网格（如基础尺寸改为1.0、0.7、0.5倍），确认监视对象（如界面平均热传达系数或固体最高温度）的变化在例如2%以内。Altair AcuSolve的CHT教程建议界面附近要素长宽比保持在5以下。

真正实现整体求解法的目前还有限。COMSOL Multiphysics用「系数偏微分方程」接口可以将流体和固体方程写成一个系统，本质上是整体的。而ANSYS Fluent的「耦合热边界条件」只在界面耦合能量方程，运动量方程分离求解，是「部分整体」。判断方法是看求解器输出日志，流体和固体温度场是否在「一次线性求解器调用」中更新，还是作为「外层/内层迭代」分别调用。

OpenFOAM的`chtMultiRegionFoam`默认采用**分割法**。它分别调用流体求解器（`buoyantSimpleFoam`等）和固体求解器（`solidTemperature`），在界面交换温度和热流的「外层迭代」循环。但开发版本或派生求解器中存在能量方程部分采用整体求解的试验性实现。如果需要完整的整体法，可能需要在`CalculiX`等基于有限元法的开源代码中自己实现自定义单元。

单独这两个都不行，需要「与外部CFD软件联动」。Abaqus/Standard有「共轭热传导分析」功能，但只处理多固体区域间的热传导，不包括流体。要做包括流体的CHT，需要用Abaqus CFD（旧的Dassault Systems `PowerFLOW`联动）或ANSYS的`System Coupling`来双向耦合Fluent（流体）和Mechanical（固体）。这属于软件间的分割法，计算成本和数据转移开销很大。

最常见原因是「初始场不匹配」。特别是固体初始温度设为室温（20°C），而流入流体温度很高（如500°C），界面会产生极端的温度梯度。第一个时间步计算出的热流会大到不合理，求解器崩溃。对策是用流入流体温度附近的值设定固体初始温度，或最初几十步用极小时间步（如1e-6秒）缓慢建立物理场，用「斜升函数」来处理。Siemens Star-CCM+的`Field Function`可以轻松设定这种初始化方案。

这是流体和固体热阻抗（热扩散率比）极端不同时发生的数值「刚度」问题。空气（热扩散率 α ≈ 2e-5 m²/s）和铜（α ≈ 1e-4 m²/s）还可以，但空气和树脂（α ≈ 1e-7 m²/s）就明显。解决办法有两个。1) 换成更稳定的时间积分格式（如从BDF1改为BDF2）。2) 加严求解器收敛条件。特别是COMSOL，将「相对误差容限」从默认0.01降到0.001以下，往往能抑制振荡。根本上，要在热扩散率低（=热易透性低）的材料一侧细化界面附近网格。

这说明求解器已经让「局部」质量、动量、能量残差变小，但还没有达到「全局」热平衡。主要原因有两个。1) **松弛系数过强**：特别是给固体热传导方程施加了太强的欠松弛，更新量太小，需要很多迭代才能达到整体热平衡。2) **边界条件遗漏**：固体外部绝热边界不完全，或流体流出边界有热流（对流流出等）。对策是监视全局热收支（流入热量 = 固体内部发热 + 系统损失），持续计算至此收敛。ANSYS Fluent可用`Surface Monitor`监视边界积分热流。

这是映射误差导致「数值热流漏」的典型现象。不能接受。用非匹配网格时要确认以下几点。