渦輪轉(zhuǎn)子內(nèi)動(dòng)態(tài)燃燒模型機(jī)理探討

鄭海飛，唐 豪

(1.上海飛機(jī)設(shè)計(jì)研究院適航工程中心，上海201210；2.南京航空航天大學(xué)能源與動(dòng)力學(xué)院，南京210016)

1 引言

長(zhǎng)期以來(lái)，為提高航空發(fā)動(dòng)機(jī)推重比，減少污染物排放，降低單位推力燃油消耗率，以及拓寬發(fā)動(dòng)機(jī)穩(wěn)定工作范圍[1-3]，經(jīng)常采用的技術(shù)手段是提高燃燒室出口溫度和縮短燃燒室長(zhǎng)度。但這種技術(shù)會(huì)縮短燃燒滯留時(shí)間，形成不完全燃燒，導(dǎo)致燃燒繼而發(fā)生于渦輪內(nèi)，這與設(shè)計(jì)的初衷相違背。為了解決這一難題，利用延續(xù)至渦輪內(nèi)的高溫高壓燃?xì)庠俅稳紵岣邿崃ρh(huán)性能，業(yè)界提出了渦輪內(nèi)增燃技術(shù)，即在渦輪內(nèi)增加燃燒室，噴油再次燃燒，以提高熱力循環(huán)性能[4-6]。

國(guó)外在上世紀(jì)末就開(kāi)始了渦輪內(nèi)增燃技術(shù)的研究。Sirignano等[4-6]首先提出了渦輪內(nèi)增燃技術(shù)的概念及射流旋流結(jié)構(gòu)方案，即渦輪葉間燃燒室(TIB)。射流旋流方案是利用周向二次射流產(chǎn)生的旋流離心力進(jìn)行離心燃燒的新技術(shù)概念，在穩(wěn)定燃燒的同時(shí)提高航空發(fā)動(dòng)機(jī)的推重比。由于離心燃燒的離心加速度與燃燒環(huán)的半徑成反比，在大型燃?xì)鉁u輪發(fā)動(dòng)機(jī)上使用TIB技術(shù)時(shí)，射流旋流方案無(wú)法達(dá)到預(yù)期的性能指標(biāo)。因此，針對(duì)大型燃?xì)廨啓C(jī)的射流渦流方案應(yīng)運(yùn)而生，其應(yīng)用駐渦凹腔穩(wěn)焰技術(shù)——即利用凹腔前后壁面的二次射流在凹腔內(nèi)形成回流區(qū)，從而達(dá)到穩(wěn)定燃燒的目的。Sekar等[7-10]首次提出將駐渦凹腔(TVC)燃燒技術(shù)應(yīng)用于TIB技術(shù)中，并分別對(duì)駐渦凹腔燃燒室、帶葉片的TIB燃燒室以及帶徑向凹槽葉片的TIB燃燒室的燃燒性能進(jìn)行研究分析。結(jié)果顯示，由于增加了燃油和空氣的質(zhì)量與動(dòng)量，燃油和空氣能夠在駐渦凹腔內(nèi)充分混合；此外，駐渦凹腔內(nèi)出現(xiàn)了多個(gè)燃燒區(qū)域，這將有利于凹腔內(nèi)的空氣和燃油的燃燒中間產(chǎn)物向主流通道內(nèi)的滲透。

但目前針對(duì)射流渦流結(jié)構(gòu)的方案，缺乏將射流渦流結(jié)構(gòu)應(yīng)用于渦輪轉(zhuǎn)子方面的研究。為此，本文建立了三個(gè)模型，即原高壓渦輪轉(zhuǎn)子模型和應(yīng)用射流渦流方案的高壓渦輪轉(zhuǎn)子模型，利用數(shù)值模擬方法，采用尺度適應(yīng)模擬湍流模型對(duì)高壓渦輪轉(zhuǎn)子內(nèi)應(yīng)用射流渦輪方案進(jìn)行了仿真研究。

2 應(yīng)用射流渦流結(jié)構(gòu)的渦輪轉(zhuǎn)子幾何模型

