甲烷/空氣燃燒簡化機(jī)理及其應(yīng)用

江　冰，季　強(qiáng)，趙昌普，武昭暉，胡建功

(1.太原學(xué)院 機(jī)電工程系，山西 太原030032；2.泛亞汽車技術(shù)中心有限公司 研發(fā)試驗(yàn)中心，上海 201201；3.天津大學(xué) 內(nèi)燃機(jī)燃燒學(xué)國家重點(diǎn)實(shí)驗(yàn)室，天津 300072)

?

甲烷/空氣燃燒簡化機(jī)理及其應(yīng)用

江冰1，季強(qiáng)2，趙昌普3，武昭暉1，胡建功1

(1.太原學(xué)院 機(jī)電工程系，山西 太原030032；2.泛亞汽車技術(shù)中心有限公司 研發(fā)試驗(yàn)中心，上海 201201；3.天津大學(xué) 內(nèi)燃機(jī)燃燒學(xué)國家重點(diǎn)實(shí)驗(yàn)室，天津 300072)

摘要：為了找到適合活塞式發(fā)動(dòng)機(jī)甲烷燃燒的化學(xué)反應(yīng)動(dòng)力學(xué)簡化機(jī)理，從甲烷燃燒詳細(xì)化學(xué)動(dòng)力學(xué)機(jī)理出發(fā)，利用敏感性分析法，分析了包括多個(gè)不同組分的22步、39步、51步和58步基元反應(yīng)的4種不同簡化機(jī)理，并根據(jù)活塞式發(fā)動(dòng)機(jī)的燃燒特點(diǎn)，將這4種簡化機(jī)理應(yīng)用于均勻攪拌反應(yīng)器模型內(nèi)甲烷/空氣預(yù)混燃燒過程的計(jì)算中。研究結(jié)果表明：與詳細(xì)反應(yīng)機(jī)理的計(jì)算結(jié)果相比，當(dāng)采用較少步數(shù)的22步基元反應(yīng)的簡化機(jī)理時(shí)，計(jì)算得到的燃燒溫度結(jié)果誤差較大，且出現(xiàn)明顯的峰值后移現(xiàn)象，計(jì)算誤差最大可達(dá)到12.5%。但隨著采用基元反應(yīng)簡化機(jī)理的步數(shù)增加，這些誤差明顯減小。當(dāng)采用58步基元反應(yīng)的簡化機(jī)理時(shí)，計(jì)算得到的出口溫度和主要組分物質(zhì)的量濃度隨當(dāng)量比變化的結(jié)果，與詳細(xì)反應(yīng)機(jī)理計(jì)算得到的結(jié)果基本吻合，最大誤差不超過1%。這說明58步的基元反應(yīng)簡化機(jī)理包括了主要的反應(yīng)式，計(jì)算結(jié)果具有足夠的精度，可以很好地預(yù)測甲烷/空氣預(yù)混燃燒現(xiàn)象。

關(guān)鍵詞：甲烷燃燒；簡化機(jī)理；敏感性分析；均勻攪拌反應(yīng)器

0引言

由于內(nèi)燃機(jī)燃用煤層氣時(shí)，可以有效降低各類廢氣的排放量，所以，煤層氣發(fā)動(dòng)機(jī)的研究與應(yīng)用已逐漸受到重視[1-5]。目前，已有學(xué)者致力于研究發(fā)動(dòng)機(jī)燃用甲烷[6]和煤層氣摻氫燃燒[7-10]時(shí)的燃燒特性。然而，要想控制好發(fā)動(dòng)機(jī)中煤層氣的燃燒過程，還應(yīng)首先對(duì)甲烷燃燒的化學(xué)動(dòng)力學(xué)機(jī)理進(jìn)行深入細(xì)致的研究。

