侯凌云 楊 縉 馬雪松 劉 巍

(1清華大學(xué)航天航空學(xué)院,北京 100084; 2北京動(dòng)力機(jī)械研究所,北京 100074)

空氣污染各組分對(duì)甲烷超聲速燃燒性能的影響

侯凌云1,*楊 縉1馬雪松2劉 巍2

(1清華大學(xué)航天航空學(xué)院,北京 100084;2北京動(dòng)力機(jī)械研究所,北京 100074)

在與甲烷詳細(xì)化學(xué)反應(yīng)機(jī)理對(duì)比驗(yàn)證基礎(chǔ)上,采用18組分24步簡(jiǎn)化反應(yīng)機(jī)理模擬甲烷超聲速燃燒過程,從化學(xué)動(dòng)力學(xué)和熱力學(xué)角度用數(shù)值方法研究了乙醇燃燒加熱空氣中的七種主要污染組分(H2O,CO2,O, OH,CO,H,H2)對(duì)甲烷超聲速燃燒性能的影響.分析結(jié)果表明：在一定條件下,進(jìn)口空氣中污染組分H2O的增加造成平均比熱容增加,總溫降低,并作為第三體抑制甲烷的燃燒過程,使超燃室的性能下降;CO2因大分子量特性使燃?xì)馄骄肿恿吭龃?降低超燃室做功能力,H2O和CO2兩組分對(duì)甲烷超燃性能都起消極作用;污染組分自由基H、O、OH和燃燒中間產(chǎn)物CO、H2使燃燒室燃燒效率上升,對(duì)甲烷超燃性能起積極作用.

超聲速燃燒;化學(xué)反應(yīng)機(jī)理;污染空氣;組分

超燃沖壓發(fā)動(dòng)機(jī)地面實(shí)驗(yàn)中要再現(xiàn)飛行時(shí)的高速高焓氣流特性,必須對(duì)進(jìn)口氣體進(jìn)行預(yù)加熱以使總焓達(dá)到高飛行馬赫數(shù)時(shí)氣流總焓.目前,主要采用電弧加熱和燃燒加熱兩種方式對(duì)空氣進(jìn)行加熱,以模擬高空高焓氣流條件.在工程應(yīng)用中,常使用比電加熱更經(jīng)濟(jì)的燃燒加熱方式來模擬高空高速、高焓進(jìn)口氣流.然而,燃燒加熱空氣與實(shí)際空氣組分組成、比熱、分子量等的差異會(huì)使地面實(shí)驗(yàn)結(jié)果不能等同于飛行實(shí)驗(yàn)結(jié)果.于是,如何使用污染空氣來模擬和復(fù)現(xiàn)超聲速燃燒飛行狀態(tài)下的實(shí)際性能,研究燃燒加熱污染空氣對(duì)超聲速燃燒的影響就成為燃燒實(shí)驗(yàn)技術(shù)需要解決的重要課題.

國外加熱空氣污染物影響超音速燃燒的研究始于1960-1970年,人們從理論、實(shí)驗(yàn)和數(shù)值模擬來研究燃燒加熱后的污染空氣對(duì)超音速燃燒過程的影響[1-3].但其發(fā)表的文獻(xiàn)大多是針對(duì)氫超燃過程影響的分析.國內(nèi)也在近幾年展開了相關(guān)研究[4-6],但大多側(cè)重在對(duì)最終產(chǎn)物影響的計(jì)算和實(shí)驗(yàn).本文針對(duì)乙醇燃燒加熱空氣中的主要污染物,包括最終產(chǎn)物H2O、CO2,自由基O、OH、H,中間產(chǎn)物CO、H2,對(duì)甲烷超聲速燃燒性能影響展開數(shù)值模擬.

1 計(jì)算模型與方法

1.1 燃燒室模型

