胡全友，劉小兵，趙 琴，安滿意，黃長(zhǎng)久

(西華大學(xué)能源與動(dòng)力工程學(xué)院，成都 610039)

空化是液體特有的物理現(xiàn)象，當(dāng)液體中局部壓力降低到一定壓力值(相應(yīng)工作溫度所對(duì)應(yīng)的飽和蒸汽壓力)時(shí)，液體汽化而引發(fā)的微氣泡爆發(fā)性的生長(zhǎng)現(xiàn)象[1]。泵在空化條件下長(zhǎng)期工作也會(huì)降低它的運(yùn)轉(zhuǎn)穩(wěn)定性和部件壽命。因此，研究流道內(nèi)空化發(fā)生、成長(zhǎng)和破滅過(guò)程對(duì)性能的影響是水力機(jī)械的重要研究課題。隨著計(jì)算流體動(dòng)力學(xué)(CFD)軟件的快速發(fā)展，計(jì)算機(jī)性能的逐步提高，數(shù)值模擬求解流體運(yùn)動(dòng)問(wèn)題已成為水力機(jī)械內(nèi)部流場(chǎng)計(jì)算的重要方法。低比轉(zhuǎn)數(shù)離心泵在各領(lǐng)域的廣泛應(yīng)用，使得其效率低，流量-揚(yáng)程特性曲線易出現(xiàn)駝峰，大流量工況電機(jī)易過(guò)載等問(wèn)題逐漸引起人們的關(guān)注。實(shí)踐表明，分流葉片的設(shè)計(jì)是改善離心泵內(nèi)部流場(chǎng)和壓力場(chǎng)的重要手段。袁建平等[2]用正交試驗(yàn)研究分析了分流葉片位置葉片數(shù)、周向位置和進(jìn)口直徑等幾個(gè)參數(shù)對(duì)離心泵性能的影響。袁壽其等[3]用正交試驗(yàn)比較確定分流葉片主要參數(shù)對(duì)水泵性能影響的主次順序。M Dular等[4]對(duì)翼型空化流動(dòng)進(jìn)行數(shù)值模擬，并與PIV試驗(yàn)結(jié)果相互對(duì)比，結(jié)果表明數(shù)值模擬能較精確反映試驗(yàn)研究結(jié)果。羅先武等[5]對(duì)低比轉(zhuǎn)數(shù)離心泵進(jìn)行空化性能實(shí)驗(yàn)，并在設(shè)計(jì)工況點(diǎn)附近，基于標(biāo)準(zhǔn)κ-ε模型和VOF空化模型對(duì)葉輪進(jìn)口幾何參數(shù)對(duì)離心泵空化性能影響進(jìn)行模擬，表明保證葉輪進(jìn)口的流動(dòng)均勻性是改善空化性能的必要措施。

由于離心泵空化過(guò)程中氣液兩相流動(dòng)的復(fù)雜性，目前雖然關(guān)于分流葉片抗空化的特點(diǎn)得到共識(shí)，但其相關(guān)研究卻未引起重視，大部分還是針對(duì)分流葉片幾何結(jié)構(gòu)對(duì)離心泵流場(chǎng)和性能的影響[6,7]。因此在離心泵內(nèi)部空化流場(chǎng)的研究中分流葉片減輕空化對(duì)水泵運(yùn)行的不利影響具有十分重要的意義。對(duì)在多個(gè)工況下對(duì)離心泵的空化流場(chǎng)進(jìn)行數(shù)值研究并預(yù)測(cè)葉片數(shù)和分流葉片對(duì)空化性能的影響是本文的研究目的。

1 幾何模型與網(wǎng)格建立

1.1 幾何模型設(shè)計(jì)

實(shí)驗(yàn)泵為比轉(zhuǎn)數(shù)ns=58 mm3/min的離心泵，流量Q=75 m3/h，揚(yáng)程H=86 m，轉(zhuǎn)速n=2 950 r/min，效率η=77%。葉片數(shù)Z=6，葉輪出口直徑D2=260 mm，葉輪出口寬度b2=12 mm，葉輪出口安放角β2=29°。

