热管的热仿真

热管热的理论基础

动作原理 — 4个热传输流程

1. 蒸发（Evaporator）: CPU等热源接触的端部，工作流体（通常是水）蒸发，吸收潜热 $h_{fg}$。
2. 蒸汽传输（Adiabatic Section）: 蒸汽由压力差驱动向冷凝部高速移动。相对于铜的热导率380 W/(m·K)，热管的实效热导率达到数千～数万 W/(m·K)。
3. 冷凝（Condenser）: 蒸汽在散热器或风扇侧冷凝，释放潜热。
4. 毛细管还流（Capillary Return）: 冷凝的液体通过吸芯（烧结金属粉末、网孔、沟槽等）的毛细管力返回蒸发部。

例如笔记本电脑的情况，从CPU（蒸发部，约95℃）到铰链附近的散热器（冷凝部，约55℃），仅凭约40℃的温度差就能传输15～30W的热量。

等效热阻网络

整体的热阻 $R_\text{hp}$ 按以下方式分解：

各项的物理意义如下：

这里 $k_l$ 是液体的热导率，$k_s$ 是吸芯固体的热导率，$\varepsilon$ 是孔隙率。烧结铜粉吸芯（$\varepsilon \approx 0.5$、$k_s = 380$ W/(m·K)、水 $k_l = 0.65$ W/(m·K)）的情况下，$k_\text{eff} \approx 30\text{--}50$ W/(m·K)。

毛细管极限（Capillary Limit）

各项展开为：

• 毛细管压力（驱动力）: $\displaystyle \Delta P_\text{cap} = \frac{2\sigma}{r_\text{eff}}$
• 液相压力损失: $\displaystyle \Delta P_l = \frac{\mu_l L_\text{eff}}{K A_w \rho_l h_{fg}} Q$
• 蒸汽相压力损失: $\displaystyle \Delta P_v = \frac{128 \mu_v L_\text{eff}}{\pi d_v^4 \rho_v h_{fg}} Q$
• 重力压头: $\displaystyle \Delta P_g = \rho_l g L \sin\phi$

这里 $L_\text{eff} = L_a + (L_e + L_c)/2$ 是有效长度。最大热传输量 $Q_\text{max}$ 从等号成立的条件求得：

Merit数（性能指数 M）和工作流体选择

M值越大，在相同的吸芯结构和尺寸下，能传输更多热量。代表性工作流体的Merit数比较：

其他动作极限

1. 毛细管极限（Capillary）: 吸芯液体还流能力的极限。常温区最严格。
2. 沸腾极限（Boiling）: 蒸发部吸芯内产生气泡并成长，阻碍液体供应。高热流密度时发生。
3. 音速极限（Sonic）: 蒸汽流速达到音速（choke）。启动或低温时注意。
4. 飞溅极限（Entrainment）: 高速蒸汽卷起对侧液膜。用Weber数判定。
5. 粘性极限（Viscous）: 蒸汽压在极低温下极低。与液态金属热管启动有关。

电子设备用（水-铜、30〜100℃）中，通常毛细管极限和沸腾极限是设计的关键。

热管热的数值计算方法

1D热阻模型的实现

其中 $R_\text{hp}$ 从数据表获取，或用刚才的热阻分解式计算。实务要点：

• $R_\text{hp}$ 对输入热量 $Q$ 非线性（低Q时高，接近毛细管极限时急剧上升）
• 需要过渡响应时，将壁、吸芯、蒸汽的热容 $C$ 并联，变成 RC 电路
• 使用温度依赖的物性值时，需要迭代计算

等效热导率模型（各向异性 $k_\text{eff}$）

CFD二相流模型（VOF / Mixture）

1. VOF法（Volume of Fluid）: 直接追踪气液界面。能看到界面形状，但计算量大且对网格分辨率敏感。用于研究蒸发部吸芯内液膜的行为。
2. Mixture / Euler二流体模型: 用体积分率表示气液。比VOF粗糙，但计算效率高。

两种方法都广泛用Lee模型做蒸发冷凝源项：

