程 明，尚守堂，劉殿春，高家春，郭瑞卿，張珊珊，李 鋒

（1.北京航空航天大學(xué)能源與動力工程學(xué)院，北京 100191；2.中航工業(yè)沈陽發(fā)動機(jī)設(shè)計研究所，沈陽 110015）

0 引言

TAPS燃燒室是GE公司為下一代商用飛機(jī)研制的低污染燃燒室，經(jīng)試驗驗證其第1代燃燒技術(shù)的NOX排放比目前常規(guī)燃燒室的降低50%，且具有使NOX排放持續(xù)降低50%～75%的潛力[1]。Mongia等人采用試驗方法研究了TAPS模型燃燒室的瞬時流場結(jié)構(gòu)、主燃級和值班級火焰的相互作用，以及周期性回火現(xiàn)象等[2-4]。國內(nèi)對于TAPS燃燒技術(shù)的研究尚處于起步階段。文獻(xiàn)[5]對單環(huán)腔中心分級燃燒室進(jìn)行了數(shù)值研究，通過改變旋流器旋向和葉片角度等幾何參數(shù)，分析其對燃燒室冷態(tài)流場的影響。文獻(xiàn)[6]采用數(shù)值模擬方法初步研究了TAPS燃燒室3維二相流場，并且研究了2種燃燒室結(jié)構(gòu)和2種噴油方式下的燃燒性能，但沒有考慮燃油分級比例的影響。由于TAPS燃燒室采用了空氣徑向分級配合燃油徑向分級燃燒技術(shù)，故本研究實際涉及到空氣流量分配和燃油流量分配2個方面的內(nèi)容。

本文從TAPS燃燒室流場結(jié)構(gòu)出發(fā)，綜合分析TAPS燃燒室在不同工況下實現(xiàn)低污染及穩(wěn)定燃燒對各子燃燒區(qū)當(dāng)量比的要求，基于中等推力航空發(fā)動機(jī)設(shè)計點(diǎn)參數(shù)提出TAPS燃燒室頭部旋流器的3級旋流流量分配方案，并在該方案下，對不同燃油分配比例對TAPS燃燒室性能的影響進(jìn)行數(shù)值模擬。

1 TAPS燃燒室流場結(jié)構(gòu)和低污染燃燒原理

TAPS旋流器由值班級和主燃級構(gòu)成，2級之間由一定高度的臺階隔開，值班級由內(nèi)、外2級旋流組成，二者之間被文氏管隔開，典型TAPS旋流器結(jié)構(gòu)如圖1所示。

文獻(xiàn)[2-3]詳細(xì)介紹了TAPS燃燒室的流場特征。在高工況下，TAPS旋流器在主燃區(qū)中產(chǎn)生3個回流區(qū)，分別為主回流（PRZ）、角回流（CRZ）和臺階回流（LPZ），如圖 2所示。相應(yīng)燃燒區(qū)域被3個回流區(qū)劃分為值班級擴(kuò)散火焰區(qū)、剪切層火焰區(qū)和主燃級預(yù)混火焰區(qū)。

圖1 典型TAPS旋流器結(jié)構(gòu)

圖2 TAPS燃燒室高工況下主燃區(qū)流場結(jié)構(gòu)

