第一次做离心泵汽蚀（空化）仿真，想预测泵的必需汽蚀余量（NPSHr）。用 Fluent 的 Schnerr-Sauer 空化模型，算出来 NPSHr = 2.5m。结果实验测出来 3.2m，差了 22%

我检查了网格、空化模型、边界条件，都没发现问题。后来请教师兄，师兄问了四个问题，每个都戳中了我的盲点。

第一问：你的不凝性气体含量设了多少？

我说：「默认值 1.5e-5。」 师兄说：「实际自来水里溶解气体含量是 1e-4 到 5e-4，是默认值的 10-30 倍。

气体含量影响空化核的数量，直接影响空化起始。」 空化核是液体里微小的气泡或气体微团，是空化的"种子"。

第二问：你的网格在叶轮进口加密了吗？

离心泵空化最先发生在叶轮进口低压区——叶片前缘背面。如果这个区域网格太粗，低压区压力分布算不准，空化起始位置就不对。

我的教训：叶轮进口网格 2mm，低压区压力梯度完全没解析出来。加密到 0.3mm 后，空化区域和高速摄影照片吻合了。

网格要求：叶片前缘附近网格 ≤ 叶片厚度的 1/5，空化区域至少 5-8 层网格。

气体含量越高，空化核越多，空化越容易发生。 我重新设成 3e-4，NPSHr 从 2.5m 变成了 3.0m，接近实测值了。 

第三问：你的空化模型参数调了吗？

Schnerr-Sauer 模型有两个关键参数：

气泡核密度 n₀：默认值 1e10 /m³，实际可能在 1e8 到 1e12 之间

气泡核半径 R₀：默认值 1e-6 m（1μm），实际可能在 0.5μm 到 10μm 之间

这两个参数对空化起始和气泡生长影响很大。

我有个项目，把 R₀ 从 1μm 调到 5μm，NPSHr 从 2.8m 变成了 3.3m。 建议：先用默认值跑一遍，然后和实验对比，反推合适的参数。

第四问：你的时间步长够小吗？

空化是瞬态过程，气泡的生长和溃灭在毫秒量级。时间步长太大，会漏掉气泡溃灭过程，空化强度算低。

经验法则：时间步长 ≤ 气泡特征时间/50 气泡特征时间 ≈ R_max × √(ρ/ΔP)

对于离心泵，典型时间步长 1e-5 到 5e-5 秒。

我之前的时间步长 1e-4 秒，气泡溃灭过程完全没解析，空化区域算小了 30%。

这些细节容易忽略

忽略 1：旋转参考系的设置

离心泵仿真用 MRF（多参考系）或滑移网格。MRF 是稳态近似，算不出空化的瞬态脉动。

如果要算空化对性能的影响，必须用滑移网格做瞬态仿真。

有个项目用 MRF 算 NPSHr 2.5m，用滑移网格算 3.1m，实验 3.2m。

忽略 2：壁面粗糙度

叶轮表面粗糙度影响边界层发展，进而影响低压区的压力分布。新泵和运行一段时间后表面粗糙度不同，NPSHr 也会变化。

有个项目，新泵仿真 NPSHr 2.8m，实验 3.0m（运行 1000 小时后粗糙度增加）。

加上粗糙度模型后，仿真 3.1m，吻合了。 

忽略 3：温度影响 

液体温度影响饱和蒸汽压和表面张力：

20°C 时水饱和蒸汽压 2.34kPa

60°C 时水饱和蒸汽压 19.9kPa（差了 8.5 倍！）

温度越高，饱和蒸汽压越高，越容易汽蚀。 

有个项目按 20°C 算 NPSHr 2.5m，实际水温 50°C，NPSHr 实测 3.5m。

验证方法

怎么知道汽蚀仿真算得对不对？

NPSHr 曲线对比：算不同流量下的 NPSHr，和实验 NPSHr 曲线对比。误差<10% 可以接受。

高速摄影对比：用高速摄影拍叶轮内的空化区域，对比仿真气相体积分数云图。

压力脉动对比：空化会产生压力脉动，对比仿真和实测压力脉动频谱。

经验总结

我现在做离心泵汽蚀仿真的流程是：

实测水中溶解气体含量（或按水质估算）

叶轮进口网格加密（≤叶片厚度 1/5）

设置合理的空化模型参数（气泡核密度和半径）

时间步长 ≤ 气泡特征时间/50 用滑移网格做瞬态仿真

设置温度相关的饱和蒸汽压

和 NPSHr 曲线 + 高速摄影对比验证