劉震雄，胡 博

(1.上海民航職業(yè)技術(shù)學(xué)院，上海 200232; 2.上海飛機(jī)設(shè)計研究院， 上海 201210)

離心壓縮機(jī)蝸殼內(nèi)部流動特性的數(shù)值研究*

劉震雄1，胡 博2

(1.上海民航職業(yè)技術(shù)學(xué)院，上海 200232; 2.上海飛機(jī)設(shè)計研究院， 上海 201210)

針對蝸殼在提高離心壓縮機(jī)整機(jī)效率中的重要性，利用數(shù)值模擬分析了非對稱圓截面、對稱圓截面和對稱正方形截面3種不同結(jié)構(gòu)形狀的蝸殼內(nèi)部的流動特性，重點研究了蝸殼結(jié)構(gòu)形狀對蝸殼內(nèi)部的流動、靜壓、流量等參數(shù)的影響規(guī)律。綜合比較3種結(jié)構(gòu)蝸殼，非對稱圓截面蝸殼最有利于流體在蝸殼流道內(nèi)流動。

離心壓縮機(jī)；蝸殼；流動特性；數(shù)值模擬

0 引 言

離心壓縮機(jī)以其單級增壓比高、穩(wěn)定工作范圍寬、結(jié)構(gòu)簡單可靠等優(yōu)點廣泛應(yīng)用于航空航天動力系統(tǒng)，能源化工、汽車渦輪增壓系統(tǒng)等方面。離心壓縮機(jī)中的蝸殼用來收集來自葉輪的增壓氣流并引入排氣管道，同時蝸殼也充分利用氣流的動能進(jìn)一步提高氣體的壓力[1-2]。蝸殼對于整個離心壓縮機(jī)的效率及穩(wěn)定工作都有很大的影響，因此有必要對蝸殼的內(nèi)部流動進(jìn)行深入研究。

針對離心壓縮系統(tǒng)蝸殼內(nèi)部的流動特性，國內(nèi)外學(xué)者采用實驗和數(shù)值模擬等手段進(jìn)行了研究。劉偉，齊維彪，王尚錦[3]利用五孔探針對離心壓縮機(jī)蝸殼內(nèi)的三維流動進(jìn)行了詳細(xì)測量，了解其三維流動的生產(chǎn)演化細(xì)節(jié)，為蝸殼設(shè)計打下基礎(chǔ)。Hagelstein D[4]針對對稱矩形截面蝸殼進(jìn)行了數(shù)值模擬研究，主要研究蝸殼內(nèi)部的流動結(jié)構(gòu)與損失，研究表明：流量的不同對蝸殼內(nèi)部漩流結(jié)構(gòu)的位置、通流速度、流動損失都產(chǎn)生一定的影響。張克松，王桂華等[5]通過CFD方法研究了蝸殼截面變化規(guī)律對蝸殼內(nèi)部流動及出口參數(shù)的影響，通過對比分析，發(fā)現(xiàn)同一蝸殼在不同流向下，其出口速度及氣流角的變化規(guī)律曲線基本相同，但絕對值隨流量的降低而下降。周莉，張鑫[6]針對整機(jī)條件下的離心壓縮機(jī)進(jìn)行了詳細(xì)的數(shù)值研究，分析了流量變化對蝸殼內(nèi)的流動及壓力的影響，計算結(jié)果表明隨著流量的增加，蝸殼各截面中心總壓顯著增加；在蝸殼的周向方向上，蝸舌附近渦量和氣流角隨流量的變化波動較大。

目前離心式壓縮機(jī)最常用的蝸殼形式從蝸殼橫截面形狀上區(qū)分主要有矩形截面與圓形截面兩種；從蝸殼結(jié)構(gòu)上區(qū)分主要有對稱性與非對稱性兩種。

雖然到目前為止對于離心壓縮機(jī)蝸殼內(nèi)部的流動特性研究較多，但大多都考慮氣流流量的改變、葉輪與蝸殼的匹配、蝸舌等方面對其影響[7-9]，專門針對不同結(jié)構(gòu)形狀的蝸殼內(nèi)部流動的研究還很少見，筆者利用數(shù)值模擬的研究方法，針對不同結(jié)構(gòu)形狀的離心壓縮機(jī)蝸殼進(jìn)行研究，重點分析蝸殼截面形狀、蝸殼對稱性對蝸殼內(nèi)部流動的影響以及變化趨勢，有助于離心壓縮機(jī)蝸殼的選型，提高整機(jī)的性能，節(jié)約成本。

1 物理模型和計算方法

1.1 物理模型

筆者計算采用3種不同結(jié)構(gòu)的蝸殼進(jìn)行內(nèi)部流動研究，即對稱圓截面蝸殼、非對稱圓截面蝸殼和對稱正方形截面蝸殼。

