徐興亞，唐 豪，鄭海飛

(南京航空航天大學(xué)能源與動力學(xué)院，南京 210016)

1 引言

渦輪葉間燃燒室(TIB)是利用美國空軍研究實驗室(AFRL)研究的超緊湊燃燒室(UCC)技術(shù)，通過在渦輪導(dǎo)向器葉片間安裝噴嘴，在渦輪葉片間的通道內(nèi)進行燃燒的一種新型燃燒室[1～3]?；贚ewis提出的高旋流產(chǎn)生的離心加速度提高火焰?zhèn)鞑ニ俣鹊膶嶒炑芯?，利用周向旋流燃燒代替?zhèn)鹘y(tǒng)的軸向組織燃燒，通過在渦輪通道外圍布置一圈燃燒環(huán)周向腔，將高速氣流以一定角度射入周向腔內(nèi)從而產(chǎn)生高旋流，使噴入的霧化燃料迅速蒸發(fā)并與空氣摻混，形成可燃混氣，發(fā)生燃燒。高溫燃?xì)怆x開燃燒環(huán)后，進入下游葉間通道，未燃盡燃料及其它可燃成分在通道內(nèi)繼續(xù)燃燒。高溫燃?xì)庠谕ǖ纼?nèi)與主流空氣摻混，發(fā)生動量、能量交換，組分?jǐn)U散，最后排出通道[4]。由于利用了高離心過載提高火焰?zhèn)鞑ニ俣?，同時延長火焰在燃燒室的滯留時間，不僅實現(xiàn)了在很短軸向長度內(nèi)火焰穩(wěn)定、高效燃燒，而且燃燒性能基本達到設(shè)計要求[5，6]。

目前，國外對TIB方案已進行了多項數(shù)值模擬及實驗研究。其中Anisko等[7]研究了3種燃燒環(huán)結(jié)構(gòu)分別在兩種工況下的流場及燃燒性能；Greenwood[8，9]設(shè)計了帶徑向槽結(jié)構(gòu)的TIB模型，即在每個噴油點下方的葉片上設(shè)計一個徑向槽，并發(fā)現(xiàn)徑向槽能提高氣流混合、改善出口溫度分布。國內(nèi)對TIB方案的研究，主要是數(shù)值模擬了不同渦輪導(dǎo)向器葉片凹槽結(jié)構(gòu)、二次氣入射角及當(dāng)量比等對渦輪葉間燃燒室的影響[10，11]。在燃燒環(huán)結(jié)構(gòu)方面，僅國外對幾種特定類型進行了研究，但缺乏系統(tǒng)、幾何參數(shù)化的討論。本文在Anisko研究的基礎(chǔ)上，保持燃燒環(huán)軸向?qū)挾纫欢?，僅改變徑向長度，通過數(shù)值模擬的方法，比較分析了不同長寬比燃燒環(huán)結(jié)構(gòu)對TIB模型內(nèi)部流場和總體性能參數(shù)的影響。

2 物理模型

圖1 TIB基本結(jié)構(gòu)示意Fig.1 Sketch of the turbine inter-vane burner structure

渦輪葉間燃燒室的幾何結(jié)構(gòu)如圖1所示。中心體模擬發(fā)動機軸，其與前后管組成一個環(huán)形通道，主氣流流動方向如圖1(a)中的箭頭所示。通道內(nèi)有6個周向均布的葉片，每個葉片的同側(cè)都開有徑向槽，葉片頂部是個由前后法蘭和環(huán)形蓋圍成的環(huán)形空腔，即燃燒環(huán)周向腔。環(huán)形蓋上均布6個燃油噴口和6×4個二次氣入射孔，每個燃油噴口底部都開有軸向槽。液體燃料垂直噴入燃燒環(huán)內(nèi)，二次氣45°傾斜射入燃燒環(huán)，做高速向心運動[4]。圖1(b)示出了TIB幾何結(jié)構(gòu)的剖面圖。設(shè)定燃燒環(huán)軸向?qū)挾萕為25 mm，燃燒環(huán)徑向長度L按如下變化：L/W=0.4、0.6、0.8、1.0、1.2、1.4。

3 計算方法

3.1 網(wǎng)格劃分

考慮到TIB幾何結(jié)構(gòu)的對稱性，為節(jié)省計算時間，選取1/6模型即只對60°的扇形域進行計算，如圖2所示。利用ANSYS ICEM CFD軟件劃分網(wǎng)格，模型網(wǎng)格數(shù)約為60萬。為保證網(wǎng)格質(zhì)量，在葉片周圍及二次氣入射孔進行O-grid網(wǎng)格劃分，還根據(jù)壓力、溫度、速度、近壁網(wǎng)格的y+，采用網(wǎng)格自適應(yīng)方法局部加密，最終網(wǎng)格數(shù)約為100萬。

圖2 TIB計算網(wǎng)格Fig.2 Computational grid for the turbine inter-vane burner

3.2 數(shù)值方法

