張 玲 黃鵬杰 鄒 永 修棟波

(東北電力大學(xué)能源與動力工程學(xué)院)

多級離心泵因結(jié)構(gòu)先進(jìn)合理、可靠性好的優(yōu)點被廣泛應(yīng)用于發(fā)電和石油化工領(lǐng)域[1]。當(dāng)離心泵流場內(nèi)某處壓力低于所輸送水的汽化壓力時，就會發(fā)生汽蝕現(xiàn)象，從而對泵產(chǎn)生破壞，影響其運行效率[2]。隨著計算流體動力學(xué)(CFD)的不斷發(fā)展及流體仿真技術(shù)的進(jìn)步，通過數(shù)值模擬對離心泵內(nèi)部流場進(jìn)行研究已被廣泛采用[3～5]。劉宜等通過CFD技術(shù)對離心泵全流道三維空化流動進(jìn)行了定常和非定常數(shù)值模擬，結(jié)果與實驗結(jié)果相一致[6]。Steinmann A 等分別采用Rayleigh- Plesset空化模型和試驗方法分析了蝸殼泵在不同汽蝕余量下的汽泡流動分布以及汽泡潰滅的過程，兩種方法得到的結(jié)果基本吻合[7]。楊孫圣等對某離心泵進(jìn)行了定常數(shù)值模擬，獲得了泵發(fā)生汽蝕時的流場汽液兩相分布規(guī)律[8]。

筆者針對某熱電廠變頻給水泵的小流量工況進(jìn)行研究。當(dāng)機組的負(fù)荷降低時，幾臺給水泵并聯(lián)運行時的總體出力超過了鍋爐的需求，造成很大的浪費。為降低電廠的廠用電率，將變頻給水泵與工頻給水泵并聯(lián)降壓運行，當(dāng)調(diào)速范圍過大頻率過低時，母管壓力降低，致使工頻泵流量過大，則變頻泵流量會相應(yīng)的減少，水泵在小流量高轉(zhuǎn)速運行時，容易發(fā)生振動，產(chǎn)生汽蝕，出力有所下降。筆者采用Fluent軟件對該變頻泵(DG85- 80×10型多級離心泵)在小流量工況時的速度及壓力分布情況進(jìn)行分析，紊流模型選用標(biāo)準(zhǔn)模型，確定易發(fā)生汽蝕的位置并提出相應(yīng)的改進(jìn)措施，為水泵的優(yōu)化設(shè)計和進(jìn)一步提高泵的效率提供依據(jù)。

1 模型的建立

1.1幾何模型

DG85- 80×10型給水泵的性能參數(shù)和首級葉輪幾何參數(shù)如下：

流量 85m3/h

級數(shù) 10級

總揚程 800m

轉(zhuǎn)速 2 950r/min

效率 65%

進(jìn)口直徑 145mm

出口直徑 252mm

出口寬度 21mm

葉片數(shù) 6

圖1為多級泵首級葉輪與導(dǎo)葉流域部分實體模型。

圖1 葉輪與導(dǎo)葉流域部分幾何模型

1.2網(wǎng)格劃分

利用Gambit軟件對模型進(jìn)行網(wǎng)格劃分，水泵的葉輪及導(dǎo)葉部分采用非結(jié)構(gòu)化網(wǎng)格，在葉輪和正、反導(dǎo)葉進(jìn)出口部分對網(wǎng)格進(jìn)行加密。如圖2所示，首級葉輪部分劃分的網(wǎng)格數(shù)為751 690，導(dǎo)葉部分劃分的網(wǎng)格數(shù)為1 412 190。

圖2 葉輪與導(dǎo)葉部分網(wǎng)格

1.3控制方程

連續(xù)性方程：

(1)

其中，ρ是水的密度，u、v、ω分別是流體在x、y、z3個方向上的速度分量。

動量方程：

(2)

其中，ui、uj(i,j=1,2,3)為流體的各速度分量；xi、xj(i,j=1,2,3)是各坐標(biāo)的分量，ρ是流體的密度，p是壓力，ρFi為作用在流體微團(tuán)上的質(zhì)量力，μ粘性系數(shù)。

文中采用標(biāo)準(zhǔn)k-ε湍流模型。k和ε對應(yīng)的輸運方程分別為：

(3)

(4)

其中p為湍動能生成項：

1.4邊界條件

采用Fluent 分離隱式求解器進(jìn)行穩(wěn)態(tài)計算，壓力和速度的耦合采用SIMPLEC算法，對動量項、湍動能項和湍動能耗散率項采用二階迎風(fēng)差分格式，壓力項采用標(biāo)準(zhǔn)差分格式[9]。進(jìn)口邊界條件為速度進(jìn)口。出口邊界為壓力出口，壁面條件為無滑移固壁條件，使用標(biāo)準(zhǔn)壁面函數(shù)確定固壁附近的流動[10]。

2 計算結(jié)果與分析

2.1性能預(yù)測

