泵空化
泵空化的理论基础
概述
泵的空化就是气泡冒出来把机器坏掉的现象吧？
大致是这样。液体的静压低于饱和蒸气压时，蒸气泡就会产生，在下游高压区域崩裂。崩裂时局部会产生数GPa的冲击压力，对叶轮表面造成侵蚀。
数GPa！？那肯定要坏了...
性能面也很严重。扬程下降、振动增大、噪音增加。因此空化回避是泵设计的最优先课题之一。
NPSH的定义
经常听到NPSH，正确的定义是什么？
NPSH（Net Positive Suction Head）表示吸入侧液体相对于蒸气压的余量。
$p_{atm}$:大气压、$p_v$:饱和蒸气压、$z_s$:液面至泵中心的高度、$h_f$:吸入管道的损失水头。这是系统侧的NPSH（Available）。
泵侧也有NPSH吧？
NPSH_r（Required）是泵自身所需的最低NPSH，定义为扬程下降3%的点。安全上要满足。
也可以用Thoma数（空化系数）表示。
Rayleigh-Plesset方程
有描述气泡成长和崩裂的公式吗？
Rayleigh-Plesset方程是基本的。
$R$:气泡半径、$p_B$:气泡内压、$p_\infty$:周围压力、$S$:表面张力。CFD的空化模型是将其简化为质量输运方程。
空化困扰潜艇的历史
泵空化的气泡力学（Rayleigh-Plesset方程）作为实际问题被重视，始于第一次世界大战期间的潜艇螺旋桨。叶尖附近气泡激烈地产生和崩裂，导致推进效率下降和金属侵蚀的双重问题。战后研究使得NPSH（正净吸入水头）的概念得以系统化，并被传承到现代泵设计中。
泵空化的数值计算手法
均质混合模型
CFD怎样计算空化？
最广泛使用的是均质混合（Homogeneous Mixture）模型。将液相和蒸气相视为单一流体，求解蒸气体积分率 $\alpha_v$ 的输运方程。
$\dot{m}^+$ 为蒸发（气泡产生），$\dot{m}^-$ 为凝聚（气泡崩裂）的源项。
源项有什么模型？
比较代表的3个模型。
| 模型 | 特点 | 使用求解器 |
|---|---|---|
| Zwart-Gerber-Belamri | 基于核生成点密度，参数易调整 | CFX（默认）、STAR-CCM+ |
| Schnerr-Sauer | 基于Rayleigh-Plesset，指定气泡数密度 | OpenFOAM（interPhaseChangeFoam）、Fluent |
| Singhal (Full Cavitation) | 考虑非凝聚气体，实用 | Fluent |
CFX以Zwart模型为标准吧。
是的。蒸发系数 $F_{vap}=50$、凝聚系数 $F_{cond}=0.01$ 为默认值。核生成点体积分率 $\alpha_{nuc}=5 \times 10^{-4}$、气泡初始半径 $R_B=10^{-6}$ m 为典型值。
数值设置的要点
空化计算收敛有什么窍门吗？
这是困难的计算，有几个要点。
- 时间步长: 必须非定常。叶片通过时间的1/20～1/50为目标
- 收敛标准: RMS残差 $10^{-5}$ 以上为目标（空化振动可能只降到$10^{-4}$）
- 初值条件: 先用非空化状态定常收敛，然后开启空化模型
- 压缩性: 由于蒸气-液体密度比大，数值上需要考虑压缩性设置
Zwart、Merkle、Singhal——三大模型之争
在CFD空化解析现场，经常听到"用哪个模型好？"的提问。Zwart模型明确处理核生成位点密度，Merkle模型以压力差直接反应的经验形式，Singhal（Full Cavitation）模型考虑最复杂的溶解气。基准研究表明不同情况排名会变，"这个最好"难以下定论。初值、网格、经验常数设置对结果影响很大，与自身实验的验证必不可少。
泵空化的实务应用
NPSH特性曲线的获取步骤
用CFD怎样描绘NPSH曲线？
逐步降低入口全压的方法最常见。
1. 基准计算：充分高的NPSHa（无空化）下定常收敛
2. 入口压力下降：入口全压以0.1～0.2 atm的步长逐步降低
3. 各点非定常计算：开启空化模型，计算数转
4. 记录时间平均扬程：扬程从基准值下降3%的点为NPSH_r
0.1 atm的步长是不是太粗了？
NPSH_r附近扬程急剧下降，所以先粗略掌握全貌，再在3%下降附近用0.02～0.05 atm的步长细化是有效率的。
可视化与评价
空化结果怎样看？
蒸气体积分率 $\alpha_v = 0.5$ 的等值面表示空化腔形状。在CFD-Post中确认以下内容。
| 可视化项目 | 确认要点 |
|---|---|
| $\alpha_v$ 等值面 | 空化腔位置、大小、对叶面的贴附 |
| 叶面压力分布 | 吸入面何处低于蒸气压 |
| 叶面侵蚀风险指标 | 气泡崩裂压力的累积（CFX的Erosion Model） |
| 径向空化分布 | 毂部 vs. 前缘的区别 |
侵蚀风险也能用CFD评估吗？
CFX内置了Cavitation Erosion Model，从气泡崩裂的能量密度生成侵蚀风险图。但定量寿命预测还需要与实验校准。
设计改进的方向
怎样改进空化性能？
汇总主要设计参数的影响。
| 参数 | 对NPSH_r的影响 | 权衡 |
|---|---|---|
| 增加诱导器 | 大幅降低 | 结构复杂，成本增加 |
| 叶片入口角优化 | 降低入口冲击角改善 | 偏离设计点性能 |
| 叶片数增加 | 叶片负荷降低而改善 | 摩擦损失增大 |
| 吸入口径扩大 | 流速降低而改善 | 泵尺寸增大 |
NPSH特性试验为何以3%落差定义
泵空化实验中经常出现"NPSH₃"一词。这是"扬程H下降3%时的NPSH值"的定义，但为什么是3%呢？实际上没有明确的物理根据，而是1940～50年代国际标准化的经验值。在CFD中预测NPSH特性时，就是"重现扬程下降3%气相体积分率的条件"，这在数值上是否能稳定检出成为实务中的重要检验点。
泵空化的软件比较
主要求解器的空化功能比较
什么求解器适合空化解析？
比较各求解器的支持情况。
| 功能 | Ansys CFX | Ansys Fluent | STAR-CCM+ | OpenFOAM |
|---|---|---|---|---|
| 均质混合空化 | Zwart（标准） | Schnerr-Sauer, Singhal | Rayleigh-Plesset | interPhaseChangeFoam |
| 侵蚀模型 | Erosion Model | UDF扩展 | Cavitation Damage | 无（可自制） |
| 涡轮专用网格生成 | TurboGrid | Turbo Meshing | 内置模板 | 无 |
| 非定常旋转计算 | Transient Rotor-Stator | Sliding Mesh | Rigid Body Motion | cyclicAMI |
CFX功能最全啊。
在泵空化解析方面，CFX实绩最多。耦合型求解器的稳定性适合密度比大的空化计算，TurboGrid的工作流也成熟了。
Fluent怎样？
Fluent 2024R1以后强化了Turbomachinery工作流，含空化泵解析也达到了实用水平。特别与多面体网格的组合，自动网格生成很方便。
OpenFOAM中的空化
OpenFOAM也能做空化吗？
interPhaseChangeFoam求解器实现了Schnerr-Sauer和Kunz模型。但与旋转机械组合时cyclicAMI设置繁琐，收敛稳定性常弱于商用版。科研用途或预算受限时的选择。
计算成本估计
空化计算要多久？
一般为非空化计算的5～10倍。必须非定常，蒸气-液体界面变动要求小时间步。
| 模型规模 | 单元数 | 核心数 | 计算时间（转数） |
|---|---|---|---|
| 单级泵（1叶列） | 200万 | 32 | 4～8小时 |
| 带诱导器 | 500万 | 64 | 12～24小时 |
| 多级泵 | 1000万 | 128 | 2～4天 |
STAR-CCM+在泵业为何强势
作为空化解析工具，STAR-CCM+在给水、化工等泵业应用广泛，背景是旋转区域（MRF、Sliding Mesh）与空化模型的一站式设置便利。特别是"VOF+空化模型"组合稳定，易于自动批量计算NPSH特性曲线的工作流。ANSYS Fluent虽有等价功能，但从GUI体验来看，实务者多数倾向STAR-CCM+。
泵空化的前沿研究
非定常空化现象
空化非定常会有什么问题？
叶面上的空化腔周期性成长、崩裂反复的"云团空化"特别严重。崩裂时的冲击波打击下游叶面，侵蚀加速。
云团空化的频率是多少？
Strouhal数 $St = fL/V \approx 0.2 \sim 0.4$ 。$f$:振动频率、$L$:空化腔长度、$V$:主流速度。叶弦长50mm、流速20m/s的话约100～160Hz。
水下噪音的预测
空化噪音能用CFD预测吗？
用Ffowcs Williams-Hawkings（FW-H）方程的混合法。CFD计算近场非定常压力，用声学类比计算远场音压。
Fluent 和 STAR-CCM+ 内置FW-H求解器。
有空化时噪音怎样变？
空化时宽带噪音增加10～20dB。特别是气泡崩裂源的高频成分（1～100kHz）显著。对噪声规制严格的住宅小区泵等，需充分确保NPSH余量的设计。
流体-结构耦合（FSI）
能分析空化导致的叶轮振动吗？
CFD非定常压力场转写到FEM结构模型的单向FSI较实用。确认空化腔崩裂的冲击激励与叶轮固有频率的共振。Ansys System Coupling或STAR-CCM+的Co-Simulation可执行。
需要完全的双向FSI吗？
泵叶轮刚性高，多数情况单向FSI足够。但树脂叶轮或薄壁翼就需要流体-结构的双向耦合了。
火箭发动机涡轮泵与空化
H-IIA火箭的涡轮泵在-253℃液态氢环境下工作，每秒升压数百升液体。在极限环境中控制空化，在泵吸入前级装有"诱导器"——低NPSH预备叶。诱导器容许部分空化同时升压，向主叶轮供应稳定流。这一设计思想是JAXA和IHI多年实验和CFD磨练的结果，是宇宙开发与泵空化研究的深度交集领域。
泵空化的故障排查
常见失误
空化解析常见错误有什么？
经验上，频繁出现的问题有。
1. 饱和蒸气压设置错误
水在25℃时3170Pa，80℃时47400Pa，随温度变化大。高温水泵用常温蒸气压计算，误判"不会空化"。
2. 用定常计算评价空化
定常也能出现空化吧？
定常计算中空化腔"冻结"显示，但真实空化本质非定常。NPSH_r附近定常扬程与非定常时间平均常差异大。NPSH_r评估必须用非定常计算。
3. 网格不足
空化界面是尖锐密度梯度，粗网格导致空化腔形状模糊，体积高估。叶面吸入侧网格特别要细。
与实验的对比验证
怎样验证CFD结果？
分阶段验证较好。
| 验证阶段 | 对比对象 | 允许误差 |
|---|---|---|
| 非空化 | QH特性曲线 | 扬程±3%以内 |
| 空化开始 | 空化腔目视观察 | 发生位置一致 |
| NPSH_r | 3%扬程下降点 | ±0.5m以内 |
| 侵蚀位置 | 实机侵蚀模式 | 定性一致 |
NPSH_r的±0.5m相当严格啊。
是的。CFD若能±0.3m精度预测NPSH_r就非常好了。更高精度则需用实验数据校准空化模型常数。
温度效应（热动力学效应）
高温液体空化会变吗？
液氮或LNG等低温流体，或高温水中，气泡产生时潜热导致周围液温下降，蒸气压降低的"热动力学效应"会出现。结果空化被抑制。高温水泵的NPSH_r比冷水低。CFX有考虑这效应的选项。
"CFD无空化但叶片被削"事件
在水力机械现场稀有地出现"CFD预测无空化，却发现叶面侵蚀"的情况。原因疑似"涡空化"——叶尖泄漏涡中心局部低于蒸气压的现象。定常RANS计算过度扩散这个涡核，难以检出。用LES或非定常计算需高分辨涡核低压，这是侵蚀预测易遗漏的典型陷阱。
价值
更多
要报告