數(shù)值計算采用Realizable k-ε湍流模型。液體燃料使用煤油(C12H23)，燃燒過程考慮C12H23、CH4、CO、CO2、H2、H2O、H2O(液)、O2、OH、N2及 C(s)等 11 種組分，其中CH4、C(s)只在富油燃燒時考慮。在TIB內(nèi)，空氣和燃油分別從不同入口進入，因此選擇非預(yù)混燃燒模型(平衡混合分?jǐn)?shù)/PDF模型)。油氣混合燃燒為多相流流動，采用離散相模型(DPM)，通過隨機顆粒軌道模型模擬兩相的相互作用，相間采用耦合計算。通過SIMPLEC算法處理壓力-速度耦合項。

3.3 邊界條件

如圖2所示，主流和燃燒環(huán)上二次氣進口均為質(zhì)量流量進口，兩側(cè)壁面為旋轉(zhuǎn)式周期性邊界條件，出口為壓力出口。燃料、壓力、溫度的設(shè)置參考文獻[7]，燃油流量 0.053 kg/min，油滴直徑 55 μm，初速度30.5 m/s，工況具體參數(shù)如表1所示。

表1 工況參數(shù)Table 1 Operating conditions

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

4.1 速度場

圖3給出了6種長寬比下燃油噴射口截面(x=72.5 mm)的速度矢量分布。對比發(fā)現(xiàn)，在長寬比較小時(如L/W=0.4)，燃燒環(huán)內(nèi)的速度場較紊亂，這是由于較小的燃燒環(huán)徑向長度會導(dǎo)致二次氣流在徑向更易打穿主流氣流，但在徑向槽處，具有較大徑向動量的二次氣流與進入徑向槽的主流氣流碰撞，主流氣流阻礙了二次氣流的下潛，從而導(dǎo)致部分氣流在燃燒環(huán)周向速度的引導(dǎo)下重新上行，上行氣流沖擊至燃燒環(huán)上部壁面后折返，隨后沿周向繼續(xù)流動。L/W=0.6時，由于徑向長度增加，燃燒環(huán)內(nèi)的速度場紊亂程度較L/W=0.4時有了明顯改善，燃燒環(huán)上部僅有周向流動，主流與二次氣流的相互擴散滲透主要在燃燒環(huán)中部與下部進行，上行氣流在流至燃燒中部時即在上部氣流的引導(dǎo)下周向流動。

在長寬比較大(如L/W=0.8、1.0、1.2、1.4)時，燃燒環(huán)內(nèi)部為周向流動，僅底部與主流有擴散滲透發(fā)生。但在徑向槽處，燃燒環(huán)底部氣流與主流相互作用形成一回流旋渦(圖中圓環(huán)位置)，隨著長寬比的增大，旋渦強度越來越大，渦也越來越飽滿，這增加了可燃混氣的駐留時間并有利于燃燒產(chǎn)物與主流的摻混。

4.2 溫度場

燃燒環(huán)內(nèi)的溫度分布如圖4所示?？梢?，不同長寬比的燃燒環(huán)內(nèi)部溫度分布存在明顯差異。L/W=0.4時，高溫區(qū)主要分布在燃油噴射口的軸向腔內(nèi)和主流通道下游帶徑向槽一側(cè)。這是由于斜射的二次氣流無法在具有極短徑向長度的燃燒環(huán)內(nèi)促成旋流燃燒，反而阻礙了燃油順利進入燃燒環(huán)內(nèi)燃燒，導(dǎo)致燃油在軸向腔內(nèi)富油燃燒，進而形成腔內(nèi)的高溫區(qū)；同時，未燃盡油滴與大量可燃?xì)怏w通過徑向槽引導(dǎo)進入主流通道下游后，繼續(xù)燃燒釋放熱量再次形成高溫區(qū)。L/W=0.6時，主要燃燒區(qū)域出現(xiàn)在燃燒環(huán)中部，軸向腔內(nèi)的高溫區(qū)已消失，主流下游帶徑向槽一側(cè)雖有高溫現(xiàn)象，但沒有出現(xiàn)大面積的極高溫度，溫度分布較均勻。

在長寬比較大(如L/W=0.8、1.0、1.2、1.4)時，由于燃燒環(huán)周向腔容積增大，旋流燃燒形成，燃料在環(huán)內(nèi)充分燃燒，且燃燒區(qū)域都在環(huán)的中部及下部，避免了燃燒環(huán)壁面出現(xiàn)高溫區(qū)；葉片兩側(cè)都出現(xiàn)了局部高溫，但面積較小且溫度分布均勻。

圖3 不同長寬比下燃油進口截面的速度矢量分布Fig.3 Velocity vectors of fuel inlet for turbine inter-vane burner with different length-width ratio

圖4 不同長寬比下TIB內(nèi)部的溫度分布Fig.4 Contours of static temperature along the turbine inter-vane burner flow path with different length-width ratio

