韓 偉 鄭 昊 王麗瓊 張振中 蘇 敏

(1.蘭州理工大學(xué)能源與動(dòng)力工程學(xué)院， 蘭州 730050； 2.甘肅省流體機(jī)械重點(diǎn)實(shí)驗(yàn)室， 蘭州 730050；3.山東雙輪股份有限公司， 威海 264203； 4.蘭州理工大學(xué)電氣工程與信息工程學(xué)院， 蘭州 730050)

引言

雙吸離心泵具有單機(jī)流量大、揚(yáng)程高、結(jié)構(gòu)簡(jiǎn)單、便于維修等特點(diǎn)[1]，應(yīng)用于諸多領(lǐng)域。在水力機(jī)械中，大部分情況下空化是一種有害的現(xiàn)象，長(zhǎng)時(shí)間在空化條件下工作，將破壞過(guò)流部件，改變流道內(nèi)的速度分布，導(dǎo)致效率下降、泵的揚(yáng)程降低，引起機(jī)器振動(dòng)，產(chǎn)生噪聲，甚至影響泵的安全運(yùn)行[2]。因此，如何進(jìn)一步改善雙吸離心泵的空化性能，是水力機(jī)械重點(diǎn)研究方向之一。

LI等[3]通過(guò)采用仿真軟件Pumplinx及完全空化模型對(duì)雙吸離心泵進(jìn)行數(shù)值模擬分析并與試驗(yàn)進(jìn)行對(duì)比，結(jié)果表明，Pumplinx軟件在分析離心泵的流動(dòng)和空化特性時(shí)具有一定的可靠性；董興華等[4]采用數(shù)值模擬和外特性試驗(yàn)相結(jié)合方法，在空化條件下對(duì)雙吸泵進(jìn)行多工況下的定常數(shù)值研究，結(jié)果表明流量、葉輪結(jié)構(gòu)和雙蝸殼結(jié)構(gòu)對(duì)雙吸泵空化發(fā)展有影響；YAO等[5]通過(guò)試驗(yàn)分析在非空化和空化條件下，雙吸離心泵在半螺旋吸水室和壓水室中的壓力波動(dòng)頻率變化，并通過(guò)分析頻譜和時(shí)間頻域研究了空化條件下的水動(dòng)力響應(yīng)；MENG等[6]在設(shè)計(jì)流量和大流量對(duì)雙吸離心泵進(jìn)行數(shù)值模擬并進(jìn)行試驗(yàn)驗(yàn)證，結(jié)果表明，數(shù)值模擬可以較準(zhǔn)確地預(yù)測(cè)在雙吸離心泵的空化現(xiàn)象且在空化條件下在葉片背面會(huì)產(chǎn)生微型旋渦；肖若富等[7]通過(guò)采用數(shù)值模擬和試驗(yàn)研究相結(jié)合的方法，研究在變流量工況下長(zhǎng)短葉片復(fù)合葉輪和常規(guī)葉輪在中比轉(zhuǎn)數(shù)雙吸離心泵下的空化性能，并通過(guò)對(duì)葉輪內(nèi)部空化流場(chǎng)分析尋求提高空化性能的方法；劉建瑞等[8]通過(guò)對(duì)5種葉片交錯(cuò)角度進(jìn)行數(shù)值分析并與試驗(yàn)進(jìn)行對(duì)比，研究了交錯(cuò)葉片對(duì)雙吸離心泵壓力脈動(dòng)特性的影響；雷明川等[9]基于RNGk-ε湍流模型對(duì)3種方案的葉輪葉片交錯(cuò)布置進(jìn)行三維非定常湍流流場(chǎng)數(shù)值分析，研究了雙吸離心泵蝸殼內(nèi)的壓力脈動(dòng)分布規(guī)律。

在雙吸離心泵葉輪葉片交錯(cuò)布置的研究過(guò)程中，有關(guān)空化性能的相關(guān)問(wèn)題研究尚不多見(jiàn)。本文為更好地研究葉輪葉片交錯(cuò)角度對(duì)泵空化性能的影響，以6葉片雙吸葉輪為研究對(duì)象，相鄰葉片夾角為60°,選取葉片交錯(cuò)角為 0°、10°、15°、20°、30°等5組角度在設(shè)計(jì)工況條件下對(duì)其分別進(jìn)行研究。

1 數(shù)值計(jì)算

1.1 建模及網(wǎng)格劃分