數(shù)值計(jì)算所使用的燃燒室模型如圖1所示,為具有中心支板和兩對(duì)凹槽的突擴(kuò)型雙模態(tài)超燃沖壓發(fā)動(dòng)機(jī)燃燒室[7].此模型的支板燃料注入系統(tǒng)能夠使燃料在高速氣流中均勻分布,并在很短的時(shí)間內(nèi)完成混合,凹槽火焰穩(wěn)定器能對(duì)點(diǎn)火和穩(wěn)定燃燒起積極作用.燃料從燃燒室壁面距凹槽15 mm處徑向垂直噴入,以及支板尾緣順向噴入主流.由于結(jié)構(gòu)的對(duì)稱性,取其一半進(jìn)行計(jì)算.采用結(jié)構(gòu)化網(wǎng)格,對(duì)支板、凹槽及壁面進(jìn)行了局部加密,網(wǎng)格數(shù)為15265.

1.2 控制方程與模型

采用包含多種組分帶化學(xué)反應(yīng)的雷諾平均、守恒型Navier-Stokes方程作為氣體湍流流動(dòng)、燃燒控制方程:

式中U為守恒變量向量,F、G為對(duì)流項(xiàng)向量,Fv、Gv為粘性項(xiàng)向量,H為源項(xiàng)向量.其中:

其中ρ為密度,u、v為軸向和徑向速度分量,E、p為內(nèi)能和壓力,Yi為i組分的質(zhì)量分?jǐn)?shù),τij和qi為粘性應(yīng)力分量和反應(yīng)吸放熱量,Sd,h為化學(xué)反應(yīng)作用源項(xiàng),ωi為組分i的質(zhì)量生成率.以SST(shear-stress transport) k-ω湍流模型、空間反應(yīng)渦耗散概念燃燒模型封閉雷諾平均方程,對(duì)超聲速流和化學(xué)反應(yīng)進(jìn)行耦合計(jì)算.

1.3 化學(xué)反應(yīng)動(dòng)力學(xué)模型

本文參考Peters等[8]的甲烷14組分18反應(yīng)模型和董剛等[9]的甲烷Y簡(jiǎn)化機(jī)理模型,在考慮氮氧化物生成的條件下,構(gòu)造而得甲烷18組分24反應(yīng)簡(jiǎn)化反應(yīng)機(jī)理模型.如下表1所示.

此甲烷模型中R1-R2為鏈引發(fā);R3-R8為鏈分支;R9-R12為鏈傳遞;R13-R19為鏈終止[9],其余為氮氧化物生成反應(yīng).利用Chemkin4.0中的一維對(duì)沖擴(kuò)散火焰模型對(duì)甲烷18組分24步簡(jiǎn)化反應(yīng)機(jī)理與甲烷53組分325步詳細(xì)反應(yīng)機(jī)理(主要取自GRIMECH2.11[10])進(jìn)行對(duì)比、驗(yàn)證.兩者的溫度分布和主要組分分布比較見圖2、圖3和圖4.

表1 甲烷燃燒化學(xué)反應(yīng)動(dòng)力學(xué)模型Table 1 Chemical reaction mechanism of methane combustion

圖2、圖3、圖4中橫坐標(biāo)d代表從燃料噴口到氧化劑噴口的距離,圖3和圖4左邊的縱坐標(biāo)代表O2和CH4的摩爾分?jǐn)?shù),右邊的縱坐標(biāo)x代表H2、H2O、HO2、CH3、CO和CO2的摩爾分?jǐn)?shù).可以看出,兩者的溫度分布和主要組分分布基本上是一致的,說明本文構(gòu)造的甲烷18組分24步簡(jiǎn)化反應(yīng)機(jī)理與53組分水嶺325步詳細(xì)反應(yīng)機(jī)理較好吻合,能夠反映出各成分的變化規(guī)律.

1.4 邊界條件與進(jìn)口空氣組分