采用建模軟件對(duì)轉(zhuǎn)輪和蝸殼等流道的計(jì)算區(qū)域三維造型。在ICEM CFD中劃分網(wǎng)格，轉(zhuǎn)輪和蝸殼采用非結(jié)構(gòu)化四面體網(wǎng)格，進(jìn)出口延長(zhǎng)段采用結(jié)構(gòu)化六面體網(wǎng)格。對(duì)網(wǎng)格質(zhì)量進(jìn)行無(wú)關(guān)性檢查，當(dāng)網(wǎng)格質(zhì)量對(duì)揚(yáng)程的影響低于0.5%，則認(rèn)為網(wǎng)格質(zhì)量滿足計(jì)算要求。最終全流道網(wǎng)格總數(shù)為1.4×106，葉輪網(wǎng)格為5.8×105，如圖1所示。

圖1 計(jì)算區(qū)域的三維建模和網(wǎng)格生成圖Fig.1 3D model and grid of computing domain

1.2 葉輪改型設(shè)計(jì)

為研究分流葉片對(duì)泵空化性能的影響，對(duì)葉輪葉片進(jìn)行重新設(shè)計(jì)，設(shè)計(jì)出3長(zhǎng)葉片，3長(zhǎng)3短分流葉片兩種葉輪方案。由于葉片數(shù)、周向位置、進(jìn)口直徑和偏轉(zhuǎn)角等幾個(gè)參數(shù)是決定分流葉片對(duì)泵性能影響的最重要因素，因此我們只考慮這幾個(gè)主要參數(shù)的設(shè)計(jì)。具體設(shè)計(jì)參數(shù)參考袁建平等正交實(shí)驗(yàn)研究出的分流葉片最佳組合參數(shù)：周向位置為0.6θ，偏轉(zhuǎn)角為-5°～-10°。進(jìn)口直徑為2/3D2，分流葉片出口角與長(zhǎng)葉片出口角相等。分流葉片偏置設(shè)計(jì)[2]如圖2所示。

圖2 分流葉片偏置設(shè)計(jì)Fig.2 Design sketch of bias splitter blades

2 數(shù)值模擬

2.1 控制方程

采用時(shí)均流動(dòng)的連續(xù)性方程、Navier-Stokes方程及基于Yakhot和Orzag提出的RNGκ-ε模型，即：

連續(xù)性方程：

(1)

動(dòng)量方程：

(2)

混合密度：

ρm=ρvαv+ρl(1-αv)

(3)

式中：ui為雷諾時(shí)均速度；μeff為等效黏性系數(shù)，μeff=μ+μt；p*為包括湍動(dòng)能κ靜壓力；ρl、ρv分別為液體和氣體的密度。

2.2 Zwart-Gerber-Belamari空化模型

兩相中相間質(zhì)量傳輸率[8]：

(4)

空泡體積變化率：

(5)

式中：F為經(jīng)驗(yàn)系數(shù)；r1汽核初始體積分?jǐn)?shù)；rg空泡體積分?jǐn)?shù)；R為氣泡半徑；pv為蒸發(fā)壓力；p為空泡周圍液體的壓力；V為空泡體積。

2.3 計(jì)算方法

采用RNGκ-ε湍流模型來(lái)考慮湍流的影響，均質(zhì)多相模型和Rayleigh-Plesset兩相空化模型來(lái)考慮空泡的初生、發(fā)展和破滅過(guò)程。設(shè)定流體的飽和蒸汽壓力為3 170 Pa，空泡的平均直徑為2×10-6m，參考?jí)毫υO(shè)置為0。轉(zhuǎn)輪進(jìn)出口邊界通過(guò)設(shè)置Frozen Rotor交界面實(shí)現(xiàn)動(dòng)靜耦合。

2.4 邊界條件

設(shè)置總壓進(jìn)口條件，方向垂直于進(jìn)口邊界。進(jìn)口邊界κ和ε由經(jīng)驗(yàn)公式確定，即：

κ=0.005u2

(6)

ε=C3/4κ3/2/l

(7)

式中：u為進(jìn)口平均速度；C=0.09；l=0.07Din；Din為進(jìn)口等效管徑。

