劉偉凱，陳林泉，楊向明

(中國(guó)航天科技集團(tuán)公司四院四十一所，西安 710025)

固體燃料超燃沖壓發(fā)動(dòng)機(jī)燃燒室摻混燃燒數(shù)值研究①

劉偉凱，陳林泉，楊向明

(中國(guó)航天科技集團(tuán)公司四院四十一所，西安 710025)

根據(jù)國(guó)外研究機(jī)構(gòu)的直連式試驗(yàn)數(shù)據(jù)，設(shè)計(jì)了固體燃料超音速燃燒室模型，建立了超音速燃燒數(shù)值計(jì)算的數(shù)學(xué)模型，通過(guò)數(shù)值模擬獲得了超音速燃燒室流場(chǎng)內(nèi)的氣體狀態(tài)參數(shù)分布。結(jié)果表明，超音速燃燒室靜壓隨軸向距離的增加而逐漸降低;流場(chǎng)中心區(qū)域?yàn)榛旌铣羲倭鲃?dòng)，而后向臺(tái)階的圓周區(qū)域?yàn)閬喴羲倭鲃?dòng);燃燒效率隨軸向距離的增加而增加。

超燃沖壓發(fā)動(dòng)機(jī);固體燃料;摻混燃燒;燃燒效率

0 引言

超燃沖壓發(fā)動(dòng)機(jī)是高超音速飛行器的理想動(dòng)力裝置，固體燃料超燃沖壓發(fā)動(dòng)機(jī)以其結(jié)構(gòu)簡(jiǎn)單、易于存儲(chǔ)、成本低等優(yōu)點(diǎn)成為各國(guó)研究機(jī)構(gòu)研究的熱點(diǎn)［1-2］。固體燃料超音速燃燒室內(nèi)的燃燒流動(dòng)過(guò)程非常復(fù)雜，存在燃料分解氣化、氣化燃料與來(lái)流空氣的混合、湍流燃燒、激波及其附面層的相互干擾、流動(dòng)分離與再附著等現(xiàn)象，流場(chǎng)數(shù)值仿真能在代價(jià)相對(duì)較小情況下充分研究燃燒室燃燒流動(dòng)過(guò)程［3］。國(guó)外Jarymowycz等通過(guò)數(shù)值模擬研究，認(rèn)為燃燒室入口溫度和壓強(qiáng)對(duì)推進(jìn)劑燃速有很強(qiáng)的影響［4］;Ben-Arosh等發(fā)展了二維軸對(duì)稱計(jì)算模型，認(rèn)為燃燒室結(jié)構(gòu)對(duì)設(shè)計(jì)參數(shù)有很強(qiáng)的敏感性，不合適的結(jié)構(gòu)會(huì)導(dǎo)致熄火或熱壅塞［5-6］;Simone和Bruno對(duì)采用LiH作為超燃沖壓發(fā)動(dòng)機(jī)固體燃料的發(fā)動(dòng)機(jī)性能進(jìn)行了分析，認(rèn)為L(zhǎng)iH是一種理想的、具有潛力的高能量密度的超燃沖壓發(fā)動(dòng)機(jī)固體燃料［7］。國(guó)內(nèi)孫波等針對(duì)文獻(xiàn)中的結(jié)構(gòu)，也進(jìn)行了燃燒室冷流和反應(yīng)的流場(chǎng)分析［8］，但在燃料超燃沖壓發(fā)動(dòng)機(jī)中，氣流通道面積隨時(shí)間是不斷變化的，之前的數(shù)值研究都未能進(jìn)行非穩(wěn)態(tài)計(jì)算。

針對(duì)上述問(wèn)題，本文通過(guò)數(shù)值方法研究了固體超燃沖壓發(fā)動(dòng)機(jī)中超音速燃燒室內(nèi)的非穩(wěn)態(tài)燃燒流動(dòng)過(guò)程，以及入口氣流參數(shù)和燃燒室結(jié)構(gòu)對(duì)燃燒室燃燒效率的影響。

1 超音速燃燒室結(jié)構(gòu)

固體燃料超燃沖壓發(fā)動(dòng)機(jī)由進(jìn)氣道、超音速燃燒室、尾噴管3部分組成，超音速燃燒室是其最核心的部件，結(jié)構(gòu)簡(jiǎn)圖如圖1所示。