以250GS80型雙吸離心泵作為研究對(duì)象。其中，葉輪葉片數(shù)Z=6，葉輪葉片交錯(cuò)角度為30°，葉輪出口直徑為490 mm。為研究葉輪葉片交錯(cuò)角度對(duì)空化性能的影響，葉輪保持原設(shè)計(jì)參數(shù)不變，改變?nèi)~輪葉片交錯(cuò)角度，選取交錯(cuò)角度分別為0°、10°、15°、20°、30°(圖1)5組葉輪方案進(jìn)行數(shù)值分析。

圖1 雙吸離心泵葉輪方案Fig.1 Impeller models of double-suction centrifugal pump

為了突出研究結(jié)論的共性，引入交錯(cuò)度的定義。雙吸葉片交錯(cuò)度計(jì)算公式為

θ=β/φ

式中φ——雙吸葉輪最大交錯(cuò)角

β——雙吸葉輪相鄰葉片交錯(cuò)角

由幾何結(jié)構(gòu)可知，6葉片雙吸葉輪最大交錯(cuò)角φ為30°，因此上述選取交錯(cuò)位置的交錯(cuò)度分別為0、1/3、1/2、2/3和1。

為使數(shù)值模擬盡可能接近于真實(shí)試驗(yàn)，通過(guò)分析CFD誤差產(chǎn)生機(jī)理及加工工藝需求得，在流動(dòng)空間建模過(guò)程中應(yīng)考慮腔體、工藝圓角等易忽視因素以消除物理模型近似誤差[10]。為使泵進(jìn)、出口水流速度保持均勻，不影響水泵效率，以泵實(shí)際管路安裝要求，在進(jìn)、出口處安裝長(zhǎng)度為水管口徑3倍以上的直管。分別取進(jìn)口延伸段管長(zhǎng)800 mm，出口延伸段管長(zhǎng)500 mm，進(jìn)行水體繪制并組裝，其組裝圖如圖2所示。

圖2 雙吸離心泵三維水體模型Fig.2 Three dimensional wave model of double-suction centrifugal pump1.進(jìn)口延伸段水體 2.吸水室水體 3.壓水室水體 4.出口延伸段水體 5.葉輪水體 6.泵腔水體

非結(jié)構(gòu)化網(wǎng)格能生成高質(zhì)量的網(wǎng)格,易控制網(wǎng)格大小和節(jié)點(diǎn)密度，有利于進(jìn)行網(wǎng)格自適應(yīng)，提高計(jì)算精度[11]。因此對(duì)吸水室、葉輪、壓水室等部件水體采用非結(jié)構(gòu)化網(wǎng)格劃分，并在扭曲較大處進(jìn)行局部網(wǎng)格加密。

在通常的計(jì)算中，離散誤差隨網(wǎng)格變細(xì)而減小，但由于網(wǎng)格變細(xì)時(shí)，離散點(diǎn)數(shù)增多，舍入誤差也隨之加大。為使數(shù)值計(jì)算更具有準(zhǔn)確性，在設(shè)計(jì)工況條件下選取6組不同數(shù)量網(wǎng)格進(jìn)行無(wú)關(guān)性驗(yàn)證，網(wǎng)格數(shù)量分別為153.055萬(wàn)、218.765萬(wàn)、326.675萬(wàn)、457.742萬(wàn)、570.803萬(wàn)、704.065萬(wàn)。根據(jù)數(shù)值模擬結(jié)果繪制出圖3所示曲線，從圖中可以看出，當(dāng)網(wǎng)格數(shù)量不小于457.742萬(wàn)時(shí)，泵揚(yáng)程、效率之間的誤差波動(dòng)均小于0.15%，因此選用網(wǎng)格數(shù)量為457.742萬(wàn)的網(wǎng)格模型作最終分析。

圖3 網(wǎng)格無(wú)關(guān)性檢查Fig.3 Mesh independence check

1.2 邊界條件

由于泵結(jié)構(gòu)中含有較多彎曲部位，因此湍流模型選用RNGk-ε模型，它既考慮了修正湍動(dòng)粘度，又考慮了平均流動(dòng)中的旋轉(zhuǎn)及旋流流動(dòng)情況。通過(guò)在大尺度運(yùn)動(dòng)和修正后的粘度項(xiàng)體現(xiàn)小尺度的影響，可以更好地處理高應(yīng)變率及流線彎曲程度較大的流動(dòng)[12]。為驗(yàn)證所選湍流模型的準(zhǔn)確性，選擇多重參考系MRF(Multiple reference frame),定義葉輪計(jì)算域?yàn)樾D(zhuǎn)坐標(biāo)系，進(jìn)口管、吸水室、壓水室、出口管和泵腔計(jì)算域?yàn)殪o止坐標(biāo)系，設(shè)置進(jìn)口條件為速度進(jìn)口，且假定進(jìn)口速度在軸向均勻分布；出口條件設(shè)置為質(zhì)量流量出口，壁面采用無(wú)滑移邊界條件，在近壁區(qū)采用標(biāo)準(zhǔn)壁面函數(shù)進(jìn)行外特性數(shù)值模擬。

