CFD仿真模拟：如何让数字风洞说出真话

做CFD仿真的人，大多在第一次拿到漂亮云图的时候兴奋过，然后被实验数据打脸——计算出来的升力系数和风洞数据差了25%，阻力系数的趋势甚至反了。

这不是CFD不行，是数字风洞和实体风洞一样，需要校准和验证。你给它的几何、网格、边界条件，每一个都在参与决定最终结果的可靠性。

几何清理：最容易被跳过的一步

CAD几何直接导入CFD软件，大概率会有问题。缝隙、重叠面、非流形边、细小特征（比如直径2mm的螺栓），这些在CAD里无所谓，在CFD网格划分里是灾难。

实际工作里，几何清理占整个CFD项目30-40%的时间，但大多数人只给它5%。后果是网格质量差，计算要么不收敛，要么收敛了解答是错的。

STAR-CCM+的几何修复工具比较强，能自动检测缝隙和重叠面。Fluent里需要用SCDM（SpaceClaim）或DesignModeler先做清理，再导入。不要跳过这步——后面所有工作都建立在几何正确的基础上。

网格策略：结构还是非结构

结构化网格（hexahedral）精度高、计算效率高，但生成难度大，复杂几何几乎做不出全结构网格。

非结构化网格（tetrahedral + prism layer）生成快，对复杂几何友好，但精度稍低，需要更多的网格单元才能达到同样的Y+分辨率。

工程实践里的折中方案：近壁面用棱柱层网格（prism layer，本质是结构化的），主流区用非结构网格。这套组合在Fluent和STAR-CCM+里都是标准配置，兼顾精度和效率。

关键参数：棱柱层的第一层厚度由目标y+决定。如果目标是y+≈1（低雷诺数模型），第一层厚度大概在10^-5量级（以特征长度为1 m计）。用Fluent的Mesh→Inflation或者STAR-CCM+的Prism Layer设置，先估算再调整。

湍流模型选型的工程判据

前面文章9已经讨论过湍流模型的技术细节，这里从工程决策角度再说一遍：

• 外部气动（汽车、建筑、桥梁）：k-ω SST，首选。
• 内部流动（管道、HVAC）：Realizable k-ε，够用且便宜。
• 强旋流、燃烧室流动：需要用RSM或者大涡模拟（LES），标准两方程模型在这里会失效。

有个简单的验证方法：先用k-ω SST跑一个基准算例，再用Realizable k-ε跑一遍，如果两个结果差5%以内，说明这个流动对湍流模型不敏感，后续可以用更便宜的模型。如果差20%以上，必须上更精细的模型。

收敛判断：残差不是唯一标准

CFD新手最容易犯的错误：看残差降到10^-5就认为计算收敛了。

残差下降只说明方程在数值上趋于稳定，不代表你要的物理量收敛了。正确的判断标准是：

1. 监控你关心的物理量（升力系数、出口平均速度、最高温度），看它们是否稳定。
2. 做网格收敛性验证（前面反复提过）。
3. 对瞬态问题，检查时间步长是否足够小——Courant数（CFL）应该控制在1-5之间。

Fluent里可以用Report Definitions定义监控量，计算过程中实时画图。STAR-CCM+里的Monitor功能同理。

结果验证的三层框架

第一层：量级检查。你算出来的阻力系数是0.3，物理上合理吗？汽车Cd通常在0.25-0.4，飞行器在0.02-0.15。如果算出来Cd=2.5，不用比实验，量级就不对。

第二层：趋势检查。改变攻角，升力系数应该线性增长然后非线性饱和。改变雷诺数，阻力系数应该按理论趋势变化。如果趋势反了，说明网格或者湍流模型有问题。

第三层：实验比对。有实验数据最好，没有的话找文献里的同类几何对比。偏差在5%以内可以接受，10%以上要警惕。

常见坑点

坑一：二维代替三维。Re>10000的流动，三维效应不可忽略。二维模拟会高估湍流混合，结果偏于危险。判据：如果几何在第三维有变化（比如弯管、带攻角的翼型），必须用三维。

坑二：边界层分辨率不够。y+≈30用壁面函数法可以，但如果流动有分离，壁面函数法会失效。更稳妥的做法是用y+≈1的低雷诺数模型，直接解析粘性底层。

坑三：时间步长太大。瞬态CFD里，时间步长由流动的最小时间尺度决定。经验法则：一个特征时间（比如来流通过翼弦的时间）至少分成100-200个时间步。

CFD的价值在哪里

CFD仿真模拟的价值，不在于算出和实验完全一致的数字，而在于让你看清流动机制——哪里有分离？哪里有涡脱落？哪个几何参数对性能影响最大？

实验能告诉你”是什么”，CFD能告诉你”为什么”。两者结合，才是工程决策的正确姿势。