應(yīng)用射流渦流結(jié)構(gòu)的渦輪轉(zhuǎn)子幾何模型，是在渦輪葉片底部耦合駐渦凹腔，并于凹腔前后壁面分別設(shè)置二次氣流射流孔，即渦輪內(nèi)增燃技術(shù)的射流渦流方案。二次氣流的作用是加固燃燒回流區(qū)穩(wěn)定燃燒及強(qiáng)化凹腔內(nèi)燃燒產(chǎn)物與主流燃?xì)獾膿交靃10]。如圖1所示，射流渦流方案中，渦輪燃燒室包括駐渦凹腔(TVC)、徑向葉片(RV)和徑向葉片凹槽(RVC)三個(gè)結(jié)構(gòu)，所以渦輪增燃技術(shù)的射流渦流方案可以簡(jiǎn)稱為T(mén)IB-TRC方案。

本文研究的幾何模型，包括原高壓渦輪轉(zhuǎn)子模型(model-B1)和應(yīng)用射流渦流方案的高壓渦輪轉(zhuǎn)子葉片有無(wú)徑向凹腔的模型(model-B2，model-B3)。圖2中渦輪轉(zhuǎn)子葉型為某發(fā)動(dòng)機(jī)高壓渦輪轉(zhuǎn)子葉片葉型，葉片高度由發(fā)動(dòng)機(jī)主燃室出口截面高度而定。

3 初始條件和數(shù)值模擬方案

為準(zhǔn)確模擬主流通道內(nèi)的燃?xì)馓匦裕瑴u輪轉(zhuǎn)子進(jìn)口截面的溫度分布采用渦輪導(dǎo)向器出口截面溫度分布。渦輪轉(zhuǎn)子的進(jìn)氣角度和轉(zhuǎn)速分別為68.5°、1.2×104r/min。渦輪轉(zhuǎn)子進(jìn)口截面的燃?xì)鈿怏w質(zhì)量分?jǐn)?shù)采用渦輪導(dǎo)向器出口截面燃?xì)鈿怏w質(zhì)量分?jǐn)?shù)，以提供較為符合實(shí)際的高溫高壓主流氣體特性。表1給出了進(jìn)口氣體組分。

采用基于壓力的隱式穩(wěn)態(tài)求解器模擬渦輪內(nèi)增燃燃燒室內(nèi)的流動(dòng)與燃燒過(guò)程。梯度差值方案采用基于單元體的最小二乘法插值，其適用于六面體網(wǎng)格，與基于節(jié)點(diǎn)的格林-高斯格式具有相同的精度且不會(huì)出現(xiàn)偽擴(kuò)散。壓力項(xiàng)的離散采用標(biāo)準(zhǔn)格式，壓力與速度耦合方程的求解算法為SIMPLEC。動(dòng)量方程、能量方程及湍流方程的離散均采用二階迎風(fēng)格式。數(shù)值模擬過(guò)程中應(yīng)用尺度適應(yīng)模擬湍流模型(SAS)，以上數(shù)值方法均參考文獻(xiàn)[11]～[13]。壁面函數(shù)為標(biāo)準(zhǔn)壁面函數(shù)，配合適當(dāng)網(wǎng)格，壁面Y+約為30，符合數(shù)值模擬要求[14]。

表1 進(jìn)口氣體組分Table 1 Gas component at inlet

對(duì)計(jì)算流體域進(jìn)行六面體網(wǎng)格劃分。葉片前緣、尾緣及二次流入射孔處采用O型網(wǎng)格劃分[11]。壁面處網(wǎng)格距離為0.1 mm，網(wǎng)格增長(zhǎng)率為1.2。邊界條件設(shè)置中，燃燒室主流進(jìn)口和二次流進(jìn)口均為質(zhì)量進(jìn)口，出口為壓力出口；周期邊界條件為旋轉(zhuǎn)周期邊界條件，其余均為壁面。表2列出了進(jìn)出口的流量、溫度、壓力等初始條件參數(shù)和燃油性能參數(shù)。為模擬航空發(fā)動(dòng)機(jī)實(shí)際情況，表中二次射流溫度取850 K，是因?yàn)閬?lái)自高壓壓氣機(jī)的二次射流引氣的溫度范圍為700～1 000 K。燃油噴嘴為錐形噴嘴，分布于駐渦凹腔后壁面，且各噴嘴相互獨(dú)立。