在分析不同葉片交錯(cuò)度對(duì)泵空化特性影響時(shí)，進(jìn)口條件設(shè)置為總壓進(jìn)口；出口條件設(shè)置為質(zhì)量流量出口，設(shè)定進(jìn)口邊界中氣相體積分?jǐn)?shù)為0，液相體積分?jǐn)?shù)為1，水的汽化壓力為3 540 Pa，收斂精度為10-6。計(jì)算時(shí)，將單相流動(dòng)的計(jì)算收斂結(jié)果作為空化流動(dòng)計(jì)算的初始值以節(jié)省空化計(jì)算時(shí)間，選用基于Rayleigh-Plesset方程的Zwart空化模型[13]，保持出口質(zhì)量流量不變，通過(guò)逐漸降低進(jìn)口總壓，以實(shí)現(xiàn)雙吸離心泵內(nèi)部的空化性能的模擬。

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

2.1 數(shù)值模擬與試驗(yàn)對(duì)比

為便于數(shù)值模擬結(jié)果與試驗(yàn)讀取值進(jìn)行對(duì)比，搭建如圖4所示的試驗(yàn)臺(tái)，對(duì)雙吸葉輪交錯(cuò)度為1的研究對(duì)象，進(jìn)行外特性和空化特性試驗(yàn)結(jié)果測(cè)試。

圖4 試驗(yàn)臺(tái)布置示意圖Fig.4 Schematic diagram of test bench1.水封閘閥 2.水箱 3.真空表 4.雙吸離心泵 5.轉(zhuǎn)速轉(zhuǎn)矩儀 6.電動(dòng)機(jī) 7.壓力表 8.流量計(jì) 9.真空泵 10.調(diào)節(jié)閥

根據(jù)葉片泵外特性計(jì)算方法[14]，泵揚(yáng)程為

(1)

式中pin——泵進(jìn)口總壓力，Pa

pout——泵出口總壓力，Pa

ρ——流體密度，kg/m3

g——重力加速度，m/s2

泵效率為

η=ρgQH/(Mω)×100%

(2)

式中M——葉輪繞軸產(chǎn)生的總力矩，N·m

ω——葉輪旋轉(zhuǎn)角速度，rad/s

Q——流體體積流量，m3/s

在數(shù)值模擬過(guò)程中選取與雙吸離心泵試驗(yàn)時(shí)流量相同的5個(gè)工況點(diǎn)進(jìn)行外特性數(shù)值模擬計(jì)算，由式(1)、(2)計(jì)算得出揚(yáng)程和效率。繪制出外特性曲線對(duì)比圖(圖5)，從圖5可以得出，泵數(shù)值模擬值與試驗(yàn)值的變化趨勢(shì)基本一致，模擬揚(yáng)程與試驗(yàn)揚(yáng)程偏差最大時(shí)達(dá)1.61%， 模擬效率與試驗(yàn)效率偏差最大時(shí)達(dá)1.27%。存在的偏差主要來(lái)自于數(shù)值模擬時(shí)未考慮試驗(yàn)時(shí)機(jī)械損失、表面粗糙度等因素的存在。

圖5 交錯(cuò)度為1的雙吸泵性能曲線Fig.5 Pump performance curves at staggered degree of 1

泵有效空化余量計(jì)算公式[15]為

(3)

式中pva——飽和蒸汽壓力，Pa

vj——進(jìn)口速度，m/s

圖6 交錯(cuò)度為1的雙吸泵空化斷裂特性曲線Fig.6 Cavitation fracture performance curves at staggered degree of 1

將單相無(wú)空化計(jì)算結(jié)果作為空化流動(dòng)計(jì)算的初始值，在最優(yōu)設(shè)計(jì)工況點(diǎn)進(jìn)行雙吸離心泵的空化流動(dòng)數(shù)值模擬，繪制出空化曲線對(duì)比圖(圖6)，從圖中發(fā)現(xiàn)，數(shù)值模擬與試驗(yàn)的變化趨勢(shì)基本一致。未空化時(shí)，數(shù)值模擬揚(yáng)程與試驗(yàn)揚(yáng)程的絕對(duì)誤差為1.23%；取揚(yáng)程下降3%時(shí)所對(duì)應(yīng)的有效空化余量作為泵的臨界空化余量，則試驗(yàn)與數(shù)值模擬的有效空化余量分別為3.030 m和2.902 m。數(shù)值模擬臨界空化余量比試驗(yàn)所測(cè)值小4.22%。偏差主要來(lái)自于數(shù)值模擬忽視了汽液相之間的相對(duì)速度和相變的熱效應(yīng)破裂等因素的變化[16-17]。

