泵CFD分析
泵CFD的理论基础
概述
离心泵的CFD分析预测什么性能？
扬程（Head）、效率和轴功率这三个是基础。H-Q特性曲线的制作是CFD的主要目的。
扬程和效率的定义
请教扬程的计算式。
泵的扬程是入口和出口的全水头差。
CFD中从全压差直接计算最为简便。$H = (p_{t2} - p_{t1})/(\rho g)$。
效率分为水力效率和全效率。
$\tau$ 是从CFD获得的叶轮转矩，$\omega$ 是角速度。
水力效率和全效率有什么区别？
水力效率仅包含流体力学损失，不包括盘摩擦和泄漏。全效率包含盘摩擦、泄漏流量和机械损失。CFD直接得到的是水力效率，盘摩擦和泄漏需要在模型中包含磨损环隙间才能得到。
Euler扬程（理论扬程）
欧拉方程可以推导出理论扬程吗？
对于离心泵，假设入口没有旋转速度（$C_{\theta 1}=0$），则 $H_{Euler} = U_2 C_{\theta 2}/g$。考虑滑移因子 $\sigma_s$ 的扬程为 $H_{th} = \sigma_s \cdot H_{Euler}$。CFD的扬程对应于此理论扬程加上水力损失。
MRF法的定常分析
泵的CFD采用MRF法是通用做法吗？
获取H-Q曲线的标准方法是MRF法（定常）。带蜗壳的情况下使用冻结转子或滑移网格。不带蜗壳、配有导叶的泵可以使用混合平面法。
比速作为指南针——斜流泵设计的出发点
涡轮泵设计的第一步是计算比速度（Specific Speed，Ns）。Ns = n√Q / H^(3/4)（n：转速，Q：流量，H：扬程）这个无量纲数决定了泵的形式（离心、斜流、轴流）。在斜流泵的设计范围内（Ns = 400～1200），从理论上讲，应该首先用1D设计（速度三角形）来确定基本尺寸，然后用CFD来验证3D流动的详细信息。这是既定的两阶段方法。比速度的理论由20世纪初的美国水力学家体系化，之后在欧美和日本出现了采用不同单位系的混乱。现代已通过IEC标准统一了采用m³/s和m的无量纲比速度作为国际标准。
泵CFD的数值计算方法
网格生成
离心泵网格生成应注意什么？
叶轮采用TurboGrid生成结构网格的品质最高。蜗壳采用非结构四面体/多面体网格。
| 区域 | 网格类型 | 单元数目安 | 工具 |
|---|---|---|---|
| 叶轮 | 结构网格（H/J/L+O-grid） | 50～150万/周期 | TurboGrid |
| 蜗壳 | 非结构四面体+棱柱 | 100～300万 | Fluent Meshing、STAR-CCM+ |
| 吸入管 | 结构或非结构 | 20～50万 | 任意 |
| 磨损环隙间 | 结构（六面体） | 10～30万 | 手动 |
磨损环隙间也要加入模型中吗？
评估泄漏流量的影响必须包含隙间。隙间为0.2～0.5mm，非常狭窄，因此径向至少10个单元，轴向建议50个以上。
湍流模型的选择
泵应采用哪种湍流模型？
SST k-omega是通用选择。对叶面逆压梯度和分离的预测优势明显。泵的情况下，叶片数少（5～7片）叶片负荷高，k-epsilon对分离预测往往不够。
壁面函数和低雷诺数哪个应该用？
边界条件
泵的典型边界条件是什么？
出口固定静压，入口流量变化，这是最稳定的设置。
轴流和离心之间——斜流泵CFD的数值难题
斜流泵的比速度（Ns）处于轴流和离心泵之间的中间范围（Ns = 400～1200），这导致两种流动机制混合，使数值分析变得困难。轴流成分中翼端涡和二次流占主导，离心成分中科里奥利力导致流体弯曲。在这种混合领域，湍流模型的选择直接影响泵效率预测，许多对比研究表明k-ω SST模型通常比常规k-ε精度好几个百分点。网格设计需要在轴向、周向、径向三个方向上达到平衡，而不是简单地在某一方向细化。
泵CFD的实务应用
H-Q特性的计算步骤
H-Q特性曲线怎样制作？
1. 在设计流量 $Q_d$ 下进行定常MRF（或冻结转子）计算直到收敛
2. 流量设定在 $0.2Q_d$ ～ $1.4Q_d$ 范围内，取7～10个点
3. 各点改变入口质量流量重新计算（使用前一个点的结果作为初值）
4. 各点计算扬程H、转矩τ、效率η并绘制
低流量下收敛不好是为什么？
低流量时叶轮入口的入射角增大，叶面产生大规模分离。定常计算试图将非定常的分离结构强制收敛到一个解，就会出现振荡。流量低于0.3$Q_d$ 时通常需要非定常计算。
与实验的对比
CFD结果与泵试验的吻合度如何？
BEP是什么意思？
Best Efficiency Point，最高效率点。泵的设计点通常设在BEP附近。按照JIS B 8301和ISO 9906标准进行泵试验并与CFD数据对比。
盘摩擦的评估
盘摩擦损失也能用CFD计算吗？
如果模型包含叶轮背面（轮缘侧、毂侧）的隙间空间，就能计算出来。如果不包含，用经验式补正。
$C_M$ 是摩擦转矩系数，是隙间比 s/r 和雷诺数的函数。
治水与斜流泵——支撑防灾基础设施的CFD
日本的地下放水路（如首都圈外围放水路等）和排水站采用的大型斜流泵需要在突发大雨时处理超过每秒200m³的流量。设计阶段的CFD分析是全流量范围（设计流量的50～120%）效率和压力特性的预测，以及吸入池内气泡卷入风险评估的必要项目。特别是吸入漩涡的发生在实机中成为致命故障，因此需要对整个吸入区域（包括吸入池形状）进行自由表面CFD（VOF法）模拟。与模型泵试验的效率相关度是国内制造商的品质标准，要求在±1%以内。
泵CFD的软件比较
求解器比较
泵CFD最适合的求解器是哪个？
CFX的实绩最多啊。
在泵行业CFX事实上是标准。耦合型求解器的稳定性适合泵的低流量运行和两相流计算。自2024年以来Fluent通过Turbo Workflow的加强正在迅速追赶。
STAR-CCM+怎么样？
多面体网格自动生成的优势明显，蜗壳等复杂形状网格工作量少。但与TurboGrid生成的结构网格相比翼间流路精度有差距，理想方法是两者的结合。
计算成本的估算
一台离心泵的分析需要多长时间？
| 分析内容 | 单元数 | 内核数 | 所需时间 |
|---|---|---|---|
| 定常MRF（1工况点） | 200万 | 16 | 2～4小时 |
| H-Q特性（8工况点） | 200万 | 16 | 16～32小时 |
| 非定常滑移网格（1工况点） | 200万 | 32 | 12～24小时 |
| 空泡评估 | 300万 | 64 | 24～48小时 |
H-Q特性1～2天就能出来啊。
用自动脚本逐次计算各流量点，夜间批处理的话，次日早上就有结果。这是CFD的一大优势。
泵CFD工具选型的现实——通用机vs专用机的角力
在斜流泵CFD分析工具选型中，"涡轮机械专用功能是否具备"是重要评估轴。ANSYS TurboGrid或Numeca AutoGrid5这样的专用叶轮网格生成工具，即使对复杂的斜流翼形也能在数小时内生成高品质的Hex-dominant网格。相比之下，通用网格划分工具（ANSYS Meshing、Pointwise）即使经验丰富的人员也要花1～2天才能达到同样品质。成本上专用工具的年度许可证达数百万日元，因此年度分析案件较少的企业倾向采用OpenFOAM的turboBladefork版本。农业、治水类泵制造商普遍重视与传统1D设计工具（如Pump-Flo等）的联动。
泵CFD的前沿研究
多级离心泵
多级泵的CFD分析困难吗？
各级用混合平面或冻结转子来连接。级间回流通道（反向导叶）的损失预测是精度的关键。
回流通道是什么？
第1级出口的旋转流去除并整理第2级入口条件的流路。这里的损失占全级效率的2～5%。虽然容易采用粗网格，但应该精细建模。
浆液泵的CFD
包含固体颗粒的浆液泵也能用CFD分析吗？
可以。用欧拉-欧拉法或欧拉-拉格朗日法计算气液两相流。
| 方法 | 适用范围 | 求解器 |
|---|---|---|
| 欧拉-欧拉法 | 高浓度浆液（体积分率>10%） | CFX多相流、Fluent混合物模型 |
| 欧拉-拉格朗日法（DPM） | 低浓度（<5%）、粒子轨迹评估 | Fluent DPM、STAR-CCM+拉格朗日 |
| DEM耦合 | 大颗粒、颗粒间接触重要 | STAR-CCM+ DEM、EDEM耦合 |
也能预测磨损吗？
能。用Finnie或Oka磨损模型估算叶面磨损量分布。从粒子碰撞速度和角度计算磨损率。Fluent的DPM Erosion Model和STAR-CCM+的Erosion Model都可用。
泵作为涡轮 (PAT)
把泵反向作为涡轮的PAT也能用CFD评估吗？
能。在小水电领域PAT备受关注，泵反流、反向运行时的性能评估。CFD中将转向和流向反向进行分析。由于泵模式的BEP位置不同，需要重新构建性能图。
斜流泵优化的最前线——位相分析捕捉不稳定流
斜流泵前沿研究聚焦于"旋转失速的早期检测"。旋转失速发生在部分流量工况，导致泵特性曲线折返，成为运行不稳定的源头。从2010年代后期开始，人们将本征正交分解（POD）和动态模态分解（DMD）应用于CFD结果，抽取流场的支配模式。日本农业用大型扬水泵（扬程20～50m、流量超10m³/s）的更新设计中已采用这种CFD + DMD分析，能够识别不稳定流动发生的条件，将可避免流量范围比传统经验法则扩大30%。
泵CFD的故障排除
低流量发散
低流量计算发散的原因是什么？
低流量时叶轮入口出现大的逆旋转（负旋转），叶前缘发生大规模分离。定常计算试图强制将这种非定常现象收敛到一个解，就会振荡。
对策:
- 切换为非定常计算（滑移网格）
- 时间步设为叶片通过时间的1/30以下
- 降低缓和系数（CFX的Timescale Factor：0.5）
反流边界条件
入口出现反流该怎么办？
低流量时叶轮入口的毂附近会出现反流。入口BC设为常规Inlet时，反流单元会产生非物理的压力，导致发散。
对策:
- 将入口边界延伸到上游（管径的5倍以上），使反流不达到边界面
- CFX中改为Opening BC
- Fluent中设为Pressure Inlet with backflow direction
蜗壳的干涉
蜗壳舌部附近压力振荡。
舌部是泵最大非定常力的发源地。冻结转子会产生过大的干涉。换成滑移网格后时间平均值会稳定。
检查清单
请总结泵CFD的最终检查项目。
斜流泵的S形特性——设计者的大敌
斜流泵的特性曲线（H-Q曲线）出现"S字形"是设计者的噩梦。曲线从右上向左下折返呈S形时，同一扬程对应多个流量，并联运转时出现不稳定的流量分配。CFD要捕捉S字特性，需要对流量低至设计流量20%的低流量工况进行分析，这个范围内RANS模型精度急剧下降，计算误差增大。经验丰富的工程师制定了"低流量范围用LES或DES验证"的规则，仅依赖RANS结果确定设计是危险的，业界多次重复这个教训。
过
更多内容
错误