流体数值模拟：从Reynolds数到网格收敛性的完整链条

做流体数值模拟的人，第一次面对计算不收敛的时候，大概率会把所有参数都调一遍——湍流模型换一个、松弛因子调小、网格加密……结果有时收敛了，有时还是不行，但始终不知道到底是哪一步起作用了。

这不是你操作不对，是流体数值模拟本身是一个链条——网格质量、湍流模型、边界条件、求解器设置，任何一环掉链子，整个计算就会以某种方式表现出来。

网格收敛性：不能跳过的一步

很多人拿到几何模型，画完网格直接跑计算。结果出来以后，发现和实验值差了30%，然后开始怀疑湍流模型选错了。

实际上，第一步应该怀疑的是网格——你的网格够不够密？关键区域的y+值对不对？

网格收敛性验证的做法很简单：做三套网格（粗、中、细），分别跑同一个算例，看关心的结果量（比如阻力系数Cd、出口平均速度）是否随网格加密趋于稳定。如果细网格和中网格的结果差小于3%，可以认为网格收敛了。

Fluent里有个实用技巧：用Adapt→Yplus功能，让边界层网格的y+值落在你选的湍流模型的要求范围内。k-ε模型要求y+>30（壁面函数法），k-ω SST可以在y+<1的时候直接解析粘性底层。

湍流模型选型逻辑

湍流模型的选择，本质上是在精度和成本之间取平衡。

实际工程里，k-ω SST是首选——它在附着流和分离流里都表现稳健，对y+的容忍度也高（1<y+<30都能用）。如果计算资源充足，可以上RSM做对比验证。

OpenFOAM里的湍流模型命名和Fluent略有不同，kOmegaSST对应的是Menter’s SST，是目前OpenFOAM用户的首选。

边界条件的设置细节

入口边界：速度入口（Velocity Inlet）最常用，但如果是内流场（比如管道流动），质量流量入口（Mass Flow Inlet）更准确——它不受入口截面网格密度的影响。

出口边界：压力出口（Pressure Outlet）是标准选择。有个细节要注意：如果出口截面有反向流动（backflow），需要给backflow湍流参数赋值，否则计算可能不稳定。

壁面边界：无滑移壁面（No-slip）是默认设置。如果壁面粗糙度不可忽略，需要在壁面设置里填等效砂粒粗糙度（Roughness Height, Roughness Constant）。

OpenFOAM vs Fluent：选型建议

Fluent：界面友好，前处理（mesh）和后处理（CFD-Post）集成度高，适合工程设计和教学。缺点是 licenses 贵，大规模并行（>128核）的扩展性不如OpenFOAM。

OpenFOAM：开源免费，自定义求解器灵活（可以改动量方程、加源项、写自己的湍流模型）。缺点是学习曲线陡，前处理需要第三方工具（blockMesh, snappyHexMesh, cfMesh），调试成本高。

实际选型有个简单的判据：如果是常规工程CFD（暖通、管道、外气动），Fluent更高效；如果要计算非标准物理（多相流+化学反应+自定义源项），OpenFOAM的灵活性优势明显。

常见坑点

坑一：y+值和湍流模型不匹配。用k-ε模型但y+<5，壁面函数失效，边界层速度剖面算错。解决方法：先跑一个初算，看壁面y+分布，再调整边界层网格层数。

坑二：收敛判断只看残差。残差降到1e-5不代表结果对。正确的做法是同时监控关心的结果量（比如升力系数、出口温度），看这些量是否稳定——有时残差还在降，但结果量已经收敛了。

坑三：二维模拟代替三维。Reynolds数较高的时候（Re>10000），三维效应不可忽略，二维模拟会高估湍流混合，结果偏于危险。判据：如果流动有强烈的二次流（比如弯管、涡轮叶片通道），必须用三维。

结果可信度怎么判断

第一步：和实验数据比对。如果有文献数据或者自己的实验测量结果，直接比。偏差在5%以内可以接受，10%以上要警惕。

第二步：做网格收敛性验证（前面提过）。

第三步：换一套湍流模型做对比。如果k-ω SST和RSM给出的结果差20%以上，说明这个流动问题对湍流模型敏感，需要更精细的模型或者DNS/LES。

流体数值模拟的价值，不在于算出和实验完全一致的结果（这几乎不可能），而在于让你看清流动机制——哪里有分离？哪里有回流？哪个几何参数对性能影响最大？这些洞察，往往比一个精确到小数点后三位的系数更有用。