基于CFturbo的泵作涡轮PAT设计及仿真分析

研究背景

全球“绿色能源”发展浪潮中，水电是可再生能源的核心组成，而泵作涡轮（PAT）因无需结构改造、成本低廉，成为无电网偏远地区小水电开发的主力设备，相比传统小型水轮机具备显著的经济优势。

但PAT的固有缺陷十分突出：最优运行区间极窄，实际工况中极易处于偏载（非设计）状态，进而引发流动失稳和压力脉动，最终导致设备振动、噪声甚至部件损坏，严重降低PAT的运行可靠性和使用寿命。而叶轮作为PAT的核心部件，其几何设计直接决定内部流场特性，因此优化叶轮设计、改善偏载工况下的压力场分布，成为提升PAT性能的关键突破口。

研究现状

目前学界对PAT叶轮的优化研究多集中于叶片前缘、叶顶间隙、叶轮切割、叶片数量等方面，比如通过前缘倒圆减少进口流动脱流、优化叶顶间隙降低泄漏损失等，这些改造均能在一定程度上降低水力损失、提升效率。

但现有研究存在明显空白：叶片尾缘轮毂位置对PAT偏载工况下压力场的影响尚未被系统探究，而尾缘设计直接影响叶轮出口流场结构，进而关联整个PAT流道的压力脉动和转子-静子相互作用（RSI）。同时，PAT偏载工况下的非定常压力场数值模拟难度较高，需精准的几何建模工具和可靠的数值计算方法，对研究手段也提出了严格要求。

实施方案

本次研究的核心亮点体现在工具选型的专业性和研究视角的创新性，其中CFturbo在几何建模阶段发挥了基石作用，为后续数值模拟奠定了精准的模型基础，具体实施方案如下。

1. CFturbo实现定制化叶轮设计

采用CFturbo这款专业的涡轮机械CAD设计软件，精准设计了一台比转速38.1的六叶片单吸离心式PAT，核心完成了三种不同叶片尾缘轮毂位置（Rh）的叶轮定制化建模——保持叶片罩壳端固定，将轮毂端沿径向分别设置为15mm（Rh15）、20mm（Rh20）、25mm（Rh25），同时保证叶轮进出口直径、叶片包角等其他关键参数一致，仅单一变量探究尾缘位置的影响，建模精度直接决定了后续研究的科学性。

图[1] 轮毂端径向15mm（Rh15）、20mm（Rh20）、25mm（Rh25）

2. 全流道非定常数值模拟

采用ANSYS ICEM CFD划分结构化六面体网格，通过网格无关性验证确定150万节点为最优网格数，并以CFX为求解器，基于k-ε湍流模型开展偏载工况（Q/Qd=0.67）下的非定常数值模拟，同时在蜗壳、叶轮、出水管全流道布置24个压力监测点，实现压力脉动的精细化捕捉。

图[2] 网格示意图

3. 多维度分析方法

结合压力云图、流线图分析流场结构。

图[3] 压力、流线变化云图

通过时间域曲线观察压力变化趋势并利用快速傅里叶变换（FFT）将压力数据转换至频率域，从流场结构、时域变化、频域特征三个维度全面解析压力脉动特性。

图[4] 压力变化频域

结果分析

研究通过对蜗壳、叶轮、出水管三大核心流道的压力场和压力脉动分析，得出一系列关键结论，各流道的特征变化均有清晰的图像支撑：

1. 蜗壳流道：压力沿程降低，舌部脉动最剧烈

蜗壳内压力从进口向叶轮进口方向持续降低，压力分布的非对称性源于蜗壳自身几何特征，且蜗壳舌部（监测点S6）是压力脉动振幅最大的区域。

图[5] 离心泵压力变化时域

频率域分析显示，三种模型的主脉动频率均为叶片通过频率（BPF，289.8Hz），由转子-静子相互作用（RSI）主导；其中Rh20模型的蜗壳各监测点主频率脉动振幅均为最大，Rh15最小。

表[1] 三种模型主脉动频率表

2. 叶轮流道：尾缘越远，低压区越大，Rh20脉动最显著

叶轮内压力从进口向出口逐渐降低，叶片压力侧前缘为高压区，吸力侧易因偏载工况出现流动分离，形成脱流涡。

图[6] 叶片流道脉动图

叶片尾缘轮毂距离（Rh）越大，叶轮出口低压区范围越广，Rh25的低压区最明显；压力脉动方面，叶轮内主频率同样为BPF，脉动振幅从叶片前缘向尾缘递减，且Rh20模型的各监测点振幅均为三者最高，Rh15仍为最低。

表[2] 叶轮流场内不同监测点处的主要频率压力脉动幅度

3. 出水管流道：压力径向递增、轴向递减，尾缘位置决定进口压力水平

出水管内压力沿径向从轴心向管壁递增，沿轴向从进口向出口递减，这与流场速度的径向、轴向变化直接相关。

图[7] 在t=0.020秒时，出水管入口处的压力分布等值线图：(a)Rh15，(b)RH20，以及(c)Rh25

Rh15模型的出水管进口压力水平最高，Rh25最低；频率域上，主频率仍为BPF，脉动振幅依旧呈现Rh20＞Rh25＞Rh15的规律。

表[3] 出口管道流道内不同监测点处的主要频率压力脉动幅度

三种模型的全流道压力脉动均由转子-静子相互作用（RSI）主导，叶片通过频率（BPF）为核心主频率，同时伴随少量由局部流动失稳引发的低频分量；叶片尾缘轮毂位置与压力脉动振幅并非线性相关，Rh20成为压力脉动的峰值点，Rh15则表现出最优的低脉动特性。

研究总结

PAT 压力脉动由流场结构主导，叶片尾缘轮毂位置变化会改变低压区、脱流涡等流场形态，进而影响脉动特征。偏载工况下转子-静子相互作用为核心诱因，叶片通过频率为全流道主频率。三种模型中 Rh15 脉动幅值最低、Rh20 最剧烈，且无明显线性规律，可为叶轮尾缘优化提供工程参考。基于CFturbo的建模与分析方法为涡轮机械研究提供了思路，对小水电 PAT 稳定运行具有重要工程价值。