穿孔叶片离心泵空化性能的数值模拟

陶成, 宋文武, 邓强, 宿科, 周月

PDF(3135 KB)
中国农村水利水电 ›› 2021 ›› (7) : 192-197.
农田水利

穿孔叶片离心泵空化性能的数值模拟

作者信息 +

Numerical Simulation of Cavitation Performance of Perforated Vane Centrifugal Pumps

Author information +
稿件信息 +

摘要

为分析进口位置附近开孔对离心泵空化性能的影响,以某低比转速离心泵作为研究对象,利用CFD仿真软件基于RNG k-ε湍流模型对模型进行全流道定常湍流空化数值模拟,研究在叶片空化初生位置,开孔面积、开孔形状对离心泵性能的影响。研究发现,对于本文的低比转数离心泵,开孔面积过大会造成能量损失,当孔型为梅花孔时,开孔面积过小会加剧叶轮空化;不同开孔面积,不同孔型对离心泵空化性能影响差异较大,当穿孔面积为1/4 π mm2,孔型为方形孔时,扬程提高了9.8%,并且显著增强离心泵的抗空化性能;当孔型和开孔面积能有效抑制叶轮空化时,其叶轮进口低压区都被截断,呈间断性分布;当开孔面积为1/4 π mm2,孔型为方形孔和开孔面积为25/4 π mm2,孔型为梅花孔时,对离心泵湍动能减小明显,减小了原模型空化带来的能量损失,有效地抑制了叶轮内空化。

Abstract

For the analysis of important position near the hole of the influence of the centrifugal pump cavitation performance, with a low specific speed centrifugal pump as the research object, the CFD simulation software is used based on RNG k-ε turbulence model to full port steady turbulent flow cavitation model numerical simulation, study in blade cavitation primary position, opening area, open hole shape affect the performance of centrifugal pumps. It is found that for the low specific revolution centrifugal pump in this paper, too large opening area will cause energy loss. When the hole type is quincunx hole, too small an opening area will aggravate impeller cavitation. The cavitation performance of the centrifugal pump is significantly affected by different hole areas and different hole types. When the hole area is 1/4 π mm2 and the hole type is square, the head increases by 9.8% and the anti-cavitation performance of the centrifugal pump is significantly enhanced. The low pressure area at the inlet of the impeller is truncated and distributed intermittently when the hole type and hole area could effectively inhibit the impeller cavitation. When the hole area is 1/4 π mm2, the hole type is square hole and the hole area is 25/4 π mm2, and the hole type is quincunx hole, the turbulent kinetic energy of the centrifugal pump is significantly reduced, the energy loss caused by the original model cavitation is reduced, and the cavitation in the impeller is effectively inhibited.

关键词

离心泵 / 穿孔面积 / 孔型 / 空化 / 低压区

Key words

centrifugal pump / perforation area / pass / cavitation / low pressure area

基金

四川省科技厅项目(2016JY0187)
四川省教育厅重大培育项目(14CZ0013)

引用本文

导出引用
陶成 , 宋文武 , 邓强 , 宿科 , 周月. 穿孔叶片离心泵空化性能的数值模拟[J].中国农村水利水电, 2021(7): 192-197
Cheng TAO , Wen-wu SONG , Qiang DENG , Ke SU , Yue ZHOU. Numerical Simulation of Cavitation Performance of Perforated Vane Centrifugal Pumps[J].China Rural Water and Hydropower, 2021(7): 192-197

0 引 言

