干了十几年流体仿真，说实话，最怕的不是模型不收敛，而是收敛了结果却离谱。你知道那种感觉吗？仿了半天，实验一测，对不上，然后从网格到边界条件一个个查，最后发现是入口段没留够。得，又白干了。
咱们今天不说教科书上的理论，就说点实际工程里的破事。流体仿真，说穿了就是CFD，但CFD不是个黑箱子，你给什么条件它就出什么结果。这玩意儿，跟做饭一样，火候差一点，味道差很多。
网格——仿真的基石，也是折磨你的源头
网格质量直接决定仿真能不能算下去，算出来准不准。很多新手上来就全局加密，然后算到天荒地老。其实没必要。你得先搞清楚你关注的是哪个区域。比如泵里的叶轮，那必须加密，但进口管道，稍微粗一点没关系。这就是局部加密的意义。
另外，边界层网格一定要画。不画边界层，壁面的传热和摩擦阻力都算不准。但画多厚？第一层高度怎么算？还是得靠经验公式估一下。Y+值，懂的自然懂。不懂的，先回去啃边界层理论。
问：网格数量越多越准吗？
答：当然不是。网格多，计算资源消耗大，而且可能因为网格质量差导致数值耗散反而更严重。关键是要在流动梯度大的地方加密，比如边界层、尾迹区。有些区域，网格太密了反而让湍流模型失效。
湍流模型——不是越复杂越好
很多商业软件默认k-epsilon，但你知道它在逆压梯度流动里表现很烂吗？分离流、旋流，用k-epsilon就是灾难。SST k-omega相对靠谱一点，但对网格要求又高。所以，你得根据你的工况来选。别一个模型用到底，那是懒人做法。
我之前做过一个三通管混合的仿真，用k-epsilon算出来混合效果极差，后来又试了RSM，结果才接近实验值。但RSM那计算量，啧啧，一个case跑了一个星期。心疼死我了！
问：怎么判断湍流模型选得对不对？
答：最好用实验数据校验。没有实验数据的话，就看流场特征是否合理。比如有没有出现不对称的涡，或者流动分离的位置对不对。再不行，拿两个模型对比一下，如果结果差很多，那就得小心了。
边界条件和求解设置——细节决定成败
入口边界给速度还是压力？出口给压力还是outflow？这些不是随便填的。给错了，结果完全两样。还有，收敛判据别只看残差。残差降到1e-4不代表物理量不变了。你得监测出口流量、壁面压力这些量，看它们稳不稳定。
常见的边界条件设置有：
- 入口：给定速度或流量；
- 出口：给定表压；
- 壁面：无滑移；
- 对称面：法向速度为零。
❗ 再说说瞬态仿真。时间步长选大了，计算发散；选小了，算到猴年马月。有些人用小的时间步长去算稳态问题，纯属浪费资源。可以用局部时间步长加速稳态收敛，但瞬态就不行了。
哦对了，还有一个常被忽视的——进口段的长度。如果你在入口直接给均匀速度，但离计算区域太近，会导致不真实的发展段。一般建议至少5倍管径。这个坑，我踩过好几次。
网格无关性验证——不是走过场
💡 其实流体仿真还有很多好用的工具，比如无量纲分析、网格无关性验证。这些不是口号，是真的能救命的。网格无关性验证，至少你要做三套网格，看关键量变不变化。变化大？继续加密。变化小？可以了。
但注意，网格无关性验证不是让你把每套网格都从头算一遍。可以先用粗网格算，再逐步细化，对比工程上关心的量，比如压降、流量。如果粗网格和细网格的差别已经小于5%，那就没必要再加密了。
仿真不是万能的
最后说一句，仿真不是万能的。有些流动现象，比如空化、多相流，模型还不成熟，有时候得靠试验。所以，老老实实做仿真的同时，也要尊重实验。两者结合，才是正道。
如果你也搞过流体仿真，一定会对上面这些有感触。欢迎留言交流，顺便点个赞再走。