3. 流体流动接口与传热接口的匹配
在COMSOL中做传热与流体联合仿真，真正决定模型能否收敛、结果是否可信的，往往不是网格有多密，而是接口选得对不对、多物理场节点耦合得顺不顺。本章围绕层流、湍流、非等温流动、共轭传热这几类核心接口，讲清它们的适用条件，以及如何通过多物理场节点把流动与传热真正“接”起来。
一、先判断流动状态，再选流动接口
流体流动接口的选择，第一步是判断流态。工程上最常用的判据是雷诺数：
Re = ρ·u·L / μ
ρ：流体密度 (kg/m³)
u：特征速度 (m/s)
L：特征长度 (m)，管内流取水力直径
μ：动力黏度 (Pa·s)以水在圆管内流动为例，取 ρ=998 kg/m³、μ=1.0×10⁻³ Pa·s、管径 d=0.02 m。当平均流速 u=0.05 m/s 时，Re≈998，属于层流；当 u=1 m/s 时，Re≈19960，已进入湍流。经验上，管内流 Re<2300 用层流接口，Re>4000 用湍流接口，2300～4000 之间属于过渡区，需要谨慎处理。
层流（Laminar Flow）接口直接求解不可压缩或可压缩的Navier–Stokes方程，计算量小、物理清晰，适合微通道、小流量、高黏度流体（如聚合物熔体、润滑油）以及需要精确解析速度剖面的场景。它的代价是：一旦流态进入湍流，层流解会严重偏离真实流动，压降和换热系数都会被低估。
湍流（Turbulent Flow）接口不直接解析所有尺度涡，而是通过雷诺平均（RANS）加湍流模型封闭方程。COMSOL常用 k-ε、k-ω、SST 等模型。k-ε 适合充分发展的管内流和自由剪切流，对壁面附近精度一般；k-ω 和 SST 对近壁区、逆压梯度和分离流更友好，是换热问题中更稳妥的选择。湍流接口必须配合壁面处理，通常用壁函数，此时第一层网格的无量纲壁面距离 y⁺ 应落在 30～300；若想解析到黏性底层，则需 y⁺≈1，网格量会显著增加。
二、非等温流动：把温度场“带进”流动方程
当温度变化会反过来影响流体物性时，就不能把流动和传热分开算。非等温流动（Nonisothermal Flow）多物理场节点正是为此设计：它把流体流动接口与传热接口双向耦合，动量方程中出现温度相关的密度（浮力项），能量方程中出现速度场带来的对流项。
典型适用条件包括：自然对流（浮力驱动）、强制对流中温差较大导致黏度或密度明显变化、以及需要同时考虑黏性耗散的场合。例如空气在温差 50 K 的腔体内自然对流，密度差产生的浮力不可忽略，必须启用非等温流动并打开重力。若温差很小、物性近似恒定，则可用“单向耦合”——先算流场，再把速度场传给传热接口，计算更省。
三、共轭传热：流体与固体一起传热
共轭传热（Conjugate Heat Transfer, CHT）处理的是流体与固体同时存在、热量在两者界面连续传递的问题，如散热器翅片、换热器管壁、电子芯片封装。它的核心特征是：不需要在流固界面人为指定换热系数，而是通过界面上的温度连续和热流连续条件自动求解。
COMSOL提供两个层次的CHT接口：Conjugate Heat Transfer 用于层流，Conjugate Heat Transfer (Turbulent) 用于湍流。选择时先看流体侧流态，再看是否需要固体域导热。若流体域是湍流，必须用湍流版CHT，否则近壁换热会被严重低估。
使用CHT时，流固界面会自动生成“共轭传热”边界条件，用户只需保证两侧网格在界面处连续或使用一致性对。若固体侧导热系数远大于流体侧（如铜翅片配空气），界面热阻主要由流体侧边界层控制，此时流体侧近壁网格质量比固体侧更关键。
四、多物理场节点的耦合关系
接口选好后，真正把物理场“接”起来的是多物理场节点。常见组合如下：
- Nonisothermal Flow：耦合 Laminar/Turbulent Flow 与 Heat Transfer in Fluids，双向传递速度、压力、温度、物性。
- Conjugate Heat Transfer：耦合流体流动与 Heat Transfer in Solids and Fluids，自动处理流固界面。
- Temperature Coupling / Flow Coupling：用于单向或分步耦合，先算流场再算温度，适合弱耦合问题。
- Thermal Expansion（可选）：当固体热膨胀影响流道尺寸时加入。
耦合方向决定求解策略：双向耦合用全耦合或分离式求解器，非线性强时建议全耦合加自动牛顿迭代；单向耦合可先解流动、再解传热，节省计算量。判断依据是物性变化幅度和浮力强弱，而不是模型大小。
五、一个可复现的选型示例
设一个水冷板模型：铝制基板 80 mm×60 mm×5 mm，内部矩形流道宽 4 mm、高 3 mm，水入口流速 0.8 m/s，入口温度 300 K，芯片热源 50 W。先算水力直径 d_h=2ab/(a+b)=2×4×3/(4+3)≈3.43 mm，Re≈998×0.8×0.00343/0.001≈2738，落在过渡区。此时若直接用层流接口，压降和换热会被低估；若直接用湍流接口，过渡区模型本身不确定度较大。稳妥做法是：先用层流接口试算，观察出口温升和压降是否合理；若与经验关联式偏差超过 20%，改用 SST 湍流模型并检查 y⁺。
固体侧用 Heat Transfer in Solids，流体侧用 Turbulent Flow，多物理场添加 Conjugate Heat Transfer (Turbulent)。求解后重点看三个量：流道进出口压降、芯片表面最高温度、流固界面热流密度。若界面热流密度在流道入口段明显偏高、沿程递减，说明入口热边界层正在发展，符合物理预期；若界面热流出现负值或跳变，通常是界面网格不连续或耦合方向设错。
六、接口匹配的检查清单
| 检查项 | 判断依据 | 处理方式 |
|---|---|---|
| 流态判据 | Re 与 2300/4000 比较 | 层流/湍流接口二选一 |
| 是否含固体域 | 存在翅片、管壁、封装 | 启用共轭传热接口 |
| 物性是否随温度变 | 温差 > 20 K 或黏度敏感 | 启用非等温流动双向耦合 |
| 浮力是否显著 | Ri = Gr/Re² ≫ 1 | 打开重力，双向耦合 |
| 壁面网格 | y⁺ 落在模型要求区间 | 调整边界层网格 |
| 耦合方向 | 弱耦合可单向，强耦合必须双向 | 选择求解器与迭代策略 |
把这张清单走一遍，接口选择和多物理场耦合基本不会出方向性错误。真正需要反复调试的，通常是网格与求解器设置，而不是接口本身选错。