由于在模擬計(jì)算中，使用包含53種組分、325個(gè)基元反應(yīng)的甲烷燃燒詳細(xì)化學(xué)動(dòng)力學(xué)機(jī)理GRI(gas research institute)-Mech3.0時(shí)，會(huì)存在“剛性”和計(jì)算效率的問題，在計(jì)算復(fù)雜流動(dòng)燃燒時(shí)，無法滿足對(duì)計(jì)算速度和存貯量的要求，給流動(dòng)/燃燒耦合模擬計(jì)算帶來困難。因此，必須對(duì)復(fù)雜的化學(xué)反應(yīng)動(dòng)力學(xué)模型進(jìn)行簡化，使其能夠應(yīng)用于湍流火焰的計(jì)算中。

目前，在甲烷化學(xué)反應(yīng)動(dòng)力學(xué)簡化模型研究中，文獻(xiàn)[11]簡化得到包括18種組分、14個(gè)基元反應(yīng)的含氮的甲烷/空氣燃燒14步簡化機(jī)理，并在寬熱力學(xué)參數(shù)變化范圍內(nèi)，很好地預(yù)測了燃燒現(xiàn)象和NO的生成量。而文獻(xiàn)[12]簡化得到包含14種組分、18個(gè)基元反應(yīng)的簡化機(jī)理，并能夠在較大熱力學(xué)參數(shù)變化范圍內(nèi)，較好地預(yù)測甲烷/空氣預(yù)混燃燒現(xiàn)象，獲得了較高精度的點(diǎn)火延遲時(shí)間估計(jì)。顯然，在簡化過程中，對(duì)于究竟應(yīng)包含多少種組分和多少個(gè)基元反應(yīng)才比較合適的問題，目前仍難以確定。為了對(duì)這個(gè)問題有比較清楚的認(rèn)識(shí)，本文從甲烷燃燒化學(xué)動(dòng)力學(xué)機(jī)理著手，對(duì)4種甲烷燃燒化學(xué)反應(yīng)動(dòng)力學(xué)的簡化機(jī)理進(jìn)行了分析和對(duì)比，找出了較適合于煤層氣發(fā)動(dòng)機(jī)燃燒特性的簡化機(jī)理。

1機(jī)理的簡化

1.1計(jì)算模型與簡化方法

目前，常用的簡化化學(xué)反應(yīng)動(dòng)力學(xué)模型的方法有準(zhǔn)穩(wěn)態(tài)假設(shè)(quasi-steady-state approximation,QSSA)法和敏感性分析法兩種。利用QSSA方法進(jìn)行簡化時(shí)，雖然不需要求解該物質(zhì)的速率微分方程，使計(jì)算量減少，但是由于平衡反應(yīng)和穩(wěn)態(tài)物種需要經(jīng)驗(yàn)來判斷，因此，該方法只適用于較小的組分和溫度范圍，且人為誤差較大。而敏感性分析通過刪去對(duì)反應(yīng)過程影響不大的基元反應(yīng)來減小化學(xué)反應(yīng)的規(guī)模，可以深入直觀地分析組分、基元反應(yīng)和反應(yīng)條件等因素對(duì)系統(tǒng)反應(yīng)參數(shù)變化的敏感程度，因而這種方法近年來得到了較廣泛的應(yīng)用。

本文采用敏感性分析法，建立起一套甲烷點(diǎn)火燃燒的化學(xué)反應(yīng)動(dòng)力學(xué)模型簡化方法，得到相應(yīng)的簡化化學(xué)反應(yīng)動(dòng)力學(xué)模型，并在均勻攪拌反應(yīng)器條件下進(jìn)行計(jì)算和檢驗(yàn)。

受限于篇幅，本文不再詳細(xì)列出模型的控制方程。

1.2機(jī)理的簡化過程

首先，將詳細(xì)機(jī)理GRI-Mech 3.0運(yùn)用到一維預(yù)混層流火焰模型中，計(jì)算典型工況下甲烷的燃燒現(xiàn)象，得出預(yù)混層流火焰反應(yīng)式；然后，通過敏感性分析，刪除計(jì)算結(jié)果中對(duì)燃燒過程影響不明顯的物質(zhì)和基元反應(yīng)，從而得到包含不同基元反應(yīng)的簡化機(jī)理；最后，再將得到的簡化機(jī)理應(yīng)用于均勻攪拌反應(yīng)器模型中，對(duì)甲烷/空氣預(yù)混燃燒過程進(jìn)行計(jì)算和檢驗(yàn)。