随着社会的快速发展,气液混输泵在国民经济中占有重要的地位,因此对离心泵在气液两相下的研究显得十分具有意义。空化空蚀破坏是离心泵运行中常见的一种现象。空化空蚀破坏会造成离心泵效率、扬程降低,水力机组振动,产生噪音,当破坏程度严重时,气泡会堵塞整个流道,致使水力机械无法正常运行1-5
气液两相中离心泵内部流动是包含多相、非定常、湍流相等结合的一种复杂的流动状态,为保证离心泵在气液混输中保持稳定、高效、安全的运行,国内外学者做了大量的研究。 Murakami等人6通过高速摄影方法得到了半开式叶轮的空泡的运动情况和含气率较小时空泡的迁移过程。Ji B等人7通过数值模拟和实验研究的方法分析得出绕NACA66周围的空穴体积变化是其周边压力脉动的主要来源。Sekoguchi等人8基于数值模拟和可视化试验的方法分析离心泵在气液两相下内部流动和气相体积分数变化情况。宋文武等人9基于高速离心泵回流旋涡及空化特性的研究,得出了诱导轮的一种优化方案用于提高高速离心泵的抗空化性能。王维军等人10采用改进的空化模型、修正RNG k-ε模型,提出空化初生判定准则,系统的划定离心泵的空化区域。刘小兵等人11运用数值模拟和试验的方法总结出混输泵在几种典型空化工况下的运输性能和叶轮内部空泡的分布规律。罗旭等人12基于RNG k-ε模型和 Rayleigh-Plesst空化模型方程分析了空化在高速离心泵叶轮内部的演变过程,得到流体诱发压力脉动的变化规律。李金琼等人13采用RNG k-ε模型和Rayleigh-Plesst方程的均相流空化模型分析了粗糙度对离心泵进口回流非定常特性的影响。胡赞熬等人14应用CFD仿真软件对分析离心泵在不同气蚀余量条件下,开孔位置和开孔直径对离心泵空化性能的影响。朱晓东等人15采用欧拉多相流模型和RNG k-ε湍流模型,分析了不同翼型结构周围的压力分布。赵伟国等人16基于SST k-ω湍流模型和Zwart空化模型,分析了叶轮进口处开孔孔径对低比转速离心泵空化性能的影响。
根据以上国内外研究现状,本文对10种不同开孔的水力模型进行空化定常数值模拟。基于减小叶片进口吸力面与压力面的压力差,减小低压区的面积,进而抑制空化在离心泵内部流场发生的原理。研究开孔面积和形状对离心泵空化性能的影响,其研究成果对解决水力机械空化空蚀问题有一定的理论基础和指导意义。

1 计算模型及网格划分

1.1 计算模型

离心泵模型参数:设计流量Q=75 m3/h,扬程H=92 m,转速n=2 950 r/min,叶轮外径D 2=290 mm,叶轮进口直径Dj =100 mm,叶轮出口宽度D 2=11 mm,轮毂直径Dh =50 mm,效率n=84%,比转数n s=52,叶片数Z=6。利用UG、CFTurbo,Pro.e等三维建模软件对离心泵建模,离心泵模型主要由进口段、叶轮、蜗壳、出口段四部分组成,离心泵的计算域如图1所示。离心泵叶轮越偏移中心,相对速度越大,且低压区主要集中在叶片头部,空化在该位置越容易发生,为研究离心泵进口位置穿孔孔型和面积对空化性能的影响,首先对该模型进行清水条件下的定常空化数值模拟,用以确认离心泵叶轮内空化的初生位置。经数值计算得到,空化初生位置距叶片头部约5 mm,本文均采用贯穿孔方式,开孔位置如图2所示。为进一步研究开孔面积、开孔形状对离心泵空化性能的影响,建立了面积为1/4 π mm2、9/4 π mm2、25/4 π mm2的圆形孔、方形孔、梅花形孔。其组合方案如表1所示,孔型展示如图3所示。
图1 离心泵三维模型

Fig.1 Three-dimensional model of centrifugal pump

Full size|PPT slide

图2 叶轮开孔位置

Fig.2 Impeller opening position

Full size|PPT slide

表1 孔型-开孔面积方案组合

Tab.1 Pass-hole area combination