揚程計算公式為：

(5)

式中pout——出口總壓，Pa；

pin——進(jìn)口總壓，Pa；

Δz——葉輪進(jìn)口與導(dǎo)葉出口在垂直方向的距離，m。

效率計算公式為：

(6)

式中M——泵軸的輸出扭矩，N·m；

n——轉(zhuǎn)速，r/min；

qv——體積流量，m3/s。

軸功率計算公式為：

(7)

利用以上公式，對DG85- 80×10型給水泵首級的能量特性曲線進(jìn)行預(yù)測。其預(yù)測值和試驗值比較結(jié)果如圖3所示。

a. 揚程-流量

b.功率-流量

c.效率-流量

由圖3可以看出:預(yù)測的揚程、軸功率和效率隨流量的變化趨勢與實際性能曲線基本一致，說明采用該紊流模型對多級泵進(jìn)行數(shù)值模擬并預(yù)測其外特性可行。同時由于計算機資源及時間的限制，網(wǎng)格數(shù)相對比較少，并簡化了幾何模型，對計算準(zhǔn)確性有一定的影響，導(dǎo)致預(yù)測值與試驗值之間也存在著一定的誤差。

2.2額定工況下泵內(nèi)部流動特性

圖4是在額定工況下離心泵首級葉輪與導(dǎo)葉中間截面內(nèi)部的壓力分布圖?？梢钥闯觯陬~定流量工況下，葉輪與導(dǎo)葉內(nèi)部的流動壓力變化均勻平滑。由于葉輪與導(dǎo)葉相對位置差異，所以壓力分布并不完全對稱。由于正導(dǎo)葉入口端對流體存在阻擋干擾，使附近流體產(chǎn)生擾動，從而在葉輪與導(dǎo)葉之間存在局部的低壓區(qū)，該低壓區(qū)葉輪與導(dǎo)葉內(nèi)部工質(zhì)流動產(chǎn)生一定的干擾，同時容易產(chǎn)生間隙汽蝕。

圖4 葉輪與導(dǎo)葉中間截面壓力分布

圖5為額定工況下葉輪與導(dǎo)葉中間截面速度分布。流體流入導(dǎo)葉后，在導(dǎo)葉各流道前段靠下端都產(chǎn)生了小范圍的渦流。這是由于導(dǎo)葉進(jìn)口角度不佳，進(jìn)入導(dǎo)葉流道內(nèi)的流體因離心力作用，使在導(dǎo)葉流道前段上端流速大于下端流速，形成小范圍內(nèi)脫流而產(chǎn)生旋渦，使此處流速降低較快，但由于導(dǎo)葉的擴(kuò)壓作用而不致發(fā)生汽蝕。在額定流量下，各流道雖有小部分渦流存在，但整體分布較均勻，流動比較順暢。

圖5 葉輪與導(dǎo)葉中間截面速度分布

圖6為額定工況導(dǎo)葉入口端壓力與流線分布圖。可以看出，工質(zhì)在導(dǎo)葉入口端擾流較平滑，但由于來流速度方向與導(dǎo)葉切線方向還存在少量偏差，導(dǎo)致在導(dǎo)葉入口端偏向葉輪方向也產(chǎn)生了小面積低壓區(qū)，從而容易產(chǎn)生局部汽蝕。

圖6 導(dǎo)葉入口端壓力與流線圖

圖7為額定工況反導(dǎo)葉內(nèi)壓力與速度分布圖。當(dāng)工質(zhì)流向反導(dǎo)葉時流速較高，葉片與工質(zhì)流動方向角度過大，形成較大的沖角，引起較大的沖擊損失。在個別反導(dǎo)葉背面產(chǎn)生小范圍邊界層分離現(xiàn)象，雖整體上能夠較好的流通過渡，但同時也存在較大的流動損失。

圖7 反導(dǎo)葉內(nèi)壓力與速度分布

2.3小流量工況下泵內(nèi)部流動特性

分別對90%、80%、70%、60%、50%、40%、30%的額定流量7個工況進(jìn)行模擬后，發(fā)現(xiàn)流量較小時離心泵內(nèi)部流動狀況趨于不穩(wěn)定，筆者取比較有代表性的40%額定流量工況進(jìn)行分析。圖8為小流量下葉輪與導(dǎo)葉速度分布圖?？梢娫谛×髁抗r下，導(dǎo)葉流道內(nèi)的流動狀態(tài)變差，由于泵內(nèi)流動分布不均勻的原因，在個別流道內(nèi)出現(xiàn)較大的旋渦流，并且存在反方向的速度分量，使流動阻力增大，流道阻塞，水頭損失增大，使工況惡化。對此應(yīng)在設(shè)計上延長動能向壓力能過渡轉(zhuǎn)換過程，減小工質(zhì)進(jìn)入的沖角，改變導(dǎo)葉型線，使工質(zhì)平滑連續(xù)的流過導(dǎo)葉，由此改善工質(zhì)的流動狀態(tài)，提高泵的穩(wěn)定性及效率。