出口恒定質(zhì)量流量出口，質(zhì)量流量出口能保證水泵運(yùn)行時(shí)的流量不變。壁面采用無(wú)滑移邊界條件。逐步設(shè)置較低的進(jìn)口總壓以達(dá)到泵內(nèi)空化的目的。空化計(jì)算時(shí)，設(shè)定進(jìn)口邊界中氣相體積分?jǐn)?shù)為0，液相體積分?jǐn)?shù)為1，收斂精度為10-5。并將單相計(jì)算結(jié)果作為空化流動(dòng)計(jì)算的初始值以節(jié)省空化計(jì)算時(shí)間。

3 離心泵無(wú)空化特性曲線

在不同工況下，模擬計(jì)算出 長(zhǎng)葉片方案在未發(fā)生空化時(shí)的外特性曲線與實(shí)驗(yàn)結(jié)果能較好吻合。說(shuō)明了數(shù)值模擬結(jié)果具有一定的精確度，如圖3所示。

圖3 無(wú)空化時(shí)能量特性預(yù)測(cè)值與試驗(yàn)值的對(duì)比Fig.3 Contrast of Energy characteristics predicted and experiment without cavitation

圖3表明預(yù)測(cè)揚(yáng)程和效率隨著流量變化的曲線與試驗(yàn)值變化趨勢(shì)相一致，預(yù)測(cè)揚(yáng)程和效率相對(duì)于試驗(yàn)值偏高，這是因?yàn)樵谟?jì)算離心泵的效率過(guò)程中只考慮了水力效率而忽略了容積效率和機(jī)械效率。預(yù)測(cè)揚(yáng)程偏高是因?yàn)榱黧w流過(guò)部件表面由于未考慮粗糙度而導(dǎo)致水里損失被忽略。

4 離心泵空化特性

4.1 不同流量下空化曲線

基于未發(fā)生空化數(shù)值計(jì)算結(jié)果上，增加空化流動(dòng)模型，逐步設(shè)置較低的進(jìn)口總壓。當(dāng)進(jìn)口總壓較大時(shí)，離心泵內(nèi)部流場(chǎng)沒(méi)有發(fā)生空化，泵的揚(yáng)程基本保持不變，其運(yùn)行特性也沒(méi)有受到影響。當(dāng)進(jìn)口總壓逐步減小時(shí)，泵內(nèi)空化程度由局部向葉輪流道擴(kuò)展，氣泡的大量產(chǎn)生使流體連續(xù)性遭到破壞，性能下降，這與理論分析基本相吻合。6長(zhǎng)葉片方案在不同流量下空化特性曲線如圖4所示。

圖4 不同流量下空化特性曲線Fig.4 Cavitation characteristics under different mass flow

離心泵有效汽蝕余量NPSHa：

(8)

式中：pv為泵工作溫度下流體的汽化壓力。

通常我們將泵揚(yáng)程下降3%時(shí)所對(duì)應(yīng)的有效汽蝕余量NPSHa作為泵的臨界汽蝕余量NPSHc[6]。離心泵在3個(gè)不同工況點(diǎn)下的NPSHc分別為1.82、3.2、4.8 m。 隨流量的變化趨勢(shì)與理論分析一致，表明模擬計(jì)算對(duì)預(yù)測(cè)空化性能具有一定的精準(zhǔn)度。

4.2 不同葉輪空化曲線

在有限葉片數(shù)葉輪流道中，除了有一個(gè)均勻的相對(duì)速度外，還有一個(gè)相對(duì)的軸向漩渦運(yùn)動(dòng)[9]。葉片數(shù)越多，流體的流動(dòng)越接近無(wú)限葉片數(shù)。葉片數(shù)越少，葉輪流道內(nèi)相對(duì)軸面漩渦運(yùn)動(dòng)越劇烈，在葉片工作面附近，兩種相對(duì)運(yùn)動(dòng)合成的相對(duì)速度較小，而葉片背面相對(duì)運(yùn)動(dòng)合成的相對(duì)速度較大。由于相對(duì)流速的不同，葉片兩面出現(xiàn)壓力差。當(dāng)葉片數(shù)發(fā)生變化時(shí)，可近似認(rèn)為其圓周速度保持不變。在設(shè)計(jì)流量相同的三組葉輪設(shè)計(jì)方案中，由于排擠系數(shù)對(duì)軸面速度影響較小，可認(rèn)為軸面速度相同。因此葉片數(shù)的減少，葉輪內(nèi)流體的絕對(duì)速度降低，造成泵壓頭的衰減，而揚(yáng)程對(duì)NPSHc的影響較敏感。圖5為不同葉片數(shù)下空化性能預(yù)測(cè)曲線。

