李 凱,曾卓雄,徐義華

(1 南昌航空大學(xué)飛行器工程學(xué)院,南昌 330063;2 上海電力學(xué)院能源與機械工程學(xué)院,上海 200090)

進氣管結(jié)構(gòu)及尺寸對旋流冷壁燃燒室性能的影響*

李 凱1,曾卓雄2,徐義華1

(1 南昌航空大學(xué)飛行器工程學(xué)院,南昌 330063;2 上海電力學(xué)院能源與機械工程學(xué)院,上海 200090)

為了考察進氣管結(jié)構(gòu)尺寸對旋流冷壁燃燒室性能的影響,采用雷諾應(yīng)力模型(RSM)對矩形進氣管、圓形進氣管和三角形進氣管的旋流冷壁燃燒室進行了數(shù)值模擬。結(jié)果表明:綜合考慮燃燒效率、總壓損失、冷壁效果以及出口溫度分布,矩形進氣管較圓形與三角形進氣管性能更優(yōu)。矩形進氣管入口長寬比對總壓損失系數(shù)的影響不大,對出口溫度分布、燃燒效率影響較大。矩形進氣管入口長寬比為0.64時,旋流冷壁室綜合性能達到最佳。

旋流冷卻;進氣管結(jié)構(gòu)尺寸;燃燒室性能;數(shù)值模擬

0 引言

燃燒室工作溫度極高,如果壁面得不到合適的冷卻,就容易變形或燒壞;同時,火焰筒壁產(chǎn)生很大的溫度梯度,引起很大的熱應(yīng)力,會使火焰筒產(chǎn)生裂紋。渦流冷卻作為一種新型燃燒室冷卻方法,能很好的實現(xiàn)冷卻壁面的目標。同時,應(yīng)用渦流冷卻技術(shù)的發(fā)動機具有成本低、維護性好以及壽命長的特點[1]。

近年來,國內(nèi)外對渦流冷卻技術(shù)在燃燒室方面的應(yīng)用開展了一些研究。唐飛等[2]開展了氧化劑噴嘴入射角優(yōu)化的仿真研究,發(fā)現(xiàn)氧化劑傾斜一定角度可以增加渦流速度與渦量,有助于提高燃燒效率。李家文等[3]研究表明,在三種不同噴嘴分布直徑的氫噴注面板下,燃燒室最高效率為97.6%,同時發(fā)現(xiàn),高溫火焰在燃燒室圓筒段59.5%半徑以內(nèi)的區(qū)域。Anderson等[4]利用示蹤粒子進行了推力室冷流實驗,得到了渦流結(jié)構(gòu)。孫得川等[5]進行了氣氧/甲烷渦流冷壁燃燒室流場與壁面耦合傳熱分析,獲得了燃燒室和噴管的壁面溫度隨時間的變化規(guī)律。李恭楠等[6]開展了渦流冷卻推力室的三維仿真研究,驗證了雙向渦旋結(jié)構(gòu),內(nèi)外渦流分界面約占燃燒室半徑的86%,燃燒區(qū)域約占推力室半徑的70%。

目前,渦流冷卻技術(shù)主要用于火箭發(fā)動機推力室,尾部多采用縮放噴管,而為了防止切向噴入的氧化劑直接向下游流向噴管,需要設(shè)計合適的噴管收斂段[7],并且大多研究只停留在內(nèi)外雙渦流結(jié)構(gòu)以及壁面冷卻效果的驗證,并未深入開展進氣管結(jié)構(gòu)尺寸對燃燒室性能影響的研究。為此,文中結(jié)合旋流分離器[8-9]提出旋流冷壁燃燒室,并針對旋流冷壁燃燒室的進氣管結(jié)構(gòu)尺寸開展了仿真研究,探究進氣管對旋流冷壁燃燒室性能的影響規(guī)律。

1 幾何模型和計算方法

1.1 幾何模型

