離心式風(fēng)機(jī)流動特性的數(shù)值分析與改型設(shè)計

唐建光, 朱懿淵, 姚 征, 蔣旭平

(上海理工大學(xué)能源與動力工程學(xué)院,上海 200093)

離心風(fēng)機(jī)中通常存在較強(qiáng)的逆壓梯度和黏性分離,加上旋轉(zhuǎn)葉輪和蝸殼通道流動狀態(tài)的非周期性組合,使其內(nèi)部流動相當(dāng)復(fù)雜,通過傳統(tǒng)的實驗手段和理論分析方法難以對其內(nèi)部流場進(jìn)行深入的研究.近十年來,隨著大容量、高速度計算機(jī)的普及,以及計算流體力學(xué)(CFD)技術(shù)的迅速發(fā)展,三維黏性流動數(shù)值模擬方法逐漸成為離心風(fēng)機(jī)設(shè)計的重要輔助手段.人們通過三維數(shù)值模擬來深入考察風(fēng)機(jī)內(nèi)部流動的結(jié)構(gòu)特點,找出改善流動特性的措施,從而提高氣動性能.

目前國內(nèi)外對離心風(fēng)機(jī)和壓縮機(jī)內(nèi)部流動的數(shù)值模擬采用的數(shù)學(xué)模型基本上都是以雷諾時均方程組為控制方程,湍流模型采用一方程或兩方程的渦黏性模型,并借助于各種CFD通用軟件.在2002年,文獻(xiàn)[1]運用CFD軟件CFX-TASCflow對一離心壓縮機(jī)進(jìn)行了數(shù)值模擬,并著重討論了蝸殼軸向流動的非均勻性及其對下游蝸殼內(nèi)流場結(jié)構(gòu)的影響.在2001年,文獻(xiàn)[2]利用三維黏性數(shù)值模擬軟件分析了離心通風(fēng)機(jī)內(nèi)部流場.在2002年,文獻(xiàn)[3]對離心通風(fēng)機(jī)在設(shè)計工況時的整機(jī)內(nèi)部流場進(jìn)行了數(shù)值模擬,捕捉到了離心通風(fēng)機(jī)內(nèi)部許多重要的流動現(xiàn)象,證實了蝸殼、葉輪間相互作用引起的流場不對稱性.在2004年,文獻(xiàn)[4]對柜式空調(diào)器用多翼離心通風(fēng)機(jī)的三維流場進(jìn)行了數(shù)值模擬.在2006年,文獻(xiàn)[5]對在輸送氣體、物料及高壓強(qiáng)制通風(fēng)等領(lǐng)域應(yīng)用廣泛的9-26型高壓離心通風(fēng)機(jī)進(jìn)行了整機(jī)三維數(shù)值模擬.在2009年,文獻(xiàn)[6]利用CFX軟件計算了蝸殼壁厚對離心風(fēng)機(jī)振動噪聲的影響,與實驗結(jié)果相吻合.

本文對某一軍用艦艇上的通風(fēng)機(jī)進(jìn)行數(shù)值模擬,獲取其性能曲線,對其流場進(jìn)行分析,在此基礎(chǔ)上提出改進(jìn)方案,以優(yōu)化其性能.該通風(fēng)機(jī)在改進(jìn)后,氣流在蝸殼內(nèi)擴(kuò)壓充分,全壓和效率有了顯著的提高.船舶機(jī)艙是船舶動力系統(tǒng)的核心所在,這里空間相對狹小,而機(jī)艙中因為機(jī)器工作產(chǎn)生的熱量、廢氣會散發(fā)到周圍空氣當(dāng)中,通風(fēng)機(jī)擔(dān)負(fù)著機(jī)艙內(nèi)排出高溫廢氣和提供新鮮空氣的任務(wù)[7-8],所以,提高船舶通風(fēng)機(jī)效率很有意義.今后在同類型通風(fēng)機(jī)的新產(chǎn)品設(shè)計中,就可以方便地參照本文的計算和分析方法來進(jìn)行性能預(yù)測和分析.本文借助CFD軟件Fluent進(jìn)行計算.

1 風(fēng)機(jī)的幾何模型與修改

離心風(fēng)機(jī)由集流器、葉輪及蝸殼等組成.葉輪包括12個長葉片和12個分流小葉片,進(jìn)、出口直徑D1和D2分別為190.6 mm和430.7 mm,葉片進(jìn)、出口安裝角分別為24°和90°,蝸殼寬度 B=148 mm,葉片寬度b=100 mm.圖1為結(jié)構(gòu)示意圖.

圖1 風(fēng)機(jī)結(jié)構(gòu)Fig.1 Fan structure

由文獻(xiàn)[9]中的公式確定葉片出口尺寸