在超音速燃燒室內(nèi)，固體燃料分解、氣化，與進(jìn)氣道進(jìn)入的超音速氣流摻混、燃燒并釋放熱量，燃料在燃燒室中的駐留時(shí)間很短(通常小于1 ms)，導(dǎo)致?lián)交烊紵屎艿汀Ｒ虼?，熱分解的氣體燃料與來(lái)流空氣的高效摻混以及火焰的穩(wěn)定燃燒是燃燒室設(shè)計(jì)中亟需解決的關(guān)鍵技術(shù)。

圖1 固體燃料超燃沖壓發(fā)動(dòng)機(jī)簡(jiǎn)圖Fig.1 Solid fuel scramjet device

在固體燃料亞燃沖壓發(fā)動(dòng)機(jī)燃燒室中，火焰穩(wěn)定是通過(guò)入口段后的后向臺(tái)階形成凹腔產(chǎn)生回流區(qū)來(lái)實(shí)現(xiàn)的，這一方法經(jīng)過(guò)大量的理論和實(shí)驗(yàn)研究證明是合理的［9-10］，回流區(qū)長(zhǎng)度與臺(tái)階高度呈線性關(guān)系，臺(tái)階越高，火焰穩(wěn)定性越好。在超音速流動(dòng)的燃燒室中，入口氣流速度很高(通常Ma＞1)，火焰穩(wěn)定更加困難，需要更高的臺(tái)階高度，通常的突擴(kuò)結(jié)構(gòu)存在一定的局限性，臺(tái)階高度越高，將進(jìn)一步減小給定燃燒室容積內(nèi)的實(shí)際裝藥量。根據(jù)超音速流動(dòng)的特點(diǎn)，Ben-Yakar等將火焰穩(wěn)定區(qū)的突擴(kuò)結(jié)構(gòu)進(jìn)行改進(jìn)，火焰穩(wěn)定區(qū)由入口后向臺(tái)階、定截面區(qū)域和傾斜的前向臺(tái)階組成，實(shí)驗(yàn)研究表明具有一定的火焰穩(wěn)定能力，并確定了火焰穩(wěn)定的幾何參數(shù)和火焰穩(wěn)定極限?；鹧娣€(wěn)定參數(shù)為dfhLfh/，表示火焰穩(wěn)定區(qū)的相對(duì)尺寸，(dfh/dcyl)2表示流速的度量［11］。

根據(jù)文獻(xiàn)［11］的實(shí)驗(yàn)結(jié)果，文中設(shè)計(jì)了可穩(wěn)定燃燒的超燃燃燒室?guī)缀螛?gòu)型，幾何構(gòu)型簡(jiǎn)圖如圖2所示。燃燒室結(jié)構(gòu)參數(shù)為din=15 mm，dfh=49 mm，dcyl=20 mm，Lih=60 mm，得到火焰穩(wěn)定參數(shù)dfhLfh/=13.1，(dfh/dcyl)2=6，處于火焰穩(wěn)定的范圍。

圖2 燃燒室?guī)缀螛?gòu)型簡(jiǎn)圖Fig.2 Sketch of chamber geometry

2 物理模型及計(jì)算方法

固體燃料超燃沖壓發(fā)動(dòng)機(jī)工作過(guò)程如圖3所示。高溫高速空氣流經(jīng)中心進(jìn)氣通道進(jìn)入燃燒室，與固體燃料的熱解氣體在邊界層內(nèi)進(jìn)行摻混，并發(fā)生擴(kuò)散燃燒，形成一薄層擴(kuò)散燃燒火焰區(qū)域，燃燒后釋放的熱量通過(guò)對(duì)流和輻射換熱等形式，反饋給固體燃料壁面，用以維持固體燃料的熱解，使燃燒室的擴(kuò)散燃燒能夠持續(xù)進(jìn)行下去。固體燃料的熱解引起燃料表面的退移，使發(fā)動(dòng)機(jī)中氣流的通道面積不斷變化，等直段藥柱表面的退移使發(fā)動(dòng)機(jī)燃燒室靜壓降低，火焰穩(wěn)定區(qū)形成亞音速的回流，有利于燃料與空氣的摻混和燃燒。本文主要研究超音速燃燒室氣流參數(shù)的分布和入口參數(shù)對(duì)燃燒效率的影響等特性，從計(jì)算工作量考慮，假設(shè)氣流為單相流動(dòng)，傳熱中不考慮輻射的影響，各組分?jǐn)U散系數(shù)相同，不考慮重力等的影響。