為了在簡化過程中，找到較為合適的基元反應(yīng)簡化機(jī)理，本文分別研究了22步、39步、51步和58步不同基元反應(yīng)的簡化機(jī)理。表1為22步基元反應(yīng)的化學(xué)式(限于篇幅的原因，39步、51步和58步基元反應(yīng)的化學(xué)式不在本文中列出。)

表1　22步基元反應(yīng)簡化機(jī)理

組分：CH4、CH3、CH2、CH、CH2O、HCO、CO、CO2、HCN、CH3O、C2H6、H2、N2、N、NO、NO2、O、O2、OH、H、HO2、H2O、H2O2

2計(jì)算分析

2.1敏感性分析

一般情況下，敏感性因數(shù)在化學(xué)反應(yīng)劇烈的著火延遲時(shí)間附近有較大的值，而在反應(yīng)完成之后或者滯燃期時(shí)，敏感性因數(shù)趨近于零。溫度敏感性因數(shù)為正，表示增大該基元反應(yīng)的指前因子等參數(shù)時(shí)，溫度增加；溫度敏感性因數(shù)為負(fù)，表示增大該基元反應(yīng)的指前因子等參數(shù)時(shí)，溫度下降。

對(duì)甲烷燃燒時(shí)溫度敏感性分析結(jié)果表明：溫度敏感性因數(shù)較大的一些基元反應(yīng)為R32、R38、R53、R98、R116、R119、R155、R156、R158和R161，表2為所對(duì)應(yīng)的基元反應(yīng)的化學(xué)反應(yīng)方程。它們對(duì)整個(gè)化學(xué)反應(yīng)的速率起著重要作用。其中，基元反應(yīng)R32、R38、R119、R155、R156和R161對(duì)溫度的敏感性因數(shù)為正，當(dāng)這幾個(gè)反應(yīng)的反應(yīng)速率增大時(shí)，燃燒的溫度也會(huì)相應(yīng)上升，即它們主要是放熱反應(yīng)?；磻?yīng)R53、R98、R116和R158對(duì)溫度的敏感性因數(shù)為負(fù)，當(dāng)這幾個(gè)反應(yīng)的反應(yīng)速率增大時(shí)，燃燒的溫度會(huì)相應(yīng)下降，即它們主要是吸熱反應(yīng)。

表2　溫度敏感性因數(shù)較大的基元反應(yīng)

表3為在過量空氣因數(shù)為1.0時(shí)，燃燒對(duì)CH4、CO2、H和OH濃度影響的敏感性分析結(jié)果。從表3可以看出：基元反應(yīng)對(duì)CH4、CO2、H和OH濃度影響的敏感性因數(shù)變化都較大，其變化分別為：-26.2～30.0、-25.7～22.4、-110.1～96.0和-98.6～86.9。相對(duì)而言，基元反應(yīng)R32、R119和R158對(duì)CH4、CO2、H和OH濃度的影響更大。這說明，在燃燒中基元反應(yīng)R32、R119和R158起著重要的作用。

表3　燃燒對(duì)CH4、CO2、H和OH濃度影響的敏感性分析結(jié)果

2.2簡化機(jī)理計(jì)算分析

由于內(nèi)燃機(jī)的燃燒過程屬于復(fù)雜的流動(dòng)燃燒現(xiàn)象，因此，選擇均勻攪拌反應(yīng)器模型作為簡化機(jī)理的應(yīng)用來進(jìn)行計(jì)算分析。

計(jì)算中均勻攪拌反應(yīng)器容積V=67.4 cm3，條件為絕熱(即反應(yīng)器散熱速率Q= 0)，反應(yīng)物進(jìn)口溫度T*= 473 K，滯留時(shí)間τ為0.002 s。

圖1　　　　在0.025 MPa條件下，均勻攪拌反應(yīng)器　　　出口溫度隨當(dāng)量比的變化