旋流冷壁燃燒室結(jié)構(gòu)如圖1所示,頭部噴注方式如圖2所示。燃燒室兩側(cè)為切向進氣管,燃料從頭部環(huán)形孔噴入。以頭部圓面圓心為坐標原點,沿Z軸正方向建立幾何模型。具體結(jié)構(gòu)尺寸參數(shù)如下:燃燒室筒體直徑D=45 mm;矩形進氣管空氣入口邊長為a,寬為b;環(huán)形孔外徑dout=32 mm,內(nèi)徑din=30 mm;燃燒室筒體段高度H=164 mm;噴管高度T=50 mm,噴管伸入長度S=30 mm,噴管直徑D1=24 mm。其他進氣管結(jié)構(gòu)為圓形(當(dāng)前常用)、三角形(等邊),均與燃燒室筒體相切(如圖3),其中保持3種情況下入口面積A=100 mm2不變。矩形進氣管長寬參數(shù)值如表1所示。

圖1 旋流冷壁燃燒室結(jié)構(gòu)示意圖

圖2 頭部噴注方式示意圖

圖3 進氣管結(jié)構(gòu)示意圖

a/mmb/mma/bCase15200.25Case25.717.50.33Caes36.7150.45Caes4812.50.64Caes510101Caes612.581.6Caes7156.72.2Caes817.55.73.1Caes92054

2.2 控制方程與邊界條件

由于旋流冷壁燃燒室內(nèi)部流動類似旋流分離器,屬于強旋流流動,較多的研究認為雷諾應(yīng)力模型(RSM)可以較準確的描述旋流器內(nèi)的流體行為[10-11],故此,湍流模型采用雷諾應(yīng)力模型(RSM),近壁面采用標準壁面函數(shù)法,擴散項采用二階中心差分,對流項采用QUICK格式。速度-壓力耦合采用SIMPLEC方法,燃燒模型為通用有限化學(xué)反應(yīng)速率模型,湍流-化學(xué)反應(yīng)模型為渦耗散模型。來流條件采用速度入口邊界條件,空氣入口速度為20m/s;燃燒室總當(dāng)量比為0.47;出口采用壓力出口邊界條件,出口壓力設(shè)為大氣壓,壁面采用絕熱壁面,速度取無滑移條件。空氣與甲烷的進口溫度均為300K。

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

2.1 算例驗證

由于旋流冷壁燃燒室內(nèi)部流動類似旋流分離器,所以按照上述的湍流模型及邊界條件對高效Stairmand型旋風(fēng)分離器的內(nèi)部流場進行了數(shù)值模擬,模擬結(jié)果與文獻[12]中的實驗數(shù)據(jù)比較吻合(如圖4),表明文中選用的模擬方法是可信的。

圖4 Stairmand旋風(fēng)分離器在Z=0.2 m處軸向速度分布與模擬結(jié)果對比

2.2 進氣管結(jié)構(gòu)的影響

2.2.1 內(nèi)部特征

頭部區(qū)域是燃燒和壁面冷卻的重點區(qū)域,所以文中選取Z/H=25%位置處的渦旋速度沿徑向分布進行分析。由圖5知,渦旋速度從壁面向內(nèi)快速上升,隨后進入速度相對平穩(wěn)的階段,最后快速下降。內(nèi)旋流核心區(qū)域渦旋速度與距離中心的長度成正比,屬于強迫渦;而外旋流內(nèi)渦旋速度與距離中心的長度成反比,屬于自由渦,這與文獻[13-15]中所述的內(nèi)外雙層旋流可簡化為外層自由渦、內(nèi)層強迫渦的結(jié)論基本一致。另外,中間渦旋速度相對穩(wěn)定的區(qū)域是內(nèi)外旋流的過渡區(qū),稱為雙層旋流的環(huán)形空間。圖5表明,矩形進氣管條件下內(nèi)旋流核心區(qū)渦旋速度最大,圓形進氣管與之接近,而三角形進氣管則渦旋速度最小。在環(huán)形空間內(nèi),三角形進氣管渦旋速度最大,矩形進氣管最小;但是圓形和三角形進氣管條件下,環(huán)形空間內(nèi)在內(nèi)外旋流分界的位置存在渦旋速度的快速上升,而矩形進氣管則過渡平穩(wěn)。

圖5 不同進氣管結(jié)構(gòu)下Z/H=25%位置處渦旋速度沿徑向分布

