李春蓄,沙 毅,胡雅文

(浙江科技學(xué)院 機械與能源工程學(xué)院,杭州 310023)

離心泵應(yīng)用領(lǐng)域極為廣泛,在生產(chǎn)制造中占據(jù)著重要位置[1]。葉輪是離心泵的核心構(gòu)件,在針對葉輪的優(yōu)化中葉片包角是一項重要的設(shè)計參數(shù),它對水泵性能有著極大的影響[2]。葛書亭等[3-4]利用數(shù)值模擬方法對不同葉片內(nèi)速度場進行大量分析,結(jié)果表明葉片存在一個最佳包角值。謝志賓等[5-6]針對離心泵包角對壓力脈動展開研究,得到一個適合水泵平穩(wěn)運行的葉片包角值。牟介剛等[7-8]對包角與水泵特性變化關(guān)系展開研究,同時加入了抗氣蝕性能的設(shè)計因素,從而完善了包角的選擇理論。但是,針對中轉(zhuǎn)速比離心泵葉片包角選值研究相對較少,故本研究以一臺轉(zhuǎn)速比為129的中轉(zhuǎn)速比離心泵為研究對象,在其他葉輪參數(shù)恒定的情況下,分別設(shè)計90°、100°、110°、120°、130° 共5種包角的葉輪。利用Fluent仿真軟件對不同包角葉輪進行非定常數(shù)值模擬分析,得到不同包角葉輪的揚程、效率數(shù)值及內(nèi)部流場分布云圖。通過在流道及蝸殼內(nèi)布置測點得到對應(yīng)點壓力的脈動特性,總結(jié)水泵內(nèi)部流場及壓力脈動變化規(guī)律。

1 計算模型的建立

研究對象為一臺轉(zhuǎn)速比ns為129的單級單吸式離心泵,額定流量qv為4.8 m3/h,設(shè)計揚程H為4.5 m,額定轉(zhuǎn)速n為3 000 r/min。葉輪幾何參數(shù)見表1。在保持葉輪其他設(shè)計參數(shù)不變的前提下,對葉片包角α分別設(shè)計為90°、100°、110°、120°、130°,根據(jù)葉輪其他參數(shù)進行三維建模。計算流體力學(xué)數(shù)值模擬需將水泵模型進行流體域劃分,且為確保仿真精度與真實性,需在水泵進出口處各增加一段圓柱形流體域以延長進出口流動路徑,使得流體在進出口處充分擴散進而將進出口對流體的影響降至最低[9]。不同包角葉輪流體域模型如圖1所示,離心泵全流道流體域模型如圖2所示。

表1 葉輪幾何參數(shù)Table 1 Geometric parameters of impeller

圖1 不同包角葉輪流體域模型Fig.1 Fluid domain model of impellers with different wrap angles

圖2 離心泵全流道流體域模型Fig.2 Fluid domain model in whole flow passage of centrifugal pump

2 數(shù)值模擬方法

2.1 模型網(wǎng)格劃分

網(wǎng)格劃分的質(zhì)量直接影響仿真精度及收斂程度的計算[10]。使用網(wǎng)格劃分軟件,利用非結(jié)構(gòu)化四面體網(wǎng)格能較好適應(yīng)葉輪及蝸殼的過渡性曲面。在蝸殼突變處、蝸舌狹小區(qū)域及葉輪較小表面進行局部加密操作能進一步合理劃分網(wǎng)格,進而減少收斂誤差以保障仿真精度。模型總網(wǎng)格數(shù)量過少會使網(wǎng)格質(zhì)量較差導(dǎo)致模擬失真甚至無法收斂的情況發(fā)生,模型網(wǎng)格密度越小越接近真實情況,但此時網(wǎng)格數(shù)量會極其多,受計算機性能限制,需進行網(wǎng)格無關(guān)性驗證以選擇最佳網(wǎng)格數(shù)量[11]。以90°包角葉輪為例進行網(wǎng)格無關(guān)性驗證(表2),當(dāng)模型網(wǎng)格總數(shù)約為50萬個時,由于網(wǎng)格質(zhì)量較低導(dǎo)致模擬計算未收斂;當(dāng)網(wǎng)格總數(shù)約為100萬個及200萬個時,揚程最大相對誤差為23.1%,未滿足無關(guān)性要求;當(dāng)網(wǎng)格總數(shù)約為250萬個及以上時,揚程最大相對誤差僅為1.68%,此時可忽略網(wǎng)格數(shù)量帶來的誤差影響。為了節(jié)省模擬仿真收斂時間,最終決定采用247萬個左右的網(wǎng)格數(shù)量為劃分標(biāo)準。離心泵各部分網(wǎng)格數(shù)見表3,離心泵全流道網(wǎng)格模型如圖3所示。