圖8 小流量下葉輪與導(dǎo)葉速度分布

圖9為小流量下導(dǎo)葉入口端壓力分布圖，可見在小流量工況下，導(dǎo)葉頭部偏向葉輪方向的小塊低壓區(qū)同樣存在，該處容易產(chǎn)生局部汽蝕，同時會對導(dǎo)葉頭部形成一定程度的水力沖擊，產(chǎn)生振動或噪聲危害，影響泵的使用壽命。

圖9 小流量下導(dǎo)葉入口端壓力分布

圖10為葉輪末端壓力分布圖。在葉輪末端邊緣處存在小范圍低壓區(qū)，這是由葉片尾部流動邊界層分離所造成。為防止發(fā)生汽蝕磨損，可將葉片幾何型線進(jìn)一步優(yōu)化。

圖10 葉輪末端壓力分布

圖11為小流量工況下反導(dǎo)葉內(nèi)壓力及速度分布圖。工質(zhì)進(jìn)入反導(dǎo)葉后，由于流動方向發(fā)生急劇改變，其流動出現(xiàn)嚴(yán)重不均勻現(xiàn)象，部分流道內(nèi)出現(xiàn)較大的渦流，使阻力增加，嚴(yán)重阻礙了工質(zhì)的流通。因此對比于額定工況，小流量工況下?lián)p失增大，可在設(shè)計上改變反導(dǎo)葉進(jìn)口角度。

圖11 小流量下反導(dǎo)葉內(nèi)壓力與速度分布

3 結(jié)論

3.1預(yù)測的揚程、軸功率和效率隨流量的性能曲線與實際性能曲線基本一致，證明了采用標(biāo)準(zhǔn)模型對多級泵進(jìn)行數(shù)值模擬研究的可靠性。

3.2在額定工況下，各導(dǎo)葉內(nèi)部流道雖有小部分渦流存在，反導(dǎo)葉入口有水力沖擊現(xiàn)象，但總體來看，工質(zhì)在設(shè)計工況下流動趨于平穩(wěn)，流動損失較小。

3.3在小流量工況下，導(dǎo)葉流道內(nèi)的流動狀態(tài)惡化，個別流道內(nèi)出現(xiàn)較大的旋渦流，導(dǎo)葉頭部偏向葉輪方向一直存在小塊低壓區(qū)，容易產(chǎn)生局部汽蝕及振動傷害。 葉輪末端邊緣處由于邊界層分離造成小范圍低壓區(qū)。在反導(dǎo)葉內(nèi)部分流道出現(xiàn)較大渦流，嚴(yán)重阻礙了工質(zhì)的流通。對比于額定工況，小流量工況下?lián)p失增大，需在設(shè)計上加以優(yōu)化。

[1] 姚睿, 李志鵬, 秦武, 等. 火電廠用泵的技術(shù)研究及發(fā)展[J]. 電站系統(tǒng)工程, 2008, 24(5): 5～8.

[2] 劉云, 盛頌恩, 白少先. 基于Labview的機械密封實驗數(shù)據(jù)采集[J]. 流體機械, 2010, 38(6): 48～50.

[3] 張學(xué)靜, 楊軍虎. 多級泵內(nèi)部流場的三維數(shù)值模擬及性能預(yù)測[J]. 流體機械, 2011, 39(8): 24～28.

[4] 陳芳芳, 李志鵬, 王昌生. 首級葉輪結(jié)構(gòu)對多級泵效率影響的數(shù)值模擬[J]. 電站系統(tǒng)工程, 2013, 29(3): 9～11.

[5] Plutecki J, Skrzypacz J. CFD Simulations of 3D Flow in a Pump Stator with a Spherical Surface [J]. World Pumps, 2003, 443: 28～31.

[6] 劉宜, 趙希楓, 蔡玲春，等. 離心泵內(nèi)部流場的三維空化湍流的數(shù)值計算[J]．蘭州理工大學(xué)學(xué)報, 2009, 35(1): 50～53.

[7] Steinmann A, Wurm H, Otto A. Numerical and Experimental Investigations of the Unsteady Cavitating Flow in a Vortex Pump[J]. Journal of Hydrodynamics, Ser. B, 2010, 22 (5): 324～329.

[8] 楊孫圣, 孔繁余, 周水清，等. 離心泵氣蝕性能的數(shù)值計算與分析[J]. 華中科技大學(xué)學(xué)報(自然科學(xué)版), 2010,38 (10): 93～95.

[9] 李艷凱, 姜樹李, 郭晨海. 離心水泵內(nèi)三維數(shù)值模擬[J]. 小型內(nèi)燃機與摩托車, 2008, 37(4): 41～43.

[10] 邵春雷, 顧伯勤, 陳曄. 離心泵內(nèi)部非定常壓力場的數(shù)值研[J]. 農(nóng)業(yè)工程學(xué)報, 2009, 25(1): 75～80.