圖6 不同進氣管結(jié)構(gòu)下燃燒室溫度分布圖

文中通過未過進氣口的中心對稱面上的溫度分布來觀察旋流冷壁燃燒室的壁面冷卻效果。由圖6可知,3種進氣管結(jié)構(gòu)下燃燒室均能起到冷壁的作用,壁面附近溫度均在常用火焰筒材料最大工作溫度要求的1 200 K[16]以下。其中三角形進氣管下壁面附近溫度最高,達到1 100 K,而其他情況下只有800 K左右。這主要受兩方面因素的影響:一方面,三角形進氣管下的內(nèi)旋流核心區(qū)域渦旋速度低,燃料摻混差,燃燒不充分,使得環(huán)形空間內(nèi)燃料濃度較高;另一方面,該情況下環(huán)形空間渦旋速度大,利于燃料摻混,使得該區(qū)域燃燒劇烈,這些高溫產(chǎn)物很快將熱量傳遞至緊挨環(huán)形空間的外旋流內(nèi),導(dǎo)致壁面附近溫度快速上升,破壞了冷壁效果。

2.2.2 總壓損失系數(shù)

總壓損失是燃燒室性能參數(shù)之一,代表了燃燒室設(shè)計的好壞。為考察旋流冷壁燃燒室總壓的損失情況,引入總壓損失系數(shù),其計算公式為:

(1)

由式(1)計算得矩形、圓形、三角形進氣管條件下總壓損失系數(shù)分別為0.6%、0.64%、0.58%,3種情況下總壓損失系數(shù)都很小,且相差不大,說明進氣管結(jié)構(gòu)對總壓損失系數(shù)的影響不大。

2.2.3 燃燒效率

燃燒效率亦為燃燒室性能參數(shù)之一,文中燃燒效率計算采用燃氣分析法,計算公式為:

(2)

式中:EI為污染排放物指數(shù)(g/kg),LHV為低位熱值,下標x表示各污染排放物,f表示燃料。

由式(2)計算得矩形、圓形、三角形進氣管下的燃燒效率分別為97.3%、97%、96.2%。由于離心力的作用,密度較小的甲烷聚集在燃燒室中心,而渦旋速度的提高有利于燃料的摻混,使燃燒更充分,因此內(nèi)旋流核心區(qū)域的渦旋速度是影響燃燒效率的關(guān)鍵因素。矩形進氣管下的內(nèi)旋流核心區(qū)域渦旋速度最大,所以燃燒效率最高;三角形進氣管下的內(nèi)旋流核心區(qū)域渦旋速度最低,所以燃燒效率最低。圓形進氣管下由于具有較高的環(huán)形空間渦旋速度,從而改善了內(nèi)旋流核心區(qū)域渦旋速度低引起的燃燒效率低的問題,使其具有與矩形進氣管接近的燃燒效率。

2.2.3 出口溫度分布系數(shù)

燃燒室出口溫度的品質(zhì)將影響渦輪葉片的壽命,文中采用出口溫度分布系數(shù)(OTDF)來衡量出口溫度分布的質(zhì)量,其計算公式為:

(3)

式中:T4max表示燃燒室出口最高溫度;T4ave表示燃燒室出口平均溫度;T3ave表示燃燒室入口平均溫度。

由式(1)計算得矩形、圓形、三角形進氣管條件下OTDF值分別為0.32、0.48、0.36。3種情況下OTDF有一定差距,其中只有矩形進氣管才滿足航空發(fā)動機燃燒室OTDF值在0.25～0.35之間的要求[11];而圓形進氣管的OTDF達到0.48,出口溫度分布較不均勻,將會影響燃燒室后續(xù)部件的壽命。綜上所述,矩形進氣管條件下,旋流冷壁燃燒室具有相對最高的燃燒效率,較低的總壓損失以及相對最佳的出口溫度分布。

2.3 矩形進氣管長寬比的影響

2.3.1 內(nèi)部特征