圖5 不同葉輪方案下空化曲線Fig.5 Cavitation characteristics under different programs

從圖5可以看出三組葉輪方案的NPSHc分別為3.2、4.18、5.8 m，泵揚(yáng)程變化曲線隨著NPSHc的變大而升高，而曲線斜率逐漸降低，當(dāng)超過(guò)NPSHc時(shí)，曲線斜率基本保持不變。從圖5可以看出3種葉輪方案在無(wú)空化時(shí)的揚(yáng)程分別為86.8、80.5、69.8 m，添加了分流葉片的葉輪相對(duì)3長(zhǎng)葉片揚(yáng)程提高了15.3%，NPSHc提高了31.4%，6長(zhǎng)葉片揚(yáng)程相對(duì)3長(zhǎng)葉片揚(yáng)程提高了24.3%，NPSHc提高了81.3%。

4.3 葉輪中間斷面空泡分布

當(dāng)進(jìn)口壓力設(shè)置到一定值后，氣泡逐漸開(kāi)始形成，當(dāng)壓力繼續(xù)降低時(shí)，氣泡將隨著液體在流動(dòng)過(guò)程中不斷長(zhǎng)大。圖6為不同葉輪方案中間斷面空泡分布，其中φ為體積分?jǐn)?shù)。從圖6通過(guò)橫向?qū)Ρ瓤梢园l(fā)現(xiàn)，隨著NPSH的減小，泵內(nèi)的氣體體積分?jǐn)?shù)逐漸增多，在隔舌流道處空泡卻明顯偏少?？张蓦S著流體的流動(dòng)逐漸向葉輪出口擴(kuò)張。這是因?yàn)槿~片對(duì)流體做功，流道內(nèi)壓力沿流體流動(dòng)方向逐漸升高，所以空泡在葉片進(jìn)口處較明顯，卻沿著流動(dòng)方向在逐漸減弱。通過(guò)縱向的比較可以知道，由于流道空間被空泡占據(jù)，過(guò)流面積減小導(dǎo)致葉輪內(nèi)流體速度和壓力發(fā)生變化。隨著葉片數(shù)的增多，葉輪流道被空泡堵塞現(xiàn)象有所緩解。

圖6 不同葉輪方案中間斷面空泡分布Fig.6 Vapor volume distribution under different program

4.4 葉片表面載荷分布

由于工作面與吸力面相對(duì)速度的不同，使得葉片兩面存在壓力差，我們將同一葉片在相同半徑處工作面與吸力面的壓力差值稱為葉片表面載荷。圖7和圖8分別為3/3分流葉片和6長(zhǎng)葉片在設(shè)計(jì)工況下長(zhǎng)葉片表面靜壓分布曲線。

圖7 分流葉片設(shè)計(jì)工況長(zhǎng)葉片表面靜壓分布Fig.7 Pressure distribution of splitter blade

圖8 6長(zhǎng)葉片設(shè)計(jì)工況葉片表面靜壓分布Fig. 8 Pressure distribution of 6 blades

