4. Fluent流体仿真基础回顾

在流固耦合与热变形联合仿真中,流体侧模型的质量直接决定耦合界面能否稳定传递压力与温度载荷。本节不重复Fluent的完整操作手册,而是围绕“为耦合仿真准备流体模型”这一目标,梳理网格、边界条件、湍流模型、求解器与收敛判据五个关键环节中真正影响耦合结果的部分。

一、网格:耦合界面处的分辨率决定载荷传递精度

流体网格的核心矛盾是:耦合界面(通常是固体壁面)上的压力与热流必须准确,而这两者都依赖近壁区的网格尺度。对于后续要传递到固体侧的壁面压力,壁面法向第一层网格高度 y+ 是关键参数。若采用壁面函数法,y+ 宜落在 30~300;若采用增强壁面处理或低雷诺数模型,则要求 y+ ≈ 1。这个选择必须在划分网格前确定,因为它决定了边界层第一层高度和层数。

以示例为例:假设管道内径 D = 50 mm,空气流速 20 m/s,密度 1.225 kg/m³,动力粘度 1.8×10⁻⁵ Pa·s,则雷诺数 Re ≈ 6.8×10⁴,属充分湍流。若目标 y+ = 1,按平板湍流边界层估算第一层高度约为 0.02 mm;若改用壁面函数、目标 y+ = 50,则第一层高度可放宽到约 1 mm。两者网格量差异可达数倍,因此必须与湍流模型选择联动决策。

耦合仿真还需注意三点:其一,界面网格不必与固体网格一一对应,但界面附近应避免尺寸突变,否则插值误差会放大;其二,若固体侧发生热变形导致界面位移,流体网格需支持动网格或网格重构,此时界面层网格质量比内部网格更敏感;其三,建议在界面附近设置至少 5~10 层边界层,并在耦合前单独检查界面处的网格正交性与长宽比。

二、边界条件:入口、出口与壁面的物理一致性

边界条件必须与真实工况一致,否则耦合传递的载荷本身就是错的。常见设置如下:

  • 入口:速度入口适合已知流量工况,质量流量入口适合可压缩或密度变化明显的工况。湍流入口需给出水力直径与湍流强度,示例中水力直径取 50 mm,湍流强度按 I = 0.16·Re^(-1/8) 估算约为 4%。
  • 出口:压力出口最常用,表压设 0 Pa;若出口可能回流,应指定回流湍流参数,避免发散。
  • 壁面:耦合界面设为无滑移壁面,并作为数据传递面;非耦合外壁可设为绝热或对流换热边界。
  • 热边界:若做热变形耦合,壁面需给出热流密度、温度或对流换热系数,且必须与固体侧热边界在界面上保持一致。
耦合界面在Fluent中通常以壁面(wall)或耦合壁面(coupled wall)形式存在。若使用System Coupling,界面需在Fluent侧标记为参与耦合的壁面,并确认其法向指向流体域内部。

三、湍流模型:按物理特征而非习惯选择

湍流模型的选择应服务于“壁面压力与热流是否准确”。对于内流、分离不明显、以壁面剪切为主的工况,k-ε Realizable 配合增强壁面处理是稳健选择;对于存在逆压梯度、分离、旋涡的工况,k-ω SST 更合适,因为它对近壁区和自由剪切层都有较好适应性。若涉及强浮力或热分层,可考虑 RNG k-ε 或雷诺应力模型。

需要强调的是:模型选择与 y+ 要求必须匹配。选择 SST 却把第一层网格做到 y+ = 100,等于浪费了模型优势;选择壁面函数却把 y+ 压到 1,则壁面函数失效。耦合仿真中,若固体热变形改变了近壁间隙,y+ 会随之变化,因此建议在耦合迭代若干步后复查 y+ 分布。

四、求解器设置:压力-速度耦合与松弛