圖7為不同入口長寬比下Z/H=25%位置處的切向速度沿徑向分布,為使圖形清晰選取了部分具有代表性的情況。由7圖可知,內(nèi)旋流核心區(qū)、環(huán)形空間的渦旋速度隨著長寬比的增加均呈下降趨勢,其中長寬比小于1時的內(nèi)旋流核心區(qū)渦旋速度相差不大,長寬比大于1時的內(nèi)旋流渦旋速度下降明顯。

不同長寬比條件下的燃燒室均具有很好的冷壁效果。因為冷壁效果受多個因素的綜合影響:長寬比小于1時,內(nèi)旋流核心區(qū)域渦旋速度高,摻混效果好,燃燒充分,使得環(huán)形空間內(nèi)燃料濃度低,即使擁有較高的環(huán)形空間渦流速度也沒有破壞冷壁效果;而長寬比大于1時,內(nèi)旋流核心區(qū)域渦旋速度低,摻混效果差,燃燒不充分,使得環(huán)形空間內(nèi)燃料濃度高,但是環(huán)形空間渦流速度低,燃燒不劇烈,所以仍保持了好的冷壁效果。

圖7 不同入口長寬比下Z/H=25%位置處的渦旋速度沿徑向分布

2.3.2 總壓損失系數(shù)

由圖8可知,總壓損失系數(shù)隨著長寬比的增加而略有下降,但總體數(shù)值較小,說明旋流冷壁燃燒室具有低阻特性。長寬比增大表明進氣管在高度方向上所占的比例越大,空氣進入燃燒室沿壁面向下旋轉(zhuǎn)的圈數(shù)就越少,經(jīng)歷的路程就越短,從而減小了氣流與壁面產(chǎn)生的摩阻損失,所以長寬比較大時總壓損失系數(shù)有所下降。

圖8 不同入口長寬比下總壓損失系數(shù)

2.3.3 燃燒效率

從圖9可以發(fā)現(xiàn),在長寬比小于1時,燃燒效率接近,在97%上下;而當(dāng)長寬比大于1時燃燒效率下降明顯,最低為95.4%。結(jié)合渦旋速度的分析可知,長寬比小于1時,影響燃燒效率的旋流核心區(qū)的渦旋速度相差不大,所以燃燒效率比較接近。然而,長寬比大于1時,內(nèi)旋流核心區(qū)渦旋速度較快下降,導(dǎo)致了燃燒效率的較快下降。

圖9 不同入口長寬比下燃燒效率

圖10 不同入口長寬比下的OTDF值

2.3.4 出口溫度分布系數(shù)

在圖10中,OTDF呈現(xiàn)“W”形狀,其最小值出現(xiàn)在長寬比為1的兩側(cè),為0.32,而且只有該兩種情況下滿足了OTDF小于0.35的要求。其他情況下出口溫度分布較不均勻,增加了燃燒室后續(xù)部件工作環(huán)境的惡劣程度。綜上分析,矩形進氣管在長寬比為0.64的情況下,具有最高的燃燒效率,較低的總壓損失,良好的壁面冷卻效果以及最佳的出口溫度分布,是較為合理的進氣管結(jié)構(gòu)與尺寸。

3 結(jié)論

文中開展了不同進氣管結(jié)構(gòu)及矩形進口尺寸對旋流冷壁燃燒室影響的熱流數(shù)值研究,得到了以下結(jié)論:

1)矩形進氣管較圓形與三角形進氣管有最高的燃燒效率,低的總壓損失,良好的冷壁效果以及最佳的出口溫度分布。

2)入口長寬比對總壓損失系數(shù)的影響不大,但對出口溫度分布影響較大;矩形進氣管入口長寬比小于1時,其對燃燒效率影響不大;矩形進氣管長寬比大于1時,其對燃燒效率的影響較為明顯。

3)矩形進氣管入口長寬比為0.64時,旋流冷壁室綜合性能達到最佳,燃燒效率為97.3%,總壓損失系數(shù)為0.6%,OTDF為0.32。

[1] Chiavertini M J, Malecki M J, Sauer J A, et al. Vortex combustion chamber development for future liquid rocket engine applications, AIAA 2002-4149 [R]. 2002.

[2] 唐飛, 李家文, 常克宇. 渦流冷卻推力室中渦流結(jié)構(gòu)的分析與優(yōu)化 [J]. 推進技術(shù), 2010, 31(2): 165-169.

