卓長(zhǎng)青,邵春雷,吉雨翔
(南京工業(yè)大學(xué)機(jī)械與動(dòng)力工程學(xué)院,江蘇 南京 211816)
離心泵是流體輸送過(guò)程中較常見(jiàn)的機(jī)械,在實(shí)際應(yīng)用中,離心泵輸送的并非是單一液體,多混有氣體、蒸汽等其他相[1-2]。離心泵吸水室的作用是將流體按一定的條件引入葉輪內(nèi),以保證葉輪進(jìn)口流體的速度大小和方向符合要求。吸水室內(nèi)部流體的流動(dòng)狀態(tài)將直接影響其在葉輪等過(guò)流部件內(nèi)的流動(dòng),從而對(duì)離心泵的效率以及汽蝕性能產(chǎn)生較大的影響[3]。因此,研究氣液混輸狀態(tài)下其在吸水室內(nèi)的流動(dòng),對(duì)提高離心泵的設(shè)計(jì)水平、保證離心泵的高效穩(wěn)定運(yùn)行起到至關(guān)重要的作用。
國(guó)內(nèi)外相關(guān)學(xué)者對(duì)離心泵內(nèi)氣液兩相流動(dòng)進(jìn)行了大量研究。Verde等[4-5]利用高速攝像機(jī)對(duì)離心泵內(nèi)部的流場(chǎng)進(jìn)行了可視化試驗(yàn),研究了進(jìn)口氣相體積分?jǐn)?shù)與葉輪內(nèi)流型之間的關(guān)系。Zhan等[6-7]通過(guò)數(shù)值研究發(fā)現(xiàn),離心泵進(jìn)口處的渦流會(huì)增加泵內(nèi)流動(dòng)的不穩(wěn)定性,從而對(duì)泵的性能產(chǎn)生不良影響。李思[8]采用數(shù)值模擬方法,對(duì)平直型、漸擴(kuò)型和漸縮型吸水室進(jìn)行了對(duì)比分析,發(fā)現(xiàn)平直型吸水室的水力性能最佳。牟介剛等[9]基于噴射原理,設(shè)計(jì)了不同噴射管數(shù)和引流管徑的引射吸水室,得到了最佳引射吸水室設(shè)計(jì)結(jié)構(gòu)。王健等[10-11]研究了吸水室內(nèi)置隔板對(duì)離心泵性能的影響。李書(shū)磊等[12-14]對(duì)水平管內(nèi)氣液兩相流的流型及渦流進(jìn)行了研究,為吸水室內(nèi)流型的分析提供了參考。相關(guān)文獻(xiàn)對(duì)離心泵輸送氣液兩相介質(zhì)時(shí)葉輪內(nèi)的流動(dòng)分析較多,但對(duì)吸水室缺乏系統(tǒng)深入的研究,尤其是吸水室內(nèi)的氣液兩相流動(dòng)規(guī)律有待進(jìn)一步揭示。
本文采用可視化試驗(yàn)結(jié)果對(duì)數(shù)值模擬結(jié)果進(jìn)行驗(yàn)證,并通過(guò)數(shù)值模擬方法研究吸水室內(nèi)氣液兩相流動(dòng)。詳細(xì)研究泵內(nèi)各種流型出現(xiàn)的條件,探討氣體螺旋長(zhǎng)度的影響因素及影響規(guī)律,闡明吸水室內(nèi)湍流強(qiáng)度的分布。通過(guò)對(duì)比分析吸水室內(nèi)流型與泵外特性,建立吸水室內(nèi)流型與泵外特性之間的關(guān)系。
1—儲(chǔ)水罐;2—進(jìn)出口調(diào)節(jié)閥;3—金屬軟管;4—渦輪流量計(jì);5—進(jìn)口壓力傳感器;6—進(jìn)氣管;7—出口壓力傳感器;8—電控柜;9—電機(jī);10—轉(zhuǎn)速傳感器;11—離心泵;12—攝像機(jī);13—轉(zhuǎn)子流量計(jì);14—?dú)怏w控制閥;15—減壓閥;16—?dú)獗茫?7—冷光源。圖1 試驗(yàn)裝置示意圖Fig.1 Schematic diagram of experimental device
圖1為氣液兩相流離心泵試驗(yàn)裝置示意圖。試驗(yàn)過(guò)程中,電控柜中變頻器控制泵的轉(zhuǎn)速(n),出口調(diào)節(jié)閥控制液體流量,氣泵提供氣體,氣體流量由氣體控制閥控制。試驗(yàn)中的模型泵由有機(jī)玻璃制作,其尺寸參考了IS80-50-250型離心泵,額定流量Qd=25 m3/h,額定揚(yáng)程Hd=20 m,額定轉(zhuǎn)速nd=1 450 r/min,比轉(zhuǎn)速ns=47。離心泵的進(jìn)口管直徑為80 mm,出口管直徑為50 mm,進(jìn)氣管直徑為10 mm。
試驗(yàn)選取清水和空氣作為試驗(yàn)介質(zhì),將攝像機(jī)拍攝方向垂直于吸水室豎直軸截面,調(diào)節(jié)攝像機(jī)與吸水室出口截面之間的軸向距離為0.2 m??紤]到有機(jī)玻璃模型泵的強(qiáng)度,試驗(yàn)驗(yàn)證時(shí)將泵的轉(zhuǎn)速調(diào)節(jié)為600 r/min,調(diào)節(jié)液體流量分別為2 m3/h、4 m3/h、6 m3/h、8 m3/h、10 m3/h、12 m3/h以及進(jìn)出口調(diào)節(jié)閥門(mén)全開(kāi)時(shí)的液體流量(18.3 m3/h),不同液體流量對(duì)應(yīng)的氣體流量均從0.5 L/min逐漸增加,直至液體流量出現(xiàn)下降趨勢(shì)(臨界狀態(tài))。待吸水室內(nèi)流動(dòng)穩(wěn)定后,采用攝像機(jī)對(duì)吸水室內(nèi)進(jìn)行拍攝,同時(shí)測(cè)量吸水室內(nèi)的氣體螺旋長(zhǎng)度,并實(shí)時(shí)記錄進(jìn)出口壓力等數(shù)據(jù)。
圖2(a)為離心泵幾何模型,其中進(jìn)氣管直徑、吸水室直徑以及出口管直徑均與試驗(yàn)設(shè)備相同,進(jìn)氣管軸線距離吸水室出口位置與試驗(yàn)裝置相同,即相距450 mm,出口管段延伸至試驗(yàn)裝置中壓力傳感器的位置,即延伸200 mm。對(duì)葉輪和蝸殼進(jìn)行四面體網(wǎng)格劃分,對(duì)吸水室部分進(jìn)行六面體結(jié)構(gòu)化網(wǎng)格劃分,在壁面設(shè)置邊界層網(wǎng)格,為了計(jì)算邊界層的黏性區(qū)域,調(diào)整邊界y+≈1[15]。圖2(b)為組合后的流體域網(wǎng)格,放大區(qū)域?yàn)槲掖怪陛S截面的網(wǎng)格。為節(jié)省計(jì)算時(shí)間,同時(shí)保證計(jì)算的精確度,對(duì)幾何模型進(jìn)行網(wǎng)格獨(dú)立性驗(yàn)證,當(dāng)網(wǎng)格數(shù)從688 514增加到1 624 907時(shí),離心泵的揚(yáng)程和效率增加幅度顯著,網(wǎng)格數(shù)從1 624 907增加到2 484 537時(shí),離心泵揚(yáng)程和效率的波動(dòng)范圍均小于0.4%,因此最終確定幾何模型的網(wǎng)格總數(shù)為1 624 907,其中吸水室部分網(wǎng)格數(shù)為365 486。
圖2 離心泵幾何模型及計(jì)算流體域網(wǎng)格Fig.2 Geometric model of centrifuge pump and grids of computational fluid domain
選用商用CFD軟件Fluent對(duì)離心泵內(nèi)流體流動(dòng)進(jìn)行數(shù)值模擬。葉輪與吸水室以及葉輪與蝸殼之間的交界面均為動(dòng)靜結(jié)合的面,數(shù)值模擬中動(dòng)靜邊界的耦合選用多重參考系模型(MRF),用靜止坐標(biāo)系對(duì)吸水室內(nèi)的流動(dòng)進(jìn)行分析。氣液兩相流模型選用Eulerian-Eulerian非均相流模型,其中水為主相,空氣為次相,氣液兩相之間的曳力模型選擇Schiller-Naumann模型。k-ε模型是一種典型的近壁湍流模型,可以對(duì)近壁黏性區(qū)域進(jìn)行求解。因吸水室內(nèi)部有渦流的存在,本文選用RNGk-ε模型。RNGk-ε模型是在標(biāo)準(zhǔn)k-ε模型基礎(chǔ)上的改進(jìn)模型,其精度得到了提高,具有更強(qiáng)的適應(yīng)性[16]。
數(shù)值模擬采用氣相和液相兩個(gè)進(jìn)口:液相進(jìn)口采用速度進(jìn)口邊界條件,水力直徑為80 mm,速度大小根據(jù)液體流量確定,氣相體積分?jǐn)?shù)為0;氣相進(jìn)口也采用速度進(jìn)口邊界條件,水力直徑為10 mm,速度大小根據(jù)氣體流量確定,氣相體積分?jǐn)?shù)為1。進(jìn)口處湍流強(qiáng)度由經(jīng)驗(yàn)公式計(jì)算求得,選擇自由出流作為出口邊界條件。采用Phase Coupled SIMPLE算法求解壓力與速度耦合的方程式,固壁面選用無(wú)滑移邊界條件,粗糙度為0,用標(biāo)準(zhǔn)壁面函數(shù)法處理近壁面區(qū)域。
流型1:如圖3(a)所示,當(dāng)液相流量Ql=4 m3/h,進(jìn)口氣相體積分?jǐn)?shù)φi=0.5%時(shí),吸水室內(nèi)出現(xiàn)形狀近似橢球形的氣泡,氣泡流動(dòng)一定距離后,開(kāi)始被液體切碎并隨旋流進(jìn)入泵內(nèi),此種流型稱之為穩(wěn)態(tài)螺旋泡狀流。流型2:如圖3(b)所示,當(dāng)Ql=4 m3/h,φi=1.5%時(shí),吸水室內(nèi)含氣量明顯增加,氣泡在非螺旋區(qū)域呈塞狀流動(dòng),在螺旋區(qū)域氣泡分布區(qū)域變大,此流型稱之為穩(wěn)態(tài)螺旋塞狀流。流型3:如圖3(c)所示,吸水室內(nèi)氣泡較小,數(shù)量較少,形狀仍近乎橢球形,有部分小氣泡附著在大氣泡附近,大氣泡之間保持一定的距離,從吸水室上部進(jìn)入葉輪內(nèi),此流型稱之為穩(wěn)態(tài)非螺旋泡狀流。流型4:如圖3(d)所示,氣泡間距減小,開(kāi)始發(fā)生聚合,變成長(zhǎng)度逐漸增大的塞狀氣泡,直至氣泡彼此相連,該流型稱之為穩(wěn)態(tài)非螺旋塞狀流。流型5:如圖3(e)所示,大量微小的氣泡分布在吸水室上部較大的區(qū)域,此時(shí)吸水室內(nèi)流型稱之為脈動(dòng)非螺旋泡狀流。流型6:如圖3(f)所示,小氣泡的數(shù)量減少,逐漸聚合成為大氣泡,且在大氣泡附近附著較多的小氣泡,此流型稱之為脈動(dòng)非螺旋塞狀流。
圖3 吸水室內(nèi)流型Fig.3 Flow patterns in suction chamber
圖4為上述各流型對(duì)應(yīng)工況下數(shù)值模擬得到的吸水室內(nèi)軸截面上氣相體積分?jǐn)?shù)φg分布情況,為了便于比較各流型,對(duì)吸水室內(nèi)φg>10%的部分用灰色表示,φg≤10%的部分按色標(biāo)顏色區(qū)分。
圖4 吸水室內(nèi)氣相體積分?jǐn)?shù)Fig.4 Gas volume fraction in suction chamber
由圖4(a)(b)可見(jiàn),吸水室內(nèi)出現(xiàn)明顯的氣相螺旋流動(dòng)現(xiàn)象,在吸水室出口位置附近,氣相不再單一地分布于吸水室上部,而是分布在整個(gè)吸水室的出口截面。保持液體流量不變,增加氣體流量至吸水室內(nèi)出現(xiàn)塞狀流時(shí),氣相在吸水室軸截面徑向的分布區(qū)域略有增加,且φg>10%的區(qū)域明顯增加。保持φi不變,增加液體流量至吸水室內(nèi)出現(xiàn)脈動(dòng)流時(shí),氣相在吸水室軸截面徑向的分布區(qū)域變大,且φg<10%的區(qū)域明顯增加。對(duì)比模擬與試驗(yàn)的氣相分布可以發(fā)現(xiàn),兩者基本吻合。數(shù)值模擬雖然不能給出氣泡大小以及氣泡間距的變化,但可以準(zhǔn)確地模擬出氣相體積分?jǐn)?shù)分布情況,因此利用該數(shù)值模擬方法可以對(duì)吸水室內(nèi)流型進(jìn)行深入研究。
圖5為試驗(yàn)與模擬所得氣體螺旋長(zhǎng)度L0的對(duì)比,工況為轉(zhuǎn)速n=600 r/min,氣相流量Qg=0.5 L/min,液相流量Ql=2 m3/h、4 m3/h、6 m3/h、8 m3/h、10 m3/h、12 m3/h,可以發(fā)現(xiàn)試驗(yàn)和數(shù)值模擬得到的L0相近,模擬得到的L0比試驗(yàn)得到的數(shù)值略小,這是因?yàn)樵囼?yàn)中模型泵葉輪與旋轉(zhuǎn)軸的銜接處為不規(guī)則結(jié)構(gòu),而數(shù)值模擬過(guò)程對(duì)模型泵的結(jié)構(gòu)進(jìn)行了適當(dāng)?shù)暮?jiǎn)化處理,但試驗(yàn)與模擬得到的L0誤差小于10%,數(shù)值模擬得到的L0在允許范圍內(nèi),因此可以通過(guò)數(shù)值模擬方法對(duì)L0進(jìn)行影響因素分析。
圖5 試驗(yàn)與模擬所得氣體螺旋長(zhǎng)度比較Fig.5 Comparison of gas spiral length between experiment and simulation
圖6為吸水室內(nèi)6種流型對(duì)應(yīng)工況下試驗(yàn)與數(shù)值模擬所得泵揚(yáng)程及效率對(duì)比。由圖6可以發(fā)現(xiàn),數(shù)值模擬所得泵的揚(yáng)程和效率均略高于試驗(yàn)所得結(jié)果,這是因?yàn)閿?shù)值模擬中只考慮了泵的水力損失,而在試驗(yàn)中除了水力損失外,機(jī)械損失和容積損失均無(wú)法避免,這些損失包括葉輪輪盤(pán)的摩擦損失、軸承的摩擦損失、密封裝置及口環(huán)處的泄漏損失等。數(shù)值模擬與試驗(yàn)得到的數(shù)據(jù)規(guī)律性一致,且兩者之間的誤差均在可以接受的范圍之內(nèi),表明本文采用的數(shù)值模擬方法是準(zhǔn)確可靠的。因此,采用數(shù)值模擬方法對(duì)其他工況下吸水室內(nèi)流體流動(dòng)進(jìn)行進(jìn)一步研究。
圖6 試驗(yàn)與模擬所得泵外特性比較Fig.6 Comparison of external characteristics of pump between experiment and simulation
當(dāng)離心泵進(jìn)口氣相體積分?jǐn)?shù)小于臨界氣相體積分?jǐn)?shù)φic,吸水室內(nèi)氣液兩相流流型可以按照脈動(dòng)流與穩(wěn)態(tài)流、泡狀流與塞狀流以及螺旋流與非螺旋流進(jìn)行區(qū)分。通過(guò)試驗(yàn)發(fā)現(xiàn):吸水室內(nèi)脈動(dòng)流與穩(wěn)態(tài)流主要由液體流量決定,定義由穩(wěn)態(tài)流過(guò)渡到脈動(dòng)流時(shí)的液體流量為Qlp;吸水室內(nèi)泡狀流與塞狀流主要由氣體流量和液體流量?jī)烧吖餐瑳Q定,定義由泡狀流轉(zhuǎn)變?yōu)槿麪盍鲿r(shí)的進(jìn)口氣相體積分?jǐn)?shù)為φip;吸水室內(nèi)的螺旋流與非螺旋流主要由泵的轉(zhuǎn)速和液體流量共同決定,泵的轉(zhuǎn)速越大,泵能夠達(dá)到的最大液體流量Qlmax以及出現(xiàn)螺旋流時(shí)的液體流量Qls也越大,定義Qls=Qlp時(shí)泵的轉(zhuǎn)速為nc。當(dāng)泵的轉(zhuǎn)速n
表1 吸水室內(nèi)出現(xiàn)的流型
試驗(yàn)可以發(fā)現(xiàn)隨著氣體流量的增大,L0基本不變。為研究液體流量和泵的轉(zhuǎn)速對(duì)L0的影響,選取5種不同的轉(zhuǎn)速,固定氣體流量為1 L/min。根據(jù)泵的比例定律計(jì)算出不同轉(zhuǎn)速下的額定流量Qd,液體流量分別選Qd、0.8Qd、0.6Qd、0.4Qd和0.2Qd。統(tǒng)計(jì)各工況對(duì)應(yīng)的L0,繪制曲線如圖7所示。由圖7可以發(fā)現(xiàn):(a)隨著液體流量的增加,L0逐漸減短,且L0越小,液體流量對(duì)L0的影響越弱。(b)隨著轉(zhuǎn)速的增大,L0逐漸增長(zhǎng),但增幅逐漸下降,這是因?yàn)槁菪魇怯筛咚傩D(zhuǎn)的葉片帶動(dòng)具有一定黏性的液體運(yùn)動(dòng)而產(chǎn)生的,由于流動(dòng)過(guò)程中能量的損失,這種帶動(dòng)作用隨著距離的增加而不斷減弱。(c)隨著轉(zhuǎn)速的增大,出現(xiàn)氣體螺旋流時(shí)的液體流量逐漸增大,不同轉(zhuǎn)速對(duì)應(yīng)的額定流量下氣體螺旋長(zhǎng)度為0,表明液體流量大于額定流量后,吸水室內(nèi)不會(huì)出現(xiàn)螺旋流。
圖7 不同轉(zhuǎn)速下氣體螺旋長(zhǎng)度隨 液體流量的變化Fig.7 Variation of gas spiral length with liquid phase flow at different rotational speeds
湍流強(qiáng)度I可以定義為速度波動(dòng)的均方根與平均速度的比值,其計(jì)算公式為
I=0.16Re-1/8
(1)
式中:Re——雷諾數(shù)。可以發(fā)現(xiàn)氣體流量對(duì)湍流強(qiáng)度的影響較小。對(duì)吸水室內(nèi)脈動(dòng)非螺旋泡狀流(圖8(a))、穩(wěn)態(tài)非螺旋泡狀流(圖8(b))、脈動(dòng)螺旋泡狀流(圖8(c))和穩(wěn)態(tài)螺旋泡狀流(圖8(d))4種流型下的湍流強(qiáng)度進(jìn)行研究。設(shè)吸水室內(nèi)不同橫截面至吸水室出口截面之間的軸向距離為z,吸水室縱截面中心高度至采樣點(diǎn)之間的距離為R。圖8為φi=0.5%時(shí)不同截面處豎直直徑上I隨z的變化趨勢(shì)。由圖8可以發(fā)現(xiàn):(a)在非螺旋流中,出現(xiàn)z較小時(shí)吸水室中心處I較大的現(xiàn)象,這是因?yàn)樵撐恢么嬖谝后w螺旋流動(dòng),但其強(qiáng)度較弱,沒(méi)有帶動(dòng)吸水室上部的氣體旋轉(zhuǎn),從而導(dǎo)致液體螺旋長(zhǎng)度大于氣體螺旋長(zhǎng)度。(b)I隨著液體流量的增加而增加,且在螺旋流區(qū)域大幅度增加,這是因?yàn)樵诼菪恢?,周向速度的增大?dǎo)致液體的速度波動(dòng)增大。(c)I在同一截面上并非對(duì)稱分布,且在螺旋區(qū)域變化幅度較大,表明螺旋流既不穩(wěn)定也不對(duì)稱。
圖8 吸水室不同截面處豎直直徑上湍流強(qiáng)度分布Fig.8 Turbulence intensity distribution on vertical diameter at different cross sections of suction chamber
圖9為額定轉(zhuǎn)速下不同工況對(duì)應(yīng)離心泵的揚(yáng)程和效率,從圖9(a)(b)可以看出,當(dāng)φi從0.5%增大至2.5%時(shí),泵的揚(yáng)程和效率均呈下降趨勢(shì),表明保持液體流量不變,增大φi使吸水室內(nèi)流型從泡狀流轉(zhuǎn)變?yōu)槿麪盍鲿r(shí),泵的性能會(huì)逐漸降低。由圖9(c)可以發(fā)現(xiàn),保持φi不變時(shí),泵的揚(yáng)程隨轉(zhuǎn)速的增大而增大,隨液體流量的增大而減小,表明減小液體流量或增大泵的轉(zhuǎn)速,使吸水室內(nèi)流型從非螺旋流轉(zhuǎn)變?yōu)槁菪鲿r(shí),泵的揚(yáng)程會(huì)大幅度增大。由圖9(d)可以發(fā)現(xiàn),額定流量對(duì)應(yīng)泵的效率最高,由額定流量下吸水室內(nèi)為非螺旋流,表明泵的效率最高時(shí),吸水室內(nèi)為非螺旋流。當(dāng)吸水室內(nèi)流型從穩(wěn)態(tài)流轉(zhuǎn)變?yōu)槊}動(dòng)流時(shí),液體流量增大,因此泵的揚(yáng)程會(huì)下降,效率在該轉(zhuǎn)速對(duì)應(yīng)的額定流量下為最大值,因此效率最高時(shí)吸水室內(nèi)流型可能為穩(wěn)態(tài)流,也可能為脈動(dòng)流。綜上所述,離心泵效率最高時(shí),吸水室內(nèi)流型為穩(wěn)態(tài)非螺旋泡狀流或脈動(dòng)非螺旋泡狀流。
圖9 不同工況下泵的外特性Fig.9 External characteristics of pump under different working conditions
a.輸送氣液兩相介質(zhì)的離心泵其吸水室內(nèi)的流型可分為穩(wěn)態(tài)螺旋泡狀流、穩(wěn)態(tài)螺旋塞狀流、穩(wěn)態(tài)非螺旋泡狀流、穩(wěn)態(tài)非螺旋塞狀流、脈動(dòng)螺旋泡狀流、脈動(dòng)螺旋塞狀流、脈動(dòng)非螺旋泡狀流和脈動(dòng)非螺旋塞狀流8種,某一轉(zhuǎn)速下吸水室內(nèi)最多出現(xiàn)6種流型。
b.氣體螺旋長(zhǎng)度隨著液體流量的增大而減短,隨著轉(zhuǎn)速的增大而增長(zhǎng)。
c.吸水室出口位置存在液體螺旋,且液體螺旋長(zhǎng)度大于氣體螺旋長(zhǎng)度;吸水室內(nèi)湍流強(qiáng)度隨著液體流量的增大而增大,且在螺旋流區(qū)域大幅度增大;螺旋流并非對(duì)稱分布。
d.輸送氣液兩相介質(zhì)的離心泵,其吸水室內(nèi)流型為穩(wěn)態(tài)螺旋泡狀流時(shí),泵的揚(yáng)程較大;離心泵效率最高時(shí),吸水室內(nèi)流型為穩(wěn)態(tài)非螺旋泡狀流或脈動(dòng)非螺旋泡狀流。