表2 網(wǎng)格無關(guān)性驗證數(shù)據(jù)

表3 離心泵各部分網(wǎng)格數(shù)

圖3 離心泵全流道網(wǎng)格模型Fig.3 Mesh model in whole flow passage of centrifugal pump

2.2 邊界條件與壓力監(jiān)測點設(shè)定

數(shù)值模擬使用Fluent軟件,選擇重整化群(renormalization group,RNG)k-ε。RNGk-ε是一種高雷諾數(shù)湍流計算方法,對求解旋轉(zhuǎn)運動中高應(yīng)變率及流線變化較大的流體運動有著較高的仿真精度,且RNGk-ε模型是對標(biāo)準k-ε模型的優(yōu)化,能較準確地模擬動靜壁面處的流動特性。流體材料為25 ℃的液體水介質(zhì),使用壓力耦合方程的半隱相容算法(semi-implicit method for pressure linked equation consistent,SIMPLEC)進行計算,求解精度會進一步提升[12]。收斂條件設(shè)置為10-4,對旋轉(zhuǎn)機械采用多重參考系模型(multireference frame,MRF),設(shè)置葉輪為轉(zhuǎn)動域,轉(zhuǎn)速為3 000 r/min。入口段、出口段及蝸殼均為固定域,進口設(shè)置為速度進口,出口采用壓力出口,其余設(shè)置保持軟件默認[13]。非定常計算中,以首次定常計算結(jié)果為初始值。取葉輪每旋轉(zhuǎn)3°的時長為一個時間步,時間步Δ=1.667 10-4s,每個時間步中最大迭代次數(shù)為20次[14]。計算5個旋轉(zhuǎn)周期,總時長t=0.1 s。為避免計算波動,選擇最后一個周期作為數(shù)據(jù)采集周期。同時為研究葉片包角對離心泵壓力脈動的影響,在葉輪流道中線處均布L1、L2、L33個測點,在水泵蝸殼內(nèi)均布P1、P2、P3、P44個測點。測點位置如圖4所示,測點極坐標(biāo)見表4。

圖4 壓力脈動監(jiān)測點分布圖Fig.4 Distribution of pressure fluctuation monitoring points

表4 測點極坐標(biāo)

3 仿真結(jié)果分析

3.1 離心泵外部特性曲線分析

用數(shù)值模擬得到的數(shù)據(jù)計算出不同包角葉輪進出口壓力差值。離心泵揚程H與水泵水力效率η[15]計算公式如下:

(1)

(2)

式(1)～(2)中:P2為葉輪出口壓力,Pa;P1為葉輪進口壓力,Pa;ρ為流體密度,kg/m3;g為重力加速度,m/s2;Δh為進出口垂直高度差,m;Q為水泵流量,m3/h;M為葉輪軸扭矩,N·m;ω為葉輪軸旋轉(zhuǎn)角速度,rad/s。