圖3 燃燒室工作過(guò)程示意圖Fig.3 Sketch of combustion process

通過(guò)上述分析，建立了固體燃料超燃沖壓發(fā)動(dòng)機(jī)工作過(guò)程模擬的數(shù)值模型?？刂品匠滩捎民詈狭藙?dòng)量、能量、連續(xù)性方程以及組分運(yùn)輸方程的雷諾平均N-S方程組;湍流的模擬采用帶壁面函數(shù)的RNGk-ε模型，相對(duì)于其他模型，該模型在超音速流動(dòng)的計(jì)算方面具有更高的可靠性;化學(xué)反應(yīng)的模擬采用渦團(tuán)耗散的化學(xué)反應(yīng)模型;固體燃料為碳?xì)湄氀跬七M(jìn)劑，主要組分為HTPB和AP，壁面熱解的氣化速度通過(guò)燃料表面的氣/固相之間的質(zhì)量守恒確定;固體燃料的動(dòng)態(tài)退移過(guò)程采用彈簧光順?lè)ê途植恐貥?gòu)的動(dòng)網(wǎng)格更新方法，用UDF程序控制每一節(jié)點(diǎn)的運(yùn)動(dòng)。計(jì)算網(wǎng)格如圖4所示。

圖4 計(jì)算網(wǎng)格Fig.4 Calculation mesh

模型中考慮了 6 種 C2H4、O2、CO2、H2O、CO 和 N2，化學(xué)反應(yīng)模型采用簡(jiǎn)化的兩方程燃燒模型:

固體燃料壁面熱解的氣化速度通過(guò)燃料表面的氣/固相之間的質(zhì)量守恒確定，采用下式計(jì)算:

式中為壁面的熱流;ρf為固相燃料速度;Hv，eff為燃料的有效汽化熱。

本文假定傳熱機(jī)理的驅(qū)動(dòng)力是火焰和壁面之間的溫度差:

固體燃料表面溫度近似恒定，設(shè)為800 K，火焰溫度由軸向位置的最大溫度確定。

為了研究固體燃料超燃沖壓發(fā)動(dòng)機(jī)燃燒室結(jié)構(gòu)對(duì)燃燒效率影響，在保持入口氣流總溫1 000 K、總壓1.5 MPa、空氣流量 0.3 kg/s、燃?xì)饬魍?1.2 kg/(m2·s)的情況下，計(jì)算了不同凹腔臺(tái)階高度和不同凹腔長(zhǎng)度的影響，計(jì)算工況如表1所示。

表1 結(jié)構(gòu)變化各工況計(jì)算參數(shù)Table 1 Calculate parameter of different structure

為了評(píng)估燃燒效率與模擬的不同飛行條件下的入口氣流參數(shù)的關(guān)系，在保持燃?xì)饬魍繛?.2 kg/(m2·s)、燃燒室結(jié)構(gòu)凹腔臺(tái)階高度為17 mm和凹腔長(zhǎng)度為60 mm的情況下，分別對(duì)3個(gè)入口參數(shù)進(jìn)行了研究，包括入口總溫、空氣流量和入口總壓，計(jì)算工況如表2所示。

表2 入口參數(shù)變化各工況計(jì)算參數(shù)Table 2 Calculate parameter of different inlet condition

3 計(jì)算結(jié)果與分析

3.1 燃燒室參數(shù)分布

圖5給出了燃燒室內(nèi)各組分質(zhì)量分?jǐn)?shù)的分布。

圖5 燃燒室內(nèi)各組分質(zhì)量分?jǐn)?shù)分布Fig.5 Quality fraction contours in the chamber

由圖5可見，燃料的氣化組分C2H4和CO的質(zhì)量分?jǐn)?shù)從壁面向火焰位置逐漸減少，超過(guò)火焰位置后，其組分逐漸消失，O2的質(zhì)量分?jǐn)?shù)從中心區(qū)最大值減少到火焰表面的近似零值，燃燒的主要產(chǎn)物是CO2和H2O，擴(kuò)散火焰位于燃料壁面附近的一薄層，此處CO2和H2O的質(zhì)量分?jǐn)?shù)最大，說(shuō)明燃燒主要發(fā)生在這一區(qū)域。