為獲得超燃室入口條件,首先需對(duì)乙醇加熱燃燒器展開其研究.對(duì)于乙醇加熱器而言,其入口條件為:空氣質(zhì)量流量為0.067 kg·s-1,氧氣的質(zhì)量流量為0.016 kg·s-1,乙醇的質(zhì)量流量為0.0077 kg·s-1,燃燒完成后,保證出口燃?xì)庵醒鯕獾哪柗謹(jǐn)?shù)仍為21%,燃?xì)饪倻貫?644 K.采用Chemkin和Lawrence實(shí)驗(yàn)室的68組分388步反應(yīng)乙醇燃燒機(jī)理[11]對(duì)乙醇燃燒加熱器燃燒過程進(jìn)行模擬,得到加熱器出口(超燃室入口)主要組分含量,出口燃?xì)庵谐薔2、O2外,含量較多的7種組分分別是H2O、CO2、O、OH、CO、H、H2.其中H2O和CO2含量占的比重非常大,摩爾分?jǐn)?shù)分別在7%和10%左右;O、OH、CO含量在同一個(gè)數(shù)量級(jí)10-3;H、H2的含量在同一個(gè)數(shù)量級(jí)10-4.其余組分如C2H4、C2H2、CH4、CH3、CH2O、HO2等,摩爾分?jǐn)?shù)數(shù)量級(jí)在10-6以下,這些組分不會(huì)對(duì)超聲速燃燒造成明顯影響,因此主要考慮7種組分即最終產(chǎn)物(H2O、CO2)、自由基(O、OH、H)和中間產(chǎn)物(CO、H2).

空氣入口采用壓力入口,總壓強(qiáng)732152.8 Pa,靜壓強(qiáng)42852.9 Pa,總溫度1644.15 K.燃料入口為壓力入口,以聲速噴入常溫(300 K)甲烷,其流量按燃燒當(dāng)量比0.723給定,第一個(gè)凹槽前沿與支板尾緣燃料噴注比為6:1.出口為壓力出口,背壓為大氣環(huán)境壓力101325 Pa,總溫300 K,壁面為無滑移絕熱壁面.表2-4為分別改變H2O、CO2、O、OH、CO、H、H2的摩爾分?jǐn)?shù),與O2、H2O和N2(剩余含量計(jì)入N2組分中)構(gòu)成不同進(jìn)口氣流工況.其中每組工況中帶“#”號(hào)的為基準(zhǔn)工況(電加熱空氣工況),帶“*”號(hào)的工況中的單個(gè)組分含量值為乙醇燃燒加熱污染空氣中該組分的含量.當(dāng)七種組分中每種組分摩爾含量改變,而其余成分不變時(shí)(指O2和H2O),研究各組分對(duì)超燃性能的影響.

1.5 性能參數(shù)定義[1,12]

(1)燃燒效率

本文選擇基于總溫變化來計(jì)算燃燒效率.定義如下:

其中質(zhì)量平均總溫Tt定義為

Tt=∫ρuTt,AdA/∫ρudA

絕對(duì)火焰溫度T0,ideal為

其中,T0為入口的質(zhì)量平均總溫,ΔH0b為定壓標(biāo)準(zhǔn)燃燒熱;nj、Cp,m,j分別為燃燒產(chǎn)物中第j組元的摩爾數(shù)及摩爾定壓熱容,pr為生成物.

(2)燃燒室參考比推力

推力由壓差和動(dòng)量差兩部分組成,定義如下:

比推力f定義為超燃室參考推力與燃?xì)赓|(zhì)量流量之比,即 f=F/m.其中F為超燃室參考推力,m為燃?xì)赓|(zhì)量流量.

(3)出口截面燃?xì)馄骄▔罕葻崛莺推骄肿恿?

其中,ni為第i個(gè)組分的摩爾分?jǐn)?shù),Cp,i為第i個(gè)組分的定壓比熱容.Mi為第i個(gè)組分的分子量.

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

2.1 污染組分最終產(chǎn)物H2O的作用

污染組分最終產(chǎn)物H2O的工作狀況設(shè)置見表2,結(jié)果分析見圖5、圖6.隨著進(jìn)口空氣中的水含量的增多,超燃室的燃燒效率、比推力均明顯下降.