其中 $r_i$ 是质量转移系数 [1/s]。$r_i$ 的值在不同文献中从0.1到 $10^6$ 变化极大，必须和实验结果对标调整。$r_i$ 太大会发散，太小会导致界面温度偏离 $T_\text{sat}$。

吸芯多孔质模型（Darcy流）

$K$ 是渗透率 [m²]，依赖吸芯结构：

热管热的实务应用

薄型热管（Vapor Chamber）的解析

1. 第1阶段 — 1D/2D筛选（1小时）: 热阻网络＋毛细管极限检查。筛选候选材料和尺寸。
2. 第2阶段 — 各向异性 $k_\text{eff}$ 模型（1天）: 系统级3D热解析。看基板的温度分布和其他部件的热影响。
3. 第3阶段 — 多孔质二相CFD（1～2周）: 预测吸芯内液分布、干涸位置。仅在开发新吸芯结构时。

实务中80%的工作只需第1～2阶段。第3阶段是R&D部门的活。

网格策略

• 厚度方向至少3层: 分开上壁、VC主体（$k_\text{eff}$）、下壁。壁是铜（$k = 380$），VC层是各向异性。
• 面内方向热源大小的1/5以下: 例如芯片10mm×10mm的话，面内网格2mm以下。
• 不要忘记接触热阻（TIM层）: 芯片↔VC间的硅脂（$k \approx 5$ W/(m·K)、厚度50μm）可能占全温差的20～40%。

边界条件的设置

• 蒸发部（芯片面）: 均匀热流或热点分布。例：SoC 5W / 8mm×8mm = 78 kW/m²
• 冷凝部（散热器面）: 对流 $h = 5\text{--}50$ W/(m²·K)（自然对流）或 $h = 50\text{--}200$ W/(m²·K)（风扇强制）
• 侧面、底面: 绝热或自然对流（$h = 5\text{--}10$ W/(m²·K)）

注意：电子设备要设最坏环境温度（例如45℃），不是室温25℃。在25℃通过的设计在车载（85℃环境）完全不行。

验证和有效性确认（V&V）

1. 等级1 — 解析验证（Verification）: 毛细管极限 $Q_\text{max}$ 和手计算对比。和公开文献基准（如Faghri 2016的钠热管实验）对标。
2. 等级2 — 和实验对比（Validation）: 用热电偶测蒸发部、冷凝部的温度。各向异性 $k_\text{eff}$ 模型误差 ±3℃以内、CFD模型 ±5℃以内为合格。
3. 等级3 — 网格收敛性: 3等级网格（粗、中、细）的Richardson外插。GCI（Grid Convergence Index）< 5%。

实务中"蒸发部温度"直接关系规格，首先要控制这里的误差。±3℃内就可以通过。

热管热的软件比较

商用工具比较

建模方法选择指南

热管热的先端研究

脉动热管（PHP / OHP）

优点是制造简单（弯曲管子而已）、易于薄型化。缺点是机制非线性，仿真极难。CFD解需用VOF+相变模型追踪每个液塞，计算量巨大。

最近Spring-Mass-Damper模型（用弹簧、质量、阻尼近似液塞）受关注，计算成本是3D CFD的千分之一，能以不错的精度预测振荡频率和平均热传输。

环路热管（LHP）

仿真需要混合方法，蒸发器内的二相流详细解，但蒸汽管和液管用1D流路处理。COMSOL的"Heat Pipe Module"或自定义代码实现。

AI替代模型的设计优化

1. 吸芯参数（粒径、孔隙率、厚度）× 工作条件（热量、倾角）组合 200～500个案例的CFD/热阻模型预先计算
2. 结果用神经网络（DNN）或高斯过程回归学习，构建替代模型
3. 用替代模型配合遗传算法或贝叶斯优化探索（数万次评估秒级完成）

例如"厚度0.5mm以下、$Q_\text{max} \geq 15$W、蒸发部温度 ≤ 85℃"约束下，找最优粒径和孔隙率。传统参数研究要数周，AI替代1小时解决。

热管热的故障处理

常见失败和对策

第二常见的是"接触热阻（TIM）被忽视"。芯片↔热管间的硅脂或焊料热阻可能是全温差的30～50%，却按理想接触估算。现实很严酷。