1. 课程导入与COMSOL多物理场工作流

传热与流体联合仿真，是工程仿真中最常见、也最容易“看起来简单、做起来踩坑”的一类多物理场问题。它的核心并不在于把两个物理场接口都点开，而在于理解温度场与流场之间是如何通过物性、浮力、对流换热和边界条件相互耦合的。本课程将以COMSOL Multiphysics为工具，围绕典型工程场景，逐步建立从几何、材料、物理场设置到网格、求解与后处理的完整工作流。

一、典型工程场景：为什么必须联合仿真

先看几个真实感很强的场景。第一个是电子散热：一块功率芯片贴装在散热器上，风扇驱动空气流过翅片。如果只做纯导热，你会得到偏高的芯片温度，因为忽略了空气流动带走的热量；如果只做纯流体，你又无法判断芯片结温是否超限。第二个是换热器：冷热流体被固体壁面隔开，热量从热流体经壁面传到冷流体，流道内的速度分布直接决定局部换热系数。第三个是自然对流：密闭机箱内发热元件加热空气，空气受热密度变小而上升，形成浮力驱动流，此时温度场和流场是双向强耦合的。

这些场景的共同点是：温度影响流体物性（尤其是密度），流体运动又反过来输运热量。在COMSOL中，这种耦合通常通过“非等温流动”接口组来实现，而不是分别建立两个独立物理场再手动传递数据。

二、COMSOL Desktop与模型开发器

打开COMSOL Desktop后，界面左侧是模型开发器（Model Builder），它是一棵按建模顺序组织的树。典型节点顺序为：全局定义 → 组件 → 几何 → 材料 → 物理场 → 网格 → 研究 → 结果。这个顺序不是装饰，而是求解逻辑的体现：几何决定计算域，材料决定物性，物理场决定方程，网格决定离散精度，研究决定求解策略。

右侧是设置窗口，选中任意节点后在此填写参数。下方是消息、进度与日志窗口，求解时的收敛信息、警告和错误都会在这里出现。很多初学者忽略日志，实际上“非线性求解器不收敛”的第一手线索往往就在日志里，例如残差下降停滞、时间步被反复截断等。

三、物理场接口选择：从问题特征出发

选择接口时，先判断流动状态和传热机制。如果流速较低、惯性效应可忽略，选层流；如果雷诺数较高、需要湍流模型，选湍流下的k-ε或SST接口。传热方面，若固体与流体都在计算域内，选共轭传热；若只关心流体侧且壁面温度已知，可用对流换热简化。

耦合方式上，COMSOL提供“非等温流动”多物理场节点，它会自动把流场速度传给传热方程的对流项，并把温度传给流动方程以计算浮力或变物性。对于自然对流，还需在流动接口中启用体积力并设置浮力参考温度。这里有一个容易出错的约定：重力方向默认为-y，若你的几何竖直方向是z，必须在重力节点中改为(0,0,-g)，否则浮力方向会完全错误。

四、求解器基本流程与收敛判断

对于稳态非等温流动，COMSOL默认使用分离式求解器，先解流场再解温度场，交替迭代。当浮力较强或物性变化剧烈时，分离式可能收敛缓慢，此时可切换为全耦合求解器，把速度和温度放在同一个大型方程组中求解，代价是内存占用更高。瞬态问题则需注意时间步进：初始阶段建议使用较小的初始步长，让浮力流场平稳建立。

判断收敛不能只看“求解完成”提示。应检查：残差曲线是否下降到设定容差以下；监测点温度是否趋于稳定；进出口能量是否平衡。一个实用的检查是计算能量不平衡率，即（输入热量−输出热量）/输入热量，若该值在1%以内，通常说明结果可信。

五、一个可复现的教学示例

假设一个二维矩形腔体，宽0.1 m、高0.1 m，左侧壁面温度350 K，右侧壁面温度300 K，上下壁面绝热，腔内充满空气，初始温度300 K，重力沿-y方向。这是一个经典的自然对流方腔问题。在COMSOL中，选择二维组件，建立矩形几何，材料选Air，物理场添加“层流”和“共轭传热”，多物理场添加“非等温流动”。层流接口中启用体积力，浮力参考温度设为300 K。网格使用较细的自由三角形网格，壁面附近添加边界层网格。研究类型选稳态，求解器可先用分离式，若残差下降缓慢再改全耦合。

求解完成后，观察温度云图：左侧高温壁面附近空气受热上升，右侧低温壁面附近空气下沉，形成顺时针环流。速度场最大值通常出现在腔体中部偏上区域。若把左侧壁温提高到400 K，浮力增强，流速和换热系数都会上升。这个示例虽然简单，却完整覆盖了本课程后续要反复使用的操作链条。

六、本课程的学习路径

后续章节将沿着“几何与材料 → 物理场与耦合 → 网格与求解 → 结果与验证”的主线展开，并逐步引入换热器、电子散热、混合对流等更复杂的案例。每一节都会给出可操作的参数和判断依据，而不是停留在界面介绍。你需要做的，是跟着示例亲手搭建模型，并在每次求解后问自己三个问题：物理场是否选对？耦合是否完整？结果是否满足守恒？

当你能独立回答这三个问题，传热与流体联合仿真就不再是“点按钮”，而是一种可预测、可验证的工程分析能力。