燃燒室出口截面(y=170 mm)徑向平均溫度隨葉片無量綱高度的變化如圖5所示。在L/W=0.8、1.0、1.2、1.4時，溫度分布趨勢類似，都是隨著高度的增加溫度逐漸升高，并在頂層達到最高溫。這是由于燃燒環(huán)內(nèi)的燃燒主要發(fā)生在燃燒環(huán)底部，二次氣沒有起到充分摻混的作用，導(dǎo)致高溫燃?xì)庠谥髁黜攲泳奂?。在L/W=0.4、0.6時，燃燒室出口溫度分布與理想的拋物線型出口溫度曲線(中間高、兩端低)相似，且L/W=0.6對應(yīng)的曲線更理想，曲線變化平滑，高低溫差較低。這是由于燃燒發(fā)生在整個燃燒環(huán)內(nèi)(圖4(b))且燃燒充分，在主流下游并未形成集中的局部高溫區(qū)。

圖5 不同長寬比下出口截面的徑向平均溫度Fig.5 Radial temperature profile at outlet of TIB with different length-width ratio

4.3 總體性能對比

表2定量給出了6種TIB模型的燃燒效率ηb、總壓損失dp/p、CO污染物及未燃碳?xì)浠衔?UHC)排放量與長寬比的關(guān)系?？梢?，不同長寬比的燃燒環(huán)的總壓損失和燃燒效率相差不大，但污染物排放量相差巨大。這是由于長寬比不同帶來不同的燃燒區(qū)域及燃燒充分程度，從而導(dǎo)致污染物排放量差異巨大。L/W=0.6時，CO排放量最低。UHC排放量隨著長寬比的增加而減少，這是由于富油狀態(tài)下燃燒環(huán)體積增加使得UHC的燃燒更加充分。

表2 不同長寬比下燃燒環(huán)的總體性能參數(shù)Table 2 Performance parameters of TIB with different length-width ratio

5 結(jié)論

對比文中6種模型的速度場、溫度場及總體性能參數(shù)發(fā)現(xiàn)，不同長寬比的燃燒環(huán)對渦輪葉間燃燒室性能有很大影響：當(dāng)燃燒環(huán)長寬比為0.6時，出口徑向溫度分布呈拋物線型，CO排放量最低，總壓損失低，燃燒效率高，最適合渦輪葉間燃燒室設(shè)計；長寬比為0.4時，燃燒環(huán)內(nèi)摻混充分，但溫度場內(nèi)易出現(xiàn)局部高溫區(qū)，導(dǎo)致燃燒效率降低；長寬比為0.8、1.0、1.2、1.4時，燃燒環(huán)溫度分布平均，但燃燒都發(fā)生在燃燒環(huán)底部且與主流摻混不劇烈，從而導(dǎo)致出口溫度分布不理想且CO排放量較高。

另外，由于TIB燃燒室屬于超緊湊燃燒室，因此為避免徑向長度過大，應(yīng)優(yōu)先考慮較低長寬比的燃燒環(huán)。同時，不同工況(如當(dāng)量比)在同一長寬比燃燒環(huán)下對TIB的影響也必不相同，因此在實際的燃燒環(huán)長寬比設(shè)計中還需綜合考慮，使TIB性能達到最優(yōu)。

[1]Sirignano A W，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.

[2]Liu F，Sirignano W A.Turbojet and Turbofan Engine Per?formance Increases through Turbine Burners[J].Journal of Propulsion and Power，2001，17(3)：695—705.

[3]Li Q，F(xiàn)an W.Parametric Cycle Analysis of Dual-Spool Mixed-Exhaust Turbofan with Interstage Turbine Burner[R].AIAA 2007-651，2007.

[4]李 明.不同結(jié)構(gòu)導(dǎo)向器對渦輪葉間補燃室性能影響的研究[D].南京：南京航空航天大學(xué)，2012.

[5]駱廣琦，孟龍，劉琨.不同葉片徑向凹槽結(jié)構(gòu)的超緊湊型渦輪級間燃燒室數(shù)值模擬[J].空軍工程大學(xué)學(xué)報(自然科學(xué)版)，2012，13(3)：6—10.

[6]李明，唐豪，張超，等.一種新型渦輪葉間燃燒室的數(shù)值模擬[J].航空動力學(xué)報，2012，26(10)：1—8.

[7]Anisko J F，Anthenien R A，Zelina J.Numerical Investiga?tion of Cavity-Vane Interactions within the Ultra Compact Combustor[R].AIAA 2006-805，2006.

[8]Greenwood R T.Numerical Analysis and Optimization of the Ultra Compact Combustor[D].Ohio：Department of the Air force Air University，2005.

[9]Greenwood R T，Anthenien R A.Computational Analysis of the Ultra Compact Combustor[R].AIAA 2005-220，2005.

[10]莫妲，唐豪，李明，等.二次氣射流角對渦輪葉間燃燒室的影響研究[J].航空發(fā)動機，2012，38(5)：18—21.

[11]李明，唐豪，莫妲，等.當(dāng)量比對渦輪葉間燃燒室性能影響的數(shù)值模擬[J].燃燒科學(xué)與技術(shù)，2012，18(2)：161—168.