开孔形状 开孔面积/mm2
1/4π 9/4π 25/4π
圆形孔 方案Ⅰ 方案Ⅱ 方案Ⅲ
方形孔 方案Ⅳ 方案Ⅴ 方案Ⅵ
梅花形孔 方案Ⅶ 方案Ⅷ 方案Ⅸ
图3 开孔形状

Fig.3 Open hole shape

Full size|PPT slide

1.2 计算网格

采用ICEM-CFD对三维模型进行网格划分,采用适应性较强的非结构性网格,取3组不同密度网格,并进行网格无关性检验,检验结果如表2所示。由于开孔尺寸较小,理论上选择更密的网格有利于提高精确度,由表2可得,方案2和方案3的扬程相对误差在0.01%内,因此可以忽略网格数对计算精度的影响。故最终选取了网格2作为计算网格,其中进口段网格数为 210 119,叶轮网格数1 401 664格,蜗壳网格数468 303格,出口段网格数181 873格,离心泵模型网格划分如图4所示。
表2 网格数无关性验证

Tab.2 Grid number independence verification

方案 网格数/万格 扬程/ m
1 158.67 91.152
2 226.20 92.557
3 352.45 92.562
图4 计算模型网格

Fig.4 Computational model grid

Full size|PPT slide

1.3 控制方法及边界条件

进口设置为总压进口,出口设置为质量流出口,在进口处设置液相体积分数为1,气相体积分数为0,参考压力设置为0 Pa,动静交界面采用固定转子法,流体介质为水,温度为25 ℃,水的饱和蒸气压设置为3 170 Pa,前后盖板和叶片设置为旋转界面,其余壁面均采用光滑无滑移边界。湍流方程以及对流项均采用高阶精度,收敛残差为10-4

1.4 数学模型

