Fluent流体工程仿真计算：湍流模型选择与边界条件的决策逻辑

Fluent 做到今天，界面越来越友好，Workbench 集成越来越紧密，新手也能在半天内跑出一个有结果的算例。但”能跑出结果”和”结果是对的”之间的距离，恰好是流体仿真的精髓所在。而在所有可能出问题的地方，湍流模型的选择和边界条件的指定，占了仿真误差来源的绝对大头。

湍流模型——不是越高级越好，是越适合越好

Fluent 提供的湍流模型列表很长：Spalart-Allmaras、Standard k-ε、RNG k-ε、Realizable k-ε、Standard k-ω、SST k-ω、Transition SST、RSM、LES、DES……每个模型都有自己的”舒适区”。

S-A 模型是一方程模型，专门为航空外流（翼型、机身绕流）设计。对于存在大分离的内部流动、旋转流、强曲率效应的场景，S-A 的性能不准——这不是模型的缺陷，是它设计之初就针对附体流。拿 S-A 算突扩管流动，分离泡尺寸偏小 30-40%。

k-ε 家族是工程界使用最广的二方程模型。Standard k-ε 在充分发展的管内湍流、无分离外流中表现稳健；Realizable k-ε 改进了圆孔射流和旋转流的预测（Standard k-ε 在圆孔射流中的扩散速率被高估 ~30%）；RNG k-ε 通过重整化群理论对低雷诺数效应有更好的处理。但所有 k-ε 模型对壁面边界层内的粘性底层都不直接求解——依赖壁面函数来桥接壁面和第一个网格节点之间的速度剖面。这是一个根本性的理论局限：在逆压梯度强或边界层分离的区域，对数律壁面函数的假设本身就不成立。

SST k-ω巧妙地在边界层内部切换为 k-ω 模型（对粘性底层和分离更优），在边界层外部切换为 k-ε 模型（对自由剪切流更优）——通过一个混合函数 F₁ 平滑过渡。这使 SST k-ω 成为今天工业 CFD 的”通用推荐”——在大多数流动场景下表现和”专门为这个场景调过的湍流模型”相当或接近。代价是 SST k-ω 对进口湍流边界条件更敏感：进口的湍流强度和湍流长度尺度设得不合理，可能导致 k 和 ω 在计算域前半段震荡发散。

换一个角度看湍流模型选型——不是从”哪个模型最新”入手，而是问自己：这个流动中粘性底层的求解重要吗？ 如果摩擦阻力或对流传热是核心目标（换热器、管道沿程损失），那么壁面分辨是必须的，SST k-ω 或低 Re k-ε + Enhaced Wall Treatment 是方向。如果远场掺混、射流贯穿、大尺度分离是核心目标，SST k-ω 或 Realizable k-ε 都能应付，壁面处理的影响退居次要。

壁面函数——y+ 的故事每个人都要走一遍

y+ 是第一个网格节点离开壁面的无量纲距离。y+ ≈ 1 意味着第一个节点在粘性底层内（速度线性分布）；y+ ≈ 30-300 意味着在 log-law 层内。不同湍流模型 + 壁面函数组合对 y+ 的要求不同：

k-ε + Standard Wall Functions：要求 y+ > 30 且最好 < 300。y+ < 15 时，壁面函数把线性速度剖面强行拟合为对数律，误差很大。 SST k-ω：要求 y+ ≈ 1（对壁面函数依赖度低，直接积分到壁面）。如果 y+ 偏离到 5-10，结果还可以接受；到 30 以上就开始偏离。 Enhanced Wall Treatment（k-ε 可用）：自适应处理——y+ 小的时候用粘性底层求解，y+ 大的时候用壁面函数。是”通用”的选择，但计算稳定性不如 SST k-ω 的全域直接积分。

网格画完之后在 Fluent 里跑一个初始解，后处理里看一下壁面 y+ 云图——如果关心换热、壁面摩擦的区域 y+ 没落到正确范围内，网格就必须调整。这件事应该在正式计算前做，而不是算完之后发现问题再回头改网格。

进口湍流条件——一个被严重低估的误差源

Fluent 的进口边界条件里，”湍流强度”和”湍流粘度比”（或”水力直径”）这两个参数，很多时候被直接填默认值（湍流强度 5%，湍流粘度比 10）。但真实的进口湍流条件——由上游管道、弯头、阀门、泵等产生——可能与这些默认值差一个数量级。

在 SST k-ω 里，进口湍流粘度比直接影响 ω 的进口值——ω_inlet ∝ k / (μ_t / ρ × ν)，湍流粘度比设得不对，ω 进口值就偏了，然后整个计算域的湍流衰减速率都偏了。对于管内流动，湍流充分发展段（距进口 > 10 倍水力直径）的湍流强度通常为 3-4%——用默认 5% 问题不大。对于从一个大空间突然收缩进一个小管道的情况，进口湍流强度可能只有 1-2%——用 5% 会高估近口段的混合和换热。

更隐蔽的问题出在多相流仿真里。进口的湍流条件不止影响主流，还影响离散相（气泡、液滴、颗粒）的扩散——如果你用了离散相模型（DPM）而且开了湍流扩散（Discrete Random Walk），湍流强度的设置直接影响离散相的扩散宽度。默认 5% 和真实值 2% 的差异，可能导致离散相的径向分布偏差 15-20%。

网格无关性——不是加密就行

流体仿真里的网格无关性验证比固体力学更麻烦——固体力学关心的是应力集中点的应力值，流体力学关心的量可能是全局的（压降、流量）也可能是局部的（壁面剪切应力、分离泡长度）。这两类量的网格收敛速度不同。

一个合理的方法是做”双层无关性”：第一层，加密网格直到全局量（如管道总压降）变化 < 1%——这个保证你的流量-压降对应关系是网格无关的。第二层，在关心区域（如分离泡附近、换热器翅片表面）继续加密，直到局部量（如壁面平均换热系数）变化 < 2%。两步都满足之后，才能说”这个网格对于当前分析目标是足够的”。

第一层验证很快——加密一次重跑，几分钟看总压降。第二层就耗时间了——需要做局部自适应网格或手动加密。很多工程计算只做到第一层，然后用”网格已收敛”来交付。在大多数工程验收标准下，这不是完全可以接受的——但需要明确声明”局部量可能仍对网格敏感”。