式中,D1m為葉片進(jìn)口平均直徑,這里取D1m=D1;μ1和μ2分別為葉輪的氣流充滿系數(shù),取μ1=μ2;Kcm為氣流離開葉片后的子午速度與葉片入口前氣流速度的子午速度的比值,這里取為1.

計算也葉片出口寬度b2=45 mm.

圖2 葉片F(xiàn)ig.2 Blade

2 數(shù)學(xué)模型和網(wǎng)格生成

本文所求解的控制方程組是三維不可壓雷諾時均納維爾-斯托克斯方程組.

湍流模型:標(biāo)準(zhǔn)κ-ε模型,κ為湍動能,ε為湍流耗散率;

壁面邊界條件:無滑移邊界條件,近壁區(qū)域采用標(biāo)準(zhǔn)壁面函數(shù)[10];

進(jìn)口邊界條件:規(guī)定進(jìn)口質(zhì)量流量;

出口邊界條件:自由出流條件.計算中添加延長段使出口處流向一致朝外,由于出口處對流強(qiáng)度較大,故可忽略下游流動的影響;

綜上所述，年齡小于60歲竇性心動過緩患者的DC值增高，迷走神經(jīng)張力增高。DC值和HRV指標(biāo)都能反映心臟自主神經(jīng)功能，DC作為一項獨立定量評估人體迷走神經(jīng)張力大小的新技術(shù)，它不易受外界環(huán)境的干擾，對竇性心動過緩患者的評估價值較高，值得進(jìn)一步探討。

動靜結(jié)合面的處理:本文采用Fluent軟件中的MRF(moving reference frame)模型,這是最簡化的處理方法,只是將運動和靜止系統(tǒng)中的速度變量在互相傳遞時進(jìn)行簡單轉(zhuǎn)換.在實際問題中,這種傳遞是非定常過程,所以,這種處理隱含一種“瞬間凍結(jié)”假設(shè),即在某一瞬間將流場“固定”,然后進(jìn)行轉(zhuǎn)換,這就帶來一個問題,“凍結(jié)”的瞬間,兩個系統(tǒng)相對位置的不同可能會導(dǎo)致不同的結(jié)果,所以,對轉(zhuǎn)換時相對位置的影響須加以考察.

葉柵通道和蝸殼通道是計算域的兩個主要部分,對網(wǎng)格質(zhì)量的要求較高.此外,還有入口段、擴(kuò)壓段、葉輪前蓋和后盤與蝸殼壁間的空隙區(qū)等,網(wǎng)格生成時分區(qū)較多,相鄰區(qū)網(wǎng)格的銜接不能錯位.為了獲得較高的網(wǎng)格質(zhì)量,并控制其數(shù)量,葉柵和蝸殼通道均采用非結(jié)構(gòu)型三梭柱和結(jié)構(gòu)型四梭柱這兩種網(wǎng)格的組合.葉柵通道內(nèi)繞葉片的邊界層區(qū)采用結(jié)構(gòu)形的四梭柱網(wǎng)格,其余部分則用非結(jié)構(gòu)的三梭柱網(wǎng)格.蝸殼通道內(nèi),彎道的大部分和擴(kuò)壓段均采用結(jié)構(gòu)形的四梭柱網(wǎng)格,但在蝸舌附近區(qū)域和彎道出口朝向擴(kuò)壓段的部分區(qū)域則采用三梭柱網(wǎng)格,并對蝸舌附近網(wǎng)格進(jìn)行加密,如圖3所示.

圖3 網(wǎng)格分布Fig.3 Gird distribution

近壁區(qū)采用了標(biāo)準(zhǔn)壁面函數(shù),按該設(shè)置要求取網(wǎng)格尺度因子y+=30,由此定出第一層網(wǎng)挌高度.在用Gambit軟件生成葉片上的邊界層網(wǎng)格時應(yīng)注意,當(dāng)在一個體的相鄰表面上各自生成的邊界層網(wǎng)格相遇時就會相互選加,形成不合理的細(xì)密網(wǎng)格.這種情況通過選用軟件中的internal continuity選項可以避免(圖4),使網(wǎng)格的質(zhì)量得到提高、數(shù)量大為減少.

圖4 邊界層網(wǎng)格Fig.4 Boundary layer grid

整個風(fēng)機(jī)的網(wǎng)格數(shù)近100萬,其中,約85%的網(wǎng)格為四梭柱與三梭柱,其網(wǎng)格角畸變度小于0.65,其余為四面體網(wǎng)格和少數(shù)金字塔網(wǎng)格,其網(wǎng)格角畸變度小于0.8,保持了較高的網(wǎng)格質(zhì)量.

3 計算結(jié)果及分析

3.1 性能曲線

風(fēng)機(jī)的設(shè)計流量為 2 500 m3/h,在1 500～3 500 m3/h這一區(qū)間內(nèi)共取8個工況點,分別對改進(jìn)前、后的風(fēng)機(jī)進(jìn)行模擬計算.圖5顯示流量-全壓和流量-效率曲線.風(fēng)機(jī)效率

