合成氣燃燒反應機理的驗證和分析

第2卷第3期新能源進(jìn)展VoL 2 No. 32014年6月ADVANCES IN NEW AND RENEWABLE ENERGYJun.2014文章編號:2095-560X(2014)03-0173-07合成氣燃燒反應機理的驗證和分析王全德1,魏賞賞2,王偉2,崔晨曉2(1.中國礦業(yè)大學(xué)低碳能源研究院,徐州221008;2.中國礦業(yè)大學(xué)電力工程學(xué)院,徐州221116)摘要:本文針對典型合成氣燃燒的詳細化學(xué)動(dòng)力學(xué)機理進(jìn)行了系統的驗證和分析。通過(guò)不同反應機理對點(diǎn)火延遲時(shí)間和層流火焰速度的預測,研究和分析了不同反應機理的區別和模擬結果的可靠性。釆用強制敏感度分析方法揭示了影響合成氣點(diǎn)火延遲時(shí)間和層流火焰速度的重要反應,并對相關(guān)反應的動(dòng)力學(xué)進(jìn)行了分析討論,為進(jìn)一步構建統一可靠的燃燒反應機理奠定了基礎。關(guān)鍵詞:合成氣;燃燒反應機理;敏感度分析;化學(xué)動(dòng)力學(xué)中圖分類(lèi)號:TK16文獻標志碼:Adoi:10.3969/1sn.2095560X.2014.03.002Validation and Analysis of Reaction Mechanisms for Syngas CombustionWANG Quan-de, WEI Shang-shang, WANG Wei, CUI Chen-xiao(1. Low Carbon Energy Institute, China University of Mining and Technology, Xuzhou 221008, China;2. School of Electric Power Engineering, China University of Mining and Technology, Xuzhou 221116, ChinaAbstract: Detailed reaction mechanisms for syngas combustion are validated and analyzed in the present work. Accordingto the predictions of ignition delay time and laminar flame speed by employing contemporary combustion mechanisms, thedifferences among these mechanisms are presented and their robustness for numerical simulations are validated. Further,on the basis of sensitivity analysis, important reactions which greatly affect the ignition and laminar flame speed areidentified, and related reaction rate constants are discussed. The present work provides fundamental information for furtherdevelopment of a universal reaction mechanism for syngas combustionKey words: syngas; combustion reaction mechanism; sensitivity analysis; chemical kinetics0前言子燃料燃燒化學(xué)動(dòng)力學(xué)模型而構建的。這些機理針對特定的燃燒性質(zhì)進(jìn)行了模擬驗證。然而,不同的反應隨著(zhù)世界對能源需求的持續增加和環(huán)境問(wèn)題的機理包括的基元反應和化學(xué)動(dòng)力學(xué)參數等存在很大日益突出,以煤基合成氣(主要成分為H2和CO)的差異,這也直接導致了不同反應機理的模擬精度存為主要燃料的整體煤氣化聯(lián)合循環(huán)系統( Integrated在差異。因此,有必要對目前的合成氣燃燒反應機理Gasification Combined Cycle,lGCC)引起了廣泛關(guān)進(jìn)行系統的驗證和分析研究,構建統一可靠的燃燒反注和重視。燃氣輪機是IGCC動(dòng)力系統的關(guān)鍵設備,應機理,這不但對合成氣燃燒的研究具有重要作用,合成氣燃燒特性是燃氣輪機先進(jìn)燃燒室設計的基同時(shí)對構建大分子燃料燃燒反應機理也具有重要意礎。構建可靠的合成氣燃燒化學(xué)動(dòng)力學(xué)機理,是燃義。本文旨在通過(guò)對已構建的典型合成氣燃燒反應機氣輪機燃燒室計算流體力學(xué)設計的前提,同時(shí)合成理進(jìn)行模擬驗證和分析,通過(guò)動(dòng)力學(xué)模擬驗證不同機氣燃燒化學(xué)動(dòng)力學(xué)機理也是構建大分子燃料燃燒化理模擬結果的精確性,采用強制敏感度分析方法揭示學(xué)動(dòng)力學(xué)模型的基礎山。目前,已有很多適用于合成影響反應機理模擬結果中的重要反應,對相關(guān)反應動(dòng)氣燃燒的化學(xué)動(dòng)力學(xué)模型被構建起來(lái),這其中有些是力學(xué)參數進(jìn)行分析討論,為進(jìn)一步構建統一可靠的專(zhuān)門(mén)用于合成氣燃燒的反應機理,有些則是為大分反應機理提供中國煤化工*收稿日期:201403-19修訂日期:201404-14CNMHG基金項目:國家中央高?；究蒲袠I(yè)務(wù)費專(zhuān)項資助基金(20130NA08)通信作者;王全德,E-mail:wqgd98686@126com174新能源進(jìn)展第2卷合成氣燃燒化學(xué)動(dòng)力學(xué)機理和模擬化的H3CO燃燒反應機理、由 Williams等構建的分析方法San diego機理、Li等人構建的C1機理和 Ranzi構建的機理,在本文中分別用 Davis機理、San1.1反應機理Diego機理、Li機理和 Ranzi機理表示。同時(shí)本文也本文選取4個(gè)典型的合成氣燃燒反應機理進(jìn)行對最近由Sun和 Keromnes等構建的兩個(gè)機理進(jìn)了系統的比較研究,4個(gè)反應機理分別為 Davis等優(yōu)行了簡(jiǎn)要的分析比較。表1列出了目前合成氣燃燒表1合成氣燃燒反應機理包含的基元反應Table 1 Reaction list in reaction mechanisms for syngas combustion反應DanzH+0,=0+OH1234OH +H,=H+H,OOH +OH=0+H,oH+h+M=H,+MO+O+M=O+M6789o+h+M=oh +MH+oh+M=H,O+O+oh+M=Ho,+MH+O2(+M)=HO2(+M)√√√√√√√√√√√√√H,+O=HO2+H√√√√√x√√√x√HO,+h=o+HOHO,+H=OH +OH15HO,+OH=O,+H,O16 OH+ OH (+M)=H,O,(+M√√√√√√√√x√√√x√√√√√√√√√√√√√√x√√√√√√√√√√√HO,+HO,=O+H,O,√√x√√√√√√√√√√√HO,+H=HO,+HH,O,+H=oh+H,oH,O,+OH= HO,+H,OCO+O(+M)=CO(+ M)√√√√√√Co+0,= CO,+0CO+HO,=CO,+OHHCO+M=co+h+M2HCO+O,=Co+HOHCO +H= Co+H,√√√√x×x√√√√√√√√√√√√√√√x√v√√√√√v√√√√√√√√√√√vv√√√√√√√HCO+O=CO,+HHCO +OH=CO +H,O√√√√√√√HCO+HO2=CO,+oH+HHCO+HCo=H,+2C0中國煤化工VCNMHGHCO + HO,=H,O,+CO注:√表示對應的反應機理包含這一反應,x表示機理不包含此反應第3期王全德等:合成氣燃燒反應機理的驗證和分析175反應機理中包含的所有反應,以及六個(gè)機理包含的延遲時(shí)間。圖1~圖3給出了不同條件下合成氣點(diǎn)對應的反應。目前,合成氣燃燒反應機理包含的物火延遲時(shí)間的實(shí)驗結果和4個(gè)反應機理的模擬結種數均為14個(gè),總的基元反應數目為35步,其中果。由圖1~圖3可以看出,4個(gè)機理對于點(diǎn)火延27步基元反應存在于所有的機理中。由表1可見(jiàn)反應機理對于氫氣部分的差異較小,主要是由于HO2+H=O+H2O、O+OH+M=HO2+M、O+Hm co: H=s0+M=OH+M三個(gè)基元反應引起的。機理包含的反空10Davis mechanism應的差別主要集中在HCO相關(guān)的反應由Smn構建自Ramzi mechanism的機理,與 Davis優(yōu)化的機理基本一致,模擬結果San diego mechanism差別較小; Keromnes機理是在Li機理的基礎上構建起來(lái)的,二者包含的反應完全一致,因此,本文重點(diǎn)針對前4個(gè)反應機理進(jìn)行驗證和分析1.2模擬方法點(diǎn)火延遲時(shí)間和層流火焰速度是燃料燃燒性質(zhì)0.700.750800.850.900.951.001.051.101000T(1/K)的重要參數。本文通過(guò)對這兩個(gè)燃燒性質(zhì)的模擬比圖1合成氣點(diǎn)火延遲時(shí)間模擬和實(shí)驗結果較不同機理的模擬精度。動(dòng)力學(xué)模擬采用 Chemkin Fig1 eling and experimental ignition delay times for2.0軟件圖,點(diǎn)火延遲時(shí)間采用 SENKIN模塊的均syngas c相反應器模型,層流火焰速度采用一維層流火焰速H/CO/Air mixtures, CO: H-20: 80度模塊模擬。層流火焰速度模擬過(guò)程中,輸運參p=llatm, (=0.5數均采用混合平均方法處理。本文采用強制敏感Davis mechanism度分析方法闡明機理中影響燃料點(diǎn)火和層流火焰速度的關(guān)鍵反應。以點(diǎn)火延遲時(shí)間為例,敏感度分析10San Diego mechanism的系數定義為121100%畫(huà)上式中,k;為第i個(gè)反應的速率常數(指前因子),0.901.05n(2k)為第i個(gè)反應的速率常數加倍后的點(diǎn)火延遲時(shí)間,x(k)為反應速率常數不變時(shí)的點(diǎn)火延遲時(shí)圖2合成氣點(diǎn)火延遲時(shí)間模擬和實(shí)驗結果21間。由此可見(jiàn),一個(gè)正的敏感度系數表明增加反應Fg2 Modeling and experimen\ignition delay times fori的速率常數會(huì )降低活性,延長(cháng)燃料的點(diǎn)火延遲時(shí)間,反之亦然。對層流火焰速度的分析采用同樣方HyCO/O / Ar mixture, H, CO=70: 30. Ar=94%法,需要注意的是,此時(shí)一個(gè)正的敏感度系數表明Davis mechanis增加反應i的速率常數會(huì )提高反應活性,提高火焰的傳播速度。103San Diego mechanis2結果與討論不同的文獻對于點(diǎn)火延遲時(shí)間的定義有多種形式,如依據溫度相對初始溫度的升高和OH自由基濃度的變化趨勢等。通過(guò)對不同定義形式的比較,0.80發(fā)現不同的定義得到的點(diǎn)火延遲時(shí)間基本沒(méi)有差v中國煤化工”1圖3合別。本文采用了 SENKIN模塊的默認點(diǎn)火延遲時(shí)間CNMHG結果2Fig 3 Modeling and experimental ignition delay times for的定義,即初始溫度升高400K時(shí)的時(shí)間即為點(diǎn)火syngas combustion 22)176新能源進(jìn)展第2卷遲時(shí)間的預測存在較大差別。在溫度高于1000K和時(shí)間主要對如下反應動(dòng)力學(xué)參數較為敏感低壓力條件下,4個(gè)機理的模擬結果基本一致;在H+O2=0+OH壓力為20atm的高壓條件下,4個(gè)機理都不能準確地模擬點(diǎn)火延遲時(shí)間的變化趨勢,其中, San DiegoH+O2(+M)=HO2(+M)機理在溫度高于1100K時(shí)的模擬結果最接近實(shí)驗OH +OH (+M)=H,O2(+M)值,而 Davis和 Ranzi機理在溫度低于1100K時(shí)的HO,+HO,=O+ho模擬結果較好為了揭示影響點(diǎn)火延遲時(shí)間的重要反應,我們H,O2+H=HO,+ H2采用第二部分給出的敏感度分析方法,采用4個(gè)機CO +HO=CO2+Oh理分別進(jìn)行了敏感度分析。圖4給出了 Davis機理的敏感度分析結果。值得注意的是,敏感度分析系Co+02= CO,+O數小的反應,不代表該反應不重要,它主要反映點(diǎn)此外,一個(gè)明顯的趨勢是隨著(zhù)溫度降低和壓力火延遲時(shí)間或層流火焰速度對于該反應速率常數的升高,OH+OH(+M)=H2O2(+M)、HO2+HO2=O2+敏感度較小。對于敏感度系數較大的反應,其速率H2O2和H2O2+H=HO2+H2這三個(gè)反應對點(diǎn)火延常數的變化對點(diǎn)火延遲時(shí)間的影響較大,因此需要遲時(shí)間的敏感度系數逐漸增大,表明這三個(gè)反應是對其動(dòng)力學(xué)參數特別關(guān)注。由圖4可見(jiàn),點(diǎn)火延遲影響低溫高壓條件下點(diǎn)火延遲時(shí)間的關(guān)鍵反應。CO+M=CO+H+MDay is mech anismCO+HO-CO+oH150KCO+0=CO+O1050KCO+OHCO +H950KCO+O(+MCO (+MHO+OH=HO +HOH O+H=ho+h1O2+HO2 =0 +H,oHO +oH=o +HoHO +0=OH+0HO+H=OH+OHOH+OH(+MH,O,(+M)H, +O=HO,+HH+O+M)HO、+M)O+H=H+OHH+0=0+OHSensitivity coefficients圖4點(diǎn)火延遲時(shí)間的敏感度分析結果(H2CO空氣混合物,當量比0.5,初始壓力20am,初始混合物摩爾分數為733%H29.71%CO、1.98%CO2、17.01%O2和63.97%N2)Fig. 4 Sensitivity analysis of ignition delay time to rate constants(H2/CO/air mixture at equivalence ratio of 0.5, P= 20 atm. Theinitial mole fractions of reactants are 7.33% H,, 9.71% co, 1.98% Co,, 17.01% O2, and 63.97%o n,)近年來(lái),隨著(zhù)人們對低溫高壓等先進(jìn)燃燒技術(shù)和O2的反應和H2與HO2的氫提取反應,其速率常和理論研究的重視,這三個(gè)反應的動(dòng)力學(xué)參數值得數存在較大差別。圖5和圖6給出了兩個(gè)反應典型詳細討論。Troe對OH+OH(+M=H2O2(+M)的速率常數隨的動(dòng)力學(xué)參數進(jìn)行了系統的分析,同時(shí)根據本文對明HO2自結小“中苗煤化工分析結果表低溫、高壓不同反應機理采用的反應速率比較表明,這一反應點(diǎn)火條件下抑」m明地這一反應速速率常數的差別較小。對于HO2自結合生成H2O2率常數隨溫度的變化趨勢,需要采用兩個(gè) Arrhenius第3期王全德等:合成氣燃燒反應機理的驗證和分析177形式的速率常數的加和,在反應機理中表現為重復HO +HO=HO+0反應。 Hippler等采用激波管實(shí)驗測量結果,提ippler et alSan Diego mechanism出了描述這一反應的速率常數,并且在Li和Sun機理中得到了應用。這一速率常數在 Davis和--.-Zhou et alKeromnes機理中應用時(shí),被減小了13%以獲得更精確的模擬結果。此后,Kapl等S擴大了實(shí)驗測量的溫度范圍,重新擬合了速率常數的表達式。最近,Zhou等采用高精度的量子化學(xué)計算結合化學(xué)動(dòng)力學(xué)理論,對這一反應的速率常數進(jìn)行了理論T(1/K)研究,并且提出了一個(gè)新的速率常數表達式。由圖圖5反應HO2+HO2=H2O2+O2的速率常數隨溫度的變化5可見(jiàn), San diego和 Ranzi機理的速率常數不能正Fig. 5 The rate constants of HO,+HO,= H,O,+0, as afunction of temperature確描述該反應速率常數隨溫度的變化趨勢。其他三個(gè)速率常數在高溫下的差別較大,而在低溫時(shí)差H,O,+ H=HO,+H,別開(kāi)始減小。通過(guò)與 Kappel等實(shí)驗測量結果比較,- GRI 3.0 mechanism我們發(fā)現 Hippler等提出的速率常數略微高估了實(shí)Sun mechanism驗結果,而理論計算的結果在700K左右也高估了IE12實(shí)驗結果。 Davis等在優(yōu)化反應機理時(shí),表明 Hippler等的速率常數減小13%后,能取得更好的模擬結果,此時(shí)的速率常數也與 Kappel等接近。H2與HO的氫提取反應是影響低溫高壓燃燒模擬結果的關(guān)0.40.60.8鍵反應之一。然而,目前這一反應的速率常數的實(shí)驗結果存在較大的差異和不確定性。圖6給出了目圖6反應HO2+H=HO2+H2的速率常數隨溫度的變化前反應機理采用的速率常數隨溫度的變化趨勢。LiFig. 6 The rate constants of H,O2+H= HO?+H as afunction of temperatureSan diego和Rani機理采用了Tang和 Hampson171在1986年提出的速率常數,Sun則采用了 Baulch圖7和圖8給出了4個(gè)機理對層流火焰速度的等最近建議的速率常數。由圖6可見(jiàn),這兩個(gè)速模擬結果和實(shí)驗結果。由圖7和8可以看出, Davis率常數差異較大。由于實(shí)驗數據較少,并且實(shí)驗結機理的模擬結果與實(shí)驗結果吻合得最好,這主要是果存在較大差異,Eing9n等19用一系列的量子由于在Dis機理優(yōu)化過(guò)程中,層流火焰速度和點(diǎn)化學(xué)計算方法結合變分過(guò)渡態(tài)理論對該反應的速火延遲時(shí)間是主要的目標函數,對應的實(shí)驗結果是率常數進(jìn)行了系統的理論研究,計算結果與 Tsang優(yōu)化的訓練集之一。和 Hampson等提出的速率常數接近。 Keromnes最H/COAir mixtures近構建的機理采用了 Ellingson等計算的速率常數。H1:CO=5050盡管如此,這一反應的速率常數仍然存在很大的不▲HCO=5:95確定性。為了提高反應機理在低溫條件點(diǎn)火模擬的精確性, Cavaliere等0人為地對該反應的速率常數進(jìn)行調整,并且對低溫點(diǎn)火延遲時(shí)間的模擬取得了更好的結果。但是,人為地對反應速率進(jìn)行調整,有可能破壞反應機理的化學(xué)合理性。因此,有必要采用實(shí)驗和理論等對這一反應進(jìn)行更深入的研究。此0.51.01.52.02.53.03.54.04.55.0外,與CO/HCO相關(guān)的反應的動(dòng)力學(xué)參數的不確定較小,對點(diǎn)火延遲時(shí)間的影響很小。由此可見(jiàn),合圖7CO和H2成氣燃燒的性質(zhì)主要受H12燃燒化學(xué)動(dòng)力學(xué)過(guò)程的結果(COHyH中國煤化工的模擬和實(shí)驗CNMHG壓力為2am)影響。Fig 7 Calculated laminar flame speeds by using variousmechanism and experimental measurements178新能源進(jìn)展第2卷H, /CO/O, /He mixture.H2CO=50:50顯而易見(jiàn),H+O+M=OH+M這一反應對于層流火焰速度有較大影響。然而,這一反應并沒(méi)有包含在 Ranzi機理中。由此可見(jiàn),對目前構建的機理進(jìn)Ranzi mechanism行系統的比較和分析,獲得反應機理的基元反應的完整列表,對構建完整統一的燃燒反應機理的必要性。此外,層流火焰速度對于CO和OH生成CO2和H的反應較為敏感,這主要是由于層流火焰速度主要受高溫燃燒化學(xué)動(dòng)力學(xué)過(guò)程控制,與絕熱火焰溫度密切相關(guān),而這一反應是燃燒過(guò)程放熱的主要圖8Co和H2在高壓條件下層流火焰速度的模擬和實(shí)驗結來(lái)源之一。與點(diǎn)火延遲時(shí)間的敏感度分析結果一致,果6(CO〃H/O/He混合物初始溫度298K,壓力10am)層流火焰速度對于HCO相關(guān)反應的敏感度系數較Fig 8 Calculated laminar flame speeds by using variousmechanism and experimental measurements小。同時(shí)由于層流火焰速度是燃料燃燒的高溫特性因此對低溫高壓條件下的點(diǎn)火延遲時(shí)間有重要影響圖9給出了對層流火焰速度采用 Davis機理的的反應,如OH+OH(+M)=H2O2(+M)、HO2+HO2敏感度分析結果。與點(diǎn)火延遲時(shí)間相比,對層流火=O2+H2O2和HO2+H=HO2+H2這三個(gè)反應對焰速度較為敏感的反應有很大的不同。層流火焰速層流火焰速度的模擬結果幾乎沒(méi)有影響。對于上述度主要受與H自由基相關(guān)的反應控制影響層流火焰速度的重要反應,我們通過(guò)對不同反H+o=0+OH應機理采用的動(dòng)力學(xué)參數的比較表明,這些反應的O+h2=h+ohoH+h=H+H,O速率常數基本一致。需要注意的是H與OH生成H2OH+oh+m=h,o+m的反應,有些機理(如Li機理)采用的是H2O分解o+h+M=oh+M的形式,由于動(dòng)力學(xué)參數不確定性較大,因此高估HO,+H=oh+oh了層流火焰速度。因此,對這一反應的速率常數需HO2+Oh=O2+H,o要進(jìn)一步研究和討論。HCO+O=CO+HOHCO+H, O=CO+H+H, OHCO+M=Co+H+MCO+OH=CO +HCO+OH=CO +HHO +oH=0+h OHO +OH=O,+H, OHO +0=0H+0COH=50:50HO +H=oH+OHH, +O, =HO,+HH+O2(+M)=HO2(+MO+H+M=OH+M=05OH+h=H+h oH+0=0+OH0050.000.050.10中國煤化工Sensitivity coef ficientCNMHG圖9層流火焰速度的敏感度分析結果( Davis機理)Fig 9 Sensitivity analysis results of laminar flame speed to rate constants by using the Davis mechanism第3期王全德等:合成氣燃燒反應機理的驗證和分析3結論9 Lutz A E, KeeR J, Miller J A Senkin: A Fortran Programfor Predicting Homogeneous Gas Phase Chemical Kinetic本文對典型合成氣燃燒的詳細化學(xué)動(dòng)力學(xué)機理with Sensitivity Analysis, SAND87-8248[R]. Livermore,CA: Sandia National laboratories. 1990.進(jìn)行了系統的分析和驗證研究,優(yōu)化的 Davis機理101 Keer, Grcar jf, Smooke md,etal. A FORTRAN對層流火焰速度的模擬結果與實(shí)驗結果吻合得最Program for Modeling Steady Laminar One-dimensionalPremixed Flames, SAND85-8240R]. San Diego, CA:好,4個(gè)機理都不能完整地模擬點(diǎn)火延遲時(shí)間在高Sandia National Laboratories. 1985壓下的點(diǎn)火延遲時(shí)間。敏感度分析結果表明支鏈反eR. Dixon-Lewis G, Warnatz J. A Fortran computer應H+O2=HO+O對點(diǎn)火延遲時(shí)間和層流火焰速multicomponent transport properties, SAND86-8246IRI度均有很大影響。隨著(zhù)溫度降低,OH+OH(+M)=Livermore. CA: Sandia National Laboratories. 1988[121 Wang Q D, Fang Y M, Wang F, et al. Systematic analysisH2O2(+M)、HO2+HO2=O2+H2O2和H2O2+H=and reduction of combustion mechanisms for ignitionHO2+H2這三個(gè)反應對點(diǎn)火延遲時(shí)間的敏感度系數lulti-component kerosene surrogate[J]. Proceedings ofthe Combustion Institute, 2013, 34(1): 187-195開(kāi)始增大。HO2+HO2=O2+H2O2和H2O2+H=[13] Troe j. The thermal dissociation/recombination reactionHO2+H2反應的速率常數存在很大的不確定性,值of hydrogen peroxide H2O2(+ M)-20H ( M)Ill得實(shí)驗和理論的進(jìn)一步研究。同時(shí),通過(guò)對不同機理Analysis and representation of the temperature andpressure dependence over wide ranges[J]. Combustion包含的重要反應的動(dòng)力學(xué)參數的比較表明影響反應and flame,2011,158(4):594601機理模擬精度的其他基元反應的動(dòng)力學(xué)參數的不確[41 Hippler H,ToeJ. Willner j. Shock wave study of thereaction HO,+ HO,-H,O,+O, Confirmation of a rate定性相對較小。此外,合成氣燃燒特性和反應機理的constant minimum near 700K[]. The Journal of模擬精度主要受氫氣燃燒化學(xué)動(dòng)力學(xué)過(guò)程的控制Chemical Physics, 1990, 93: 1755-1760[15 Kappel C, Luther K, Troe J. Shock wave study of theunimolecular dissociation of H,O2 in its falloff range and參考文獻Chemical Physics, 2002, 4(18): 4392-4398[1 Marcos C, Dryer F L. Syngas combustion kinetics and16 Zhou D Y, Han k, Zhang P, et al. Theoreticalapplications[J]. Combustion Science and Technology,Determination of the rate Coefficient for the ho,+ho2008,180(6):1053-1096-+H,O,+O, Reaction Adiabatic Treatment of2] Davis S G Joshi A V, Wang H, et al. An optimized kineticAnharmonic Torsional Effects [J]. The Journal of Physicalmodel of H2/CO combustion!J]. Proceedings of theChemistry A,2012,116(9):20892100Priyank S, Williams F A. Testing a small detailed [171 Tsang W, Hampson R F Chemical kinetic data base forCombustion Institute, 2005, 30(1): 1283-1292combustion chemistry. Part I. Methane and relatedhemical-kinetic mechanism for the combustion ofcompounds[J]. Journal of Physical and Chemicalhydrogen and carbon monoxide[J]. Combustion andReference Data. 1986. 15: 1087-1279Flame,2006,145(1):316-323[18 Baulch D L, Bowman C T, Cobos C J, et al. Evaluated[4] Li J, Zhao Z, Kazakov A, et al. A comprehensive kinetickinetic data for combustion modeling: supplement I[]mechanism for CO, CH,O, and CH3OH combustion[JIJournal of Physical and Chemical Reference Data, 200International Journal of Chemical Kinetics, 2007, 39(3)34:757-1397[19 Ellingson B A, Theis d P, Tishchenko O, et al[5 Frassoldati A, Faravelli T, Ranzi E. The ignitionReactions of hydrogen atom with hydrogen peroxide[J]combustion and flame structure of carbon monoxideThe Journal of Physical Chemistry A, 2007, 111(51)hydrogen mixtures. Note 1: Detailed kinetic modeling of13554-13566syngas combustion also in presence of nitrogen [20] Cavaliere D E, Joannon M D, Sabia P, et al.Acompounds[J]. International Journal of Hydrogencomprehensive kinetic modeling of ignition of syngas-airEnergy,2007,32(15):3471-3485mixtures at low temperatures and high pressures[J][6] Sun H, Yang S I, Jomaas G et al. High-pressure laminarCombustion Science and Technology, 2010, 182(4-6)flame speeds and kinetic modeling of carbon692-701monoxide/hydrogen combustion[J]. Proceedings of the [21] Kalitan D M, Mertens J D, Crofton M W, et alCombustion Institute, 2007, 31(1): 439-446and oxidation of lean CO/H, fuel blends in air[J] Journal[7 Keromnes A, Metcalfe W K, Heufer K A, et al. Anof propulsion and power, 2007, 23(6): 1291-1301experimental and detailed chemical kinetic modeling study [22] Thi L D, Zhang Y J, Fu J, et al. Study on ignition delayof hydrogen and syngas mixture oxidation at elevatedof multi-component syngas using shock tube[J]. Thepressures[J]. Combustion and Flame, 2013, 160(6): 995-1011Canadian Journal of Chemical Engineering, 2014, 92(5)[8] KeeR J, Rupley F M, Miller J A. Chemkin-ll: A Fortran861-870metics Package for the Analysis ofGas-Phase Chemical Kinetics, SAND89-8009IRIAlbequerque, NM: Sandia National Laboratories, 1989作者簡(jiǎn)介:山中國煤化工Lutz a e, Kee R J, Miller J A Senkin: A Fortran Programfor Predicting Homogeneous Gas Phase Chemical Kinetics王全德(1986-),CNMHG:要從事燃燒化ith Sensitivity Analysis,SAND87-8248[R]. Livermore,學(xué)、計算化學(xué)和分子模擬方向的研究CA: Sandia National laboratories, 1990