通過(guò)在雙吸葉輪葉片交錯(cuò)度為1下進(jìn)行的試驗(yàn)與數(shù)值模擬對(duì)比發(fā)現(xiàn)，兩者誤差較小，變化趨勢(shì)基本一致，表明所選數(shù)值計(jì)算方法能較為準(zhǔn)確地進(jìn)行泵的外特性和空化特性的預(yù)測(cè)。

2.2 葉片交錯(cuò)布置時(shí)空泡體積分?jǐn)?shù)

圖7為不同有效空化余量下葉輪中間流面內(nèi)的空泡分布，從圖中可以看出，空泡首先在葉輪葉片進(jìn)口邊背面低壓區(qū)附近產(chǎn)生，這是因?yàn)檫M(jìn)口邊背面的圓周速度大于進(jìn)口其他位置，由速度三角形可知，其進(jìn)口壓力損失及進(jìn)口繞流引起的壓降相應(yīng)變大；葉輪流道為光滑曲面，介質(zhì)由于葉輪高速運(yùn)轉(zhuǎn)易產(chǎn)生離心力加快流動(dòng)速度，在葉片進(jìn)口邊，由于葉輪未對(duì)液體做功，根據(jù)伯努利方程守恒，則壓力會(huì)隨著速度的增大而減小，在低壓區(qū)液體易發(fā)生空化。并且葉輪不同流道內(nèi)的空泡分布不均勻，這是因?yàn)殡p吸離心泵中螺旋形壓水室就是為了入流均勻，但實(shí)際上即使螺旋形壓水室也沒(méi)有達(dá)到入流的絕對(duì)均勻一致，從而使得葉輪葉片表面的壓力分布不對(duì)稱。隨著有效空化余量的減小，泵內(nèi)部的空泡含量逐漸增大，沿葉片背面向出口方向擴(kuò)散，并向葉片工作面擴(kuò)展。通過(guò)對(duì)比，當(dāng)有效空化余量為5.552 m時(shí)，雙吸葉片交錯(cuò)度為1/3和1/2的部分流道中葉輪葉片背面已發(fā)現(xiàn)少許空泡，在交錯(cuò)度為0、2/3、1中未發(fā)現(xiàn)空泡；隨著有效空化余量下降不同交錯(cuò)度葉輪中開(kāi)始有空泡產(chǎn)生且在不同流道區(qū)域中分布不均，當(dāng)有效空化余量為2.799 m時(shí)，從圖7中可知，氣體體積分?jǐn)?shù)0.9～1.0隨無(wú)量綱徑向位置r/R(r表示徑向位置，R表示葉輪外徑)的分布位置：交錯(cuò)度為0時(shí)分布在0.482～ 0.818之間且不同葉輪流道分布程度不同；交錯(cuò)度為1/3時(shí)分布在0.482～0.836之間；交錯(cuò)度為1/2時(shí)分布在0.473～0.873之間；交錯(cuò)度為2/3時(shí)分布在0.463～0.909之間；交錯(cuò)度為1時(shí)分布在0.473～0.727之間且不同葉輪流道分布程度不同。與此同時(shí)，交錯(cuò)度為1/3、1/2和2/3時(shí)空泡占據(jù)葉輪流道大約80%的區(qū)域，完全發(fā)生空化，破壞泵內(nèi)的能量轉(zhuǎn)換過(guò)程。

圖7 葉輪中間流面內(nèi)空泡分布Fig.7 Vapor distributions on middle stream surface of impeller

同一交錯(cuò)度下，隨著有效空化余量的降低葉輪流道中的空泡含量逐漸增加；同一有效空化余量下的不同交錯(cuò)度中空泡含量不同，當(dāng)有效空化余量為2.799 m，低于臨界空化余量時(shí)，交錯(cuò)度為2/3中空泡含量相對(duì)最多，其次為1/2、1/3、0，交錯(cuò)度為1時(shí)相對(duì)最少。

2.3 不同交錯(cuò)度下泵揚(yáng)程變化

從圖8中可以看到：隨著有效空化余量的減小，揚(yáng)程先基本保持不變，當(dāng)有效空化余量減小到一定值時(shí)，不同交錯(cuò)度下?lián)P程會(huì)不同程度減小。這是因?yàn)殡S著有效空化余量逐漸減小，葉輪葉片背面進(jìn)口靠近輪緣處有空泡產(chǎn)生，雙吸離心泵開(kāi)始發(fā)生空化，但此時(shí)有效空化余量仍較大于臨界空化余量，泵的揚(yáng)程變化波動(dòng)不大；隨著有效空化余量繼續(xù)減小，當(dāng)有效空化余量小于臨界空化余量后，離心泵空化加劇，產(chǎn)生了大量的空泡，堵塞了流道，破壞流動(dòng)連續(xù)性，導(dǎo)致離心泵的揚(yáng)程急劇下降。

