CFD仿真模拟的网格收敛性与湍流模型选择——从工程实践出发

CFD仿真模拟的门槛，不在软件操作，而在物理理解和数值验证。一个典型的工程项目里，前处理（几何清理+网格划分）通常占60-70%的时间，求解器设置占20%，后处理占10%。但决定CFD结果能不能用的，往往不是时间花得最多的前处理，而是几个关键物理模型的选择——湍流模型、近壁处理、对流格式、求解器算法。

湍流模型：为什么标准k-ε不是万能的

工业CFD仿真中最常用的湍流模型是标准k-ε模型。它简单、收敛性好、对高雷诺数充分发展湍流的描述在工程设计精度范围内是够用的。但它的局限也很明确：对分离流、强曲率流、旋流的描述较差，近壁区域需要用壁面函数来连接对数律区，如果网格分辨率不够（y+太大或太小），壁面函数的适用前提被破坏，近壁速度剖面的预测就会偏离。

对于搅拌釜、旋风分离器、弯管流动这类涉及强旋流和分离流的工程问题，标准k-ε的预测精度有明显局限。RNG k-ε和Realizable k-ε是对标准版本的改进，其中Realizable k-ε对分离流的预测通常比标准版本好一些，但仍未根本解决旋流问题。

更进一步的模型是雷诺应力模型（RSM）和大涡模拟（LES）。RSM放弃涡黏性假设，直接求解雷诺应力输运方程，对强各向异性湍流的描述更准确，但计算量比k-ε模型高约30-50%，收敛性也更差。LES的计算量又是RSM的5-10倍，对大多数工程设计来说不现实，通常只用于高精度研究或关键局部流动的精细分析。

在搅拌釜项目中，对比了标准k-ε、Realizable k-ε和RSM三种模型对桨叶尖端涡流强度的预测。RSM给出的涡心位置和实验PIV（粒子图像测速）测量的结果最接近，标准k-ε低估了涡强度约20%。但RSM的收敛性差了很多——同样的网格和求解器设置，k-ε模型500步收敛，RSM跑了3000步还在小幅振荡。最终项目选择了Realizable k-ε作为折中——预测精度比标准k-ε好，计算稳定性又比RSM好。

网格收敛性：不验证就交付的CFD结果是危险的

网格收敛性验证是CFD仿真中最容易被”赶进度”跳过的步骤。原理很简单：把网格逐层加密（比如粗网格150万、中网格300万、细网格600万），看关键输出量（如阻力系数、出口平均速度、涡脱落频率）是否随网格加密趋于稳定。

但在工程项目中，跑三组不同网格的CFD计算意味着3-5倍的计算量。如果设计迭代需要跑10种不同几何的方案，每组做网格收敛性验证意味着30-50倍的计算量——这在时间和算力上通常不可行。

实际操作中更实用的做法是：用标准算例（比如零压力梯度平板边界层、后向台阶分离流）先做一次网格收敛性验证，确定在这类流动中多少网格密度能达到工程精度（比如关键量变化<3%），后续的新几何方案用这个网格密度经验来指导，不做完整的三级网格验证。这个”标准算例标定法”在团队内部已经成了常规操作——代价是第一次做标定需要跑几组网格，但后续效率提升是显著的。

近壁处理：y+的选择不只是网格的事

近壁区域的处理是CFD仿真中和湍流模型紧密耦合的部分。标准壁面函数要求第一个网格节点的y+在30-300之间（对数律区），如果y+<5，壁面函数的前提不成立，需要用低雷诺数湍流模型或直接求解黏性底层。

但在复杂几何中，同一个壁面不同位置的y+可能差一个数量级——比如搅拌釜的桨叶表面，前缘的y+可能<5（因为当地速度高，壁面距离小），尾缘的y+可能>100。用统一的壁面处理策略就无法同时满足前缘和尾缘的精度要求。

处理方案有两种：一是用自动壁面处理（Automatic Wall Treatment），很多现代CFD软件（如ANSYS Fluent的增强壁面处理、OpenFOAM的omega-based wall function）在低y+和高y+之间做平滑过渡，对y+分布不均匀的复杂几何更友好。二是用很细的近壁网格把整个壁面都推到y+<1，然后用低雷诺数模型直接求解黏性底层——这个方案精度最高，但网格量也最大。

求解器稳定性：从发散到收敛的路径

CFD求解器的收敛性差，很多时候不是求解器算法的问题，而是初场和松弛因子的设置问题。冷流场的初场（速度、压力全设为零或常数）对复杂几何来说太远了，需要先用一阶迎风格式跑几百步得到一个粗糙的收敛解，再用二阶格式继续求精。

松弛因子（Relaxation Factor）控制每次迭代中变量的更新幅度。压力修正的松弛因子通常设在0.3-0.7之间，动量和湍流量的松弛因子设在0.5-0.8之间。设得太小，收敛慢；设得太大，容易发散。对于不稳定流动（如涡脱落、过渡流），用较小的松弛因子（0.2-0.4）配合双精度求解器，通常比用大松弛因子加单精度求解器更稳定。

在科研学术网首页上能看到更多关于CFD仿真参数设置和结果验证的工程实践经验。