(武漢理工大學(xué) 能源與動(dòng)力工程學(xué)院，武漢 430063)

天然氣的燃點(diǎn)較高，在壓燃式發(fā)動(dòng)機(jī)中不易被壓燃。通?？紤]向氣缸中噴入微量的柴油以引燃天然氣的方式來實(shí)現(xiàn)天然氣在船用壓燃式機(jī)的應(yīng)用[1]。噴入氣缸的柴油蒸發(fā)后會(huì)與缸內(nèi)的天然氣及空氣混合，其預(yù)混物的燃燒是天然氣-柴油雙燃料發(fā)動(dòng)機(jī)燃燒過程中的重要燃燒階段之一[2]，對發(fā)動(dòng)機(jī)的運(yùn)行參數(shù)有一定的影響[3]。

甲烷作為天然氣的主要成分，是常見的天然氣燃燒特性研究的替代物,通常選用高碳烷烴作為其單組分替代物[4]。正十二烷(n-C12H26)的碳值、碳?xì)浔?、放熱率、熱解性等化學(xué)性質(zhì)與柴油相近，且其蒸發(fā)與噴霧特性也與柴油相似[5-6]。因此,選擇甲烷-正十二烷作為天然氣-柴油研究的替代物燃料,通過試驗(yàn)及模擬手段研究正十二烷的含量對甲烷-正十二烷-空氣混合物預(yù)混層流燃燒特性的影響。

1 試驗(yàn)方法

1.1 試驗(yàn)裝置及數(shù)據(jù)處理

試驗(yàn)裝置主要由圓柱形定容燃燒彈、溫度控制系統(tǒng)、配氣供液系統(tǒng)、點(diǎn)火控制系統(tǒng)、高速攝像紋影系統(tǒng)和數(shù)據(jù)采集處理系統(tǒng)等組成。容彈的容積為12.8 L，側(cè)面開有直徑為120 mm的石英玻璃視窗孔，供高速攝像機(jī)拍攝記錄火焰?zhèn)鞑サ募y影圖像。試驗(yàn)相機(jī)拍攝幀率為10 000幀/s。

試驗(yàn)初始溫度選擇為423 K以保證足夠的過熱度，使正十二烷在試驗(yàn)條件下均能完全蒸發(fā)可作為理想氣體。初始?jí)毫?.1 MPa，當(dāng)量比φ從0.7～1.4。研究中空氣為79%的氮?dú)馀c21%的氧氣的混合物，正十二烷含量指正十二烷在甲烷-正十二烷混合燃料中的摩爾分?jǐn)?shù)，其值在0%～100%的范圍內(nèi)變動(dòng)。

試驗(yàn)開始前，先通過加熱系統(tǒng)將容彈內(nèi)溫度加熱并保持在設(shè)定溫度[7]。容彈被加熱到設(shè)定溫度后容彈內(nèi)徑向距離上的溫度分布見圖1，其中橫坐標(biāo)中0點(diǎn)為容彈中心。

圖1 定容燃燒彈內(nèi)沿重力軸向不同位置的溫度分布

由圖1可知，彈內(nèi)溫度沿著重力相反方向略有上升，在本試驗(yàn)裝置的觀察視窗范圍(±60 mm)內(nèi)溫差在±2 K以內(nèi)。

在將各組分氣液按比例充入容彈后，為保證彈內(nèi)氣體均勻混合，需靜置等待至少10 min。之后通過點(diǎn)火裝置點(diǎn)燃容彈內(nèi)的可燃混合物，同時(shí)觸發(fā)高速相機(jī)記錄下火焰的傳播過程。為保證試驗(yàn)結(jié)果的準(zhǔn)確性和可重復(fù)性，每組試驗(yàn)進(jìn)行至少3次，將試驗(yàn)結(jié)果的平均值作為最終結(jié)果[8]。