從化學(xué)反應(yīng)動(dòng)力學(xué)模型可以看出:H2O在鏈引發(fā)階段通過反應(yīng)CH3+H2O→CH4+OH發(fā)揮重要作用;在鏈分支階段,H2O在反應(yīng)HCO+H2O→CH2O+OH中也非常關(guān)鍵;在這兩個(gè)反應(yīng)中,H2O的增加促進(jìn)了自由基生成,對(duì)燃燒著火過程起到了一定促進(jìn)作用,但作用有限;在一定的條件下,H2O作為第三體在反應(yīng)H+O2+M→HO2+M和H+H+M→H2+M中碰撞頻率因子較高,造成自由基復(fù)合,因此對(duì)燃燒起負(fù)面作用.另外,由圖7可見,當(dāng)水含量增加,平均比熱容會(huì)增加,各截面的質(zhì)量平均溫度下降,也是導(dǎo)致燃燒效率下降的一個(gè)重要原因.當(dāng)燃燒效率下降,會(huì)使燃?xì)庾龉δ芰ο陆?從而使燃燒室比推力下降.

表2 含H2O污染空氣組分Table 2 Vitiation cases containing H2O

表3 含CO2污染空氣組分Table 3 Vitiation cases containing CO2

根據(jù)熱力學(xué)分析,由圖7可見,平均分子量隨H2O含量的增加而減小,一定溫度條件下,燃?xì)馄骄肿恿繙p小應(yīng)該促使燃?xì)庾龉δ芰υ鰪?qiáng),也就是說從熱力學(xué)角度分析進(jìn)口空氣中H2O含量增加會(huì)使比推力增大.但由于分子量的增加有限,該影響較小.

綜合化學(xué)動(dòng)力學(xué)和熱力學(xué)分析,定性地說明H2O對(duì)甲烷超燃過程的影響在化學(xué)能轉(zhuǎn)化為熱能階段起的消極作用強(qiáng)于在熱能轉(zhuǎn)化為動(dòng)能階段起的積極作用,從而降低了甲烷超燃性能.

2.2 污染組分最終組分CO2的作用

污染組分最終產(chǎn)物CO2的工況設(shè)置見表3,結(jié)果分析見圖8、圖9.從圖8、圖9、圖10可以看出,燃燒效率未隨CO2含量增加而明顯變化,而比推力明顯下降了,同時(shí)燃?xì)馄骄肿恿吭龃罅?從化學(xué)反應(yīng)動(dòng)力學(xué)角度看,CO2屬于比較惰性的最終產(chǎn)物,它對(duì)燃燒過程的影響是有限的.污染物CO2因其大分子量特性使燃?xì)馄骄肿恿吭龃?從而使燃?xì)庾龉δ芰ο陆?這是CO2對(duì)超燃室燃燒性能起消極作用的主要原因.

2.3 污染組分自由基H、O、OH的作用

污染組分自由基H、O、OH的工況設(shè)置如表4,結(jié)果分析如圖11,隨著進(jìn)口空氣中的自由基H、O和OH含量的增多,超燃室的燃燒效率均上升.

從化學(xué)反應(yīng)動(dòng)力學(xué)角度,根據(jù)本文化學(xué)反應(yīng)動(dòng)力學(xué)模型,說明自由基O主要是通過鏈分支反應(yīng)CH3+O?CH2O+H和鏈傳遞反應(yīng)O+H2O?OH+OH來促進(jìn)燃燒,從而提升超燃室性能;自由基OH是非常活躍的組分,它主要通過促進(jìn)鏈引發(fā)反應(yīng)CH4+ OH?CH3+H2O和鏈分支反應(yīng)CH2O+OH?HCO+H2O、CO+OH?CO2+H來促進(jìn)燃料燃燒過程,從而使燃燒效率、比推力上升;自由基H主要通過鏈終止反應(yīng)H+OH+M?H2O+M和CH3+H+M?CH4+M發(fā)生作用,它與甲烷燃燒過程形成了競(jìng)爭(zhēng),抑制了甲烷的燃燒,但由于自由基H本身作為相當(dāng)活躍的燃料燃燒而使超燃室性能提升.