圖6給出了2、5、8、10 s不同時(shí)刻燃燒室內(nèi)靜壓壓強(qiáng)分布。在2 s時(shí)，燃燒室凹腔靜壓約為0.66 MPa，沿著軸向靜壓逐漸降低，擴(kuò)張段出口處?kù)o壓最低，隨著燃面沿徑向不斷退移，等直段氣流通道面積不斷增大，燃燒室凹腔靜壓逐漸降低，工作至10 s時(shí)，燃燒室靜壓降低為 0.17 MPa。

圖6 不同時(shí)刻燃燒室靜壓分布Fig.6 Static-pressure contours at different times

圖7給出了2、5、8、10 s不同時(shí)刻燃燒室馬赫數(shù)分布。2 s時(shí)，燃燒室入口部位和出口部位的馬赫數(shù)為超音速流動(dòng)(Ma=1.9)，凹腔部分中心區(qū)域?yàn)榛旌统羲倭鲃?dòng)，流速約為Ma=1.1。隨工作時(shí)間增加，工作至10 s時(shí)，中心區(qū)域逐漸全部成為超音速流動(dòng)，凹腔部分流速達(dá)到Ma=1.8，出口馬赫數(shù)達(dá)到2.2。

圖7 不同時(shí)刻燃燒室馬赫數(shù)分布Fig.7 Mach number contours at different times

圖8給出了2、5、8、10 s不同時(shí)刻燃燒室內(nèi)靜溫分布。由圖8可見，隨工作時(shí)間增加，固體燃料表面附近的反應(yīng)區(qū)域逐漸加大。2 s時(shí)，燃燒室后向臺(tái)階處的燃燒溫度約2 059 K，燃燒室出口溫度約1 843 K，工作至10 s時(shí)，后向臺(tái)階處的燃燒溫度略有上升，達(dá)到2 139 K，燃燒室出口溫度升至1 991 K。

3.2 燃燒室結(jié)構(gòu)對(duì)燃燒效率的影響

評(píng)價(jià)燃燒室性能的主要方法之一是計(jì)算燃燒效率。在此，燃燒效率定義為已反應(yīng)的燃料質(zhì)量流率和能反應(yīng)的燃料最大質(zhì)量流率之比:

式中mf為燃料的質(zhì)量分?jǐn)?shù);mfR為已反應(yīng)的質(zhì)量分?jǐn)?shù)。

圖8 不同時(shí)刻燃燒室靜溫分布Fig.8 Temperature contours at different times

圖9給出了燃燒效率隨軸向距離的分布。從圖9可看出，燃燒效率隨軸向距離的增加而增加，由入口的65%增加到出口的83%，這是由于隨軸向距離增加燃料的穿透能力得到增強(qiáng)。

圖10(a)給出了凹腔臺(tái)階高度為分別為13、17、21 mm時(shí)燃燒效率的計(jì)算結(jié)果;圖10(b)給出了凹腔長(zhǎng)度為50、60、70 mm時(shí)對(duì)燃燒效率的影響。

圖9 燃燒效率沿軸向距離的分布Fig.9 Combustion efficiency distribution in the axis distance

圖10 凹腔高度和長(zhǎng)度對(duì)燃燒效率的影響Fig.10 Effects of height and length of cavity on combustion efficiency

由圖10(a)可見，隨臺(tái)階高度增加，燃燒效率有所增加。在凹腔部位燃燒效率增加較明顯，主要是由于增加臺(tái)階高度，將有利于渦流的形成，增強(qiáng)固體燃料熱解氣體與入口空氣的摻混作用。由圖10(b)可見，隨凹腔長(zhǎng)度增加，燃燒效率逐漸增加。這種影響在凹腔部位較明顯。隨凹腔長(zhǎng)度增加，流動(dòng)在凹腔形成較強(qiáng)的渦流。

3.3 入口氣流參數(shù)對(duì)燃燒效率的影響

圖11給出了不同入口總溫、空氣流量和入口總壓對(duì)燃燒效率影響曲線。由圖11可看出，隨入口總溫、空氣流量和入口總壓的增加，燃燒效率呈現(xiàn)下降趨勢(shì)。對(duì)于燃燒效率，隨著這些參數(shù)的增加，燃燒效率有所降低，這種現(xiàn)象是由于氣化燃料和中心氣流在很短的滯留時(shí)間內(nèi)混合不充分所致。

圖11 入口總溫、空氣流量及總壓對(duì)燃燒效率的影響Fig.11 Effects of inlet temperature ，inlet flux and inlet pressure on combustion efficiency