表2 初始條件和燃油性能參數(shù)Table 2 Initial conditions and fuel performance parameters

4 結(jié)果與分析

4.1 動(dòng)態(tài)燃燒對(duì)流場(chǎng)的影響

航空發(fā)動(dòng)機(jī)高壓渦輪轉(zhuǎn)子葉間應(yīng)用射流渦流方案，勢(shì)必對(duì)渦輪轉(zhuǎn)子葉間流場(chǎng)產(chǎn)生影響。圖3為主流通道內(nèi)渦輪轉(zhuǎn)子葉片徑向截面(Radial=300 mm)處，原高壓渦輪轉(zhuǎn)子葉間和應(yīng)用射流渦流方案后高壓渦輪轉(zhuǎn)子葉間流場(chǎng)的相對(duì)速度矢量分布。由圖可知，在葉片中間部位，model-B3的相對(duì)速度矢量場(chǎng)形態(tài)與model-B1和model-B2的不同。一部分流體被卷吸入model-B3的徑向凹腔后，在凹腔內(nèi)形成低速區(qū)，造成一定的流動(dòng)損失。這是因?yàn)閺较虬记怀隹谔幋嬖诩饨牵偈沽黧w在此急速轉(zhuǎn)向形成瞬間高速，進(jìn)而造成流動(dòng)損失。比較分析model-B1和mod?el-B2可發(fā)現(xiàn)，無(wú)徑向凹腔時(shí)，駐渦凹腔的增加與否對(duì)渦輪轉(zhuǎn)子葉間的相對(duì)速度矢量場(chǎng)形態(tài)影響不大，射流渦流方案的應(yīng)用只是增加了渦輪轉(zhuǎn)子葉間的氣流速度，以及渦輪轉(zhuǎn)子葉間氣流對(duì)渦輪轉(zhuǎn)子做功的能量值。

4.2 溫度場(chǎng)比較及渦輪轉(zhuǎn)子內(nèi)等溫燃燒可行性分析

圖4展示了主流通道內(nèi)渦輪轉(zhuǎn)子葉片徑向截面(Radial=300 mm)處，原高壓渦輪轉(zhuǎn)子葉間和應(yīng)用射流渦流方案后，高壓渦輪轉(zhuǎn)子葉間溫度場(chǎng)分布。由圖可看出，是否采用射流渦流方案溫度場(chǎng)的分布存在較大差異。model-B1中，主流通道內(nèi)渦輪轉(zhuǎn)子葉片進(jìn)口截面處的平均溫度值較高，且沿主流通道流動(dòng)方向單調(diào)遞減；出口截面處的溫度分布較均勻，其平均溫度低于進(jìn)口截面處的平均溫度，這是由于主流燃?xì)鈱?duì)渦輪轉(zhuǎn)子葉片做功的結(jié)果。model-B2中，主流通道內(nèi)渦輪轉(zhuǎn)子葉片進(jìn)口截面處的平均溫度略高于出口截面處的平均溫度，且主流通道內(nèi)的溫度分布基本趨于均勻化；在渦輪轉(zhuǎn)子葉片無(wú)徑向凹腔的情況下，應(yīng)用射流渦流方案可實(shí)現(xiàn)渦輪內(nèi)的近似等溫燃燒過(guò)程。model-B3中，主流通道內(nèi)的溫度分布十分均勻，渦輪轉(zhuǎn)子葉片進(jìn)出口截面處的平均溫度基本相等；在渦輪轉(zhuǎn)子葉片帶徑向凹腔的情況下，應(yīng)用射流渦流方案可基本實(shí)現(xiàn)渦輪內(nèi)的等溫燃燒過(guò)程。這是因?yàn)閺较虬记淮龠M(jìn)了主流燃?xì)馀c駐渦凹腔內(nèi)二次氣流和燃燒中間產(chǎn)物的摻混，使燃燒較為充分，二次噴油燃燒所得能量對(duì)渦輪轉(zhuǎn)子葉片做功；在部分能量消耗的前提下，原主流燃?xì)獾臏囟容^為均勻地分布于主流通道內(nèi)。