自由基O、OH、H促進(jìn)超聲速燃燒的能力是不一樣的.O的摩爾分?jǐn)?shù)從0到0.00726,使燃燒效率提升5.9%;OH的摩爾分?jǐn)?shù)從0到0.01224,使燃燒效率提升5.6%;H的摩爾分?jǐn)?shù)從0到0.0004,使燃燒效率提升6.1%.于是可以定性地看出,這三種自由基對(duì)甲烷超聲速燃燒效率提升能力強(qiáng)弱應(yīng)該依次為H、O、OH.

從熱力學(xué)角度分析,由于進(jìn)口污染空氣中自由基O、OH、H含量非常少,數(shù)量級(jí)在10-3左右.因此自由基O、OH、H含量引起燃?xì)獾钠骄▔罕葻崛莺推骄肿恿康母淖儗?duì)超燃室做功性能的影響是非常有限的,自由基主要通過促進(jìn)化學(xué)反應(yīng)動(dòng)力學(xué)過程來提升超燃室的性能.

2.4 污染組分中間產(chǎn)物CO和H2的作用

污染組分中間產(chǎn)物CO和H2的工況設(shè)置如表4,結(jié)果分析如圖12,隨著進(jìn)口空氣中的中間產(chǎn)物CO和H2含量的增多,超燃室的燃燒效率、比推力均明顯上升了.

因?yàn)镃O本身是性質(zhì)比較活躍的燃料,通過OH+CO?H+CO2作為額外添加的燃料提升了甲烷超燃性能;H2對(duì)超聲速燃燒的作用與H、CO對(duì)超聲速燃燒的作用是類似的.H2的加入會(huì)抑制甲烷的消耗,但H2本身通過O+H2?H+OH和OH+H2?H+ H2O,作為添加燃料使超燃室的性能得到提升.

3 結(jié)論

根據(jù)乙醇燃燒加熱器燃燒化學(xué)反應(yīng)得到的主要七種組分,采用甲烷18組分24步反應(yīng)簡(jiǎn)化機(jī)理,得到污染空氣中最終組分H2O、CO2,自由基H、O、OH和中間產(chǎn)物CO、H2對(duì)甲烷超聲速燃燒性能的主要影響如下:

(1)進(jìn)口空氣中污染組分H2O含量增多使燃燒效率、比推力均下降.H2O的增加造成平均比熱容增加,總溫度降低,并作為第三體抑制甲烷的燃燒過程,使超燃室的性能下降.

(2)進(jìn)口空氣中污染組分CO2含量增多對(duì)燃燒效率沒有明顯影響,但使超燃室參考比推力下降.CO2主要因?yàn)槠浯蠓肿恿刻匦允谷細(xì)馄骄肿恿吭龃?從而降低超燃室做功能力,對(duì)甲烷超燃性能起消極作用.

(3)進(jìn)口空氣中自由基O、OH、H含量增多均使超燃室燃燒效率上升.自由基O、OH主要通過化學(xué)反應(yīng)動(dòng)力學(xué)過程,加快燃燒過程,促進(jìn)甲烷和煤油消耗而提高燃燒效率,從而提升超燃室性能.自由基H作為額外添加的活躍燃料對(duì)甲烷超燃性能起積極作用.

(4)進(jìn)口空氣中污染組分CO和H2均對(duì)甲烷超燃性能起積極作用.它們的作用原理與自由基H的作用類似.

1 Edelman,R.B.;Spadaccini,L.J.J.Spacecraft,1969,6(12):1442

2 Mattick,S.J.;Frankel,S.H.Numerical modeling of supersonic combustion_validation and vitiation studies using FLUENT// 41stAIAA/ASME/SAE/ASEE Joint Propulsion Conference& Exhibit.Tucson,Arizona,2005:10-13

