直接耦合（单片式）共轭热传导
直接耦合（单片式）共轭热传导 ― 用1个联立方程组同时求解固体与流体的温度场
直接耦合（单片式）共轭热传导的理论基础
直接耦合与分离型的区别
先生，"直接耦合"与分离型到底有什么区别？哪个更好？
简单地说，将流体与固体的控制方程作为1个联立方程组同时求解就是直接耦合（单片式法）。分离型是固体求解器与流体求解器分别运行，在界面交换温度与热流量进行反复，这是两个方式的区别。
反复…？就是说多次往返交换吗？
对，分离型需要"固体端计算温度 → 把界面热流量传给流体 → 流体计算温度 → 把界面温度传回固体"这样的外部反复。这个过程需要反复进行到收敛为止。收敛要花很多时间，有时候还会发散。
直接耦合则是界面的往返交流本身就不需要，1次求解联立方程组就能同时确定固体·流体的温度场。界面的温度连续性与热流量保存在代数上被保证，不用担心收敛。
哇，那直接耦合是全面优越啊？
有得必有失。直接耦合需要专用求解器，固体与流体的网格必须一体化。已有的结构求解器与CFD求解器组合使用的话，分离型更灵活。不过，电子基板上的芯片冷却这样的含有薄固体层的问题中，Biot数会很大，直接耦合在精度·效率上都更占优势。
| 特性 | 直接耦合（单片式） | 分离型（分割式） |
|---|---|---|
| 方程组 | 1个联立方程组 | 固体·流体分别求解 |
| 界面处理 | 代数上自动保证 | 外部反复到收敛 |
| 稳定性 | 无条件稳定（隐式） | 条件稳定 |
| 求解器灵活性 | 需要专用求解器 | 可以组合既有求解器 |
| 网格 | 一体网格（conforming） | 可独立网格 |
| 内存使用量 | 大（矩阵巨大） | 小（分割求解） |
控制方程与单片式矩阵
从数式的角度怎么理解"1个联立方程组"呢？具体长什么样？
先从固体领域 $\Omega_s$ 与流体领域 $\Omega_f$ 各自的热方程开始。
$$ \rho_s c_{p,s} \frac{\partial T_s}{\partial t} = \nabla \cdot (k_s \nabla T_s) + Q_s $$
流体侧（对流+扩散）:
$$ \rho_f c_{p,f} \frac{\partial T_f}{\partial t} + \rho_f c_{p,f}\,\mathbf{u} \cdot \nabla T_f = \nabla \cdot (k_f \nabla T_f) + Q_f $$
直接耦合是把这2个方程用有限元法（FEM）或有限体积法（FVM）离散化后，界面节点共享的1个块矩阵系统中组立。像这样:
$$ \begin{bmatrix} \mathbf{K}_s & \mathbf{K}_{s\Gamma} & \mathbf{0} \\ \mathbf{K}_{\Gamma s} & \mathbf{K}_{\Gamma\Gamma} & \mathbf{K}_{\Gamma f} \\ \mathbf{0} & \mathbf{K}_{f\Gamma} & \mathbf{K}_f \end{bmatrix} \begin{bmatrix} \mathbf{T}_s \\ \mathbf{T}_\Gamma \\ \mathbf{T}_f \end{bmatrix} = \begin{bmatrix} \mathbf{Q}_s \\ \mathbf{Q}_\Gamma \\ \mathbf{Q}_f \end{bmatrix} $$
哎，3行3列的块矩阵！中间的 $\Gamma$ 是什么意思？
$\Gamma$ 是固体与流体的界面的意思。界面上的节点既属于固体也属于流体，温度 $\mathbf{T}_\Gamma$ 被共享。这就是"温度连续性"代数上被保证的机制。
实际上为了简洁，也常常写成2行2列：
$$ \begin{bmatrix} \mathbf{K}_s & \mathbf{K}_{sf} \\ \mathbf{K}_{fs} & \mathbf{K}_f \end{bmatrix} \begin{bmatrix} \mathbf{T}_s \\ \mathbf{T}_f \end{bmatrix} = \begin{bmatrix} \mathbf{Q}_s \\ \mathbf{Q}_f \end{bmatrix} $$
这里 $\mathbf{K}_{sf}$ 和 $\mathbf{K}_{fs}$ 是界面耦合矩阵（off-diagonal coupling block），表达界面通过的热交换。分离型是没有这个耦合项的，代替是外部反复进行数据交换。
界面一致性条件
界面上"温度与热流量连续"的条件，用数式怎么写？
界面 $\Gamma$ 上要满足的条件有2个。首先是温度连续性（Dirichlet条件）:
$$ T_s\big|_\Gamma = T_f\big|_\Gamma $$
然后是热流量的保存（Neumann条件）:
$$ k_s \frac{\partial T_s}{\partial n}\bigg|_\Gamma = k_f \frac{\partial T_f}{\partial n}\bigg|_\Gamma $$
$n$ 是界面的法向向量。直接耦合因为界面节点共享，温度连续性通过节点匹配就自动满足。热流量保存是在弱形式（变分定式化）中，界面积分在相邻单元间相消而代数上被保证。
好比说，分离型是"两个房间之间通过电话传温度"，直接耦合是"拆掉墙把成一个房间"。用电话转达信息会有"游戏中传话"的信息偏差，没有了墙也就没有偏差的余地。
稳定性的理论根据
听说直接耦合"稳定性高"，为什么呢？
分离型中"固体→流体→固体→…"的反复是Gauss-Seidel型的逐次更新。固体与流体的热物性比（特别是 $k_s/k_f$ 或 $\rho c_p$ 的比）很极端时，这个反复的谱半径会超过1而发散。
直接耦合是把全自由度同时求解，用隐式时间积分的话就无条件稳定。特别是:
- $k_s \gg k_f$（金属固体 + 空气流体）这样热导率差别大的情况
- 薄固体层Biot数 $\text{Bi} = hL/k_s$ 是 $O(1)$ 以上的情况
- 非定常分析中时间步长很大的情况
这些条件下分离型反复收敛困难，但直接耦合就没问题。
原来如此...物性值差别大的话分离型就吃不消啊。
Giles（1997）的理论结果有个界面条件数（interface condition number）：
$$ \kappa_\Gamma \approx \frac{k_s}{k_f} \cdot \frac{\Delta x_f}{\Delta x_s} $$
$\kappa_\Gamma$ 越大收敛越慢。钢（$k_s = 50$ W/(m$\cdot$K)）与空气（$k_f = 0.026$）的组合比率达到约2000倍。直接耦合就不受这个 $\kappa_\Gamma$ 的影响，这是最大的优势。
Coffee Break 知识角
单片式CHT的历史
直接耦合CHT的定式化最初是1990年由Saxena与Launder提出的。但那时的计算机性能无法应对矩阵的巨大规模，实际问题几乎用不了。进入2000年代，多前沿法与AMG前处理的反复法实用化，才总算能处理百万自由度以上的直接耦合CHT。现在COMSOL和STAR-CCM+可以通过GUI操作设置单片式CHT，从研究者到实务工程师都广泛使用。
直接耦合（单片式）共轭热传导的数值计算方法
离散化与矩阵组立
单片式矩阵实际上怎样在计算机上组立？
基本和普通的FEM/FVM一样。固体和流体领域各自离散化生成单元矩阵，然后整体组立。关键是界面节点的编号。
FEM情况下，固体单元的弱形式是:
$$ \int_{\Omega_s} k_s \nabla N_i \cdot \nabla N_j \, d\Omega + \int_{\Omega_s} \rho_s c_{p,s} N_i N_j \, d\Omega \cdot \frac{1}{\Delta t} $$
流体单元还要加对流项 $\int \rho_f c_{p,f} N_i (\mathbf{u} \cdot \nabla N_j) \, d\Omega$。
组立时的窍门是，界面单元要对固体·流体两个领域归属。界面节点通过编号让其在总体矩阵中只出现1次，温度连续性就自动被保证。实际求解器要注意不让共享节点的自由度（DOF）重复计算。
FVM的话怎么样？CFD求解器一般都用FVM吧？
好眼光。STAR-CCM+、Fluent等CFD求解器是FVM型，但CHT功能中固体领域也用同样的FVM框架处理。界面单元面上热流量保存被保证。FVM中相邻单元间的"面流量"相当于界面耦合。
具体来说界面面上的热流量:
$$ q_\Gamma = \frac{k_s k_f}{k_s \delta_f + k_f \delta_s} (T_{s,P} - T_{f,P}) $$
用这样的调和平均计算。$\delta_s$、$\delta_f$ 是从界面到固体·流体单元中心的距离。把这个分配到固体矩阵与流体矩阵对应的行就完成了单片式矩阵。
求解器策略
这么大的矩阵怎么求解？直接法的话内存爆炸…
完全同意。单片式矩阵包含固体+流体的全自由度非常大。求解法按问题规模区分使用：
| 求解法 | 适用规模 | 特征 |
|---|---|---|
| 直接法（LU分解） | ~50万DOF | 稳健但内存O(n²)。小规模基准问题向 |
| GMRES + ILU前处理 | ~500万DOF | 非对称矩阵对应。流体Advection项导致非对称时有效 |
| GMRES + AMG前处理 | 1000万DOF~ | 大规模CHT的标准选择。STAR-CCM+的默认 |
| 块前处理付GMRES | 数千万DOF | 固体块与流体块分别前处理。研究级 |
最近常听AMG。直接耦合CHT也有效吗？
超有效。AMG（代数多重网格）自动层级化"粗近似"与"细近似"，同时衰减低频·高频误差成分。CHT的单片式矩阵因为固体流体物性值相差巨大，条件数容易恶化，但AMG能大幅缓解。
实务上AMG的收敛判定设为残差 $10^{-6}$ 以下是标准。
网格要求
直接耦合的话网格有特殊限制吗？
最重要的是界面网格一致（conforming网格）。界面的固体端与流体端节点（FEM）或面（FVM）必须一致。有些求解器可对应非conforming界面，但精度容易下降。
实务中的留意点:
- 固体壁面边界层网格: 流体侧界面附近需充分的prism层（inflation layer）。根据$y^+$ 的网格设计不可或缺
- 固体内部: 温度梯度急的薄壁部分至少3层单元
- 网格尺寸比: 界面处固体端与流体端单元尺寸比5倍内为目安
分离型用独立网格就行，这里是直接耦合的不便之处啊…
是的。不过STAR-CCM+这样"在1个软件内一括生成固体+流体网格"的工具，conforming网格制作相对容易。反而分离型要"2个求解器间对齐网格"费劲的情况也有。
直接耦合（单片式）共轭热传导的实务应用
分析工作流程
直接耦合CHT的分析，从开始到结束做什么？
基本流程以下5个步骤：
- 形状准备: CAD数据中定义固体领域与流体领域。明确界面面
- 网格生成: 生成界面conforming的一体网格。流体侧要有prism层，固体薄壁部分最少3层单元
- 物性值·边界条件: 定义固体与流体的热物性（$k$、$\rho$、$c_p$）。外部边界设定温度固定（Dirichlet）或热流量（Neumann）
- 求解器运行: 选择单片式CHT求法。AMG前处理+GMRES。残差目标 $10^{-6}$
- 后处理·验证: 确认界面温度分布与热流量平衡。与理论解、基准（例：平板上的强制对流CHT）比较
通过Biot数进行方法选择
直接耦合与分离型，迷茫的时候怎样判断？
判断基准中最有用的是Biot数 Bi：
$$ \text{Bi} = \frac{h \cdot L}{k_s} $$
$h$ 是对流热传导系数，$L$ 是固体代表长度（肉厚），$k_s$ 是固体热导率。
- $\text{Bi} \ll 1$: 固体内温度梯度小 → 分离型就足够（收敛也快）
- $\text{Bi} \sim O(1)$: 固体表面与内部温度差不可忽视 → 直接耦合有利
- $\text{Bi} \gg 1$: 固体内温度梯度急 → 直接耦合基本必须。分离型界面收敛困难
例如，电子基板（厚1.6mm、$k_s = 0.3$ W/(m·K)）吹风机强制对流（$h = 50$）的话 $\text{Bi} = 50 \times 0.0016 / 0.3 \approx 0.27$，这已经是推荐直接耦合的范围了。
Bi 0.27 也推荐直接耦合？意外地低啊。
嗯，约 $\text{Bi} > 0.1$ 就符合教科书的"固体内温度分布不能忽视"标准。反过来，厚钢制罐（$L = 20$ mm、$k_s = 50$）自然对流（$h = 10$）的话 $\text{Bi} = 0.004$，这种情况分离型完全没问题。
应用事例
实际现场直接耦合怎样用的？
代表应用例子来看：
| 分野 | 对象 | 直接耦合有利的理由 |
|---|---|---|
| 电子机器冷却 | 基板上IC芯片+冷却风扇 | 薄基板+低k环氧树脂 → 高Bi。温度精度关系到产品寿命 |
| 燃气涡轮 | 涡轮叶片内部冷却 | 薄壁（0.5~2mm）+高温气体（1500K超）。壁面温度分布决定寿命 |
| 微通道 | CPU/GPU水冷板 | 微细流路（100~500μm宽）。固体壁与流体热交互极强 |
| EV电池 | 电池单体+液冷板 | 单体间温度均一性（±2K内）对劣化防止必需 |
| LED照明 | 封装+散热器 | 微小区域散热设计。结温的指数函数式决定寿命 |
EV电池冷却也用直接耦合？Tesla也这么做吗？
Tesla不知道，但多数EV生产商都用直接耦合CHT。液冷板的蛇形流路与电池单体（18650或4680型）的发热体同时建模，模拟单体温度均一性。1模块（96单体）的计算花8~12小时，分离型要花3~5倍时间的情况也有报告。
Coffee Break 知识角
电子基板CHT的实务TIP
电子基板CHT分析最常见的失败是"基板用各向同性计算"。FR-4基板是面内（x-y方向）与板厚方向（z方向）热导率差10倍以上的异向材料（面内：约2~3 W/(m·K)，板厚方向：0.2~0.3 W/(m·K)）。忽视这个异向性，即使用直接耦合正确组立单片式矩阵，芯片温度也会被低估10~20K。物性值输入时要特别小心。
直接耦合（单片式）共轭热传导的直接耦合（单片式）共轭热传导软件与求解器比较
主要工具的直接耦合CHT对应
直接耦合CHT能做的软件具体是哪个？
主要的比较来看：
| 软件 | CHT方法 | 离散化 | 特征 |
|---|---|---|---|
| STAR-CCM+ | 单片式（Coupled Energy） | FVM | 固体·流体统一求解器。多面体网格对应。大规模CHT强 |
| COMSOL Multiphysics | 完全单片式 | FEM | GUI操作多物理耦合设置。学术分野的事实标准。参数sweep容易 |
| Ansys Fluent | Coupled Energy Solver | FVM | 2022R2从大幅强化。1亿单元超对应。Workbench连携构造连成可 |
| OpenFOAM (chtMultiRegionFoam) | 准单片式 | FVM | OSS。内部反复最小但完全单片式不是。自定义自由度高 |
| Abaqus+CFD连携 | 分离型主流 | FEM | 热-构造耦合优。流体是外部求解器连携基本 |
COMSOL学术向？企业STAR-CCM+或Fluent多？
是。汽车·航空宇宙STAR-CCM+和Fluent主流。STAR-CCM+是"网格生成到后处理在1个软件完结"强点，直接耦合CHT设定也相对直观。Fluent历史长，UDF（用户定义函数）通过扩展性高。
COMSOL数万~数十万DOF程度的精密模型（MEMS、传感器、生物设备等）强。百万DOF超的产业规模CHT是STAR-CCM+/Fluent分。
选择指南
最后选哪个，判断的Point教我！
3个轴来考虑就好：
- 问题规模: 100万DOF以下COMSOL，超过那个STAR-CCM+/Fluent
- 多物理耦合: 电磁+热+流体的3物理以上 → COMSOL。流体+固体热2物理 → STAR-CCM+
- 成本: OSS派OpenFOAM。但直接耦合品质chtMultiRegionFoam实现有限界
Coffee Break 知识角
COMSOL vs STAR-CCM+ 直接耦合CHT比较
COMSOL年间许可（Floating版）约200万日元。STAR-CCM+是Power-on-Demand方式算计时间従量课金也可（1核时数十日圆~）。学生·研究者向COMSOL的许可圧倒的良。一方、企业的量产设计流程组込、STAR-CCM+的 Design Manager 自動化强力。
直接耦合（单片式）共轭热传导的前沿研究
研究动向
直接耦合CHT分野今后怎样进化？
注目的话题几个介绍：
- GPU加速求解器: 单片式矩阵的GMRES反复在GPU上运行。NVIDIA AmgX资源库活用AMG前处理GPU实装進、CPU的5~10倍高速化報告
- 适应网格再分割（AMR）: 界面附近自动细分。固体壁温度勾配応网格最适化。STAR-CCM+ 2024以降実装
- 机械学习代理模型: 直接耦合CHT计算结果PINN（Physics-Informed Neural Network）学習、设计参数探索1000倍以上高速化。汎化性能未解决课题
- 等几何分析（IGA）: CAD的NURBS表示的FEM的形状関数使、网格生成过程省略。界面的精度格段向上
- 数字孪生连携: 实机数据CHT模型、温度场常時更新「CHT」的概念。的寿命管理実用化始
GPU加速5~10倍？！ 行列非対称GPU実装？
鋭指摘。対称行列的CG法GPU实装、CHT的非对称行列对GMRES不规则GPU性能引出。最近的研究処理mixed-precision（FP32前处理+FP64本计算）効率上。AMG前処理的（粗行列的構築）CPU行、反復的GPU実行方式現时点的。
直接耦合（单片式）共轭热传导的故障应对
常见问题与对策
直接耦合CHT困、何？
故障整理：
| 症状 | 原因 | 对策 |
|---|---|---|
| 界面温度不連続 | 网格、界面的定义 | 界面的再定义。网格整合的確认（STAR-CCM+: Interface > In-place） |
| 残差振動收敛 | 固体流体的 $k$ 比大AMG苦戦 | AMG的回数增加（3→5）。前处理切替。緩和係数0.7下 |
| 界面的熱流量合 | 网格品質不良（界面附近的） | 界面附近的细分化。層的成长率1.2以下抑制 |
| 计算時間分離型遅 | 問題本来分離型十分（低Bi）的使用 | Bi数確認。Bi < 0.1分离型切替検討 |
| 不足 (Out of Memory) | 直接法（LU分解）解 | 反復法（GMRES+AMG）变更。設定的有効化 |
| 固体側的温度非物理的高 | 固体的熱伝導率流体的値入（単位·材料的取違） | 物性值的確认。特SI单位系的一貫性 |
「界面温度不連続」、直接耦合的起？
理論上起、実装上的罠。例STAR-CCM+固体流体的界面「In-place Interface」正定义、内部的2的独立扱。見目的、数値的「壁断熱材挟」状態。
界面横切温度的線描、不連続目视確認的基本的手順。
…线確认、覚！ 直接耦合CHT奥深。先生的說明的整理。
直接耦合「正使最正确安定」方法、使见極強力武器。COMSOL的無料簡単2D問題（平板上的强制对流CHT）試。行列的结构自分的手確认、理解一気深。
本文章的评价
感谢您的回答！
参考
了
了
想要
详情
详情
报告
错误
错误