式中,Ne為有效功率;N為軸功率;P為風(fēng)機(jī)全壓; Qs為風(fēng)機(jī)流量;I為風(fēng)機(jī)葉輪扭矩;ω為風(fēng)機(jī)轉(zhuǎn)速.

改進(jìn)風(fēng)機(jī)的全壓在設(shè)計流量處與原風(fēng)機(jī)相近,大流量的全壓增加了,小流量的全壓則減小了,從而使全壓隨流量的變化趨于平緩.改進(jìn)后全壓效率顯著提高,設(shè)計流量時提高達(dá)9.47%,由于設(shè)計工況下全壓變化不大,所以,效率提高主要是軸功率降低導(dǎo)致的.改型后風(fēng)機(jī)的全壓效率曲線也比較平坦,有利于風(fēng)機(jī)的運行[11].另外,最大全壓效率更接近設(shè)計工況(風(fēng)機(jī)流量Q=2 500 m3/h),符合風(fēng)機(jī)設(shè)計的要求.

圖5 性能曲線Fig.5 Performance curve

3.2 流場分析

3.2.1 壓力分布

現(xiàn)以流量為2 500 m3/h的設(shè)計工況為例,對離心風(fēng)機(jī)的壓力分布進(jìn)行分析.在距輪蓋和輪盤20%葉高處各取一個剖面,又在半葉高處取一剖面,得到軸向不同位置的3個橫剖面,在這3個剖面上對原始模型及改進(jìn)模型進(jìn)行分析比較,如圖6所示.由于改型風(fēng)機(jī)的輪蓋側(cè)和中葉高的剖面都剖切了錐形輪蓋,從而在壓力云圖的相應(yīng)位置上留下了白圈分隔葉輪的內(nèi)外空間.

圖6 3個橫剖面的示意圖Fig.6 Illustrations of cross sections

圖7～9分別為輪盤側(cè)、中間面和輪蓋側(cè)這3個橫剖面的壓力云圖.在葉輪內(nèi),靜壓沿半徑方向逐漸增大.靜壓最小值出現(xiàn)在朝向蝸殼出口的葉片進(jìn)口處,如圖7(a)中箭頭所指,由于這一方位的葉片通道出口阻力小,使其進(jìn)口速度增大,靜壓降低.在環(huán)形蝸殼通道中,靜壓隨蝸殼半徑的增大而增大,表明速度在迅速減小,動能轉(zhuǎn)化為壓力能.在蝸舌附近由于存在較強(qiáng)的渦流損失,使靜壓最小,接近外壁面處氣流趨于滯止,靜壓最大.通過比較這3個軸向位置的橫剖面和子午面的壓力云圖可以看出,除了葉輪入口處以外,靜壓在軸向上的差別不是很大.

圖7 輪盤側(cè)Fig.7 Hub side

圖8 半葉高處Fig.8 Half of the blade

圖9 輪蓋側(cè)Fig.9 Hubcap side

各剖面在蝸殼通道的出口處,如圖7(a)中的L處,均有一壓力較低的區(qū)域,其形成原因有兩個,一是從葉輪出來的高速氣流由于沒有蝸殼環(huán)壁的阻擋,所以,不必馬上減速擴(kuò)壓,而是將通過一個緩慢過程減速升壓,從而該處壓力相對較低;二是從蝸殼通道中沖出來的氣流由于前方阻力減小而加速,將壓力能再轉(zhuǎn)化為動能,從而降低了靜壓,如圖10所示.由于原設(shè)計蝸殼通道面積較窄小,這一效應(yīng)相當(dāng)明顯.而在改進(jìn)設(shè)計后,蝸殼通道面積變寬,氣流在蝸殼通道出口處加速不明顯,低壓區(qū)L也就縮小了(見圖11).

由于葉輪對氣流做功,全壓從葉輪進(jìn)口沿徑向不斷增加,在葉輪出口處達(dá)到最大值,然后,在氣流滲混過程中漸趨均勻,并因渦流、摩擦等損耗而有所下降,蝸舌附近氣流最紊亂,損耗較大.軸向3個剖面相比較,對于原風(fēng)機(jī),輪盤側(cè)較低,輪蓋側(cè)較高,這是因為輪盤側(cè)氣流速度高,導(dǎo)致流動損失(包括沖擊損失、摩擦損失及渦流損失等)也高.而對于改型風(fēng)機(jī),由于錐形輪蓋的阻檔,機(jī)械能都聚集在輪盤側(cè)剖面附近區(qū)域,該處全壓很高,輪蓋側(cè)剖面上大部分區(qū)域處在輪蓋外面,全壓就很低.

3.2.2 速度分布