3 Tomioka,S.;Hiraiwa,T.;Kobayashi,K.;Izumikawa,M.;Kishida, T.;Yamasaki H.J.Propulsion Power,2007,23(4):789

4 Li,W.Q.;Song,W.Y.JournalofAirForceEngineering University, 2006,5(7):10 [李衛(wèi)強(qiáng),宋文艷.空軍工程大學(xué)學(xué)報(bào),2006,5(7): 10]

5 Liu,W.X.;He,W.;Li,H.B.;Li,X.Y.;Le,J.L.Chinese Science Bulletin,2009,54(8):1317 [劉偉雄,賀 偉,李宏斌,李象遠(yuǎn),樂嘉陵.科學(xué)通報(bào),2008,53(8):2257]

6 Shao,J.X.;Tan,N.X.;Liu,W.X.;Li,X.Y.Acta Phys.-Chim. Sin.,2010,26(2):270 [邵菊香,談寧馨,劉偉雄,李象遠(yuǎn).物理化學(xué)學(xué)報(bào),2010,26(2):270]

7 Liu,O.Z.;Cai,Y.H.;Hu,Y.L.;Liu,J.H.;Ling,W.H.Journal of Propulsion Technology,2004,25(5):463 [劉歐子,蔡元虎,胡欲立,劉敬華,凌文輝.推進(jìn)技術(shù),2004,25(5):463]

8 Peters,N.;Kee,R.J.Combust.Flame,1987,68:17

9 Tong,G.;Huang,Y.;Chen,Y.L.Journal of Fuel Chemistry and Technology,2000,28(1):49 [董 剛,黃 鷹,陳義良.燃料化學(xué)學(xué)報(bào),2000,28(1):49]

10 Bowman,C.T.;Hanson,R.K.;Davidson,D.F.;Gardiner Jr.,W. C.;Lissianski,V.;Smith,G.P.;Golden,D.M.;Frenklach,M.; Goldenberg,M.1994,http://www.me.berkeley.edu/gri_mech/

11 Marinov,N.M.Int.J.Chem.Kinet.,1999,31(3):183

12 Davidenko,D.M.;G?kalp,I.;Dufour,E.;Magre,P.Systematic numericalstudyofthesupersoniccombustioninanexperimental combustion chamber.14thAIAA/AHI Space Planes and Hypersonic Systems and Technologies Conference,AIAA 2006-7913

Effects of Species in Vitiation Air on Methane-Fueled Supersonic Combustion

HOU Ling-Yun1,*YANG Jin1MA Xue-Song2LIU Wei2
(1School of Aerospace,Tsinghua University,Beijing 100084,P.R.China;2Beijing Power Machinery Research Institute,Beijing 100074,P.R.China)

Based on a detailed chemical reaction mechanism,a reduced reaction mechanism with 18 species and 24 steps was used to simulate the supersonic combustion of methane.Heated air calculations showed that seven main vitiated species,i.e.,H2O,CO2,O,OH,CO,H,and H2,were present in ethanolfueled heated air.We analyzed the effects of these species on methane-fueled supersonic combustion using chemical kinetics and thermodynamics.H2O inhibits the combustion process,decreases the combustion efficiency,and decreases the specific thrust.The relatively large molecular weight of CO2contributes to an increase in the mean molecular weight of the fuel gas,which is a negative factor in the mechanism of specific thrust.Free radicals O,OH,H can effectively promote the combustion process and thus increase the combustion efficiency.Intermediate products CO and H2increase the combustion efficiency,and this is a function of the additional fuel.

Supersonic combustion;Chemical reaction mechanism;Vitiation air;Species

O643

Received:May 10,2010;Revised:August 13,2010;Published on Web:October 26,2010.

*Corresponding author.Email:lyhou@tsinghua.edu.cn;Tel:+86-10-62772157.

The project was supported by the National Natural Science Foundation of China(50306011).

國家自然科學(xué)基金(50306011)資助項(xiàng)目

?Editorial office ofActa Physico?Chimica Sinica