稳态耦合常用压力基求解器,压力-速度耦合推荐 COUPLED,它对强耦合和较大网格较稳健;若计算资源有限或工况温和,也可用 SIMPLEC。离散格式方面,压力用二阶,动量与能量用二阶迎风,湍流方程可用一阶迎风起步、收敛后转二阶。

松弛因子是耦合稳定性的关键。默认设置下,压力 0.3、动量 0.7、湍流 0.8 通常可用;若耦合迭代出现压力振荡,可先将动量和压力松弛降至 0.5 和 0.2,待稳定后再逐步恢复。对于热变形耦合,能量方程松弛建议不超过 0.9,避免温度场剧烈跳变导致固体侧热应变突变。

示例求解设置(稳态,压力基):
压力-速度耦合:COUPLED
压力离散:Second Order
动量离散:Second Order Upwind
能量离散:Second Order Upwind
湍流离散:First Order Upwind(起步)
松弛因子:压力 0.3,动量 0.7,湍流 0.8,能量 0.9

五、收敛判据:残差不是唯一标准

残差降到 1×10⁻⁴(能量 1×10⁻⁶)只是基本门槛,不能作为唯一判据。耦合仿真中更应关注三类监测量:

  1. 界面物理量:耦合界面上的平均压力、平均温度、总热流是否随迭代趋于平稳,波动幅度应小于 1%。
  2. 全局守恒:入口与出口的质量流量差应小于 0.5%,能量不平衡应小于 1%。
  3. 关键点监测:在固体侧关注点或流体出口布置监测点,观察其温度或压力是否稳定。
若残差已达标但界面压力仍在缓慢漂移,说明耦合尚未真正收敛。此时应继续迭代或收紧收敛标准,而不是直接进入固体求解。耦合仿真的收敛是“流体+固体+界面传递”三者共同稳定,而非单侧残差达标。

六、为耦合仿真准备的检查清单

在把流体模型交给System Coupling之前,建议逐项确认:界面壁面已正确标记且法向一致;y+ 与湍流模型匹配;入口湍流参数已给定;出口回流参数已设置;能量方程已开启(若做热耦合);松弛因子已按耦合稳定性调整;界面监测已布置。完成这些检查后,流体模型才具备与固体模型联合求解的基础。

流体模型准备:从网格到耦合界面的关键链路 网格分辨率 → 边界条件 → 湍流模型 → 求解器 → 收敛判据 → 界面传递 ① 网格划分 • 界面边界层 5~10 层 • y+ 与壁面处理匹配 • 界面网格避免突变 • 支持动网格/重构 ② 边界条件 • 入口:速度/质量流量 • 出口:压力出口 + 回流 • 壁面:无滑移耦合面 • 热边界与固体一致 ③ 湍流模型 • k-ε Realizable • k-ω SST(分离流) • 与 y+ 要求匹配 • 耦合后复查 y+ ④ 求解器 • COUPLED 压力-速度 • 二阶迎风离散 • 松弛因子调稳 • 能量方程开启 ⑤ 收敛判据:残差 + 界面物理量 + 全局守恒 残差:连续 1e-4,能量 1e-6(基本门槛) 界面:平均压力/温度/热流波动 < 1% 守恒:进出口质量差 < 0.5%,能量不平衡 < 1% 耦合界面数据传递(System Coupling) 流体 → 固体:壁面压力 p、壁面热流 q、温度 T 固体 → 流体:位移 d、温度 T(热变形后界面位置更新) 检查清单:界面标记正确 · y+ 匹配 · 入口湍流参数 · 出口回流 · 能量开启 · 松弛调稳 · 监测布置

把上述五个环节串起来看,流体模型准备并不是孤立的“前处理”,而是耦合仿真的第一道精度关口。网格决定载荷的空间分辨率,边界条件决定载荷的物理含义,湍流模型决定近壁预测的可信度,求解器设置决定迭代能否稳定推进,收敛判据决定何时可以安全地把载荷交给固体侧。只有这五者协同一致,后续的热变形联合仿真才有可靠的基础。