2. 传热与流体耦合的物理基础

在COMSOL中做传热与流体联合仿真，第一步不是急着拖拽接口，而是先判断“热”和“流”之间到底通过什么机制交换能量。共轭传热、自然对流、强制对流、辐射换热和流固热耦合，本质上对应五类不同的耦合路径：有的通过固体壁面导热与流体对流串联，有的由温度差驱动浮升力，有的由外部压差驱动流动，有的通过表面辐射交换能量，还有的通过固体变形改变流道与换热面积。理解这些路径，才能正确选择接口、设置边界条件并判断结果是否合理。

共轭传热：固体导热与流体对流的串联

共轭传热（Conjugate Heat Transfer, CHT）的核心是固体域与流体域在交界面上同时满足温度连续和热流密度连续。固体域求解导热方程：

ρ_s C_p,s ∂T/∂t = ∇·(k_s ∇T) + Q_s

流体域求解能量方程，若忽略粘性耗散和压力功，可写为：

ρ_f C_p,f (∂T/∂t + u·∇T) = ∇·(k_f ∇T) + Q_f

在流固交界面上，COMSOL自动施加：

T_s = T_f
n·(k_s ∇T_s) = n·(k_f ∇T_f)

这意味着不需要手动指定界面热流，软件会通过“非局部耦合”或“共轭传热”多物理场节点完成能量守恒。实际建模时，若使用“共轭传热”接口，层流或湍流流动与固体传热会自动耦合；若分别使用“层流”和“固体传热”，则必须添加“热通量”边界或“温度”边界，并确保界面成对出现。

自然对流：浮升力驱动的流动

自然对流不依赖外部风机或泵，而是由温度差引起的密度差产生浮升力。在Boussinesq近似下，动量方程中的体积力项为：

F_b = -ρ_0 β (T - T_ref) g

其中β为热膨胀系数，T_ref为参考温度，g为重力矢量。COMSOL中启用自然对流通常需要：开启重力、设置参考温度、在流体材料中定义密度随温度变化或使用Boussinesq近似。判断自然对流强弱可用瑞利数：

Ra = g β ΔT L³ / (ν α)

当Ra较小时流动为层流，Ra超过约10⁹后可能转捩为湍流。教学示例中，一个高0.2 m的竖直热板，ΔT=20 K，空气β≈1/300 K⁻¹，ν≈1.5×10⁻⁵ m²/s，α≈2.1×10⁻⁵ m²/s，可估算Ra≈1.3×10⁷，属于层流自然对流，适合用层流接口配合Boussinesq近似。

强制对流：外部压差主导的换热

强制对流由外部驱动（如入口速度、压力差或风扇）主导，能量输运以对流项u·∇T为主。其换热强度常用努塞尔数衡量：

Nu = h L / k_f

在COMSOL中，强制对流通常通过“入口”边界指定速度或质量流量，通过“出口”边界指定压力。若流速较高、Re数超过临界值，应选用湍流接口（如k-ε或SST）。强制对流与共轭传热结合时，重点检查入口温度、壁面热通量和出口温度是否满足能量平衡：

Q_in = ṁ C_p (T_out - T_in)

若仿真得到的出口温升与输入功率不匹配，往往说明边界条件或网格在热边界层处分辨率不足。

辐射换热：表面间的能量交换

辐射换热在高温或真空场景中不可忽略。COMSOL中常用“表面对表面辐射”接口，基于角系数和表面发射率计算净辐射热流：

q_rad = ε σ (T_s⁴ - T_amb⁴)

对于参与性介质，还需考虑气体吸收与散射。实际设置时，需要为每个表面指定发射率ε、环境温度T_amb，并确保辐射接口与传热接口在同一几何上耦合。若忽略辐射，高温壁面的热损失会被低估，导致温度预测偏高。

流固热耦合：变形与换热的双向作用

流固热耦合（Thermo-Fluid-Structure Interaction）在流固耦合基础上加入温度场。固体受热膨胀产生热应变：

ε_th = α (T - T_ref)

变形会改变流道尺寸和换热面积，进而影响对流换热；流体温度变化又反过来改变固体温度。COMSOL中通常使用“流固耦合”多物理场节点，并同时耦合“固体传热”和“流体传热”。求解时建议先做稳态传热分析获得初始温度场，再进行瞬态或稳态流固热耦合，避免直接强耦合导致收敛困难。

耦合机制的选择与验证

实际建模时，先判断主导机制：若流速很低且温差明显，自然对流主导；若入口速度较高，强制对流主导；若温度超过数百度，辐射不可忽略；若固体变形量级接近流道尺寸的百分之几，流固热耦合必须考虑。一个实用的验证方法是做能量平衡检查：输入功率、出口焓增、壁面热损失三者应大致相等。若偏差超过5%，应检查网格在热边界层和流固界面处是否足够细，以及边界条件是否与物理场景一致。