圖10為兩個模型靠近蝸舌與蝸殼出口的部分葉柵通道,顯示輪盤側(cè)剖面的相對速度分布.在所有的葉柵通道中,這部分通道的流態(tài)通常是最差的.通過比較可以看出,改型后左邊幾個通道內(nèi)原有的渦消除了,右邊緊靠蝸舌的通道內(nèi)仍有渦,但強(qiáng)度也減弱了.

圖10 速度分布Fig.10 Relative velocity

圖11 原始模型Fig.11 Prototype model

圖12 改進(jìn)模型Fig.12 Improved model

圖11和圖12分別為兩種模型在蝸殼內(nèi)輪盤側(cè)的絕對速度分布云圖和局部矢量圖,由于葉片斜切后葉輪出口面積小,改進(jìn)模型在這個區(qū)域的速度普遍大于原始模型,但較為均勻.原模型在蝸舌附近速度的大小、方向變化很大,還有向葉輪倒流的現(xiàn)象,流態(tài)十分紊亂.而改進(jìn)模型在該處的流速分布均勻、變化平緩,沒有明顯旋渦,流態(tài)己有很大改善.

4 瞬間凍結(jié)法中葉輪方位的影響

由于蝸殼通道周向不對稱,用MRF(多重參考坐標(biāo))方法進(jìn)行動靜轉(zhuǎn)換時,要考慮葉輪方位的影響.每個主葉片通道有30°周期性變動范圍.為了考察方位角的影響,在設(shè)計流量下,以某一方位的基準(zhǔn),按0°,10°和20°這3種方位分別生成網(wǎng)格進(jìn)行模擬計算,其全壓為3 674,3 680,3 774 Pa,10°和20°方位的全壓相對于0°方位的誤差分別為1.10%和3.68%.這說明葉片的方位對計算結(jié)果存在一定的影響,但導(dǎo)致的誤差比較小.由于動、靜系統(tǒng)目前只能用近似的方法處理,所以,這樣的誤差還是可以接受的.

5 結(jié) 論

a.由于葉片寬度沿徑向逐漸減小,各半徑處的截面積接近相等,氣體在流道內(nèi)接近等速前進(jìn),從而不會產(chǎn)生明顯的逆壓梯度,抑制了氣流分離.

b.改進(jìn)風(fēng)機(jī)的全壓在設(shè)計流量處與原風(fēng)機(jī)相近,而大流量的全壓增加了,小流量的全壓則減小了,全壓隨流量的變化趨于平緩,有利于風(fēng)機(jī)運行中的工況變動.

c.改進(jìn)后全壓效率顯著提高,在設(shè)計流量時提高達(dá)9.47%,說明改進(jìn)方案是很成功的.由于設(shè)計工況下全壓變化不大,所以,效率提高主要是軸功率降低導(dǎo)致的.

[1] GU Fahua,ABRAHAM E.A numerical investigation on the volute/diffuser interaction due to the axial distortion at the impeller exit[J].ASME Journal of Fluid Engineering,2001,123(4):475-483.

[2] 朱之墀,唐旭東,李嵩,等.離心通風(fēng)機(jī)葉輪氣動設(shè)計工程方法的改進(jìn)[J].風(fēng)機(jī)技術(shù),2001(3):3-8.

[3] 李新宏,何慧偉,宮武旗,等.離心通風(fēng)機(jī)整機(jī)定常流動數(shù)值模擬[J].工程熱物理學(xué)報,2002,23(4): 453-456.

[4] 王嘉冰,區(qū)穎達(dá).柜式空調(diào)用多翼離心通風(fēng)機(jī)內(nèi)流場的數(shù)值分析[J].風(fēng)機(jī)技術(shù),2003(4):23-29.

[5] 李建鋒,呂俊復(fù).風(fēng)機(jī)流場的數(shù)值模擬[J].流體機(jī)械, 2006,34(4):10-13.

[6] 唐囡,任剛,蔡建程,等.蝸殼壁厚對離心風(fēng)機(jī)振動噪聲影響的數(shù)值研究[J].流體機(jī)械,2009,37(2): 15-18.

[7] 繆云仙,張言才.柴油機(jī)船舶機(jī)艙通風(fēng)量計算方法探討[J].江蘇船舶,2007,24(1):17-19.

[8] 寧寶東.船舶機(jī)艙通風(fēng)和風(fēng)機(jī)的選型[J].船舶, 2001,10(5):62-63.

[9] 商景泰.通風(fēng)機(jī)實用技術(shù)手冊[M].北京:機(jī)械工業(yè)出版社,2004.

[10] 王福軍.計算流體動力學(xué)分析[M].北京:清華大學(xué)出版社,2004.

[11] 成心德.離心通風(fēng)機(jī)[M].北京:化學(xué)工業(yè)出版社,2006.