從靜壓變化趨勢(shì)可以看出，葉片表面靜壓分布與理論相符合。隨著徑向距離的增大，靜壓先降低到低谷，然后緩慢上升，到葉片出口處達(dá)到峰值。由于流體在葉輪吸入口轉(zhuǎn)彎時(shí)在離心力的作用下，吸入口流線在葉輪進(jìn)口存在一定的速度漩渦，靜壓最低點(diǎn)位于葉片進(jìn)口吸力面稍后處，此處也是最容易發(fā)生空化的區(qū)域。從曲線變化分析知道，葉輪內(nèi)靜壓分布均勻，即使葉片數(shù)不同，但沿流動(dòng)方向靜壓分布趨勢(shì)相近。在葉片吸力面和壓力面均存在壓力最小值，但是壓力面最小壓力值卻出現(xiàn)的比吸力面晚，這段徑向位置差可能是由于葉片厚度的存在和葉片進(jìn)口安放角的影響。隨著葉片數(shù)的增加，壓力面和吸力面靜壓都升高，與前面得出的葉片數(shù)增加可以改善空化性能結(jié)果相一致。但是葉片數(shù)的增加卻引起葉片進(jìn)口低壓區(qū)面積的增加，在大流量工況下運(yùn)行會(huì)增加空化風(fēng)險(xiǎn)。從圖7和圖8中還可以看出，工作面附近流體壓力升高的徑向距離明顯短于背面，而且采用了分流葉片的葉輪進(jìn)口處背面靜壓在較長(zhǎng)徑向距離內(nèi)壓力變化不大，導(dǎo)致葉片載荷較大。根據(jù)上節(jié)我們知道，分流葉片葉輪進(jìn)口背面產(chǎn)生的空泡數(shù)明顯多于6葉片葉輪，因此進(jìn)口背面靜壓的差異可能與葉片背面生成的空泡數(shù)有關(guān)。在設(shè)計(jì)工況下，隨著葉片數(shù)的增加，葉片載荷逐漸減小，對(duì)于采用分流葉片的葉輪，長(zhǎng)葉片載荷在徑向位置大約為2/3D2附近下降斜率增大，此處剛好在分流葉片進(jìn)口直徑位置圓周上。

5 結(jié) 語(yǔ)

(1)本文通過(guò)對(duì)3長(zhǎng)葉片,3/3分流葉片和6長(zhǎng)葉片三組葉輪方案進(jìn)行空化數(shù)值模擬計(jì)算，較直觀地得到空泡體積分布和葉片表面靜壓分布曲線。

(2)通過(guò)對(duì)大流量，設(shè)計(jì)流量和小流量下的空化性能曲線分析可以看出，隨著流量的增大，臨界汽蝕余量逐漸降低。由三組葉輪方案的空化曲線和中間斷面空泡分布可以看出，空化性能的改善與葉片形式緊密相關(guān)。分流葉片的采用，既減輕了葉輪進(jìn)口處的排擠嚴(yán)重現(xiàn)象，葉片進(jìn)口壓力也有了升高，空泡體積分?jǐn)?shù)減小，離心泵空化性能也有了明顯的改善。

(3)從不同葉片數(shù)葉片表面靜壓分布圖分析，可以觀察到隨著葉片數(shù)的增加，葉片表面載荷逐漸減小，對(duì)于采用分流葉片的葉輪，長(zhǎng)葉片表面載荷在分流葉片進(jìn)口直徑位置下降曲線變陡。葉片進(jìn)口背面空泡數(shù)與葉片進(jìn)口載荷密切相關(guān)。

□

[1] 關(guān)醒凡．現(xiàn)代泵技術(shù)手冊(cè)[M]．北京：宇航出版社，1995.

[2] 袁建平，李淑娟，付躍登．用正交試驗(yàn)研究分流葉片對(duì)離心泵性能的影響[J]．排灌機(jī)械，2009，27(5)：306-309．

[3] 袁壽其，張金鳳，袁建平，等．用正交試驗(yàn)研究分流葉片主要參數(shù)對(duì)性能影響[J]．排灌機(jī)械，2008，26(2)：1-5．

[4] Dular M，Bachert R，Bernd S，et al．Experimental evaluation of numerical simulation of cavitating flow around hydrofoil[J]．European Journal of Mechanics B(Fluids)，2005，24：522-538．

[5] 羅先武，張 瑤，彭俊奇，等．葉輪進(jìn)口幾何參數(shù)對(duì)離心泵空化性能的影響[J]．清華大學(xué)學(xué)報(bào)，2008，48(5)：836-839．

[6] 王 勇，劉厚林，袁壽其，等．離心泵內(nèi)部空化特性的CFD模擬[J]．排灌機(jī)械，2011，29(2)：99-103．

[7] 童啟忠．高汽蝕性能礦用排水泵葉輪設(shè)計(jì)的探討[J]．煤礦機(jī)械，1989,(3)：11-24．

[8] 施衛(wèi)東，吳蘇青，張德勝，等．葉片數(shù)對(duì)高比轉(zhuǎn)數(shù)軸流泵空化特性的影響[J]．農(nóng)業(yè)機(jī)械學(xué)報(bào)，2013,44(11)：72-77．

[9] 丁成偉．離心泵與軸流泵[M]．北京：機(jī)械工業(yè)出版社，1991.