4. 流体流动模块核心功能
在 COMSOL 多物理场仿真中，流体流动模块（CFD Module）承担着动量与质量守恒的求解任务，是热流固耦合的“动力源”。本节围绕三个核心决策展开：选层流还是湍流、按不可压缩还是可压缩处理、如何设置流固界面与壁面条件。这三者共同决定了方程形式、边界约束和耦合方式，一旦选错，后续的温度场与结构场都会跟着失真。
一、层流与湍流模型的选择依据
层流与湍流的判断首先看雷诺数 Re = ρuL/μ。工程上通常以管道 Re≈2300、外掠平板 Re≈5×10⁵ 作为转捩参考。COMSOL 中“层流（Laminar Flow）”接口直接求解 Navier–Stokes 方程，适用于 Re 较低、流动有序的场景，例如微通道散热、小流量阀门内流、黏性流体搅拌。当 Re 明显超过转捩值、流动出现随机脉动时，应改用湍流接口。
COMSOL 提供的湍流建模路径主要有三类：RANS 类（k-ε、k-ω、SST）、大涡模拟（LES）和分离涡模拟（DES）。工程热流固耦合中，SST k-ω 因兼顾近壁低雷诺数区与远场自由剪切层，是首选。RANS 求解的是时均方程，需要额外求解湍动能 k 与比耗散率 ω，并引入壁面函数处理近壁面速度梯度。若使用壁面函数，第一层网格的 y⁺ 应落在 30～300；若采用低雷诺数解析到壁面，则要求 y⁺≈1。这一点直接决定网格划分策略，是层流与湍流切换时最容易被忽略的差异。
二、不可压缩与可压缩流动的取舍
是否可压缩由马赫数 Ma = u/c 决定。Ma < 0.3 时密度变化通常小于 5%，可按不可压缩处理，此时 COMSOL 默认使用“不可压缩流动”接口，密度为常数，压力只以梯度形式出现，求解稳定、收敛快。当 Ma > 0.3，或存在明显温度变化导致密度显著改变（如热浮力、气体加热），应启用“可压缩流动”，密度由理想气体定律 ρ = p/(RT) 或用户自定义状态方程给出。
需要注意：热流固耦合中，即使流速不高，若温差很大（例如电子器件从 25 ℃ 升到 150 ℃），空气密度变化可达 30% 以上，此时仍建议用可压缩或至少用 Boussinesq 近似。不可压缩模型无法捕捉热膨胀引起的自然对流，会低估散热能力。因此“低速”不等于“一定不可压缩”，要结合温度场判断。
三、流固界面与壁面条件设置
流固耦合的核心是界面传递：流体侧把压力与黏性应力传给固体，固体把位移与速度反馈给流体。在 COMSOL 中，流固耦合通过“流固耦合（FSI）”多物理场节点实现，流体域使用“流固耦合边界”，固体域使用“流体载荷”与“移动网格”。关键设置包括：
- 壁面条件：无滑移（u=0）是默认且最常用；滑移（Navier 滑移）用于稀薄气体或超疏水表面；移动壁面用于旋转机械或活塞。
- 入口/出口：入口可给定速度、质量流量或压力；出口通常用压力出口（p=0）并勾选“抑制回流”，防止数值发散。
- 对称与周期性：用于缩减计算域，需保证几何与流动对称。
- 网格变形：FSI 中流体网格随固体移动，需设置“变形域”并监控网格质量，避免负体积。
下面用一张图梳理从“判断流态”到“设置边界”的完整决策链，帮助你在建模时快速定位当前该选哪个接口、哪个边界条件。
四、实操示例：二维通道内热流固耦合参数设置
以一个“空气流过加热铜块”的二维模型为例，说明关键参数如何落地。几何为长 0.2 m、高 0.02 m 的通道，中间放置 0.05 m × 0.01 m 的铜块。空气入口速度 2 m/s，入口温度 293.15 K，铜块底部热流密度 5000 W/m²。
先估算：空气 ρ≈1.2 kg/m³，μ≈1.8×10⁻⁵ Pa·s，特征长度取通道高度 0.02 m，则 Re = 1.2×2×0.02/1.8×10⁻⁵ ≈ 2667，略高于 2300，处于转捩区。工程上可先用层流试算，若残差振荡或出口速度剖面出现非对称脉动，再切换 SST k-ω。Ma = 2/340 ≈ 0.006，远小于 0.3，按不可压缩处理即可；但铜块附近空气被加热到约 400 K，密度下降约 25%，若关注浮力效应，可开启“体积力”并采用 Boussinesq 近似。
边界设置：入口用“速度入口”，u = 2 m/s，方向沿 x；出口用“压力出口”，p = 0，勾选“抑制回流”；上下壁面用“无滑移壁面”；铜块与空气接触面用“流固耦合边界”。固体力学中铜块底部施加“边界热源”5000 W/m²，其余表面与流体耦合。网格方面，近壁面第一层高度取 0.1 mm，保证 y⁺ 在层流下小于 1；若切换湍流并采用壁面函数，可放宽到 0.5 mm 使 y⁺ 落在 30～300。
总结本节决策顺序：先算 Re 定层流/湍流，再算 Ma 与温差定可压缩性，最后按物理场景选壁面条件与耦合方式。三者相互关联——湍流壁面函数影响网格，可压缩性影响密度与浮力，FSI 影响网格变形。把这条决策链走顺，后续传热与结构模块的耦合才有可靠基础。