利用軟件對高速相機(jī)拍攝的球形傳播火焰過程紋影圖進(jìn)行處理，獲得火焰半徑R隨時(shí)間t的變化規(guī)律，拉伸火焰?zhèn)鞑ニ俣萐n可由求導(dǎo)Sn=dR/dt獲得。將球形火焰表面上一個(gè)無限小面積的數(shù)值對時(shí)間的求導(dǎo)可得到球形傳播火焰的拉升率α[21]，即α=2Sn/R。通過對Sn和α的關(guān)系進(jìn)行擬合后可推算出α為0時(shí)的無拉伸火焰?zhèn)鞑ニ俣萐l，非線性關(guān)系式可以更好地與試驗(yàn)數(shù)據(jù)點(diǎn)擬合[9]。

(Sn/Sl)2ln(Sn/Sl)2=-2Lbα/Sl

(1)

式中:Lb為馬克斯坦長度，用于評(píng)價(jià)球形火焰?zhèn)鞑サ牟环€(wěn)定性。

0.1 MPa、423 K條件下正十二烷-空氣混合物的拉伸火焰?zhèn)鞑ニ俣入S拉伸率變化的非線性擬合結(jié)果見圖2。

圖2 不同當(dāng)量比下正十二烷-空氣混合物拉伸火焰?zhèn)鞑ニ俣扰c拉升率的關(guān)系

可見擬合線與試驗(yàn)結(jié)果貼合較好，因此，后續(xù)工況的數(shù)據(jù)處理均采用非線性擬合。層流燃燒速度如下。

ul=Slρb/ρu

(2)

式中：ρu為未燃?xì)怏w的密度，由理想氣體的初始條件確定；ρb是已燃?xì)怏w的密度，其比值可以通過化學(xué)平衡法計(jì)算求得。

1.2 模擬方法及機(jī)理選擇

通過CHEMKIN PRO[10]軟件模擬甲烷-正十二烷-空氣的一維平面層流預(yù)混火焰，可以獲得其層流燃燒速度和火焰結(jié)構(gòu)等參數(shù)。模擬設(shè)置數(shù)值解梯度為0.03，曲率為0.03，網(wǎng)格點(diǎn)數(shù)在600以上，確保模擬計(jì)算有足夠的精度。

采用You等[11]提出的正十二烷的詳細(xì)化學(xué)動(dòng)力學(xué)機(jī)理(簡稱You機(jī)理)，該機(jī)理是在USC-Mech II[12]的基礎(chǔ)上添加了正十二烷等高碳烷烴的熱解機(jī)理所得，包含了171個(gè)組分和1 306步反應(yīng)。USC-Mech II是用于計(jì)算C1—C4低碳烷烴及合成氣燃燒的反應(yīng)機(jī)理，包含了甲烷燃燒詳細(xì)的基礎(chǔ)反應(yīng)，所選機(jī)理能較好地預(yù)測甲烷的層流燃燒速度。

1.3 試驗(yàn)裝置及機(jī)理驗(yàn)證

將試驗(yàn)所測得的甲烷及正十二烷的層流燃燒速度與現(xiàn)有文獻(xiàn)的試驗(yàn)結(jié)果及You機(jī)理的模擬結(jié)果進(jìn)行對比，結(jié)果見圖3?？梢?，試驗(yàn)裝置所獲得的數(shù)據(jù)與文獻(xiàn)中試驗(yàn)數(shù)據(jù)重合性較好。此外，You機(jī)理的預(yù)測結(jié)果與本文試驗(yàn)結(jié)果也具有較好的一致性。

圖3 正十二烷-空氣混合物及甲烷-空氣混合物層流燃燒速度的試驗(yàn)及模擬結(jié)果的對比

2 結(jié)果與討論

2.1 試驗(yàn)結(jié)果

2.1.1 層流燃燒速度