圖1為在0.025 MPa壓力條件下，將采用22步、39步、51步和58步不同基元反應(yīng)簡化機(jī)理應(yīng)用于均勻攪拌反應(yīng)器時(shí)，出口溫度隨當(dāng)量比φ(實(shí)際的燃料空氣比與當(dāng)量條件下燃料空氣比的比值)變化的計(jì)算分析結(jié)果。表4為在0.025 MPa壓力條件下各簡化模型的計(jì)算誤差。

從圖1可以看出：采用不同的基元反應(yīng)簡化機(jī)理時(shí)，均勻攪拌反應(yīng)器出口溫度的計(jì)算結(jié)果是有差別的。當(dāng)采用22步基元反應(yīng)簡化機(jī)理時(shí)，無論在低壓還是在高壓條件下，均勻攪拌反應(yīng)器出口溫度的計(jì)算結(jié)果都與詳細(xì)基元反應(yīng)機(jī)理(GRI-Mech 3.0)的計(jì)算結(jié)果有明顯的差別，誤差較大。在當(dāng)量比φ較小情況下，22步基元反應(yīng)簡化機(jī)理的計(jì)算結(jié)果明顯小于詳細(xì)基元反應(yīng)機(jī)理的計(jì)算結(jié)果，如在0.025 MPa壓力條件下，計(jì)算誤差最大達(dá)7.30%(見表4)；而在當(dāng)量比φ較大情況下，22步基元反應(yīng)簡化機(jī)理的計(jì)算結(jié)果明顯大于詳細(xì)基元反應(yīng)機(jī)理的計(jì)算結(jié)果，且出現(xiàn)明顯的峰值后移現(xiàn)象，計(jì)算誤差在當(dāng)量比φ為1.6時(shí)已達(dá)-12.50%。但隨著所采用基元反應(yīng)簡化機(jī)理的步數(shù)增加時(shí)，這些誤差明顯減小。當(dāng)采用51步和58步基元反應(yīng)簡化機(jī)理時(shí)，均勻攪拌反應(yīng)器出口溫度的計(jì)算結(jié)果與詳細(xì)基元反應(yīng)機(jī)理的計(jì)算結(jié)果的誤差均明顯減小，如在0.025 MPa壓力條件下，在當(dāng)量比φ為0.6～1.6的范圍內(nèi)，采用51步基元反應(yīng)簡化機(jī)理計(jì)算時(shí)的最大計(jì)算誤差為1.31%(見表4)，而采用58步基元反應(yīng)簡化機(jī)理計(jì)算時(shí)的最大計(jì)算誤差已降到0.77%。特別是隨著壓力的增大，51步和58步基元反應(yīng)簡化機(jī)理的計(jì)算結(jié)果已十分接近詳細(xì)基元反應(yīng)機(jī)理的計(jì)算結(jié)果(具體的計(jì)算誤差本文不再一一列出)，只是在當(dāng)量比φ大于1.6后，51步和58步基元反應(yīng)簡化機(jī)理的計(jì)算結(jié)果比詳細(xì)基元反應(yīng)機(jī)理的計(jì)算結(jié)果偏大一些，其主要原因是混合氣濃度過稀，即對(duì)稀薄的混合氣燃燒，需要更多步數(shù)的簡化機(jī)理來處理才合適。

表4　在0.025 MPa條件下各簡化模型的計(jì)算誤差　%

圖2　　　　不同壓力下，均勻攪拌反應(yīng)器反應(yīng)物和主　　　產(chǎn)物物質(zhì)的量濃度隨當(dāng)量比的變化

圖2為將采用22步、39步、51步和58步的不同基元反應(yīng)簡化機(jī)理應(yīng)用于均勻攪拌反應(yīng)器時(shí)，反應(yīng)物(O2、CH4)和主產(chǎn)物(CO2、H2O)物質(zhì)的量濃度隨當(dāng)量比φ變化的計(jì)算結(jié)果，壓力分別為0.025 MPa、0.100 MPa和2.000 MPa。

