共役传热（CHT）解析手法的全体概况
共役传热（CHT）解析手法的全体概况的理论基础
共役传热的基本概念
共役传热（CHT）解析是不是就是「将流体和固体的热解析一起做」以上的特别含义吗？
这种认识是半对半错的。重要的是「在共享边界的界面处直接满足热流密度和温度的连续性进行求解」这一点。在传统的逐次解析中，从流体侧求得的热传达系数（h）是作为固体侧的边界条件给定的，但在CHT中这种假设就不再必要了。例如，在汽车制动盘中，当盘表面温度超过400°C时，会发生热膨胀，流体（空气）的流道形状会改变。捕捉这样的强耦合现象就是CHT的本质。
支配方程是分别求解流体的Navier-Stokes方程和固体的热传导方程吗？
在数学上将其作为单个联立方程组处理。将区域分成流体（f）和固体（s）离散化后，可以用以下形式的块矩阵表示。
界面处的条件具体用什么公式表达？
最基本的条件有两个。首先是温度连续性：
共役传热（CHT）解析手法的全体概况的数值计算手法
离散化和求解器的选择
流体和固体的支配方程和离散化手法都不同，怎样才能合并成一个联立方程组？
在软件内部使用「多物理场耦合」框架。例如，COMSOL Multiphysics将各物理场的弱形式集成成单个雅可比矩阵。而Ansys Fluent这样的分离型求解器则采用「分区（Partitioned）耦合」手法，分别求解流体和固体领域，同时在每个迭代步骤更新界面条件。在后一种情况下，固体侧的热传导方程通常用FVM（有限体积法）在流体求解器内离散化。
求解器的迭代控制是怎样的？流体和固体各求解一遍就结束了吗？
不是的，需要反复迭代直到收敛。典型的工作流程如下：1. 求解流体方程多次（例如5次）。2. 将得到的界面温度、热流密度传递给固体领域。3. 完全求解固体的热传导方程。4. 将固体侧的界面条件传递给流体侧。将1～4重复进行，直到界面处的热流密度残差低于例如1e-06 W/m²为止。为了提高这个「外循环」的收敛性，可以设置类似SIMPLE算法的松弛系数（0.5～0.8）。
流体和固体的网格可以独立生成吗？还是要在界面处让节点一致？
这取决于软件。在Abaqus/Standard和CFD（Fluent/Star-CCM+）的Co-simulation中，即使网格不一致，也可以通过插值来转移数据。但在Ansys Fluent的「固体领域」功能或Siemens Star-CCM+的「共役热传递」模型中，推荐在界面处使网格一致（一致网格）。使用不一致网格时，需要特殊的补间手法（GGI: General Grid Interface等）来保证热流密度的保存性，计算成本约增加5～10%。
共役传热（CHT）解析手法的全体概况的实务应用
解析工作流程和验证
实际开始CHT解析时，首先要确认什么？
首先要判断「是否需要CHT」。计算无量纲数Biot数
网格生成时有什么特别要注意的地方？
界面附近的网格分辨率是关键。在固体侧，从陡峭温度梯度出现的界面起，至少贴3～5层边界层网格（膨胀层）。其厚度的参考标准是热渗透长度
收敛判定怎么做？只看残差就行吗？
仅监视残差是不够的。必须监视的是「界面处的净热流密度」。流体流向固体的流入热流密度与固体流向流体的流出热流密度（计算上符号相反）的差，要在计算区域整体能量收支误差（例如0.1%以下）范围内。另外，绘制代表点（固体内部的热点、界面一点）的温度时间历程，确认其是否达到稳定状态。Ansys使用「Surface Report」的「Heat Transfer Rate」，Star-CCM+使用「Interface Flux Balance」报告。
共役传热（CHT）解析手法的全体概况的软件比较
各求解器的特点和适用例
Ansys、Siemens、COMSOL在CHT的处理方式上有什么区别？
架构上有根本区别。Ansys主要采用「Fluent」内FVM的流体和固体统一求解方式。另一方面，Abaqus（FEM）和Fluent（FVM）联动的「Co-simulation」通过MPI通信并行运行各求解器，定期在界面处交换数据。这对于发动机支架等复杂结构热应力解析很有优势。Siemens Star-CCM+采用统一的FVM，将流体、固体、电磁场全部用同一FVM框架求解，号称「多物理场耦合」。COMSOL Multiphysics
计算速度和大规模并行计算的适应性会因软件而异吗？
差异很大。一般来说，单一代码库内存中耦合的Fluent和Star-CCM+，数据转移开销较少，大规模并行（1000核心以上）的扩展效率较高。特别是Star-CCM+在独有的并行化技术上很强。Co-simulation（Abaqus+Fluent）中，两个求解器间的数据交换（通常每个时间步一次）会成为瓶颈，并行效率下降的趋势。COMSOL的完全耦合精度较高，但雅可比矩阵容易变得巨大、密化，需要直接法求解器，内存消耗和计算时间成为课题。根据基准测试，在1000万单元的问题中，Fluent的CHT比Star-CCM+快约15%（依赖硬件）。
免费或低价软件无法进行CHT解析吗？
开源软件可以。OpenFOAM的`chtMultiRegionFoam`求解器基于FVM，可以进行CHT解析。但前处理（网格生成、边界条件设置）和后处理比商用软件费时。另外，CalculiX（FEM结构求解器）和Elmer（多物理场FEM求解器）组合也可行。商用软件中，Autodesk CFD以相对低价提供CHT功能，但乱流模型和复杂物理模型的丰富度不如高端软件。从成本效益来看，首先用OpenFOAM做学术基准测试（例如T型管的CHT）是不错的选择。
共役传热（CHT）解析手法的全体概况的故障对应
常见错误和对策
解析完全不收敛。界面温度振荡并发散。为什么？
最常见原因是流体和固体的「热容量不匹配」和「时间尺度差异」。固体（例如钢铁）的热扩散率
出现「界面热流密度失衡」警告。这是不满足能量守恒定律吗？
这个警告表示数值误差超过了允许值。原因主要有3个。1. 网格不一致：界面处网格未对齐，补间误差较大。采用GGI等保守补间或使网格对齐。2. 收敛不足：外循环（流体←→固体迭代）在充分收敛前停止。检查界面热流密度残差历程。3. 边界条件遗漏：固体侧设置发热（焦耳热等），但流体侧能量方程缺少对应项（如辐射），等等。首先应该在没有发热和外部热流的简单定常热传导问题中，与解析解相比验证界面平衡。
非定常CHT中，减小时间步长反而导致结果变得奇怪。为什么？
这是「数值振荡」或「条件稳定性」问题。特别是在显式方法或Co-simulation中数据交换间隔（Coupling Frequency）粗于时间步长的情况下容易发生。例如，流体时间步长为1e-4秒，数据交换间隔为1e-3秒，则流体推进10步间固体仍用旧的边界条件。这会产生不稳定。对策：1. 数据交换间隔与流体时间步长同步（「每步交换」）。2. 使用隐式方法（Implicit）。3. 尝试关闭固体和流体使用不同时间步长的「子循环」。Star-CCM+中启用「Implicit Coupling」选项可以缓解这个问题。
并行计算中，界面跨越的区域分割导致计算停止或结果异常。
这是重大的并行化陷阱。在CHT中，构成界面的单元/面不能被分割到不同处理器（CPU核心）。因为界面条件的计算需要同一内存空间中的两侧（流体和固体）信息。对策：在预处理器中设置「不跨界面分割」的区域分割。Ansys Fluent中，在「Mesh/Partition/Domain」中将界面指定为`boundary`后再分割。Star-CCM+使用「Partition Surface」来保护界面。使用自动分割时必须进行此检查。这是导致界面热流密度局部变成2倍或0的误差原因。
有用
更详细
错误