甲烷-正十二烷-空氣火焰的層流燃燒速度在正十二烷含量較小時(shí)變化比較明顯[13]，因此設(shè)置正十二烷含量步長在0%～30%的范圍內(nèi)取5%，在正十二烷含量較大時(shí)選取較大步長。在423 K、0.1 MPa，當(dāng)量比φ為0.8、1.0、1.3條件下甲烷-正十二烷-空氣混合物的層流燃燒速度隨正十二烷含量變化的試驗(yàn)及模擬值見圖4。

圖4 不同當(dāng)量比下甲烷-正十二烷-空氣混合物的層流燃燒速度隨正十二烷含量的變化

由圖4可知，在當(dāng)量比φ=1.0時(shí)，純甲烷與純正十二烷的層流燃燒速度相近，隨著正十二烷含量的增加，甲烷-正十二烷-空氣的層流燃燒速度提升甚微。而在當(dāng)量比為0.8和1.3的條件下，甲烷-正十二烷-空氣的層流燃燒速度隨著正十二烷含量的增加表現(xiàn)出了明顯的非線性增長，其增長幅度隨著正十二烷含量的增加而減少；當(dāng)正十二烷含量達(dá)到25%后，混合物層流燃燒速度的增長趨于平緩。此外，還可以看出，You機(jī)理模擬結(jié)果所顯現(xiàn)的混合物層流燃燒速度變化趨勢與試驗(yàn)值基本保持一致，表明其可以較好地預(yù)測甲烷-正十二烷混合物的層流燃燒速度。

2.1.2 馬克斯坦長度

用馬克斯坦長度評(píng)價(jià)火焰的不穩(wěn)定性。在當(dāng)量比φ為0.8、1.0和1.3的情況下，隨著正十二烷含量的增加，甲烷-正十二烷-空氣火焰的馬克斯坦長度的變化見圖5。

圖5 不同當(dāng)量比下甲烷-正十二烷-空氣混合物火焰的馬克斯坦長度隨正十二烷含量的變化

在φ=1.0時(shí)，馬克斯坦長度隨正十二烷含量的增加而有略微的增加；在φ=0.8時(shí)，隨著正十二烷含量的增加，馬克斯坦長度在初期表現(xiàn)了出了較大幅度的增長，表明其火焰的不穩(wěn)定性減弱，拉伸對火焰的傳播的抑制作用增強(qiáng)；當(dāng)含量超過約25%后，其增加幅度在慢慢變小并逐漸趨于穩(wěn)定；在φ=1.3時(shí)，甲烷-正十二烷-空氣火焰的馬克斯坦長度隨正十二烷含量的變化趨勢與φ=0.8時(shí)的趨勢相反。在φ為0.8和1.3條件下火焰的馬克斯坦長度均在正十二烷含量較小的條件下有大幅度的變化，在正十二烷含量達(dá)到約為25%后趨于穩(wěn)定。

2.2 化學(xué)動(dòng)力學(xué)分析

2.2.1 層流燃燒速度的敏感性分析

選取φ=0.8，正十二烷含量分別為0%、25%、50%及100%的情況進(jìn)行層流燃燒速度的敏感性分析。計(jì)算獲得正敏感性系數(shù)最大的10個(gè)基元反應(yīng)及負(fù)敏感性系數(shù)[14]最小的10個(gè)基元反應(yīng)，結(jié)果見圖6。

圖6 甲烷-正十二烷-空氣混合物層流燃燒速度的敏感性分析

由圖6可見，正十二烷含量從0%增加到25%后主要基元反應(yīng)對層流燃燒速度的敏感性系數(shù)有了很大的變化，而正十二烷含量提升到50%后變化則不是很明顯。

2.2.2 自由基及基元反應(yīng)速率分析

423 K、0.1 MPa下φ=0.8的甲烷-正十二烷-空氣混合物燃燒過程中五種自由基濃度的峰值與正十二烷含量的關(guān)系見圖7。

圖7 φ=0.8時(shí)甲烷-正十二烷-空氣火焰中主要自由基濃度峰值與正十二烷含量關(guān)系