圖5為不同葉片包角離心泵外部特性曲線圖。由圖5流量-揚程曲線可知,隨著葉片包角的增大,離心泵呈現(xiàn)小流量時揚程上升而大流量時揚程下降的趨勢,且包角越大下降趨勢越明顯,并以額定流量qv=4.8 m3/h為大小流量分界點。隨葉片包角增大,額定流量下?lián)P程H降低,且包角增至130°時揚程比包角90°時降低13%。這是由于包角較大時,葉輪喉部的平均安放角減小,導(dǎo)致葉輪的相對液流角逐漸變小,從而使葉輪揚程逐漸降低。由圖5流量-效率曲線可知,隨著葉片包角增大,離心泵在小流量工況下效率有所提升,而在大流量工況下效率下降較為明顯,且包角越大效率下降幅度越大,同時效率最高點也逐漸向小流量方向偏移。額定流量下,130°包角時的效率比90°包角時下降5.2%。這是由于隨著葉片包角的增大,葉輪流道逐漸變得狹長,葉片對流體的約束能力增強,流體運動更符合葉片形狀流動,也更趨于平緩。但同時葉片表面對流體的摩擦損失也逐漸增大,從而導(dǎo)致水泵效率降低。分析圖5曲線趨勢可知,較小的葉片包角可得到較高的水泵揚程及效率,但過小的葉片包角會使流道對流體的約束力不足,減少葉片的有效部分,導(dǎo)致水泵效率降低,因此必存在一個較合適的葉片包角使水泵性能達到最佳狀態(tài)。

圖5 不同包角離心泵外部特性曲線Fig.5 External characteristic curve of centrifugalpump with different wrap angles

3.2 離心泵內(nèi)部流場分析

3.2.1 葉輪靜壓力分布云圖

圖6為在額定流量下不同包角葉輪內(nèi)部靜壓力分布云圖。由圖6可知,各葉輪靜壓力均從進口至出口逐漸增大,這是由于葉輪旋轉(zhuǎn)時對流體做功的結(jié)果。隨著葉片包角的增大,葉輪進口低壓區(qū)域變大而出口高壓區(qū)域減少。其原因是包角增大使得葉片型線變平順,更貼合流體運動,從而減少了脫流現(xiàn)象的發(fā)生;但過大的包角會產(chǎn)生較大的摩擦,進而降低水泵的揚程與效率。

圖6 額定流量下不同包角葉輪內(nèi)部靜壓力分布云圖Fig.6 Static pressure distribution of impellers with different wrap angles under rated flow

3.2.2 葉輪速度分布云圖

圖7為額定流量下不同包角葉輪速度分布云圖。由圖7可知,葉輪內(nèi)的速度分布是不完全均勻的。葉輪內(nèi)低速區(qū)主要集中在葉片的工作面,隨著包角的增大低速區(qū)逐漸減少,葉片背面及葉輪出口的速度逐漸增大。這是由于包角的增加,減少了流道內(nèi)的當(dāng)量擴散角,葉片對流體的束縛力增強,但同時大包角會使流道變窄,進而增加了葉輪的摩擦損失。

圖7 額定流量下不同包角葉輪速度分布云圖Fig.7 Velocity distribution of impellers with different wrap angles under rated flow

3.3 壓力脈動特性分析

由于流體的黏性、慣性及葉輪與蝸殼間的動靜干涉造成水泵內(nèi)出現(xiàn)一種周期性的壓力脈動,這種規(guī)律性的脈動沖擊會造成水泵運行不穩(wěn)定[16]。為研究中轉(zhuǎn)速比離心泵壓力脈動特性,對圖4布置的各監(jiān)測點的壓力數(shù)據(jù)采用無量綱處理,引入壓力系數(shù)Cp[17],其計算公式如下:

(3)

式(3)中:Δp為各點實測壓力與平均壓力之差,Pa;ρ為流體密度,kg/m3;u2為葉輪出口圓周速度,m/s。