從圖2可以看出：在不同的壓力下，采用不同的基元反應(yīng)簡化機(jī)理時(shí)，反應(yīng)物(CH4、O2)和主產(chǎn)物(CO2、H2O)的物質(zhì)的量濃度隨當(dāng)量比φ的變化情況也是有差別的。當(dāng)采用22步基元反應(yīng)簡化機(jī)理時(shí)，無論在低壓還是在高壓條件下，反應(yīng)物(CH4、O2)和主產(chǎn)物(CO2、H2O)的物質(zhì)的量濃度都與詳細(xì)基元反應(yīng)機(jī)理的計(jì)算結(jié)果有明顯的差別，且誤差較大，如0.025 MPa壓力條件下，在當(dāng)量比φ為0.6～1.6時(shí)，反應(yīng)物(CH4、O2)和主產(chǎn)物(CO2、H2O)的最大計(jì)算誤差分別為-1 300%(φ為1.6)(此處出現(xiàn)這么大的相對(duì)誤差是由于詳細(xì)基元反應(yīng)機(jī)理的計(jì)算結(jié)果接近于0)、-39.1%(φ為1.0)、-25.7%(φ為1.6)和25.6%(φ為0.6)。產(chǎn)生這些較大誤差的原因是，采用22步基元反應(yīng)簡化機(jī)理計(jì)算時(shí)，由于機(jī)理過于簡化，有一些重要的基元反應(yīng)沒有被考慮，所以造成了計(jì)算較大誤差的出現(xiàn)。然而，隨著采用基元反應(yīng)簡化機(jī)理步數(shù)的增加，這些誤差明顯減小。當(dāng)采用51步和58步基元反應(yīng)簡化機(jī)理時(shí)，特別是隨著壓力的增大，其計(jì)算結(jié)果已十分接近詳細(xì)基元反應(yīng)機(jī)理的計(jì)算結(jié)果，其最大誤差均不超過1%。這說明此時(shí)對(duì)燃燒起主要作用的重要基元反應(yīng)已基本上得到了考慮。

以上的計(jì)算結(jié)果表明：采用22步基元反應(yīng)簡化機(jī)理時(shí)，其計(jì)算結(jié)果與詳細(xì)基元反應(yīng)機(jī)理的計(jì)算結(jié)果相比，誤差較大，只有采用較高步數(shù)的基元反應(yīng)簡化機(jī)理時(shí)，如采用58步基元反應(yīng)簡化機(jī)理，其計(jì)算結(jié)果才會(huì)接近詳細(xì)基元反應(yīng)機(jī)理的計(jì)算結(jié)果。

3結(jié)論

(1)對(duì)甲烷/空氣燃燒采用不同步數(shù)的基元反應(yīng)簡化機(jī)理時(shí)，出口溫度、反應(yīng)物和主產(chǎn)物物質(zhì)的量濃度的計(jì)算結(jié)果與采用詳細(xì)反應(yīng)機(jī)理計(jì)算所得結(jié)果相比，存在明顯的差別。當(dāng)采用較少步數(shù)基元反應(yīng)的簡化機(jī)理，計(jì)算誤差較大，但采用較多步數(shù)的基元反應(yīng)簡化機(jī)理時(shí)，計(jì)算誤差逐步減小。

(2)在均勻攪拌反應(yīng)器的燃燒條件下，應(yīng)用22步甲烷/空氣基元反應(yīng)簡化機(jī)理時(shí)，計(jì)算得到的出口溫度和主要組分物質(zhì)的量濃度隨當(dāng)量比的結(jié)果與詳細(xì)機(jī)理的計(jì)算結(jié)果相比，存在著較大的差別。在當(dāng)量比較小的情況下，22步基元反應(yīng)簡化機(jī)理的計(jì)算結(jié)果明顯小于詳細(xì)基元反應(yīng)機(jī)理的計(jì)算結(jié)果，如在0.025 MPa壓力條件下，計(jì)算誤差最大達(dá)7.3%；而在當(dāng)量比較大情況下，計(jì)算結(jié)果明顯大于詳細(xì)基元反應(yīng)機(jī)理的計(jì)算結(jié)果，且出現(xiàn)明顯的峰值后移現(xiàn)象，如在0.025 MPa壓力條件下，計(jì)算誤差在當(dāng)量比φ為1.6時(shí)已達(dá)12.5%。