表3列出了原渦輪轉(zhuǎn)子葉片和高壓渦輪轉(zhuǎn)子采用射流渦流方案情況下，渦輪轉(zhuǎn)子葉片進(jìn)出口截面處總溫平均值和溫降值。可看出，原渦輪轉(zhuǎn)子葉片進(jìn)出口的溫降值近200 K，這是由于主流燃?xì)鈱?duì)渦輪轉(zhuǎn)子葉片做功的結(jié)果。應(yīng)用射流渦流方案后，渦輪轉(zhuǎn)子葉片進(jìn)出口的溫降值分別為60 K和25 K，基本實(shí)現(xiàn)了渦輪內(nèi)等溫燃燒過(guò)程。

表3 渦輪轉(zhuǎn)子進(jìn)出口截面總溫比較Table 3 Comparison between total temperature at inlet and outlet of turbine rotor

4.3 動(dòng)態(tài)燃燒對(duì)壓力場(chǎng)的影響

表4列出了原渦輪轉(zhuǎn)子葉片和高壓渦輪轉(zhuǎn)子采用射流渦流方案情況下，渦輪轉(zhuǎn)子葉片進(jìn)出口截面總壓平均值和落壓比。從表中可看出，射流渦流方案應(yīng)用于渦輪轉(zhuǎn)子葉片通道內(nèi)時(shí)，高壓渦輪轉(zhuǎn)子葉片的落壓比與原渦輪轉(zhuǎn)子葉片的落壓比基本相等，這說(shuō)明射流渦流方案的應(yīng)用不會(huì)對(duì)原渦輪轉(zhuǎn)子葉片的做功能力和做功效果造成十分明顯的影響。

表4 渦輪轉(zhuǎn)子進(jìn)出口截面總壓比較Table 4 Comparison between total pressure at inlet and outlet of turbine rotor

圖5為主流通道內(nèi)渦輪轉(zhuǎn)子葉片徑向截面(Ra?dial=300 mm)處，原高壓渦輪轉(zhuǎn)子葉片葉間和應(yīng)用射流渦流方案后，高壓渦輪轉(zhuǎn)子葉間的靜壓力場(chǎng)分布。從圖中可看出，是否采用射流渦流方案對(duì)靜壓力分布形態(tài)無(wú)較大差別，均是渦輪葉片壓力側(cè)為高壓區(qū)，吸力側(cè)為低壓區(qū)。不同之處在于，model-B2和model-B3渦輪葉片通道內(nèi)的靜壓力要高于mod?el-B1渦輪葉片通道內(nèi)的靜壓力，這是因?yàn)閙od?el-B2和model-B3耦合了射流渦流結(jié)構(gòu)，對(duì)主流燃?xì)獾脑俅螄娪腿紵偈沽薽odel-B2和model-B3燃燒室內(nèi)壓力升高。對(duì)比model-B2和model-B3可發(fā)現(xiàn)：射流渦流方案應(yīng)用于渦輪轉(zhuǎn)子葉片通道時(shí)，無(wú)論是否帶徑向凹腔，對(duì)渦輪轉(zhuǎn)子葉間靜壓力場(chǎng)的形態(tài)及分布均無(wú)較大影響。

5 結(jié)論

(1)是否帶駐渦凹腔對(duì)渦輪轉(zhuǎn)子葉間的相對(duì)速度矢量場(chǎng)形態(tài)影響不明顯，射流渦流方案的應(yīng)用只是增加了渦輪轉(zhuǎn)子葉間的氣流速度和渦輪轉(zhuǎn)子葉間氣流對(duì)渦輪轉(zhuǎn)子做功的能量值。