圖8 不同交錯(cuò)度下泵揚(yáng)程變化曲線Fig.8 Head curves of pump at different staggered degrees

通過(guò)圖8中對(duì)比可以發(fā)現(xiàn)：不同交錯(cuò)度下泵的空化曲線變化趨勢(shì)與試驗(yàn)和理論基本一致。有效空化余量大于5.552 m時(shí)，雙吸離心泵揚(yáng)程變化基本趨于穩(wěn)定，交錯(cuò)度為1的揚(yáng)程分別比交錯(cuò)度為0、2/3、1/3、1/2的揚(yáng)程高0.72%、1.57%、2.19%、2.26%。從圖中可以看出，不同交錯(cuò)度下雙吸離心泵的臨界空化余量不同，交錯(cuò)度為1的臨界空化余量分別比交錯(cuò)度為0、1/3、1/2、2/3小14.94%、16.19%、18.89%、19.68%。當(dāng)有效空化余量繼續(xù)下降，由于空泡含量的增多導(dǎo)致流動(dòng)連續(xù)性遭到破壞，泵揚(yáng)程急劇下降。當(dāng)有效空化余量為2.789 m時(shí)交錯(cuò)度為1的葉輪揚(yáng)程也開(kāi)始急劇下降。這表明2.789 m為各種交錯(cuò)度下改型雙吸葉輪的極限空化余量。

2.4 不同交錯(cuò)度下葉輪所受的扭矩變化

圖9為不同葉輪葉片交錯(cuò)度在設(shè)計(jì)工況下發(fā)生空化時(shí)雙吸離心泵葉輪所受扭矩的變化曲線。由理論分析知，泵在運(yùn)轉(zhuǎn)中，葉輪內(nèi)部的壓力是不同的，進(jìn)口處壓力較低，出口處壓力較高。若其過(guò)流部件的局部區(qū)域因?yàn)槟撤N原因，抽送液體的絕對(duì)壓力降至當(dāng)時(shí)溫度下的液體汽化壓力時(shí)，液體便在該處開(kāi)始汽化，形成空泡，當(dāng)含有大量空泡的液體向前流經(jīng)葉輪內(nèi)的高壓區(qū)時(shí)，空泡周圍的高壓液體致使空泡急劇地縮小以至破裂。在空泡凝結(jié)的同時(shí)，液體質(zhì)點(diǎn)將以高速填充空穴，發(fā)生互相撞擊而產(chǎn)生水擊，使得葉輪葉片受力突然增大，從而使得葉輪所受扭矩變大。因此隨著有效空化余量的降低，葉輪所受扭矩先基本保持穩(wěn)定，當(dāng)開(kāi)始發(fā)生空化時(shí)葉輪所受扭矩逐漸增大，隨著空化現(xiàn)象的繼續(xù)發(fā)展，葉輪載荷突降，導(dǎo)致葉輪所受扭矩先增大然后降低。當(dāng)空化發(fā)展至嚴(yán)重空化時(shí)葉輪破損嚴(yán)重扭矩發(fā)生急劇下降。

圖9 不同交錯(cuò)度下葉輪的扭矩變化曲線Fig.9 Torque curves of impeller at different staggered degrees

