第4章：CFD仿真基础——控制方程、湍流模型与边界条件
各位工程师朋友，咱们今天聊聊CFD仿真的根基。说实话，我见过太多人上来就点“求解”，结果算出来的流场乱七八糟，还找不到原因。其实啊，根子往往就在控制方程、湍流模型和边界条件这三个地方没吃透。
我个人习惯，拿到一个新叶轮设计，第一件事不是画网格，而是先想清楚：我要解什么方程？用哪个湍流模型？边界怎么给？这三件事定下来，后面的工作才有意义。
4.1 控制方程：从N-S到RANS
流体运动的根本规律，就是三个守恒：质量、动量、能量。数学上，它们被写成了N-S方程。你想想看，这组方程几乎能描述所有流动现象——从叶片表面的层流，到尾迹里的涡脱落。
但问题来了：直接求解N-S方程（DNS）太贵了。我在做离心叶轮优化时试过，一个简单的算例就要算几周，工业设计根本等不起。
所以，工程上我们几乎都用RANS（雷诺平均N-S方程）。说白了，就是把瞬时量分解成“平均量”和“脉动量”，然后只求解平均量。代价是——多了一堆雷诺应力项，需要湍流模型来封闭。
核心要点：
- N-S方程：精确但计算量巨大，适合学术研究
- RANS方程：工程标配，计算量适中，精度够用
- LES（大涡模拟）：介于两者之间，我偶尔用于验证关键工况
4.2 湍流模型选择：没有万能药
说到湍流模型，我得先泼盆冷水：没有哪个模型是万能的。我在项目中踩过不少坑，比如用标准k-ε算一个带分离的扩压器，结果压力恢复系数差了20%。
下面这张图是我自己总结的选型逻辑，你照着走，基本不会跑偏：
我个人的经验是：
- 叶片表面流动：首选SST k-ω。它能很好地捕捉逆压梯度下的分离，我在轴流风机优化中一直用它。
- 内部通道流动：如果流动简单、无分离，标准k-ε就够了，计算快、收敛稳。
- 过渡流（层流→湍流）：用Transition SST。我记得有一次算低压涡轮叶片，用标准模型算出的损失偏小，换成Transition模型才和实验对上。
小技巧：不确定选哪个模型？先跑一个粗网格的稳态RANS，对比两三个模型的出口总压。哪个和设计值接近，就用哪个做精细计算。这招帮我省了不少时间。
4.3 边界条件设置：细节决定成败
边界条件，说白了就是告诉求解器“外面发生了什么”。设置错了，方程解得再准也没用。
我见过最典型的错误——进口给速度，出口也给速度。你想想看，质量守恒怎么满足？求解器只能硬算，结果要么发散，要么收敛到错误解。
4.3.1 进口边界
| 类型 | 适用场景 | 注意事项 |
|---|---|---|
| 速度进口 | 不可压缩流动 | 需给定速度大小和方向；湍流强度建议5%左右 |
| 质量流量进口 | 可压缩流动（如压气机） | 给定质量流量，总温总压可选 |
| 压力进口 | 进口压力已知，流量未知 | 常用于风机进口，给定总压和方向 |
4.3.2 出口边界
| 类型 | 适用场景 | 注意事项 |
|---|---|---|
| 压力出口 | 最常用，出口压力已知 | 设置合理的回流湍流强度，避免回流发散 |
| 自由出流 | 出口速度/压力未知 | 仅用于不可压缩、充分发展流动 |
| 质量流量出口 | 与质量流量进口配对 | 确保进出口质量守恒 |
避坑指南：我曾经犯过一个低级错误——出口边界离叶片太近。结果尾迹还没充分发展就被截断了，导致出口总压分布严重失真。后来我养成了习惯：出口至少延长5倍叶片弦长。
4.3.3 壁面边界
壁面处理是CFD里最容易出问题的地方。嗯，这里要注意：
- 无滑移壁面：默认设置，流体在壁面速度为零。适合大多数叶片表面。
- 滑移壁面：切向速度不为零，用于对称面或远场。
- 壁面粗糙度：如果叶片表面有涂层或铸造痕迹，记得设置粗糙度高度。我算过一台旧风机，加上0.1mm粗糙度后效率降了2%。
4.4 实战建议：从设置到收敛
最后，分享几个我自己的实操习惯：
- 先稳态，后瞬态：新模型先用稳态RANS跑通，再考虑瞬态。稳态都发散，瞬态更没戏。
- 监控关键量：我习惯监控进出口流量差、叶片表面压力、残差曲线。流量差超过1%就说明没收敛。
- 网格无关性验证：至少算三套网格——粗、中、细。如果关键结果（如效率）变化小于0.5%，就用中等网格。
好了，这一章的内容就到这里。控制方程是骨架，湍流模型是血肉，边界条件是灵魂。三者缺一不可。下次你打开CFD软件时，不妨先花10分钟想清楚这三个问题——相信我，这10分钟能帮你省下后面10小时的调试时间。
公众号:蓝海资料掘金营，微信deep3321