文獻(xiàn)[7]對TAPS燃燒室低污染燃燒原理進(jìn)行了詳細(xì)闡述。在低工況下，僅值班級噴嘴噴射燃油。氣流在值班級2級旋流、文氏管和套筒的共同作用下，在套筒內(nèi)及其下游形成1個強(qiáng)大的值班級主流區(qū)，燃油在該區(qū)域內(nèi)進(jìn)行富油充分燃燒，其富油燃燒產(chǎn)物在剪切層及主燃區(qū)下游與主燃級新鮮空氣進(jìn)行快速摻混燃燒，以此實現(xiàn)低工況下的高效燃燒，同時獲得較低的污染物排放。TAPS燃燒技術(shù)通過增加文氏管和臺階裝置來實現(xiàn)其目的。TAPS旋流器的徑向分級結(jié)構(gòu)決定了主燃級冷的新鮮空氣射流會對值班級火焰存在焠熄作用，如果在值班級燃燒不完全的情況下，加入大量主燃級新鮮空氣，焠熄作用會導(dǎo)致燃燒室內(nèi)CO排放的大幅增加，且燃燒效率大幅下降。主燃級和值班級射流之間的臺階將2級射流隔開一段距離，延緩了主燃級射流與值班級火焰的混合，為值班級的充分燃燒提供時間；文氏管結(jié)構(gòu)將值班級分為2級旋流，燃油與值班級內(nèi)旋流混合燃燒，外旋流則被包裹在內(nèi)旋流混氣的外側(cè)，延緩內(nèi)外2級旋流的混合，以此調(diào)節(jié)值班級套筒內(nèi)燃燒區(qū)的局部當(dāng)量比，避開污染物排放量較大的當(dāng)量比區(qū)間，以實現(xiàn)低污染燃燒。在高工況下，主燃級分配了大部分的燃油和空氣，在主燃級預(yù)混腔體中形成均勻貧油混氣，以實現(xiàn)低污染燃燒；值班級分配較少的燃油和空氣，形成的擴(kuò)散火焰用于穩(wěn)定主燃級火焰。由于大部分燃油進(jìn)行貧油燃燒，故在高工況下，TAPS燃燒室能夠大幅度降低NOX的排放。

由于TAPS燃燒室的流場呈現(xiàn)出清晰的分區(qū)燃燒結(jié)構(gòu)，故在進(jìn)行TAPS燃燒室流量分配與燃油分配方案設(shè)計時，需要綜合考慮在不同工況下各子燃燒區(qū)低污染和穩(wěn)定燃燒對局部當(dāng)量比的要求。

2 TAPS燃燒室流量分配方案設(shè)計

以中等推力發(fā)動機(jī)燃燒室設(shè)計點(diǎn)參數(shù)作為設(shè)計依據(jù)（見表1），應(yīng)用TAPS燃燒技術(shù)對主燃區(qū)各子區(qū)域按低污染燃燒原理和燃燒室穩(wěn)定運(yùn)行要求進(jìn)行針對性的流量分配。

參考目前燃燒室火焰筒和渦輪冷卻水平，燃燒室流量分配方案如下：按總進(jìn)氣量的13%用作渦輪冷卻，20%用作火焰筒冷卻，10%用作頭部冷卻，余下57%的氣流為TAPS旋流器3級旋流流量。在對3級旋流進(jìn)行流量分配時，將值班級內(nèi)、外旋流看作值班級旋流處理。按設(shè)計點(diǎn)參數(shù)，可以作出2級當(dāng)量比隨燃油分級比例的變化曲線趨勢，如圖3所示。在高工況下，考慮NOX、CO排放的要求時[8]，希望主燃級和值班級當(dāng)量比位于處于0.6～0.8的貧油燃燒區(qū)間[9]；而從燃燒穩(wěn)定性來考慮，則希望貧油預(yù)混燃燒的當(dāng)量比大于0.5。從圖3中可見，當(dāng)主燃級當(dāng)量比取為0.6～0.8時，值班級當(dāng)量比在主燃級燃油比例從0%～100%變化時大致位于處于0.2～0.6區(qū)間。

表1 某中等推力航空發(fā)動機(jī)設(shè)計點(diǎn)參數(shù)

圖3 在高工況下主燃級和值班級當(dāng)量比隨燃油分級比例的變化

按照圖3中主燃級當(dāng)量比和燃油分級設(shè)計流量分配方案時，在低工況下值班級的當(dāng)量比隨燃油分級比例的變化如圖4所示。由此可見，在高工況下，按主燃級當(dāng)量比為0.6～0.8，燃油分級比例從0%～100%變化設(shè)計流量分配方案時，會造成在低工況下值班級當(dāng)量比大幅度變化。顯然，值班級當(dāng)量比過大或過小都不利于污染物排放的控制，并且不利于發(fā)動機(jī)工況的轉(zhuǎn)換。

綜合考慮圖3、4中的比例變化，在初步設(shè)計時，按在高工況下主燃級當(dāng)量比為0.65、主燃級燃油比例為90%進(jìn)行流量分配?？梢钥闯?，該設(shè)計能夠保證在高工況下低排放和穩(wěn)定燃燒要求，同時能夠適應(yīng)燃油的大幅度改變，有利于發(fā)動機(jī)不同工況的轉(zhuǎn)換。并且在低工況下值班級當(dāng)量比為3左右，不至于過大。