從圖9中可以看出，有效空化余量大于5.552 m時(shí)，雙吸離心泵扭矩的變化基本趨于穩(wěn)定，交錯(cuò)度為1的扭矩分別比交錯(cuò)度為0、2/3、1/3、1/2的扭矩高0.31%、0.65%、0.76%、1.33%。隨著有效空化余量的降低，不同交錯(cuò)度的葉輪扭矩出現(xiàn)不同的變化趨勢(shì)，交錯(cuò)度為1、2/3和1/3的變化幅度不大，其中，交錯(cuò)度為1的葉輪所受的扭矩變化較為穩(wěn)定，交錯(cuò)度為1/2時(shí)，葉輪的扭矩變化幅度較大，扭矩下降最明顯。這是由于在臨界空化狀態(tài)時(shí)發(fā)生完全空化時(shí)，交錯(cuò)度為1/2發(fā)生完全空化區(qū)域占據(jù)葉輪流道最多，交錯(cuò)度為1發(fā)生完全空化區(qū)域最少，完全空化的發(fā)生使得葉輪葉片受力發(fā)生不同程度的改變；當(dāng)葉輪完全空化時(shí)，葉輪受損嚴(yán)重，流體流動(dòng)性被破壞，導(dǎo)致葉輪所受扭矩急劇下降，從圖9中可知，交錯(cuò)度為1/2的葉輪所受扭矩最早急劇下降，而交錯(cuò)度為1的葉輪最遲，這表明在空化狀態(tài)下，扭矩急劇下降時(shí)交錯(cuò)度為1/2的葉輪中的損失較大，交錯(cuò)度為1的葉輪中損失最小，交錯(cuò)度為1時(shí)葉輪扭矩急劇下降點(diǎn)的有效空化余量與交錯(cuò)度為0、2/3、1/3、1/2的葉輪扭矩急劇下降點(diǎn)絕對(duì)誤差分別為0.17%、3.89%、7.15%、10.21%。當(dāng)有效空化余量為2.789 m，交錯(cuò)度為1時(shí)葉輪所受扭矩開(kāi)始急劇下降，因此，這5種不同交錯(cuò)角度下葉輪扭矩產(chǎn)生陡降對(duì)應(yīng)的有效空化余量均不同。

2.5 不同有效空化余量下靜壓變化分析

由式(3)可知，當(dāng)葉輪進(jìn)口壓力小于當(dāng)?shù)販囟人柡驼羝麎毫r(shí)，靜壓低于當(dāng)?shù)販囟人柡驼羝麎毫λ趨^(qū)域會(huì)發(fā)生完全空化。圖10給出了不同有效空化余量下葉輪流道內(nèi)的靜壓分布，其中紅色區(qū)域?yàn)殪o壓高于當(dāng)?shù)販囟人柡驼羝麎毫Φ牟糠郑溆酁樾∮诘扔诋?dāng)?shù)仄瘔毫Φ牟糠?。從圖中可以看出，相同有效空化余量下葉輪中間流面發(fā)生完全空化的面積在不同交錯(cuò)度中是不相同的。當(dāng)有效空化余量為3.921 m時(shí)，不同交錯(cuò)度葉輪流道中存在完全空化區(qū)域且在不同流道區(qū)域中分布不均，交錯(cuò)度為1時(shí)汽相體積分?jǐn)?shù)約占葉輪流道5%的區(qū)域，交錯(cuò)度為2/3中含量最少。隨著有效空化余量的減小，當(dāng)有效空化余量為2.799 m時(shí)，交錯(cuò)度為0和1時(shí)發(fā)生完全空化區(qū)域較小且各流道間分布不均勻，其中交錯(cuò)度為1發(fā)生完全空化區(qū)域最少，約占葉輪流道15%；交錯(cuò)度為1/3、1/2和2/3中完全發(fā)生空化區(qū)域基本填充各個(gè)葉輪流道，其中，交錯(cuò)度為2/3發(fā)生完全空化區(qū)域占據(jù)葉輪流道最多，約占80%。

圖10 不同有效空化余量下靜壓分布區(qū)域Fig.10 Static pressure distributions at different NPSH

圖11 不同交錯(cuò)度葉片表面和葉輪出口面靜壓分布Fig.11 Static pressure distributions of blade surface and impeller outlet surface at different staggered degrees

通過(guò)分析當(dāng)有效空化余量為2.799 m時(shí)不同交錯(cuò)度下葉輪葉片表面和葉輪出口面靜壓云圖(圖11)可以看出：葉輪出口處從背面到工作面存在著明顯的壓力梯度，方向指向工作面。當(dāng)葉輪交錯(cuò)度為0時(shí)，雙吸葉輪兩側(cè)工作面與工作面，背面與背面完全重合，雙吸葉輪兩側(cè)出口處壓力沿葉輪隔板呈對(duì)稱分布，葉輪出口兩側(cè)基本無(wú)壓差存在。雙吸葉輪兩側(cè)出口壓力無(wú)法進(jìn)行相互補(bǔ)償，導(dǎo)致葉輪流道內(nèi)壓力以較大梯度變化；當(dāng)交錯(cuò)度為1/3、1/2和2/3時(shí)，雙吸葉輪一側(cè)工作面附近出流為高壓區(qū)，由于葉輪葉片交錯(cuò)布置，該側(cè)葉輪高壓區(qū)出流將對(duì)應(yīng)另一側(cè)葉輪背面低壓區(qū)，由于壓差的存在，一側(cè)葉輪出口高壓區(qū)與另一側(cè)葉輪低壓區(qū)進(jìn)行非對(duì)稱性壓力補(bǔ)償，從而影響葉輪流道內(nèi)局部區(qū)域壓力變化。相對(duì)于交錯(cuò)度為0的對(duì)稱布置雙吸葉輪，交錯(cuò)度為1/3、1/2和2/3的流道內(nèi)壓力梯度更大，使得葉輪流道內(nèi)的流動(dòng)形態(tài)變得更加紊亂，能量損失增加，低于臨界空化壓力的區(qū)域相應(yīng)增加，這樣反而進(jìn)一步加劇了空化的發(fā)生。其中交錯(cuò)度為2/3時(shí)，在已研究的交錯(cuò)度范圍內(nèi)，葉輪出口處的非對(duì)稱性壓力補(bǔ)償對(duì)葉輪內(nèi)壓力梯度的負(fù)面影響達(dá)到最大，空化最嚴(yán)重，這也與空化斷裂特性一致；當(dāng)交錯(cuò)度為1時(shí)，雙吸葉輪兩側(cè)葉片為對(duì)稱交錯(cuò)布置，一側(cè)高壓區(qū)與另一側(cè)低壓區(qū)形成對(duì)稱性壓力補(bǔ)償，有效降低了兩側(cè)流道內(nèi)的壓力梯度，降低了葉輪內(nèi)流動(dòng)損失，也在一定程度上減少了低于臨界空化壓力的區(qū)域，從而有助于減小葉輪內(nèi)空化的嚴(yán)重程度。