蝸殼3維示意圖如圖1所示。蝸殼2維示意圖如圖2所示，其中θ為蝸殼周向角。

計算中采用的蝸殼截面面積沿周向角的變化規(guī)律基本相同，具體如圖3所示。

圖1 蝸殼3維示意圖

圖2 蝸殼二維示意圖

圖3 蝸殼截面面積沿周向角的變化規(guī)律

1．2 網(wǎng)格劃分和計算模型

數(shù)值模擬過程中，計算結(jié)果的準(zhǔn)確與否同計算模型網(wǎng)格的劃分和生成直接相關(guān)。筆者蝸殼模型網(wǎng)格劃分采用分塊技術(shù)，整個計算區(qū)域生成非結(jié)構(gòu)化六面體網(wǎng)格，并且對蝸殼壁面進(jìn)行網(wǎng)格加密處理。在計算之前，為保證計算結(jié)果的可靠性，進(jìn)行了網(wǎng)格無關(guān)性計算。網(wǎng)格數(shù)量從36萬增加到85萬，從圖4可看出，當(dāng)網(wǎng)格數(shù)目增加到70萬時，以蝸殼270°周向角截面處的面平均總壓為參考量，基本沒有變化，此時可認(rèn)為數(shù)值計算結(jié)果與網(wǎng)格數(shù)目無關(guān)。

筆者計算采用商用軟件ANSYS CFX，湍流模型采用標(biāo)準(zhǔn)k-ε雙方程湍流模型進(jìn)行計算，文獻(xiàn)[6]和[10]的研究結(jié)果驗證了采用標(biāo)準(zhǔn)k-ε湍流模型對離心蝸殼內(nèi)流動模擬的可靠性。計算結(jié)果收斂的標(biāo)準(zhǔn)是各變量殘差小于1×10- 5。

實際蝸殼內(nèi)的流動受到上游葉輪、擴(kuò)壓器的影響，蝸殼進(jìn)口呈現(xiàn)出周向非均勻性的特點?？紤]到該項重點研究蝸殼的截面形狀、對稱性對蝸殼內(nèi)部流場的影響，因此暫不引入周向非均勻的作用。蝸殼進(jìn)口邊界條件設(shè)為總壓進(jìn)口，出口邊界條件為靜壓出口，壁面采用絕熱、無滑移條件，流體采用理想氣體。

圖4 網(wǎng)格無關(guān)性實驗圖

2 計算結(jié)果分析

2.1 蝸殼結(jié)構(gòu)對流場的影響

圖5所示給出了3種蝸殼在90°、180°、270°、360°周向角截面處的速度流線圖。從圖中可看出，三種蝸殼內(nèi)都存在明顯的漩渦流動。在對稱蝸殼截面內(nèi)，如圖 5(a)、(b)所示，在流場中存在一對旋轉(zhuǎn)方向相反的對轉(zhuǎn)漩渦結(jié)構(gòu)，而在圖5(c)中，在蝸殼截面內(nèi)只存在一個明顯的漩渦結(jié)構(gòu)。

圖5 蝸殼不同位置截面流線圖

蝸殼內(nèi)部不同的漩渦結(jié)構(gòu)主要受到蝸殼進(jìn)口徑向速度和蝸殼截面形狀的影響。在對稱蝸殼內(nèi)，流體徑向流入蝸殼，受到蝸殼壁面約束速度降低，同時蝸殼為對稱形狀，流體沿蝸殼壁面均勻分流，對稱的分成兩部分。受壁面剪切力的影響，速度降低，靜壓升高，在蝸殼截面中心兩側(cè)形成低壓區(qū)，流體在壓力梯度的作用下，不斷的從高壓區(qū)域運動到低壓區(qū)域形成一對旋轉(zhuǎn)方向相反的漩渦結(jié)構(gòu)，蝸殼截面的對稱性導(dǎo)致蝸殼內(nèi)的流場基本為對稱的。

在非對稱蝸殼中，氣流切向流入蝸殼，在蝸殼的起始位置(周向角較小)，在接近蝸舌附近，截面半徑較小，流體充滿流道，蝸殼內(nèi)漩渦結(jié)構(gòu)的形成主要受到壁面剪切力的影響，隨著蝸殼周向角的增大，蝸殼徑向速度增加，當(dāng)下游徑向速度更大的流體流入蝸殼流道包圍在蝸殼內(nèi)部已存在的流體外，漩渦結(jié)構(gòu)自動形成。同時從圖5也可發(fā)現(xiàn)，蝸殼截面內(nèi)的漩渦結(jié)構(gòu)隨著通流面積的增大，漩渦的尺寸也逐漸擴(kuò)大。