隨著正十二烷含量的提高，HO2的濃度峰值沒有較大的變化，CH3有略微的下降，OH有略微的上升，而O與H的濃度峰值先有顯著的上升，而在正十二烷含量達(dá)到25%后趨于穩(wěn)定。

自由基會(huì)影響基元反應(yīng)進(jìn)行，進(jìn)而影響火焰的層流燃燒速度。從敏感性分析的結(jié)果得到主要影響混合物層流燃燒速度的基元反應(yīng)，速率峰值與正十二烷含量的關(guān)系見圖8。

圖8 φ=0.8時(shí)甲烷-正十二烷-空氣火焰中重要反應(yīng)的反應(yīng)速率峰值與正十二烷含量的關(guān)系

可見具有正敏感性系數(shù)的基元反應(yīng)(R1、R31、R40、R92)的反應(yīng)速率峰值明顯遠(yuǎn)高于具有負(fù)敏感性系數(shù)的反應(yīng)(R88、R12、R35)，這保證了自由基濃度的持續(xù)增加。

對甲烷-正十二烷-空氣混合物燃燒過程影響較大的自由基H、O的濃度以及基元反應(yīng)R1、R31的反應(yīng)速率峰值都隨著正十二烷含量的提升而迅速上升，并在達(dá)到約25%以后開始趨于平緩，這也是導(dǎo)致甲烷-正十二烷-空氣混合物層流燃燒速度隨正十二烷含量而發(fā)生變化的內(nèi)在原因。

2.2.3 甲烷反應(yīng)路徑分析

使用CHEMKIN PRO軟件計(jì)算所得的正十二烷含量分別為0%、25%及50%時(shí)甲烷的反應(yīng)路徑，見圖9。

圖9 不同正十二烷含量下甲烷的反應(yīng)路徑

可以看出，正十二烷含量為25%時(shí)甲烷的反應(yīng)路徑與正十二烷含量為0%時(shí)有很大的差別，而與正十二烷含量為50%時(shí)甲烷的反應(yīng)路徑相差較小。隨著正十二烷含量的增加，甲烷形成碳數(shù)更高烷烴的反應(yīng)占甲烷消耗路徑中的比重在增加，這與上述層流燃燒速度敏感性分析所獲得的與C2—C3相關(guān)的基元反應(yīng)對層流燃燒速度的影響在加入正十二烷后逐漸增強(qiáng)的趨勢相符；此外，隨著正十二烷的添加，在CH4、CH3、C2H6等多數(shù)物質(zhì)的脫氫反應(yīng)中，消耗H的反應(yīng)的比重在減少，消耗OH的反應(yīng)的比重在增加，這促使了R1(H+O2=O+OH)往正反應(yīng)方向的進(jìn)行，從而增加了活性自由基的生成，提高了混合物的層流燃燒速度。

3 結(jié)論

1)隨著正十二烷含量的增加，甲烷-正十二烷-空氣混合物的層流燃燒速度和火焰的馬克斯坦長度快速變化并向正十二烷-空氣的層流燃燒特性接近，當(dāng)正十二烷含量增加到約25%后變化趨于平緩。

2)化學(xué)動(dòng)力學(xué)分析的結(jié)果顯示，混合物燃燒中主要的自由基H、O的濃度，以及對層流燃燒速度正敏感性系數(shù)最大的基元反應(yīng)R1(H+O2=O+OH)、R31(CO+OH=CO2+H)的反應(yīng)速率峰值隨著正十二烷含量的增加表現(xiàn)出了較大的升高，其提升量在正十二烷含量達(dá)到約25%后不再明顯。這也是混合物層流燃燒速度發(fā)生變化的內(nèi)在原因。

3)對燃燒過程中甲烷的反應(yīng)路徑分析結(jié)果表示，正十二烷的加入對甲烷自身的燃燒反應(yīng)也有一定的促進(jìn)作用。