值班級2級旋流流量分配按氣液比為2來確定。故燃燒室頭部流量分配方案見表2。

表2 TAPS燃燒室頭部流量分配方案

3 數(shù)值模擬

3.1 物理模型的建立

采用FLUENT軟件中的Realizable k-ε湍流模型、DDM噴霧模型、PDF燃燒模型、DO輻射模型、熱力模型NO排放模型進(jìn)行計算，文獻(xiàn)[5]將該模型的單環(huán)腔燃燒室數(shù)值模擬結(jié)果與試驗結(jié)果進(jìn)行對比研究，驗證了其計算結(jié)果的可信度。為此，本文采用與其相同的數(shù)理模型對不同結(jié)構(gòu)的燃燒室性能進(jìn)行對比，以研究在相同數(shù)理模型下，選擇合適的燃燒室結(jié)構(gòu)。

3.2 計算域

TAPS燃燒室計算域如圖5所示。在研究旋流杯霧化噴嘴時采用Thru-the-Vane方法[10-11]，計算域包括旋流器上游、旋流器及火焰筒區(qū)域，不考慮冷卻的影響。進(jìn)口段直徑為122 mm的圓柱形進(jìn)氣段，主燃區(qū)為100 mm×70 mm的矩形域，并帶有收斂段。

圖5 TAPS燃燒室計算域

3.3 求解及邊界條件設(shè)置

對TAPS模型燃燒室進(jìn)行穩(wěn)態(tài)二相燃燒流場計算。采用基于壓力的隱式求解器，近壁面采用標(biāo)準(zhǔn)壁面函數(shù)，湍流模型采用標(biāo)準(zhǔn)k-epsilon模型，壓力和動量方程采用2階精度離散，采用SIMPLE算法求解。

取燃燒室高工況設(shè)計點(diǎn)參數(shù)進(jìn)行計算。進(jìn)口空氣流量為2.72 kg/s，溫度為860 K，工作壓力為3300000 Pa，總?cè)加土髁繛?.073 kg/s，在不改變幾何及計算參數(shù)的條件下，調(diào)節(jié)主燃級燃油比例從0%～100%，共11個算例。燃料為煤油，應(yīng)用非預(yù)混二相噴霧燃燒模型，采用離散相模型追蹤油珠軌跡，通過簡化PDF模型模擬湍流與燃燒的相互作用。按熱力和瞬發(fā)NOX計算污染物排放。

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

3.4.1 燃油分級比例對燃燒室速度場的影響

對比11個算例后發(fā)現(xiàn)，TAPS燃燒室的流場在整體流型上差別不大，基本與圖6類似。從圖6中可見，值班級外旋流包裹在值班級內(nèi)旋流外側(cè)，且內(nèi)旋流與外旋流共同在燃燒室內(nèi)形成較強(qiáng)的主回流區(qū)，外旋流在套筒出口臺階處，形成臺階回流。主燃級旋流在到達(dá)滯止點(diǎn)后，部分形成角回流，其余的則包裹在值班級旋流外側(cè)。

圖6 TAPS燃燒室流場

雖然燃油分級比例對燃燒室整體流型的改變不大，但是從通過3級旋流器的流量來看，燃油分級比例對流場的改變非常大。燃油分級比例對3級旋流流量的影響如圖7（b）～（d）所示。由此可見，隨著主燃級燃油比例的增大，值班級內(nèi)、外旋流流量迅速增加，主燃級流量迅速減少，在燃油比例達(dá)到60%后，變化趨勢均變緩。說明在旋流器結(jié)構(gòu)參數(shù)不變的情況下，值班級旋流的平均速度會大幅增加，主燃級速度將會減小，產(chǎn)生這一現(xiàn)象的原因應(yīng)該與燃燒熱釋放所增加的阻力有關(guān)。

3.4.2 燃油分級比例對燃燒室溫度場的影響