在数值计算中,选择RNG k-ε湍流模型,该模型能较为准确的模拟不同孔型对离心泵内部流场的影响,并在一定程度上考虑了各项异性效应和湍流旋涡特性,对预测复杂湍流的精度有较高的改善。
液体相连续方程:
tρfϕf¯ +xjρfϕf¯ Vfj¯=xjρfνtσfϕf¯xj
气体相连续方程:
tρaϕa¯+xjρaϕa¯ Vaj¯=xjρqνtσaϕa¯xj
液体相动量方程:
tρfϕf¯ Vfi¯+xjρfϕf¯ Vfi¯ Vfj¯=xjρfϕf¯νf+νtVfi¯xj+Sf
气体相动量方程:
tρaϕa¯ Vai¯+xjρaϕa¯ Vai¯ Vaj¯=xjρaϕa¯νa+νtVai¯xj+Sa
湍动能k方程:
tρfϕf¯k+xjρfϕf¯k Vfj¯=xjρfϕf¯νf+νtσkkxj+GkGb-ϕf¯ρfε-YM+Sk
湍动能扩散率 ε方程:
tρfϕf¯ε+ xjρfϕf¯ε Vfj¯ =xjρfϕf¯νf+νtσεεxj
+C1εGk+ϕf¯G3εGbεk-ρfϕf¯C2εε2k-Re+Sε
式中:下脚标fa分别表示液相和气相;下脚标ij为张量坐标;V为速度;ρ为材质密度;v为材质黏性系数;φ为体积分数; ¯为某变量的平均值;C 1 ε =1.42;C 1 ε =1.68;G 3 ε =tanh(Vx /Vz );Vx 是平行于重力矢量的速度分量;Vz 是垂直于重力矢量的速度分量;σk =σε =0.75;vt 为湍流涡流运动黏性系数;vt =Cμk 2/ε;Cμ =0.0845;Gk 是平均速度梯度产生的湍流动能;Gb 是浮力产生的湍流动能;YM 是由于过渡的扩散对耗散率的贡献;Re =Cμ ρ{[η 3(1-ηη 0]/(1+β η 3)}(ε 2/k);η 0=4.38;β=0.012;η=Sk /ε S=2SijSij Sij=0.5Vj/xi+Vi/xjSfSaSkSε,为定义的源项;φf +φa =1。
本文空化模型17选用均相流模型中的质量输送空化模型:
ρvαvt+(ρvαvui)xi=m˙+-m˙-
m˙+=Cvapo3rg(1-αv)Rb23max (pv-p,0)ρl
m˙-=Ccond3αvpv  Rb23max (pv-p,0)ρl
式中:αv 为空泡体积分数;ρv 为气体密度,kg/m3 m˙+ m˙-为质量蒸发速率和质量凝结速率;rg 为单位液体中所含气核体积分数,5×10-4pv 为该温度下水的饱和蒸汽压,Pa;p为空泡周围液体的压力; Cvapo为蒸发项经验系数,50; Ccond为凝结经验系数,0.01。

2 计算结果分析

2.1 原模型空化特性曲线

通过数值模拟定常计算,利用NPSHa计算公式18,可以得到不同进口压力下的NPSHa值和扬程,在设计工况下,预测得到了未穿孔下的H-NPSHa曲线,如图5所示。经数值计算得到,未穿孔模型空化初生点NPSHa=3.758 m。
NPSHa=pin-pvρ g
式中:pin 为离心泵的进口静压力;pv 为该温度下水的饱和蒸汽压,Pa,这里取pv =3 170 Pa;ρ为流体密度,ρ=1 000 kg/m3g为重力加速度,g=9.81 m/s2
图5 空化性能曲线

Fig.5 Cavitation performance curve

Full size|PPT slide

图5空化性能曲线可知,当空化余量较大时,离心泵内无空化产生,泵内能量特性不受影响,离心泵扬程基本维持不变。随着进口压力的减小,空化余量逐步降低,当进口压力为40 kPa、空化余量为3.758 m时,得到离心泵空化初生点。开孔在空化不同阶段对叶轮的空化性能影响不同,在空化初生点时穿孔效果明显优于其他空化点。故本文以空化初生点为基础研究孔型、孔面积对离心泵空化性能的影响。

2.2 特性曲线分析

NPSH=3.758 m时,对各组合方案进行定常计算,分析开孔面积、开孔形状对离心泵扬程和效率的影响,结果如图6所示。
图6 开孔对离心泵外特性影响

Fig.6 Influence of opening hole on external characteristics of centrifugal pump

Full size|PPT slide

图6可得,随着开孔面积的增大,各孔型对离心泵扬程、效率影响的趋势基本一致。当开孔面积为1/4 π mm2,孔型为方形孔时对离心泵扬程提高最为显著,增加了8.6 m,达到峰值。方案Ⅱ、方案Ⅲ、方案Ⅳ、方案Ⅸ对离心泵外特性均有提高。开孔面积对梅花形孔离心泵影响最为敏感。

2.3 叶轮进口位置压力分布分析

在离心泵设计工况下进行定常数值计算,将进口压力设置为40 kPa,进口截面为叶片穿孔中心所在的平面,为方便比较穿孔前后叶轮进口位置压力的变化情况,将压力云图压力变化范围设置为0~4.7×105 Pa。叶轮进口位置压力分布如图7~图9所示。
图7 圆形孔穿孔截面压力分布

Fig.7 Pressure distribution of circular hole perforation section

Full size|PPT slide

图8 方形孔穿孔截面压力分布

Fig.8 Pressure distribution of square hole perforation section

Full size|PPT slide

图9 梅花形孔穿孔截面压力分布

Fig.9 Pressure distribution of the cross section of the quinted-shaped hole

Full size|PPT slide

图7图9可见,穿孔对叶轮进口位置的压力分布有较大的影响,当开孔面积S=1/4 π mm2,方形孔对叶轮进口位置低压区改善明显,大大地减少低压区面积分布并将叶片进口低压区区域截断,说明该开孔模型能有效地抑制空化。然而梅花孔形却加剧了低压区的分布,致使模型的效率、扬程下降了约60%,表明叶轮内大部分区域发生了空化,说明当开孔面积过小(单孔的直径0.408 mm)时,流体在流经过流断面时,流体并未通过孔,反而造成叶轮内的能量损失,加剧低压区的面积分布,诱导叶轮空化发生。开孔面积S=9/4 π mm2时,圆形孔对减少叶轮进口低压区区域分布最为明显,同时也截断了低压区的连续分布。开孔面积S=25/4 π mm2时,梅花孔对叶片进口低压分布区改善效果最好,能有效减小低压区面积,对抑制空化有明显的作用,表明此时开孔造成的局部阻力较小,叶轮内能量损耗较小。

2.4 叶轮中截面气相体积分析

通过研究叶栅切片上的空泡分布进一步分析设计工况下的叶轮流场,通过观察原型叶轮与穿孔叶轮上的空泡分布对比,分析不同穿孔面积、不同形状对叶轮的影响。叶轮内设计流量下的叶栅切片空泡云图如图10-图12所示。
图10 圆形穿孔叶轮汽泡分布

Fig.10 Circular perforated impeller bubble distribution

Full size|PPT slide

图11 方形穿孔叶轮汽泡分布

Fig.11 Square perforated impeller bubble distribution

Full size|PPT slide

图12 梅花形穿孔叶轮汽泡分布

Fig.12 The bubble distribution of quinted-shaped perforated impeller

Full size|PPT slide

图10~图12可以看出,当穿孔形状为圆形孔时,随着穿孔面积的不断增大,叶轮内部气泡分布区域呈先减小后增大的趋势,在S=1/4π mm2时,叶轮内空泡分布区域最小。当孔型为方形孔时,叶轮内部分布区域在S=1/4π mm2时最小,且此时空泡体积变小,随着开孔面积的不断增大,叶轮空泡区域分布变大,这是由于穿孔能量损失过大所造成的。当孔型为梅花孔时,在开孔面积S=25/4π mm2时,叶轮内气泡分布区域大大减少,说明能有效地抑制叶轮内空化情况,在开孔面积为1/4π mm2和9/4π mm2时,气泡分布区域较原模型反而增加,这是由于叶片工作面的流体不能通过狭小的圆孔进入到叶片的背面,反而造成了局部能量损失,加剧了叶轮内部空化。

2.5 叶轮流道湍动能分析

图13图14所示为设计流量下进口压力40 kPa时穿孔前后叶轮流道内湍动能分布图,分别选取本文中抗空化性能最好的各孔型进行分析,其中图13(a)、(b)、(c)、(d)分别是开孔面积为0 mm2、9/4π mm2、1/4π mm2、25/4π mm2的原模型、圆形孔、方形孔、梅花孔叶轮中间截面上湍动能分布图。其中图14(a)、(b)、(c)、(d)分别为开孔面积为0π mm2、9/4π mm2、1/4π mm2、25/4π mm2的原模型、圆形孔、方形孔、梅花孔叶轮穿孔位置处湍动能分布图。
图13 叶轮中间截面湍动能分布

Fig.13 The turbulent kinetic energy distribution in the middle section of the impeller

Full size|PPT slide

图14 叶轮穿孔截面湍动能分布

Fig.14 Turbulent kinetic energy distribution in perforated section of impeller

Full size|PPT slide

湍动能主要反映了叶轮内部能量的变化情况。通过对比原模型与穿孔模型穿孔位置处湍动能分布和叶轮中间截面上湍动能分布图。原模型叶轮总的湍动能均大于穿孔叶轮,且梅花形孔模型内叶轮总的湍动能最小。产生这种现象是由于穿孔面积为25/4π mm2时,孔型为梅花孔型时,开孔带来的能量损失最小,且对抑制空化有明显的效果,降低了空化造成的能量损失,致使湍动能最小。

3 结 论

(1)在气泡的初生位置对叶片进行穿孔后,不同穿孔面积、形状对离心泵的扬程和效率有较大的影响。当开孔形状为方形孔,开孔面积为1/4π mm2时对离心泵外特性改善最佳,扬程提高了9.8%。
(2)穿孔会造成能量损失,不同孔型的最佳开孔面积不同,开孔面积对梅花形孔影响最为敏感。当穿孔面积、形状能有效抑制叶轮空化时,叶轮流道内的低压区都会被截断,低压区呈不连续分布状态。
(3)综合本文的十种穿孔模型,穿孔形状为方形孔,开孔面积为1/4π mm2和开孔面积为25/4π mm2的穿孔模型对抑制叶轮内空化效果最佳。且叶轮内湍动能低于原型叶片,减小了离心泵运行时的振动。

参考文献

1
关醒凡. 现代泵技术手册[M]. 北京:中国宇航出版社,1995.
2
黄继汤. 空化与空蚀的原理及应用[M]. 北京:清华大学出版社,1991.
3
袁丹青,陈向阳,白滨,等. 水力机械空化空蚀问题的研究进展[J]. 排灌机械200927(4):269-272.
4
PLESSET M S PROSPERETTI A. bubble dynamics and cavitation[J]. Annual Review of Fluid Mechanics19779(1):145-185.
5
RAFAEL, Campos-Amezcua, Farid, et al. Numerical analysis of unsteady cavitating flow in an axial inducer[J]. Applied Thermal Engineering201575:1 302-1 310.
6
MURAKAMI M MINEMURA K. Effects of entrained air on the performance of a centrifugal pump: 1st report, performance and flow conditions[J]. Bulletin of JSME197417(110):1 047-1 055.
7
JI B LUO X W AMDT R E A, et al. Large Eddy simulation and theoretical investigations of the transient cavitating vortical flow structure around a NACA66 hydrofoil[J]. International Journal of Multiphase Flow201568(68):121-134.
8
SEKOGUCHI K TAKADA S KANEMORI Y. Study of air-water two-phase centrifugal pump by means of electric resistivity probe technique for void fraction measurement: 1 st report, measurement of void fraction distribution in a radial flow impeller[J]. Bulletin of JSME198427(227):931-938.
9
宋文武,石建伟,魏立超,等. 高速离心泵回流漩涡及空化特性[J]. 机械工程学报202056(4):95-103.
10
王维军,王洋,刘瑞华,等. 离心泵空化流动数值计算[J]. 农业机械学报201445(3):37-44.
11
刘小兵,胡全友,史广泰,等. 混输泵小流量工况的空化特性[J]. 排灌机械工程学报201836(1):15-20.
12
罗旭,宋文武,万伦,等. 高速离心泵空化特性研究[J]. 水力发电201945(11):74-78,94.
13
李金琼,宋文武,万丽佳. 粗糙度对离心泵进口回流非定常特性影响的研究[J]. 热能动力工程202035(2):70-77.
14
胡赞熬,王俊雄,祝宝山,等. 离心泵叶轮穿孔对空化性能的影响[J]. 热能动力工程201833(10):44-51,106.
15
朱晓东,刘小兵,田长安. 翼型表面开孔对其空化性能的影响[J]. 中国农村水利水电2020(2):83-87.
16
赵伟国,潘绪伟,宋启策,等. 叶片进口边穿孔对离心泵空化性能的影响[J]. 排灌机械工程学报201937(6):461-468.
17
曹玉良,贺国,明廷锋,等. 水泵空化数值模拟研究进展[J]. 武汉理工大学学报(交通科学与工程版)201640(1):55-59,65.
18
王勇. 离心泵空化理论与技术[M]. 北京:科学出版社,2018.
PDF(3135 KB)

访问

引用

详细情况

段落导航
相关文章

/