2.2 蝸殼結(jié)構(gòu)對壓力的影響

圖6所示給出了計算工況下3種蝸殼內(nèi)壁面靜壓分布云圖。從圖中可看出，沿蝸殼徑向方向，靜壓隨著徑向距離的增大而增大，并在蝸殼外壁面處靜壓達(dá)到較大值；在蝸殼起始段，靠近蝸舌附近區(qū)域，流體外偏流入蝸殼，靜壓顯著升高。比較圖6(a)、(b)、(c)，沿蝸殼壁面，非對稱蝸殼的靜壓梯度變化小于對稱蝸殼，在蝸舌附近非對稱蝸殼的靜壓變化相對也較平緩，說明在蝸殼內(nèi)部采用非對稱蝸殼氣流分布更加均勻，流動更容易。

圖6 計算工況下蝸殼內(nèi)壁面靜壓分布云圖

圖7給出了計算工況下3種蝸殼分別在90°、180°、270°、360°周向角處截面的靜壓分布云圖。從靜壓云圖發(fā)現(xiàn)，蝸殼外側(cè)壁面壓力較高，由外側(cè)壁面向蝸殼內(nèi)側(cè)方向，壓力呈降低趨勢。

由于蝸殼通流速度和蝸殼環(huán)向半徑在蝸殼內(nèi)外壁面之間也形成了一個壓力梯度，由式(1)確定。

(1)[11]

式中：VT為蝸殼截面通流速度，R為蝸殼環(huán)向半徑。

因此從蝸殼流道內(nèi)壁面向外壁面逐漸升高，高壓區(qū)一直靠近蝸殼流道外壁面。

另一方面，在非對稱蝸殼中，由漩流速度引起的離心力被從蝸殼截面中心到蝸殼壁面的壓力梯度所平衡，由式(2)確定。

(2)[11]

式中：VS為蝸殼截面漩流速度，rc為渦流速度的曲率半徑。

在非對稱蝸殼截面中，壓力從蝸殼流道外壁面向蝸殼流道截面中心逐漸降低，低壓區(qū)在蝸殼流道中間，而高壓區(qū)位于蝸殼流道外壁面。

在對稱蝸殼中，氣流正向進(jìn)入蝸殼流道后直接沖擊到蝸殼壁面上，氣流速度急劇下降，動能轉(zhuǎn)換為壓力能，相對非對稱蝸殼，蝸殼外壁面壓力更高。

圖7 蝸殼不同位置截面靜壓云圖

圖8給出了蝸殼不同周向角處截面的平均靜壓值曲線圖。圖中曲線表明，對稱蝸殼的靜壓值明顯高于非對稱蝸殼的靜壓值，對稱圓截面蝸殼與對稱正方形截面蝸殼相比，靜壓值相差不大；在計算工況下，隨著蝸殼周向角的增大，蝸殼截面靜壓都隨之降低。

圖8 蝸殼不同周向角截面平均靜壓曲線圖

在本文計算工況中，3種蝸殼進(jìn)、出口采用相同的邊界條件，非對稱蝸殼截面靜壓最小，結(jié)合蝸殼內(nèi)部速度分析，表明在非對稱蝸殼內(nèi)部動能損失較小，相應(yīng)的動壓相對最大，流體在蝸殼流道中流動的動能較大，更加有利于流體流出蝸殼，而在對稱蝸殼中，蝸殼內(nèi)部的一對漩渦結(jié)構(gòu)流動損失相對較大，流體流動動能相對較小，不利于流體流出蝸殼流道。

2.3 蝸殼結(jié)構(gòu)對流量的影響

圖9所示給出了計算工況下蝸殼截面單位質(zhì)量流量沿蝸殼周向角θ的變化規(guī)律曲線圖。

從圖中可發(fā)現(xiàn)，蝸殼截面流量隨著蝸殼周向角的增大而增大，且基本呈線性變化規(guī)律，說明在蝸殼模型計算中，流體在蝸殼流道中沿周向角的流動較均勻。從圖中曲線也可發(fā)現(xiàn)，除了在蝸殼起始段位置，在蝸殼其他相同的截面位置，非對稱圓截面蝸殼的單位質(zhì)量流量最大，對稱圓截面蝸殼次之，對稱正方形截面蝸殼最小，結(jié)合前文關(guān)于蝸殼內(nèi)速度和壓力分析，說明采用切向進(jìn)氣的非對稱圓截面蝸殼，流體在流道中流動阻力小，損耗低，流體更加容易流出蝸殼。

圖9 蝸殼截面單位質(zhì)量流量隨周向角的分布曲線

