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;若出口可能回流,应指定回流湍流参数,避免发散。
- 壁面:耦合界面设为无滑移壁面,并作为数据传递面;非耦合外壁可设为绝热或对流换热边界。
- 热边界:若做热变形耦合,壁面需给出热流密度、温度或对流换热系数,且必须与固体侧热边界在界面上保持一致。
三、湍流模型:按物理特征而非习惯选择
湍流模型的选择应服务于“壁面压力与热流是否准确”。对于内流、分离不明显、以壁面剪切为主的工况,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%。
- 全局守恒:入口与出口的质量流量差应小于 0.5%,能量不平衡应小于 1%。
- 关键点监测:在固体侧关注点或流体出口布置监测点,观察其温度或压力是否稳定。
六、为耦合仿真准备的检查清单
在把流体模型交给System Coupling之前,建议逐项确认:界面壁面已正确标记且法向一致;y+ 与湍流模型匹配;入口湍流参数已给定;出口回流参数已设置;能量方程已开启(若做热耦合);松弛因子已按耦合稳定性调整;界面监测已布置。完成这些检查后,流体模型才具备与固体模型联合求解的基础。
把上述五个环节串起来看,流体模型准备并不是孤立的“前处理”,而是耦合仿真的第一道精度关口。网格决定载荷的空间分辨率,边界条件决定载荷的物理含义,湍流模型决定近壁预测的可信度,求解器设置决定迭代能否稳定推进,收敛判据决定何时可以安全地把载荷交给固体侧。只有这五者协同一致,后续的热变形联合仿真才有可靠的基础。