(2)渦輪轉(zhuǎn)子內(nèi)應(yīng)用射流渦流方案，主流通道內(nèi)的溫度分布較均勻(model-B2)和十分均勻(mod?el-B3)，渦輪轉(zhuǎn)子葉片進(jìn)出口截面處的平均溫度基本相等。帶徑向凹腔的情況下，應(yīng)用射流渦流方案渦輪轉(zhuǎn)子葉片進(jìn)出口的溫降值分別為60 K和25 K，基本實(shí)現(xiàn)了渦輪內(nèi)等溫燃燒過(guò)程，驗(yàn)證了渦輪內(nèi)增燃技術(shù)可實(shí)現(xiàn)渦輪內(nèi)等溫燃燒的可行性。

(3)射流渦流方案的應(yīng)用不會(huì)對(duì)原渦輪轉(zhuǎn)子葉片的做功能力和做功效果造成影響。此外，射流渦流方案應(yīng)用于渦輪轉(zhuǎn)子葉片通道時(shí)，無(wú)論高壓渦輪葉片是否帶徑向凹腔，對(duì)渦輪轉(zhuǎn)子葉間靜壓力場(chǎng)的形態(tài)及分布均無(wú)較大影響。

[1]洪 杰，陳 光.航空發(fā)動(dòng)機(jī)半個(gè)世紀(jì)來(lái)發(fā)展的回顧[J].燃?xì)鉁u輪試驗(yàn)與研究，1998，11(4)：56—60.

[2]趙堅(jiān)行.民用發(fā)動(dòng)機(jī)污染排放及低污染燃燒技術(shù)發(fā)展趨勢(shì)[J].航空動(dòng)力學(xué)報(bào)，2008，23(6)：986—996.

[3]徐華勝，鄧遠(yuǎn)灝，馬存祥.民用航空發(fā)動(dòng)機(jī)低排放燃燒室技術(shù)[J].航空科學(xué)技術(shù)，2012，(4)：5—10.

[4]Sirignano W A，Delplanque J P，Liu F.Selected challeng?es in jet and rocket engine combustion research[R].AIAA 1997-2701，1997.

[5]Sirignano W A，Liu F.Performance increases for gas-tur?bine engines through combustion inside the turbine[J].Journal of Propulsion and Power，1999，15(1)：111—118.

[6]Liu F，Sirignano W A.Turbojet and turbofan engine per?formance increases through turbine burners[R].AIAA 2000-0741，2000.

[7]Sekar B，Thornburg H J，Briones A M，et al.Effect of trapped vortex combustion with radial vane cavity ar?rangements on predicted inter-turbine burner perfor?mance[R].AIAA 2009-4603，2009.

[8]Thornburg H J，Briones A M，Sekar B.Enhanced mixing in trapped vortex combustor with protuberances part 1:single-phase nonreacting flow[R].AIAA 2011-3421，2011.

[9]Briones A M，Sekar B，Thornburg H J.Enhanced mixing in trapped vortex combustor with protuberances part 2:two-phase reacting flow[R].AIAA 2011-3422，2011.

[10]Thornburg H J，Sekar B，Zelina J，et al.Numerical study of an inter-turbine burner(ITB)concept with curved radi?al vane[R].AIAA 2007-649，2007.

[11]Menter F R，Egorov Y.Re-visiting the turbulent scale equation[C]//.Proceeding of IUTAM Symposium on One Hundred Years of Boundary Layer Research，2004.

[12]Menter F R，Egorov Y.The scale-adaptive simulation method for unsteady turbulent flow predictions.part 1:theory and model description[J].Flow Turbulence Com?bust，2010，85：113—138.

[13]Egorov Y，Menter F R，Lechner R，et al.The scale-adap?tive simulation method for unsteady turbulent flow predic?tions.part 2:application to complex flows[J].Flow Turbu?lence Combust，2010，85：139—165.

[14]ANSYS FLUENT theory guide[M].Canonsburg：ANSYS Inc.，2013.