因此可以得出，該型雙吸離心泵當(dāng)葉輪葉片交錯(cuò)度為1時(shí)，空化性能相對(duì)最好。

3 結(jié)論

(1) 通過(guò)數(shù)值模擬與試驗(yàn)對(duì)比知，兩者誤差較小，模擬揚(yáng)程與試驗(yàn)揚(yáng)程偏差最大時(shí)達(dá)1.61%， 模擬效率與試驗(yàn)效率偏差最大時(shí)達(dá)1.27%；臨界空化余量數(shù)值模擬比試驗(yàn)值小4.22%，數(shù)值模擬值與試驗(yàn)數(shù)據(jù)的變化趨勢(shì)基本一致，所選數(shù)值計(jì)算方法能較為準(zhǔn)確地進(jìn)行泵外特性和空化特性的預(yù)測(cè)。

(2) 同一交錯(cuò)度下，隨著有效空化余量的降低葉輪流道中的空泡含量逐漸增加；同一有效空化余量下的不同交錯(cuò)度中空泡含量不同，當(dāng)有效空化余量為2.799 m，小于臨界空化余量時(shí)，交錯(cuò)度為2/3時(shí)空泡含量最多，其次依次為1/2、1/3、0，交錯(cuò)度為1時(shí)相對(duì)最少。

(3) 發(fā)生空化斷裂時(shí)，5種不同交錯(cuò)度下雙吸葉輪扭矩和揚(yáng)程急劇下降時(shí)對(duì)應(yīng)的有效空化余量均不同，表明雙吸葉輪的交錯(cuò)度對(duì)其空化性能有不同程度的影響。

(4) 當(dāng)有效空化余量低于臨界空化余量、雙吸葉輪葉片交錯(cuò)度為2/3時(shí)，發(fā)生完全空化區(qū)域占據(jù)葉輪流道最多，交錯(cuò)度為1時(shí)發(fā)生完全空化區(qū)域最少。因此對(duì)該型雙吸葉輪而言，交錯(cuò)度為1時(shí)雙吸離心泵的空化特性相對(duì)最好。

1 丁成偉.離心泵與軸流泵原理及水力設(shè)計(jì)[M].北京:機(jī)械工業(yè)出版社,1981:5-6.

2 張克危.流體機(jī)械原理[M].北京:機(jī)械工業(yè)出版社,2012:137-170.

3 LI P, HUANG Y F, LI J. Cavitation simulation and NPSH prediction of a double suction centrifugal pump[C]∥IOP Conference Series: Earth and Environmental Science, 2012,15(6):2025.

4 董興華，畢禎，周棟，等.雙吸離心泵葉輪空化流動(dòng)預(yù)測(cè)的數(shù)值分析 [J]. 西華大學(xué)學(xué)報(bào),2016,35(6)：17-22.

DONG Xinghua, BI Zhen, ZHOU Dong, et al. The prediction and numerical simulation of impeller cavitation performance in double suction pump [J].Journal of Xihua University, 2016,35(6): 17-22. (in Chinese)

5 YAO Z F, WANG F J, XIAO R F, et al. Experimental investigation of pressure instabilities affected by cavitation for a double-suction centrifugal pump[C]∥IOP Conference Series: Earth and Environmental Science, 2012, 15(6): 2040.