3 結(jié) 論

利用數(shù)值計算的方法研究了3種不同壓縮機(jī)蝸殼內(nèi)部的流動特性。主要研究了不同蝸殼流量、內(nèi)部速度、壓力等變化的特點。得出如下結(jié)論。

(1) 受蝸殼進(jìn)口徑向速度和蝸殼截面形狀的影響，在對稱蝸殼截面內(nèi)，存在一對旋轉(zhuǎn)方向相反的漩渦結(jié)構(gòu)；在非對稱蝸殼截面內(nèi)，只存在一個漩渦結(jié)構(gòu)；

隨著蝸殼通流截面積的增大漩渦尺寸增大。

(2) 蝸殼截面靜壓由外側(cè)向內(nèi)側(cè)，呈降低趨勢。在蝸殼靠近蝸舌附近，受流體外偏流進(jìn)蝸殼的影響，靜壓顯著升高。在計算工況下，隨著蝸殼周向角的增大，蝸殼截面靜壓都隨之降低且非對稱蝸殼的靜壓值最?。粚τ趯ΨQ結(jié)構(gòu)的蝸殼，截面形狀對靜壓的影響較小。

(3) 蝸殼截面流量隨著蝸殼周向角的增大而增大，且基本呈線性變化規(guī)律；除蝸殼起始段位置，在蝸殼其他相同的截面位置，非對稱圓截面蝸殼的單位質(zhì)量流量最大。

(4) 就計算結(jié)果相比，非對稱圓截面蝸殼對流動最有利。

[1] 冀春俊,王雅君,王學(xué)軍,等.離心壓縮機(jī)排氣蝸殼內(nèi)部流動分析與優(yōu)化[J].機(jī)械工程學(xué)報,2009,45(5):312-317.

[2] 李 杜,楊 策,王一棣,等.渦輪增壓器用離心壓氣機(jī)蝸殼流場計算研究[J].內(nèi)燃機(jī)學(xué)報,2010,28(3):254-262.

[3] 劉 偉，齊維彪，王尚錦.離心壓縮機(jī)圓形截面蝸殼內(nèi)部三維流場的實驗研究[J].流體機(jī)械,2006,34(11):1-5.

[4] Hagelstein D,VandenBraembussche R,Keiper R,et al.Experimental and Numerical Investigation of the Flow in a Centrifual Compressor Volute[J].Transactions of the ASME,Journal of Turbomachinery,2000(122):28-31.

[5] 張克松,王桂華,李國祥,等.蝸殼截面變化規(guī)律對蝸殼內(nèi)流動及出口參數(shù)的影響[J].內(nèi)燃機(jī)與動力裝置,2009,112(4):20-24.

[6] 周 莉,張 鑫.離心壓縮機(jī)蝸殼內(nèi)部流動的研究[J].風(fēng)機(jī)技術(shù),2010(5):5-14.

[7] Steglich T, Kitzinger J, Seume J R. Improved Diffuser/Volute Combinations for Centrifugal Compressors[J].Journal of Turbomachinery,2008,130(011014):1-10.

[8] Semi Kim, Junyoung Park, Kukyoung Ahn.Improvement of the Performance of a Centrifugal Compressor by Modifying the Volute Inlet[J]. Fluids Eng,2010(132): 91-101.

[9] Chen Hua, Guo Strong. Numerical Simulations of Onset of Volute Stall Inside a Centrifugal Compressor . ASME Conf. Proc. 2008.

[10] 李安虎，李 紅，劉煒巍，離心泵蝸殼內(nèi)部流場的數(shù)值模擬[J].技術(shù)與應(yīng)用,2004(11):50-53.

[11] 郭 強(qiáng).帶無葉擴(kuò)壓器的離心壓縮機(jī)失速現(xiàn)象的實驗和數(shù)值研究[D].上海:上海交通大學(xué),2007.

Numerical Simulation of the Internal Flow in the Volute of Centrifugal Compressor

LIU Zhen-xiong， HU Bo

(1.ShanghaiCivilAviationCollege,Shanghai200240,China;2.ShanghaiAircraftDesignandResearchInstitute,Shanghai201210,China)

Numerical simulation was carried on three different centrifugal compressor volutes to study the effect of volute structure on the internal flow characteristics. The studies focus on the volute flow visualization, static pressure of the volute cross-section, mass flow in the volute and so on. Above all, the asymmetric circle cross-section volute was the most useful for fluid flowing in it.

centrifugal compressor; volute; flow characteristics; numerical simulation

2014-02-20

劉震雄(1981-)，男，上海人，講師，在讀博士，主要從事葉輪機(jī)械方面的研究工作。

TK474.82

A

1007-4414(2014)02-0020-04