在主燃級燃油比例為0%、30%、70%和100%時，TAPS燃燒室Z=0截面溫度分布情況如圖8所示。由于計算的溫度場變化趨勢與燃油分級比例相關(guān)性一致，故沒有顯示全部云圖。對比圖 8（a）～（d），隨著主燃級燃油量的增加，燃燒室高溫區(qū)縮小，溫度場分布更加均勻。

燃燒室出口溫度分布系數(shù)（OTDF）如圖9所示。從圖中可見，隨著主燃級燃油分配比例的增大，分布系數(shù)快速降低，當(dāng)主燃級燃油比例增大到60%后，出口溫度分布變化趨勢減緩。值得一提的是，文獻(xiàn)[2]指出，在沒有值班級擴(kuò)散火焰時，主燃級火焰無法穩(wěn)定存在，而本文基于簡化快速PDF燃燒模型并沒有體現(xiàn)這一點(diǎn)。

圖9 燃油分級比例對出口溫度分布系數(shù)的影響

3.4.3 燃油分級比例對燃燒室排放的影響

污染物排放與燃燒室溫度存在強(qiáng)烈的相關(guān)性。高溫區(qū)減小、溫度場更加均勻必然會導(dǎo)致污染物排放大幅降低。燃燒室NOX、CO質(zhì)量分?jǐn)?shù)隨主燃級燃油比例的變化如圖10所示。從圖中可見，燃油分級比例的變化對減少污染物排放作用非常大，NOX排放降低了90%以上，而CO則降低95%。隨著主燃級燃油比例的增大，NOX、CO的排放呈遞減趨勢，且趨勢漸漸變緩。但是,在燃油比例為90%時，NOX排放最低，這與之前流量分配時設(shè)計點(diǎn)的燃油分配比例一致。再增大主燃級燃油比例時，NOX排放略有增加，而CO的排放依然為減少趨勢。

圖10 燃油分級比例對污染物排放的影響

4 結(jié)論

（1）燃油分級比例對TAPS燃燒室流場總體結(jié)構(gòu)影響不大，但是可以通過影響3級旋流流量分配，從而改變?nèi)紵业木植慨?dāng)量比，這需要對燃燒室的氣動性能進(jìn)行優(yōu)化設(shè)計。

（2）燃油分級比例對燃燒時溫度場影響巨大，增加主燃級供油能夠顯著減小高溫區(qū)，并改善溫度分布的均勻性。

（3）NOX和CO的排放隨主燃級燃油比例的增大而大幅減少，其減少幅度漸漸變緩，對于NOX存在1個最佳排放燃油比例。

[1]Mongia H C.TAPS-A 4th generation propulsion combustor technology for low emissions[R].AIAA-2003-22657.

[2]Dhanuka S K.Unsteady aspects of lean premixed-prevaporized(LPP)gas turbine combustors:flame-flameInteractions[R].AIAA-2010-1148.

[3]Dhanuka S K，Temmea J E,Driscoll J F,et al.Vortex-shedding and mixing layer effects on periodic flashback in a lean premixed prevaporized gas turbine comustor[J].Proceedings of the Combustion Institute，2009,32（2）：2901-2908.

[4]Dhanuka S K.Instantaneous flow structures in a reacting gas turbine combustor[R].AIAA-2008-4683.

[5]李鋒，尚守堂，程明，等.雙環(huán)腔燃燒室置換單環(huán)腔燃燒室可行性研究[J].航空動力學(xué)報，2008，23（1）：145-149.

[6]劉殿春，董玉璽，尚守堂，等.單環(huán)腔中心分級燃燒室流場數(shù)值模擬[J].航空動力學(xué)報，2010，25（6）：1251-1257.

[7]Hura H S.Methods and apparatus for decrease combustor emissions:USA,US6354072[P].2002-03-12.

[8]Mohammad B S,Jeng S M.Design procedures and a developed computer code for preliminary single annular combustor design[R].AIAA-2009-5208.

[9]林宇震，許全宏，劉高恩，等.燃?xì)廨啓C(jī)燃燒室[M].北京：國防工業(yè)出版社，2008：30-60.

[10]Mongia H C.Swirl cup modeling part II:inlet boundary conditions[R].AIAA-2003-1350.

[11]Mongia H C.Swirl cup modeling part III:grid indipendent solution with differentturbulentmodels[R].AIAA-2003-1349.