[3] 李家文, 王化余, 葉漢玉, 等. 渦流冷卻推力室燃燒效率分析 [J]. 推進技術(shù), 2013, 34(11): 1507-1512.

[4] Anderson M, Rom C, Bonazza R, et al. Vortex chamber flow field characterization using laser diagnostics [C]∥ Las Vegas:52nd JANNAF propulsion meeting,2004.

[5] 孫得川, 楊建文, 白榮博. 氣氧/甲烷渦流冷壁燃燒室流場與壁面耦合傳熱分析 [J]. 推進技術(shù), 2011, 32(3): 401-406.

[6] 李恭楠, 俞南嘉, 路強. 渦流冷卻推力室流場特征與性能仿真 [J]. 航空動力學(xué)報, 2014, 29(2): 420-426.

[7] 白榮博. GOX/CH4渦流冷壁燃燒室點火過程數(shù)值模擬 [D]. 西安: 西北工業(yè)大學(xué), 2011.

[8] 郭雪巖, 王斌杰, 楊帆. 水力旋流器流場大渦模擬及其結(jié)構(gòu)改進 [J]. 排灌機械工程學(xué)報, 2013, 31(8): 696-701.

[9] Martinez L F, Lavin A G, Mahamud M M, et al. Vortex finder optimum length in hydrocyclone separation [J]. Chemical Engineering and Processing, 2008, 47(2): 192-199.

[10] 孫得川, 白榮博, 劉上. 渦流燃燒發(fā)動機燃燒室數(shù)值模擬 [J]. 彈箭與制導(dǎo)學(xué)報, 2011, 31(2): 111-113.

[11] 王海剛, 劉石. 不同湍流模型在旋風(fēng)分離器三維數(shù)值模擬中的應(yīng)用和比較 [J]. 熱能動力工程, 2003, 18(4): 337-342.

[12] Slack M D, Prasad R O, Bakker A, et al. Advances in cyclone modeling using unstructured grids [J]. Chemical Engineering Research and Design, 2000, 78(8): 1098-1104.

[13] Vyas A B, Majdalani J, Chiaverini M J. The bidirectional vortex: Part 1 an exact inviscid solution, AIAA 2003-5052 [R]. 2003.

[14] Vyas A B, Majdalani J, Chiaverini M J. The bidirectional vortex: Part 2 viscous core corrections, AIAA 2003-5053 [R]. 2003.

[15] Vyas A B, Majdalani J, Chiaverini M J. The bidirectional vortex: Part 3 multiple solution, AIAA 2003-5054 [R]. 2003.

[16] 林宇震, 許全宏, 劉高恩. 燃氣輪機燃燒室 [M]. 北京: 國防工業(yè)出版社, 2008.

Effect of Inlet Duct on Performance of Swirl-cooled Combustor Chamber

LI Kai1,ZENG Zhuoxiong2,XU Yihua1

(1 School of Aircraft Engineering, Nanchang Hangkong University, Nanchang 330063, China;2 College of Energy and Mechanical Engineering, Shanghai University of Electris Power，Shanghai 200090, China)

In order to study effect of inlet duct on performance of swirl-cooled combustor chamber, Reynolds stress model (RSM) was used to carry out simulations of swirl-cooled combustor chamber with different inlet duct structures (rectangular inlet, circular inlet and triangular inlet). The results show that considering combustion efficiency, total pressure loss, cooling effect of wall outlet temperature distribution, the performance of rectangular inlet duct is better than circular or triangular inlet duct. Entrance aspect ratio has little effect on total pressure loss coefficient, but has greater impact on outlet temperature distribution and combustion efficiency. Under conditions of this paper, when the aspect ratio of the rectangular inlet duct gets to 0.64, the comprehensive performance of swirl-cooled combustor reaches optimal.

swirl-cooled; structures and size of inlet duct; performance of combustor chamber; numerical simulation

2014-07-02

國家自然科學(xué)基金(51066006;51266013);航空科學(xué)基金(2013ZB56002;2013ZB56004)資助

李凱(1991-),男,浙江嘉興人,研究生,研究方向:航空工程。

V434

A