當(dāng)入口總溫由800 K增加到1 100 K時(shí)，燃料退移速度隨之增加，主要是由于中心氣流與燃料壁面的換熱加強(qiáng)，質(zhì)量流量增加而引起，但同時(shí)燃燒效率隨之降低。與空氣流量和入口總壓等其他參數(shù)相比，入口總溫對(duì)燃料的退移速度影響較大。同樣，與上述變化趨勢(shì)相同，隨空氣流量和入口總壓的增加，燃料的燃燒效率降低。

4 結(jié)論

(1)超音速燃燒室靜壓隨軸向距離的增加而逐漸降低，火焰穩(wěn)定區(qū)有明顯的回流。靜溫和總溫由入口段向擴(kuò)張段逐漸增加，在燃燒室中心軸線上，溫度幾乎保持不變。

(2)流場(chǎng)中心區(qū)域?yàn)榛旌铣羲倭鲃?dòng)，而后向臺(tái)階的圓周區(qū)域?yàn)閬喴羲倭鲃?dòng)，出口馬赫數(shù)達(dá)到最大值。

(3)燃燒效率隨軸向距離的增加而增加，隨入口總溫、空氣流量和入口總壓的增加，燃燒效率呈現(xiàn)下降趨勢(shì)。

［1］李巖芳，鄭凱斌，陳林泉.固體燃料超燃沖壓發(fā)動(dòng)機(jī)技術(shù)研究進(jìn)展［C］//第二屆沖壓發(fā)動(dòng)機(jī)學(xué)術(shù)會(huì)議文集，2007.

［2］劉小勇.超燃沖壓發(fā)動(dòng)機(jī)技術(shù)［J］.飛航導(dǎo)彈，2003，(2):38-42.

［3］Ben-Yakar A.Investigation of the combustion of solid fuel at supersonic conditions in a ramjet engine［D］.M.Sc Thesis，Dec.1995.

［4］Jarymowycz T A，et al.Numerical study of solid-fuel combustion under supersonic crossflows［J］.Journal of Propulsion and Power，1992，8(2):346-353.

［5］Ben-Arosh R，et al.Mixing of supersonic airflow with fuel added along the wall in a sudden expansion chamber［R］.AIAA 97-3241.

［6］Ben-Arosh R，et al.Theoretical study of a solid fuel scramjet combustor［J］.Acta Astronautica，1999，45(3):155-166.

［7］Simone D，Bruno C.Preliminary investigation on lithium hydride as fuel for solid-fueled scramjet engines［J］.Journal of Propulsion and Power，2009，25(4):875-884.

［8］孫波，武曉松，夏強(qiáng).固體燃料超燃燃燒室數(shù)值模擬［J］.氣體物理，2010，5(1).

［9］Ben-Yakar A，et al.Experimental study of a solid fuel scramjet［R］.AIAA 94-2815.

［10］Cohen B，et al.Experimental investigation of a supersonic combustion solid fuel ramjet［R］.AIAA 97-3237.

［11］Ben-Yakar A，et al.Investigation of a solid fuel scramjet combustor［J］.Journal of Propulsion and Power，1998，14(1):447-455.

Numerical study of mixing flows in a solid fuel scramjet combustor

LIU Wei-kai，CHEN Lin-quan，YANG Xiang-ming
(The 41st Institute of the Academy of China Aerospace Science and Technology Corporation，Xi＇an 710025，China)

Based on direct-connected experimental data obtained at abroad，the model of solid fuel supersonic combustion was designed，the mathematical model for reaction flow field of supersonic combustion was established，the distributions of gas parameters in the flow field were obtained through numerical calculation.Results show that supersonic combustion static-pressure reduces gradually with increase of axial distance;central area of flow-field was mixing supersonic flow，circle area of back sidestep was subsonic flow;with increase of axial distance，combustion efficiency increases.

scramjet;solid fuel;mixing combustion;combustion efficiency

V435

A

1006-2793(2012)04-0457-06

2011-10-17;

2012-03-06。

航天科技集團(tuán)科技創(chuàng)新研發(fā)項(xiàng)目。

劉偉凱(1973—)，男，高級(jí)工程師，研究方向?yàn)楣腆w火箭發(fā)動(dòng)機(jī)。E-mail:wakenliu@gmail.com

(編輯:崔賢彬)