本研究的中轉(zhuǎn)速比離心泵額定轉(zhuǎn)速n為3 000 r/min,故泵轉(zhuǎn)動頻率為50 Hz。泵葉片數(shù)Z為6,故葉頻為300 Hz。將計算得到的壓力系數(shù)Cp進行快速傅里葉變換,從而得到不同包角下各點的壓力脈動頻域圖。圖8為額定流量下葉輪流道監(jiān)測點壓力脈動頻域特性圖,圖9為額定流量下蝸殼內(nèi)監(jiān)測點壓力脈動頻域特性圖。

圖8 額定流量下流道監(jiān)測點壓力脈動頻域特性Fig.8 Frequency domain characteristics of pressure fluctuation at monitoring points of flow-path under rated flow

圖9 額定流量下蝸殼監(jiān)測點壓力脈動頻域特性Fig.9 Frequency domain characteristics of pressure fluctuation at monitoring points of volute under rated flow

由圖8可知,在額定工況下流道內(nèi)各點脈動主頻以一倍轉(zhuǎn)頻為主,次頻以一倍葉頻為主;隨著葉片包角的增大,各點脈動幅值均有所提高;相同葉片包角情況下沿流體流動方向,各測點主頻與次主頻幅值均明顯提高且低頻段波動增多,出現(xiàn)葉頻代替轉(zhuǎn)頻成為主頻的趨勢。這是由于流道外側(cè)監(jiān)測點L3最靠近蝸殼,受蝸殼與葉輪間的動靜干涉作用影響最大,并且由于前段脈動在蝸舌處消散不及時,第二波脈動在此重疊,使得L3點經(jīng)過蝸舌處時再次受到脈動影響,這是導(dǎo)致L3點在低頻段波動增加的主要原因。

由圖9可知,額定工況下蝸殼內(nèi)監(jiān)測點頻率集中在0～900 Hz,脈動主頻約為一倍葉頻,次主頻約為一倍轉(zhuǎn)頻。隨著葉片包角的增大,主頻與次頻幅值均整體呈上升趨勢,但120°包角時壓力脈動幅值略有所降低,尤其在接近水泵出口P4的點降低較為明顯;相同葉片包角情況下沿P1、P2、P3、P4方向,脈動幅值逐漸減弱,其中P1點脈動最為劇烈,P4點脈動最弱。這是由于P1最靠近蝸舌,蝸舌處徑向P4間隙較小,蝸殼與葉片動靜干涉最為強烈,隨著監(jiān)測點向P4點移動,逐漸遠離蝸舌,徑向間隙變大,故壓力脈動逐漸減弱。

4 結(jié) 論

以一臺比轉(zhuǎn)速為129的中轉(zhuǎn)速比離心泵為研究對象,通過數(shù)值模擬對比分析5種不同包角葉輪的水泵外部特性、內(nèi)部壓力與速度分布及壓力脈動特性,得出以下結(jié)論:

1)隨著葉片包角的增大,離心泵的揚程和效率均有所下降且效率最高點向小流量方向偏移。130°葉片包角相比90°包角時揚程及效率分別降低13%與5.2%。

2)額定工況下,葉片包角的增大導(dǎo)致葉輪進口低壓區(qū)面積增大,而出口高壓區(qū)面積減少,流體更貼合葉片流動,進而減少了脫流現(xiàn)象的發(fā)生。

3)額定工況下,低速區(qū)主要分布在葉片工作面,高速區(qū)集中在葉輪出口及葉片背面處。隨著葉片包角增大,低速區(qū)域逐步減少,出口速度明顯上升。

4)離心泵壓力脈動特性受蝸舌影響較大。額定流量下蝸殼壓力脈動主頻約為一倍葉頻,葉輪流道內(nèi)壓力脈動主頻約為一倍轉(zhuǎn)頻。隨著葉片包角的增大,各點主次頻均整體呈現(xiàn)上升趨勢,其中130°包角時各點壓力脈動幅值最大。相比其他大包角,120°包角時壓力脈動幅值略微降低。綜合分析離心泵外部特性、流場特性及脈動特性,于是可取120°作為中轉(zhuǎn)速比離心泵最優(yōu)包角。