(3)在均勻攪拌反應(yīng)器的燃燒條件下，應(yīng)用58步甲烷/空氣基元反應(yīng)簡化機(jī)理時(shí)，計(jì)算得到的出口溫度和主要組分物質(zhì)的量濃度隨當(dāng)量比變化的結(jié)果與詳細(xì)機(jī)理的計(jì)算結(jié)果吻合較好，最大誤差不超過1%，說明在均勻攪拌反應(yīng)器的燃燒條件下，采用58步甲烷/空氣基元反應(yīng)簡化機(jī)理來模擬是合適的。

(4)對(duì)屬于復(fù)雜流動(dòng)燃燒的發(fā)動(dòng)機(jī)煤層氣燃燒問題，采用如58步的多步數(shù)甲烷/空氣基元反應(yīng)簡化機(jī)理來進(jìn)行模擬計(jì)算較為合適。

參考文獻(xiàn)：

[1]李從心,張欣,劉建華.發(fā)動(dòng)機(jī)燃用煤層氣燃料時(shí)燃燒特性的試驗(yàn)研究[J].汽車工程,2010,32(1):951-955.

[2]李從心,張欣,胡準(zhǔn)慶,等.發(fā)動(dòng)機(jī)燃用煤層氣燃料燃燒和排放性能試驗(yàn)研究[J].內(nèi)燃機(jī)工程,2009,30(5):31-35.

[3]左承基,郭威,錢葉劍,等.火花點(diǎn)火式變組分煤層氣發(fā)動(dòng)機(jī)的工作穩(wěn)定性和排放特性[J].內(nèi)燃機(jī)學(xué)報(bào),2003,21(5):329-332.

[4]左承基,錢葉劍,顧紅蘭.煤層氣發(fā)動(dòng)機(jī)HC和NO排放模型及試驗(yàn)[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2005,36(9):1-4.

[5]顧紅蘭.火花點(diǎn)火式發(fā)動(dòng)機(jī)燃用變組分煤層氣的燃燒模型研究[J].內(nèi)燃機(jī)工程,2004,25(5):72-75.

[6]TAIBLSKANDAR M.Combustion characteristics of methane in a direct injection engine using spark plug fuel injector[J].Jurnal kejuruteraan,2010,22:43-52.

[7]SILVANA D I,PAOLO S,BIANCA M V.Experimental investigation on the combustion process in a spark ignition optically accessible engine fueled with methane/hydrogen blends[J].International journal of hydrogen energy,2014,39(18):9809-9823.

[8]MANFRED K,HELMUT E,MARKUS S.Mixtures of hydrogen and methane in the internal combustion engine-synergies,potential and regulations[J].International journal of hydrogen energy,2012,37(15):11531-11540.

[9]MCTAGGART C G P,ROGAK S,MUNSHI S R,et al.Combustion in a heavy-duty direct-injection engine using hydrogen-methane blend fuels[J].International journal of engine research,2009,10(1):1-13.

[10]MANFRED K,HELMUT E,MARKUS S.Variable mixtures of hydrogen and methane in the internal combustion engine of a prototype vehicle:regulations,safety and potential[J].International journal of vehicle design,2010,54(2):137-155.

[11]董剛,任祝寅,陳義良.CH4/空氣含氮燃燒的14步簡化機(jī)理[J].燃料化學(xué)學(xué)報(bào),2003,31(1):39-43.

[12]劉合,陳方,劉洪,等.甲烷/空氣預(yù)混超聲速燃燒的18步簡化機(jī)理[J].燃燒科學(xué)與技術(shù),2012,18(5):467-472.

中圖分類號(hào)：TK432

文獻(xiàn)標(biāo)志碼：A

收稿日期：2015-10-09

作者簡介：江冰(1963-),女,河南正陽人,副教授,碩士,主要從事代用燃料發(fā)動(dòng)機(jī)工作過程控制方面的研究.

基金項(xiàng)目：山西省科技攻關(guān)基金項(xiàng)目(20130322002)

文章編號(hào)：1672-6871(2016)03-0024-05

DOI:10.15926/j.cnki.issn1672-6871.2016.03.006