6 MENG G, TAN L, CAO S L, et al. Numerical simulation and analysis of cavitation flows in a double suction centrifugal pump[C]∥IOP Conference Series: Materials Science and Engineering,2015,72(3):2020.

7 肖若富，王娜，楊魏，等.復(fù)合葉輪改善雙吸式離心泵空化性能研究[J/OL].農(nóng)業(yè)機(jī)械學(xué)報(bào), 2013, 44(9): 35-39. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20130907&flag=1.DOI:10.6041/j.issn.1000-1298.2013.09.007.

XIAO Ruofu, WANG Na, YANG Wei, et al. Improvement of cavitation characteristics of double suction centrifugal pump with compound impeller [J/OL].Transactions of the Chinese Society for Agricultural Machinery, 2013,44(9):35-39.(in Chinese)

8 劉建瑞,付威,高振軍,等.交錯(cuò)葉片對(duì)雙吸離心泵性能影響的數(shù)值分析[J].排灌機(jī)械工程學(xué)報(bào),2015, 33(3):196-202.

LIU Jianrui, FU Wei, GAO Zhenjun, et al. Numerical analysis of effects of staggered blades on characteristics of double- suction centrifugal pump[J].Journal of Drainage and Irrigation Machinery Engineering, 2015, 33(3):196-202. (in Chinese)

9 雷明川，賴喜德，宋冬梅，等.交錯(cuò)葉片葉輪對(duì)雙吸離心泵蝸殼內(nèi)壓力脈動(dòng)的影響研究[J].中國(guó)農(nóng)村水利水電,2014(7):177-181.

LEI Mingchuan, LAI Xide, SONG Dongmei, et al. A study of staggered arrangement blades on the pressure pulsation of a double-entry pump with volute casing [J].China Rural Water and Hydropower,2014(7): 177-181. (in Chinese)

10 JOHN D DENTON. Some limitations of turbo machinery CFD [C]∥ASME Turbo Expo 2010: Power for Land, Sea and Air, 2010,7:735-745.

11 劉厚林,董亮,王勇,等.流體機(jī)械CFD中的網(wǎng)格生成方法進(jìn)展[J].流體機(jī)械,2010,38(4):32-37.

LIU Houlin, DONG Liang, WANG Yong, et al. Overview on mesh generation methods in CFD of fluid machinery [J]. Fluid Machinery,2010,38(4)：32-37. (in Chinese)

12 王福軍.計(jì)算流體動(dòng)力學(xué)分析——CFD軟件原理與應(yīng)用[M].北京:清華大學(xué)出版社,2010:124-125.

13 曹玉良,賀國(guó),明廷鋒,等.水泵空化數(shù)值模擬研究進(jìn)展[J].武漢理工大學(xué)學(xué)報(bào),2016,40(1):55-65.

CAO Yuliang, HE Guo, MING Tingfeng, et al. Review on numerical simulation of pump cavitation [J]. Journal of Wuhan University of Technology, 2016, 40(1):55-65. (in Chinese)

14 董亮，劉厚林.葉片泵CFD數(shù)值計(jì)算實(shí)例詳解[M].北京:機(jī)械工業(yè)出版社，2015:4.

15 施衛(wèi)東,李通通,張德勝,等. 軸流泵葉輪區(qū)域空化特性數(shù)值模擬[J].農(nóng)業(yè)工程學(xué)報(bào),2012,28(13):88-93.

SHI Weidong, LI Tongtong, ZHANG Desheng, et al. Numerical simulation on cavitation characteristic in impeller of axial-flow pump [J]. Transactions of the CSAE, 2012, 28(13): 88-93. (in Chinese)

16 王秀禮,袁壽其,朱榮生,等.離心泵汽蝕過(guò)渡過(guò)程瞬態(tài)特性分析[J].西安交通大學(xué)學(xué)報(bào),2012,46(7)：38-43.

WANG Xiuli, YUAN Shouqi, ZHU Rongsheng, et al. Analysis on transient hydraulic characteristics of cavitation process in centrifugal pumps[J].Journal of Xi’an Jiaotong University,2012,46(7)：38-43. (in Chinese)

17 施衛(wèi)東,張光建,張德勝,等.水氣最大密度比對(duì)軸流泵空化計(jì)算的影響[J/OL].農(nóng)業(yè)機(jī)械學(xué)報(bào),2014,45(5):65-70. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20140510&flag=1.DOI:10.6041/j.issn.1000-1298.2014.05.010.

SHI Weidong, ZHANG Guangjian, ZHANG Desheng, et al. Effect of the maximum density ratio between liquid and vapor on cavitation simulation of axial-flow pump [J/OL].Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(5):65-70. (in Chinese)