汽油多組分表征燃料簡化動力學(xué)模型的系統(tǒng)化構(gòu)建及驗證

肖　干　張煜盛,*　姜光軍,(華中科技大學(xué)能源與動力工程學(xué)院，武漢430074；南昌大學(xué)機(jī)電工程學(xué)院，南昌33003)

汽油多組分表征燃料簡化動力學(xué)模型的系統(tǒng)化構(gòu)建及驗證

肖干1張煜盛1,*姜光軍1,2
(1華中科技大學(xué)能源與動力工程學(xué)院，武漢430074；2南昌大學(xué)機(jī)電工程學(xué)院，南昌330031)

提出了一套系統(tǒng)化多級機(jī)理簡化策略，包含基于誤差傳播的直接關(guān)系圖法、峰值濃度分析法、線性同分異構(gòu)體集總法、主組分分析法、溫度敏感性分析和產(chǎn)率分析法，并將其應(yīng)用于汽油四組分表征燃料詳細(xì)反應(yīng)機(jī)理的簡化，構(gòu)建了適用于HCCI發(fā)動機(jī)燃燒邊界下的簡化機(jī)理模型，包含149個物種、414個反應(yīng)。通過與激波管、快速壓縮機(jī)、增壓HCCI發(fā)動機(jī)實驗數(shù)據(jù)的對比驗證表明，新機(jī)理可以準(zhǔn)確地預(yù)測較寬范圍條件下的著火滯燃期，在HCCI發(fā)動機(jī)的單區(qū)模型計算中，該機(jī)理對缸內(nèi)燃燒和排放的預(yù)測結(jié)果是令人滿意的。放熱率分析表明，R+O2反應(yīng)是控制中間溫度區(qū)放熱的關(guān)鍵基元反應(yīng)，在高壓低溫下，異辛烷的放熱起到?jīng)Q定性作用。添加2-戊烯之后，使得四組分模型相較于三組分模型更為準(zhǔn)確，尤其是對于第一階段著火滯燃期有顯著影響，為進(jìn)一步探索調(diào)和燃料組分比例控制HCCI燃燒提供了一條新思路。

汽油表征燃料；多級機(jī)理簡化；HCCI發(fā)動機(jī)；中間溫度放熱；化學(xué)動力學(xué)模型

[Article]

www.whxb.pku.edu.cn

1　引言

面對日益嚴(yán)格的排放法規(guī)及有限的石油資源存儲量，下一代“低溫、預(yù)混合、壓燃”發(fā)動機(jī)成為內(nèi)燃機(jī)界的研究熱點(diǎn)1。與傳統(tǒng)燃燒模式相比，均質(zhì)充量壓燃(HCCI)、預(yù)混合充量壓燃(PCCI)、火花輔助壓燃(SACI)、分層充量壓燃(SCCI)、活性可控壓燃(RCCI)以及低溫燃燒(LTC)等新型燃燒理念中，燃料自燃著火化學(xué)反應(yīng)動力學(xué)扮演著重要的角色2。因此，在內(nèi)燃機(jī)燃燒數(shù)值模擬中，旨在精確預(yù)測發(fā)動機(jī)著火、燃燒和排放特點(diǎn)的化學(xué)反應(yīng)機(jī)理的開發(fā)，必然處于重中之重的位置。這一問題的解決不僅有利于洞悉上述新型燃燒模式的燃燒細(xì)節(jié)，而且對控制實際燃料的HCCI燃燒過程及污染物排放，進(jìn)而豐富內(nèi)燃機(jī)的燃燒理論均具有重要意義3。

汽油是由成百上千個連續(xù)光譜變化的碳?xì)淙剂辖M成的復(fù)雜混合物，包括直連和支鏈烷烴、環(huán)烷烴、烯烴、芳香烴及含氧劑，同時，其組分變化也受到原油來源、煉制工藝、季節(jié)等因素的影響4,5。假使以汽油的實際組分分布來描述其燃燒過程，將會導(dǎo)致化學(xué)動力模型中的基元反應(yīng)數(shù)目的爆炸性增長，其數(shù)值求解將大大超出現(xiàn)有計算機(jī)的計算能力。因此，從工程應(yīng)用的角度出發(fā)，一般采用表征燃料來描述實際汽油的燃燒特性，即用有限的幾個代表性單體碳?xì)浠衔锏幕旌先剂蟻肀碚鲗嶋H燃料的性能6。異辛烷作為最簡單的汽油表征燃料，通常用于描述火焰?zhèn)鞑ニ俣?。但是，由于其高辛烷值特?研究法(RON)或馬達(dá)法(MON)辛烷值均為100)，RON=MON=100)，并不能代表實際汽油的自燃著火特性?；A(chǔ)參比燃料(PRF)，即正庚烷和異辛烷的混合物，被用于調(diào)配實際汽油的辛烷值。一般來說，PRF的辛烷值敏感度(S=RON－MON)為0，而實際汽油由于包含其它類別的烴類組分(烯烴、芳香烴、含氧劑等)，辛烷值敏感度為10左右，在不同的發(fā)動機(jī)運(yùn)行工況下，需要使用不同辛烷值的PRF來描述實際汽油的自燃著火特性。另一方面，燃料的氫碳(H/C)比是決定新鮮充量空燃比和絕熱火焰溫度的主要指標(biāo)，PRF燃料的H/C比范圍局限于2.25－2.30之間，而實際汽油的氫碳比小于2。因此，PRF并不是理想的汽油表征燃料，在其基礎(chǔ)上添加敏感性燃料組分來更加真實地逼近實際汽油的辛烷值、辛烷值敏感度、H/C等燃燒化學(xué)特性指標(biāo)，成為一種迫切的需求。

當(dāng)前，國際上公認(rèn)的汽油表征燃料是正庚烷、異辛烷、甲苯(TRF)的混合物，大量的激波管及發(fā)動機(jī)實驗證明7－11，TRF可以準(zhǔn)確地描述實際汽油的自燃著火特性。隨著詳細(xì)化學(xué)反應(yīng)機(jī)理研究的不斷深入，一些研究者提出更為復(fù)雜的四組分乃至七組分的汽油表征燃料模型。Mehl等12依據(jù)RD387汽油的特性，提出了包含正庚烷、異辛烷、甲苯、2-戊烯的四組分表征模型；Perez和Boehman13進(jìn)一步添加了甲基環(huán)己烷，形成了五組分模型，并在定容燃燒彈上進(jìn)行了試驗研究。Naik等14使用1-戊烯代替上述2-戊烯，構(gòu)建了適用于汽油HCCI發(fā)動機(jī)的五組分表征模型。Fikri等15在TRF的基礎(chǔ)上，添加了二異丁烯，構(gòu)建了四組分表征模型。Andrae16提出了正庚烷、異辛烷、甲苯、二異丁烯和乙醇的五組分表征模型。Yahyaoui等17進(jìn)行了異辛烷、甲苯、1-己烯、ETBE混合燃料射流攪拌反應(yīng)器實驗研究。Sarathy18和Ahmed19等提出了適用于富烷烴汽油的五組分和六組分表征燃料模型。

盡管上述多組分表征燃料模型較之PRF可以更好地描述實際汽油的燃燒特性，但其詳細(xì)化學(xué)反應(yīng)動力學(xué)模型的規(guī)模十分龐大，往往包含數(shù)以千計的物種和反應(yīng)。在將其與發(fā)動機(jī)計算流體力學(xué)(CFD)模型相耦合進(jìn)行發(fā)動機(jī)性能的數(shù)值模擬時，會使計算工作量過于繁重。同時，由于各物種及其相關(guān)的反應(yīng)時間尺度相差很大，極易引起嚴(yán)重的計算“剛性”問題。因此，必須對詳細(xì)反應(yīng)機(jī)理進(jìn)行大規(guī)模簡化，以便在保證模擬計算精度的條件下，提高計算效率20。此外，使用簡化反應(yīng)機(jī)理可以方便地研究對燃燒過程起到關(guān)鍵作用的基元反應(yīng)和重要物種，有助于加深人們對燃料設(shè)計和污染物生成機(jī)理的認(rèn)識。

傳統(tǒng)的機(jī)理簡化方法嚴(yán)重依賴于人們對于詳細(xì)機(jī)理的認(rèn)識和“經(jīng)驗性”理解，由此獲得的簡化機(jī)理其適用范圍十分有限。其原因在于，首先，簡化機(jī)理的構(gòu)建大多是通過組裝不同來源的子機(jī)理而得到的，這將會造成動力學(xué)參數(shù)、熱力學(xué)參數(shù)和輸運(yùn)參數(shù)的不一致性。其次，為了達(dá)到簡化的極限，簡化機(jī)理往往使用高度集總的總包反應(yīng)步，通過經(jīng)驗性的擬合動力學(xué)參數(shù)來吻合實驗數(shù)據(jù)，缺失了詳細(xì)機(jī)理固有的層級拓展性特征。此外，簡化機(jī)理不包含對預(yù)測污染物很重要的基元反應(yīng)步驟，比如，多環(huán)芳香烴(PAH)被認(rèn)為是生成碳煙的重要前驅(qū)體，而C2－C4的中間物種對PAH生成路徑十分重要，這在簡化機(jī)理中無法得到充分體現(xiàn)。有鑒于此，一些理論化、自動化、可控化的機(jī)理簡化方法被陸續(xù)提出，其中，Lu和Law21于2005年提出的直接關(guān)系圖法(DRG)獲得高度的關(guān)注和認(rèn)可，進(jìn)而衍生出多種DRG方法的變形算法。比如，基于誤差傳播的直接關(guān)系圖法(DRGEP)22、路徑通量法(PFA)23、修正的DRG方法(DRG-max)24等等。但大多數(shù)機(jī)理簡化方法多應(yīng)用于單組分燃料的直接簡化，特別是正庚烷、異辛烷、正十二烷等單體碳?xì)淙剂系难趸磻?yīng)機(jī)理的簡化25,26，而對多組分表征燃料機(jī)理的簡化目前受到的關(guān)注較少。通常的做法是將每一個組分的詳細(xì)機(jī)理進(jìn)行簡化之后再組合起來，形成多組分混合燃料機(jī)理27,28。這種做法不僅忽略了組分間強(qiáng)烈的非線性耦合作用，而且也忽略了不同工況、不同組分比例條件下，多組分燃料間的反應(yīng)路徑會發(fā)生較大變化等實際情況。綜合考慮上述因素對多組分簡化機(jī)理構(gòu)建的影響勢必給相關(guān)研究工作提出了新的挑戰(zhàn)。

鑒于單一機(jī)理簡化方法的局限性，通過對多組分詳細(xì)機(jī)理固有特征的分析，本文從全方位多角度考慮，以基于誤差傳播的直接關(guān)系圖法DRGEP為起點(diǎn)，結(jié)合峰值物種濃度分析、線性同分異構(gòu)化集總、主組分分析、溫度敏感性分析和產(chǎn)率分析等方法，提出了一套新的系統(tǒng)化多級機(jī)理簡化策略。該策略摒棄了人為經(jīng)驗因素的干擾，保證每一步的簡化過程都在可控誤差精度的條件下進(jìn)行，盡量實現(xiàn)計算機(jī)自動簡化，以達(dá)到簡化極限。通過將其直接應(yīng)用于大規(guī)模汽油多組分機(jī)理的簡化過程，構(gòu)建出適用于HCCI發(fā)動機(jī)燃燒邊界條件，可用于CFD耦合計算的簡化機(jī)理模型，并藉此開展了激波管、快速壓縮機(jī)及增壓HCCI發(fā)動機(jī)的實驗驗證，闡明了簡化模型的有效性。最后，從化學(xué)反應(yīng)動力學(xué)角度，分析了控制中間溫度放熱的關(guān)鍵基元反應(yīng)，以及添加烯烴組分的優(yōu)勢所在。

2　構(gòu)建方法與策略

2.1汽油表征燃料詳細(xì)反應(yīng)機(jī)理的選擇

為了排除外界因素(原油來源、煉制工藝、季節(jié)等)對實際汽油組分不確定性的影響，本文采用美國RD387(research grade)汽油為研究對象，其辛烷值為87，H/C摩爾比為1.85。它代表一種通用的不含氧研究級汽油，近幾年被廣泛用于燃燒科學(xué)研究。Gauthier等8使用激波管測量了RD387汽油及其表征燃料在高壓高溫下的滯燃期；Kukkadapu等9進(jìn)一步完善了數(shù)據(jù)，利用快速壓縮機(jī)測量了中低溫和高壓下的滯燃期；Dec等10,11以RD387汽油為燃料，進(jìn)行了增壓HCCI發(fā)動機(jī)和部分燃料分層(PFS)發(fā)動機(jī)的實驗研究。這些實驗數(shù)據(jù)都為開發(fā)和驗證化學(xué)反應(yīng)動力學(xué)模型提供了寶貴的數(shù)據(jù)資料。根據(jù)Pitz和Dryer的建議4,5，表征燃料的組配應(yīng)該根據(jù)其具體的應(yīng)用目標(biāo)來進(jìn)行。以RD387汽油為目標(biāo)汽油，Mehl等12按照辛烷值等級、辛烷值敏感度、分子類別、H/C等配比指標(biāo)，提出了正庚烷、異辛烷、甲苯、2-戊烯的四組分機(jī)理模型，并確定其最佳摩爾比為0.153:0.488:0.306:0.053。其詳細(xì)機(jī)理29共包含1389物種、5935個基元反應(yīng)。在較寬廣的實驗條件下和各種基礎(chǔ)反應(yīng)器中，分別就單組分、兩組分及各種混合物進(jìn)行了系統(tǒng)的驗證，均表現(xiàn)出了較高的預(yù)測精度。因此，本文以該詳細(xì)機(jī)理為基礎(chǔ)，運(yùn)用多級化系統(tǒng)機(jī)理簡化策略，構(gòu)建適用于HCCI發(fā)動機(jī)燃燒條件的汽油多組分表征燃料的簡化化學(xué)動力學(xué)模型。

2.2機(jī)理簡化方法

在此提出的一套系統(tǒng)化多級機(jī)理簡化策略，包含了基于誤差的直接關(guān)系圖法、峰值濃度分析法、線性同分異構(gòu)體集總法、主組分分析法、溫度敏感性分析和產(chǎn)率分析法等、它們的具體應(yīng)用分述如下。

2.2.1基于誤差傳播的直接關(guān)系圖法

Pepiot-Desjardins和Pitsch22于2008年提出了帶有誤差傳播的直接關(guān)系圖法(DRGEP)，引入了搜索路徑及誤差傳播的概念，重新定義了物種之間直接作用系數(shù)。其基本思路是：通過目標(biāo)物種A進(jìn)行最優(yōu)路徑搜索可以得到所有其他組分對物種A的重要性系數(shù)，然后設(shè)定一個閾值，當(dāng)重要性系數(shù)大于閾值時，認(rèn)為此組分是重要的，予以保留。其物種之間直接相互作用系數(shù)為：

式中，ωi為第i個基元反應(yīng)的凈反應(yīng)速率，νA,i為物種A在第i個基元反應(yīng)中的化學(xué)計量系數(shù)，nR為總的基元反應(yīng)數(shù)目，如果第i個基元反應(yīng)包含組分B，則δ=1，否則為0。路徑依賴性系數(shù)為：

式中，s1=A，sn=B，RAB表示物種之間重要性系數(shù)的最終表達(dá)式。DRGEP方法的快速高效性，使其成為近些年最為流行的簡化方法之一，常常應(yīng)用于多級機(jī)理簡化方法的第一步，直接移除大規(guī)模冗余物種和反應(yīng)。

2.2.2峰值濃度分析法

在碳?xì)淙剂系淖匀贾疬^程中，不同的反應(yīng)階段起關(guān)鍵作用的基元反應(yīng)不同，各物種的時間尺度差異很大，進(jìn)而造成其濃度水平也不一樣。尤其是在多組分詳細(xì)反應(yīng)機(jī)理中，存在許多大分子低濃度物種，隨著反應(yīng)鏈的斷裂，它們往往對于著火點(diǎn)的預(yù)測影響很小，甚至可以忽略。據(jù)此，可采用類似于傳統(tǒng)的CSP理論與方法30，按照時間尺度的劃分，來選取準(zhǔn)穩(wěn)態(tài)物種。本文采用峰值濃度分析法，在寬廣的取值范圍內(nèi)，以物種的峰值濃度為依據(jù)，將每一個物種的峰值濃度按照升序大小排列，不斷地迭代移除低濃度物種，保留燃料、氧化劑、主要產(chǎn)物等高濃度物種。每一次移除后所形成的簡化機(jī)理，被用來計算著火滯燃期，并將其與詳細(xì)機(jī)理計算所得的滯燃期進(jìn)行比較，如果小于初始定義的誤差精度(一般為5%)，則該物種可以安全移除；如果超過誤差精度，則重新移回，針對每一個低濃度物種(小于1×10－7，摩爾分?jǐn)?shù))進(jìn)行測試，直到最終機(jī)理超出用戶定義的誤差精度。峰值濃度分析方法是一種系統(tǒng)化、自動化的簡化方法，可以直接衡量單一物種及其反應(yīng)對于著火點(diǎn)的影響，并且兼顧了簡化精度和簡化效率這一矛盾的兩個方面。其誤差精度判定式為：

2.2.3線性同分異構(gòu)體集總法

詳細(xì)機(jī)理中包含有大量的同分異構(gòu)體，這一情況在低溫反應(yīng)階段尤甚。由于同分異構(gòu)體物種具有相同的摩爾質(zhì)量、相似的熱力學(xué)特性和輸運(yùn)參數(shù)，可以考慮將其合并為一個代表性物種，將其所參與的基元反應(yīng)合并為一條主干反應(yīng)路徑，以進(jìn)一步減小機(jī)理的規(guī)模，同時也保留了詳細(xì)機(jī)理的固有層級性。本文主要采用Lu和Law31的線性集總處理方法，按照同分異構(gòu)體簇的統(tǒng)計平均質(zhì)量分布，修正了集總后反應(yīng)的動力學(xué)參數(shù)。

物種的比例系數(shù)為：

式中，χi為物種i的濃度，χL為物種i所在同分異構(gòu)體簇L的濃度。

集總后反應(yīng)的動力學(xué)參數(shù)為：

式中，kfj為反應(yīng)j的速率常數(shù)，v′ij為反應(yīng)j中物種i的正反應(yīng)化學(xué)計量系數(shù)。

2.2.4主組分分析法

主組分分析方法由Turanyi32提出，可以針對濃度敏感性矩陣進(jìn)行分析，也可以針對速率敏感性矩陣進(jìn)行分析。本文采用速率敏感性矩陣分析法，通過對矩陣進(jìn)行特征值和特征向量分析，特征值所對應(yīng)的反應(yīng)組中特征向量元素較大的值就代表了重要的基元反應(yīng)，據(jù)此刪除次要反應(yīng)而使反應(yīng)機(jī)理得到簡化。速率敏感性矩陣的元素定義為：

式中，fi代表物種i的凈生成速率，kj為反應(yīng)j的速率常數(shù)，ωj為第j個基元反應(yīng)的凈反應(yīng)速率，νi,j為物種i在第j個基元反應(yīng)中的化學(xué)計量系數(shù)。主組分分析法可以有效減少冗余反應(yīng)，進(jìn)一步縮減機(jī)理尺寸。

2.2.5溫度敏感性分析和產(chǎn)率分析法

對于燃燒系統(tǒng)而言，敏感性分析法是在給定的溫度、壓力等條件下，考察每一個基元反應(yīng)的速率常數(shù)的微小擾動所引起的系統(tǒng)特征量(如濃度、溫度、反應(yīng)速率)的變化，其值稱為敏感度。通過比較各敏感性度的量級，即可確定各基元反應(yīng)的相對重要性，由此忽略那些相對次要的反應(yīng)，從而得出簡化的化學(xué)反應(yīng)動力學(xué)模型。由于溫度對HCCI自燃著火過程最為敏感，本文采用溫度敏感性分析的方法33來衡量基元反應(yīng)的重要程度，歸一化對數(shù)敏感性系數(shù)定義為：

式中，ki為反應(yīng)i的速率常數(shù)，T為系統(tǒng)溫度。通過求解一系列非線性常微分方程組，可方便得出上述敏感性系數(shù)值。值得注意的是，敏感性分析法得到的只是在某一個時間點(diǎn)或者空間點(diǎn)上的局部結(jié)果，所得簡化機(jī)理也僅僅是詳細(xì)機(jī)理的一個子集，無法得到最大形式上的簡化。因此，有必要輔以物種產(chǎn)率分析來彌補(bǔ)上述缺陷。產(chǎn)率分析直接描述單一基元反應(yīng)對物種生成率或消耗率的貢獻(xiàn)值，因此可以迅速得出鏈傳播過程所經(jīng)歷的反應(yīng)。生成率和消耗率分別定義為：

通過溫度敏感性分析和產(chǎn)率分析兩種方法的耦合使用，作為最后一步機(jī)理簡化方法，計算代價雖然較小，卻可以最大限度地使反應(yīng)機(jī)理得以簡化。

2.3機(jī)理簡化策略

基于以上機(jī)理簡化方法的特點(diǎn)，本文提出了一套系統(tǒng)化多級機(jī)理簡化策略，如圖1所示。該策略包含了六種簡化方法和七個簡化階段。由于DRGEP方法具有簡單、高效的特點(diǎn)，于是被重復(fù)使用一次。圖1中陰影部分為自動機(jī)理簡化方法，其開發(fā)基于CHEMKIN-SENKIN模塊34,35，程序采用Fortran語言編寫。其余簡化方法也簡單易行，手動簡化效果十分顯著。Pepiot-Desjardins和Pitsch22指出，簡化方法的順序?qū)τ谧罱K簡化機(jī)理的規(guī)模影響可以忽略?；谶@一考慮，采用本文提出的簡化策略，直接應(yīng)用于汽油多組分表征燃料詳細(xì)反應(yīng)動力學(xué)模型的簡化。在簡化過程中，首先給定重要物種、簡化初始條件的取值范圍及閾值等，選取著火滯燃期、最高燃燒溫度和最大燃燒壓力為每一次簡化的基本標(biāo)準(zhǔn)。不同的簡化階段，選擇不同的誤差精度界限，反復(fù)比較階段性簡化機(jī)理與詳細(xì)機(jī)理的吻合程度，直到最終機(jī)理的最大相對誤差大于人為設(shè)定的著火滯燃期、最高燃燒溫度和最高燃燒壓力的誤差為止。具體誤差精度判定準(zhǔn)則為：

圖1　多級機(jī)理簡化流程圖Fig.1　Multi-stage mechanism reduction flowchart grey block:automatic mechanism reduction,white block:manual mechanism reduction

式中，τ為著火滯燃期，T為最高燃燒溫度，p為最高燃燒壓力，ετ、εT、εp分別為著火滯燃期、最高燃燒溫度和最高燃燒壓力的相對誤差。

3　結(jié)果與討論

3.1簡化機(jī)理的構(gòu)建

本文的目標(biāo)是構(gòu)建適用于HCCI發(fā)動機(jī)燃燒邊界條件，可用于CFD耦合計算的汽油多組分表征燃料簡化機(jī)理模型，所設(shè)置的參數(shù)取值范圍為溫度650－1400 K，壓力1×105－60×105Pa，化學(xué)計量比0.2－1.0，基本覆蓋整個HCCI發(fā)動機(jī)的運(yùn)行工況范圍。詳細(xì)反應(yīng)機(jī)理的特征，見2.1節(jié)。

DRGEP方法作為第一步，旨在考察組分間的強(qiáng)耦合關(guān)系和去除非重要的物種。在此，選取IC8H18、NC7H16、C5H10-2、C6H5CH3、O2、N2、H、CO、HO2、CO2和H2O為重要物種，在自燃著火點(diǎn)附近稠密取樣。著火滯燃期定義為初始溫度升高400 K所經(jīng)歷的時間，諸多研究者采用此方法來判斷著火滯燃期21,25。最大滯燃期相對誤差定義為40%，最高燃燒溫度和最高燃燒壓力相對誤差為20%，算法閾值按照對數(shù)尺度線性增加。通過不斷地循環(huán)迭代，直到獲得最小的滿足誤差精度的簡化機(jī)理。DRGEP方法所形成的最終機(jī)理包含509個物種、2625個基元反應(yīng)，顯見機(jī)理規(guī)模縮小了一半以上。

圖2　物種摩爾分?jǐn)?shù)分布圖Fig.2　Species mole fraction distribution

然后，采用峰值濃度分析方法，使用SENKIN模塊34，在定容狀態(tài)下，計算整個參數(shù)取值范圍中，每一個物種的峰值濃度分布情況，如圖2所示。設(shè)置1×10－7，1×10－8，1×10－9，1×10－10(摩爾分?jǐn)?shù))四個濃度值為閾值，依次將小于閾值的物種，按照峰值濃度大小升序排列，逐個移除低濃度物種。每一次移除之后，計算所形成的簡化機(jī)理與詳細(xì)機(jī)理的滯燃期相對誤差，如果相對誤差小于5%，或者有利于提高滯燃期預(yù)測的精度，則表明該物種對于著火點(diǎn)的預(yù)測無明顯影響，可安全移除；反之，當(dāng)移除某一低濃度物種，造成滯燃期相對誤差大于5%，則保留該物種。如此不斷地迭代計算，測試每一個低濃度物種的重要性，形成最終的簡化機(jī)理，包含352個物種和1802個基元反應(yīng)。

此后，采用線性同分異構(gòu)體集總方法，針對大分子物種的低溫反應(yīng)路徑，進(jìn)行合并處理。在此共劃分25個同分異構(gòu)體群簇，主要集總處理結(jié)果見表1。Luo等36的研究表明，在低溫反應(yīng)區(qū)域，同分異構(gòu)體物種之間的線性摩爾質(zhì)量分布并不是很明顯。本文采用統(tǒng)計平均方法，獲得每一個同分異構(gòu)體群簇中物種的摩爾質(zhì)量分?jǐn)?shù)，隨后盡可能細(xì)分同分異構(gòu)體群簇的大小，任意兩個同分異構(gòu)體物種只要存在強(qiáng)烈的線性關(guān)系，即可劃分一個群簇，以達(dá)到最大的集總處理能力而不破壞詳細(xì)機(jī)理的層級拓展性。每一次集總處理之后，保證滯燃期相對誤差在5%之內(nèi)，最高燃燒溫度和最高燃燒壓力相對誤差在10%之內(nèi)，最終簡化機(jī)理包含290個物種，1644個基元反應(yīng)。

表1　主要同分異構(gòu)體物種簇Table 1　Main lumped isomer groups

經(jīng)過峰值濃度分析和線性集總處理之后，物種之間耦合關(guān)系發(fā)生急劇的變化，衍生出一批冗余物種和反應(yīng)。因此，再次使用DRGEP方法，快速高效地移除72個物種、334個反應(yīng)。為了進(jìn)一步縮減機(jī)理的規(guī)模，依照作者此前的工作37，使用主組分分析法，在包含218個物種和1310個反應(yīng)的簡化機(jī)理基礎(chǔ)上，又可刪除113個反應(yīng)。

最后，采用溫度敏感性分析和產(chǎn)率分析相結(jié)合的方法，盡可能使機(jī)理規(guī)模達(dá)到可以直接和CFD耦合計算的程度，同時保持應(yīng)有的計算精度。溫度敏感性分析主要是在著火點(diǎn)附近存在劇烈溫度變化梯度的時刻進(jìn)行，覆蓋整個參數(shù)的取值范圍。首先設(shè)定一個敏感度閾值，當(dāng)基元反應(yīng)的敏感性系數(shù)值小于此閾值時，表明該反應(yīng)相對不重要，可以移除，否則，予以保留。在簡化過程中，滯燃期相對誤差設(shè)為5%，最高燃燒溫度和最高燃燒壓力相對誤差設(shè)為10%，每一次移除都要保證所形成的簡化機(jī)理在誤差精度之內(nèi)，結(jié)果共移除626個反應(yīng)，簡化效果較為顯著。產(chǎn)率分析亦是如此，但計算代價更小，通過物種生成率和消耗率的分析，可以清晰地得到整個機(jī)理的主干反應(yīng)路徑，如圖3所示。藉此進(jìn)而不斷篩選移除非主干邊緣物種，得到最終簡化機(jī)理，包含149個物種和414個基元反應(yīng)步驟。

3.2簡化機(jī)理的驗證

為了驗證上述簡化機(jī)理的有效性，本文以著火滯燃期和發(fā)動機(jī)缸壓、污染物排放以及放熱率為目標(biāo)，選擇激波管、快速壓縮機(jī)和增壓HCCI發(fā)動機(jī)的實驗數(shù)據(jù)，在隨后的章節(jié)中，進(jìn)行模擬計算和實驗數(shù)據(jù)的比較分析。

3.2.1著火滯燃期的實驗驗證

選取Gauthier等8在激波管中測量RD387汽油高溫滯燃期，以及Kukkadapu等9在快速壓縮機(jī)上進(jìn)行中低溫滯燃期測量的實驗條件，作為模擬計算的輸入?yún)?shù)，利用CHEMKIN-SENKIN模塊34,35進(jìn)行著火滯燃期的預(yù)測計算，其實驗結(jié)果和計算結(jié)果的對比見圖4(a,b)。由圖可見，表征燃料簡化機(jī)理的模擬計算結(jié)果與試驗結(jié)果符合較好。實際上，模擬計算不僅很好地預(yù)測了化學(xué)計量比對滯燃期的影響，而且還成功地揭示出氧化歷程中所經(jīng)歷的負(fù)溫度系數(shù)區(qū)(760－920 K)以及整個滯燃期的變化趨勢。在化學(xué)計量比分別為0.3和0.5，初始壓力20×105Pa，初始溫度小于800 K時，著火滯燃期計算值小于實驗值，兩者之間存在一定的偏差。這是因為在快速壓縮機(jī)的定容模擬計算過程中，忽略了壓縮沖程中存在的少量化學(xué)反應(yīng)和傳熱損失，致使著火滯燃期縮短，而機(jī)理本身的精度滿足要求。

圖4(c)所示為EGR(exhaust gas recirculation)條件下表征燃料計算的滯燃期與實驗值的對比。EGR率定義為：燃料/空氣的混合氣共占(100－X)%摩爾，其中X%為燃料/空氣完全燃燒的產(chǎn)物(CO2、H2O、N2和O2)的摩爾分?jǐn)?shù)。由圖可見，表征燃料簡化機(jī)理仍然給出了與實驗值較為吻合的結(jié)果，在化學(xué)計量比0.5、EGR率為20%的條件下，預(yù)測值略顯偏高，但仍處于工程應(yīng)用允許的誤差范圍內(nèi)。

3.2.2增壓HCCI發(fā)動機(jī)的實驗驗證

Dec和Yang10以RD387汽油為燃料進(jìn)行了一系列增壓HCCI發(fā)動機(jī)的燃燒及排放性能實驗，為研究表征燃料簡化機(jī)理模型在HCCI發(fā)動機(jī)中的適應(yīng)性，選取其中一組實驗數(shù)據(jù)進(jìn)行綜合比較。模擬計算主要集中于燃料化學(xué)因素的研究，盡可能的減小物理因素干擾。采用零維單區(qū)模型，假設(shè)缸內(nèi)混合氣的熱力學(xué)狀態(tài)，溫度、壓力、組分分布完全均勻，忽略缸內(nèi)殘余廢氣，缸壁及余隙效應(yīng)的影響。按照Sjoberg等38的建議，單區(qū)模型代表缸內(nèi)混合氣的“高溫核心區(qū)”，約占整個充量20%的質(zhì)量，由于該區(qū)首先發(fā)生自燃著火，使用單區(qū)模型可以精確預(yù)測HCCI發(fā)動機(jī)著火時刻，同時，由于其簡單、計算代價小，可以方便的進(jìn)行變參數(shù)匹配的研究。

圖4　不同的初始溫度、壓力、化學(xué)計量比和EGR率下，簡化機(jī)理與實驗數(shù)據(jù)對于著火滯燃期的比較Fig.4　Comparison of experimental and calculated ignition delay times under different temperatures, pressures,equivalence ratios,and EGR ratios conditions

Dec和Yang10為達(dá)到HCCI發(fā)動機(jī)的負(fù)荷上限及最佳燃燒相位，同時保持較高的熱效率，不斷調(diào)整進(jìn)氣壓力、進(jìn)氣溫度和EGR率協(xié)同控制缸內(nèi)充量的狀態(tài)，以實現(xiàn)穩(wěn)定燃燒。在模擬過程中，發(fā)動機(jī)結(jié)構(gòu)參數(shù)來自實驗發(fā)動機(jī)，降低壓縮比至13.3，來彌補(bǔ)壓縮沖程的傳熱損失，轉(zhuǎn)速保持1200 r?min－1不變，EGR氣體為燃料/空氣完全燃燒的產(chǎn)物CO2、H2O、N2和O2，隨著進(jìn)氣壓力的增加，調(diào)整進(jìn)氣溫度和EGR率來匹配實驗中的燃燒相位。

圖5(a,b)分別顯示了實驗和計算下所得到的缸內(nèi)壓力在不同進(jìn)氣壓力下隨時間的變化曲線。如圖可以看到，計算曲線與實驗曲線變化趨勢基本一致，燃燒相位(10%累積放熱時刻)均由自然進(jìn)氣狀態(tài)下的9°CAATDC(after top dead centre)，先推遲至180×105Pa下的15°CA ATDC，再提前至325×105Pa的9°CAATDC，這主要是由于隨著進(jìn)氣壓力的增加，著火時刻提前，缸內(nèi)最大爆發(fā)壓力升高，在(0－180)×105Pa范圍內(nèi)，通過不斷增加循環(huán)噴油量，降低進(jìn)氣溫度，來推遲燃燒相位，達(dá)到負(fù)荷上限；在(180－325)×105Pa范圍內(nèi)，保持最低進(jìn)氣溫度不變，進(jìn)一步添加冷EGR來稀釋缸內(nèi)的充量，減緩燃燒速率，延長燃燒持續(xù)期，來抵消增加進(jìn)氣壓力，提前著火時刻的效應(yīng)，同時保持上止點(diǎn)附近放熱，盡可能實現(xiàn)最大負(fù)荷上限，實質(zhì)上在高增壓(大于260×105Pa)下，缸內(nèi)混合氣接近于化學(xué)計量比狀態(tài)，因此，燃燒相位是進(jìn)氣壓力、進(jìn)氣溫度、冷EGR率三者相互競爭的結(jié)果，呈現(xiàn)出先推遲后提前的趨勢，簡化機(jī)理很好的捕捉到這一現(xiàn)象。值得注意的是，在計算中使用Mehl等29的詳細(xì)機(jī)理與本文所開發(fā)的簡化機(jī)理在不同進(jìn)氣壓力下所計算的結(jié)果完全相同，說明簡化機(jī)理可以很好地重現(xiàn)詳細(xì)機(jī)理的燃燒特性，在預(yù)測著火點(diǎn)和燃燒持續(xù)期方面，兩者效果一樣。

圖5　不同的進(jìn)氣壓力下，簡化機(jī)理和實驗數(shù)據(jù)缸內(nèi)壓力曲線的比較Fig.5　Cylinder pressure for different intake pressure with reduced mechanism and experimental data

圖6　不同的進(jìn)氣壓力下，簡化機(jī)理和實驗數(shù)據(jù)CO和UHC排放曲線的比較Fig.6　CO and UHC emissions for different intake pressures with reduced mechanism and experimental data

為進(jìn)一步驗證機(jī)理預(yù)測排放指標(biāo)的能力，圖6 (a,b)對比了實驗和計算所得到的CO和UHC濃度隨進(jìn)氣壓力的變化曲線。實驗和計算中均采用體積濃度單位10－6，即在一定的排氣容積中，有害排放物所占到的容積比例。計算中UHC排放指的是排放物中所有未燃碳?xì)浠衔锏捏w積分?jǐn)?shù)乘以該化合物中含碳原子的數(shù)目之后，加在一起的總和，與實驗測量方法相同。CO排放是碳?xì)淙剂显谌紵瘜W(xué)反應(yīng)過程中的中間產(chǎn)物和不完全燃燒產(chǎn)物之一，其生成機(jī)理十分復(fù)雜，主要受缸內(nèi)混合氣濃度和溫度的影響。UHC排放也是一種不完全燃燒的產(chǎn)物，其來源主要有不完全燃燒、壁面淬熄效應(yīng)、壁面油膜和積碳的吸附效應(yīng)等等，與發(fā)動機(jī)的結(jié)構(gòu)設(shè)計、物理化學(xué)因素都有較大的關(guān)系。從圖中可以看到，兩者變化的趨勢完全相同，但是在數(shù)量級上有一定差別，這是因為單區(qū)模型僅僅代表“高溫核心區(qū)”，忽略了許多因素的影響，缸內(nèi)燃燒較充分，CO和UHC排放要明顯降低。同時可以看到，對于CO來說，隨著進(jìn)氣壓力的升高而排放增大，這是因為增壓后，降低進(jìn)氣溫度和較大EGR率，使得缸內(nèi)平均溫度降低，燃燒效率下降，不完全燃燒而引起的；UHC濃度隨著進(jìn)氣壓力升高而降低，是由于增壓后較高的壁面溫度，促進(jìn)了冷壁面及縫隙區(qū)的氧化反應(yīng)。但是由于模型的局限性，單區(qū)模型只能定性給出污染物的排放趨勢，并不能完全衡量那種機(jī)理所預(yù)測的CO和UHC更為準(zhǔn)確，只有求助于多維模型和更多的實驗數(shù)據(jù)才能做進(jìn)一步的比較分析。

3.3簡化機(jī)理的分析

3.3.1中間溫度放熱的化學(xué)反應(yīng)動力學(xué)分析

HCCI燃燒方式的難點(diǎn)之一在于運(yùn)行工況向高負(fù)荷拓展，Dec和Yang10在實驗中使用高增壓和推遲燃燒相位的策略來達(dá)到最大負(fù)荷上限，避免了壓力升高率過大直接引起的敲缸現(xiàn)象，同時，研究發(fā)現(xiàn)，在主放熱階段之前，也就是整體燃燒過程的中后期，存在一個平緩的放熱階段，不同于低溫放熱階段的“雙峰”現(xiàn)象，呈現(xiàn)出單階段放熱，由于它介于低溫放熱和高溫放熱階段之間，被稱為中間溫度放熱(ITHR)39。

中間溫度放熱一般發(fā)生于850－1000 K之間，允許缸內(nèi)溫度在上止點(diǎn)之后持續(xù)升高，抵消了由于活塞下行而引起的局部充量冷卻效應(yīng)。因此，在增壓條件下，可以更多地推遲燃燒相位，增加循環(huán)噴油量，從而有效地拓展負(fù)荷上限，并保持穩(wěn)定燃燒，另一方面，中間溫度放熱可以減小傳熱損失，降低敲缸強(qiáng)度，從而提高熱效率。Dec等40的研究進(jìn)一步表明，中間溫度放熱具有強(qiáng)烈的燃料和壓力依賴性。

為了揭示中間溫度放熱的化學(xué)本質(zhì)，闡明不同組分對于中間溫度放熱的貢獻(xiàn)，本文選取Dec和Yang10的一組典型放熱率實驗數(shù)據(jù)，進(jìn)一步驗證了簡化機(jī)理的預(yù)測能力，并進(jìn)行了化學(xué)反應(yīng)動力學(xué)的分析。如圖7(a,b)所示，在典型的進(jìn)氣壓力和溫度條件下，機(jī)理計算結(jié)果與實驗數(shù)據(jù)基本上一致，但稍有不同的是，模擬計算的放熱率更加集中快速，而實驗值則相對平緩，這主要是由于單區(qū)模型的假設(shè)而引起的。同時由實驗和計算可以看到，在保持進(jìn)氣溫度不變，提高進(jìn)氣壓力時，出現(xiàn)了明顯了中間溫度放熱，進(jìn)一步降低進(jìn)氣溫度，中間溫度放熱更加顯著，在較早的時刻，放熱率曲線就會分離。

圖7　典型的進(jìn)氣溫度和壓力下，簡化機(jī)理和實驗數(shù)據(jù)放熱率曲線的比較Fig.7　Normalized heat release rate under the specified intake temperature and pressure conditions with reduced mechanism and experimental dataThe traces have been shifted to match the same ignition timing.

圖8　典型的進(jìn)氣溫度和壓力下，對整體放熱率貢獻(xiàn)最大的基元反應(yīng)分析Fig.8　Contributions of important reactions towards total heat release rate under the specified intake temperature and pressure conditionsITHR:intermediate temperature heat release,BTDC:before top dead centre

圖9　典型的進(jìn)氣溫度和壓力下，不同類別組分對整體放熱率影響的分析Fig.9　Contributions of different reaction groups to total heat release rate under the specified intake temperature and pressure conditions

表2　汽油表征燃料組分分布(摩爾分?jǐn)?shù))Table 2　Composition distribution(mole fraction)of different surrogate model

圖10　不同的初始溫度、壓力和化學(xué)計量比下，三組分和四組分表征燃料及實驗數(shù)據(jù)的壓力曲線比較Fig.10　Comparisons of experimental and calculated pressure traces using three-component,four-component surrogate fuel model under different temperature,pressures,equivalence ratios conditions

由于簡化機(jī)理表現(xiàn)出較好的預(yù)測能力，選擇在中間溫度放熱區(qū)早期(CA50前10°CA)，即圖7 (b)中桔色曲線的位置，進(jìn)行了關(guān)鍵基元反應(yīng)的放熱率分析。如圖8所示，分別給出了高壓低溫和高溫低壓下，對于放熱率影響最為重要的20個基元反應(yīng)?？梢钥吹?，R+O2反應(yīng)(R為烷基自由基)，反應(yīng)物脫氫反應(yīng)，HO2生成反應(yīng)，CH3氧化反應(yīng)起到?jīng)Q定性作用，值得注意的是，C8H17O2=C8H17+O2表現(xiàn)出極大的放熱能力。在高壓低溫工況下，基元反應(yīng)活性整體提高，加快了整個反應(yīng)系統(tǒng)的進(jìn)程，同時，由于較低的進(jìn)氣溫度，使得上止點(diǎn)溫度也相應(yīng)較低，有利于HO2的大量積累，通過再結(jié)合反應(yīng)生成H2O2，而反應(yīng)H2O2+M=OH+OH+M是吸熱反應(yīng)，此時反應(yīng)速率較慢，因此，會形成一個緩慢放熱的中間溫度區(qū)，大約占到總體累積放熱量的5%左右。

對于反應(yīng)進(jìn)一步歸類處理，分別探討2-C5H10、NC7H16、IC8H18、C6H5CH3、C0－C4對于中間溫度放熱率的貢獻(xiàn)，如圖9(a,b)所示。可以看到，隨著曲軸轉(zhuǎn)角的變化，不同類別的碳?xì)浠衔锓艧崧手饾u升高。在低壓高溫工況下，C0－C4小分子碳?xì)浠衔锓艧崞鸬街饕淖饔茫藭r整體放熱量較少，而在高壓低溫工況下，由于中間溫度放熱的影響，整體放熱時刻開始較早，放熱量也較大，這里更應(yīng)該注意的是，IC8H18的放熱會占據(jù)主導(dǎo)位置，其它分子類別的放熱也增大明顯，這與圖8的基元反應(yīng)放熱率分析結(jié)果一致，是由于高壓低溫下，反應(yīng)C8H17O2=C8H17+O2的活性顯著增加而引起的。

3.3.2烯烴對于四組分機(jī)理模型的影響分析

為了判明在汽油三組分表征燃料基礎(chǔ)上添加烯烴的必要性及其影響情況，借助于Kukkadapu等9利用快速壓縮機(jī)測得的RD387汽油低溫自燃著火特性以及壓力變化數(shù)據(jù)，對上述兩種表征燃料的低溫自燃著火特性進(jìn)行了模擬計算和比較研究。計算從壓縮沖程始點(diǎn)開始，計入了壓縮沖程中少量的化學(xué)反應(yīng)以及整個過程的傳熱損失，然后選取其中幾組典型的壓力數(shù)據(jù)進(jìn)行對比分析。燃料組分分布見表2，三組分模型surrogate A比例值來自于Gauthier等8的建議，四組分模型比例值來自于本文的工作，兩者均采用本文所構(gòu)建的同一簡化動力學(xué)模型。圖10分別示出在不同初始溫度、初始壓力和化學(xué)計量比條件下，三組分、四組分表征燃料自燃著火過程中的壓力變化計算值與RD387汽油自燃著火特性實測值的對比。由圖可見，三組分和四組分機(jī)理模型的計算曲線和實驗曲線變化趨勢較為一致，均能較好揭示出單階段或兩階段著火燃燒現(xiàn)象。值得注意的是，在相同的工況條件下，三組分表征燃料顯得化學(xué)活性更強(qiáng)，著火時刻提前，滯燃期變短；與之相比，四組分模型則與實驗數(shù)據(jù)吻合較好，尤其是在化學(xué)計量比等于0.3和1時，四組分模型的計算值與實驗曲線基本重合。顯然，添加2-戊烯使表征燃料模型的預(yù)測精度有所提高。究其原因，主要在于2-戊烯分子中含有雙鍵，其低溫反應(yīng)十分活躍，有利于增加OH基的數(shù)量，導(dǎo)致反應(yīng)沿著Waddington路徑進(jìn)行。此時，部分烷烴的鏈分支反應(yīng)被替代，致使整個反應(yīng)系統(tǒng)的活性有所降低，從而推遲了自燃著火時刻41。

圖11　10%組分?jǐn)z動引起的第一階段和整體著火滯燃期敏感度變化Fig.11　Percent sensitivity analysis on first-stage ignition and total ignition delays based on 10%composition variationinitial conditions of 760 K,20×105Pa and equivalence ratio 0.5 under constant volume conditons

為了進(jìn)一步說明單個組分對于著火滯燃期的影響，每一次額外添加10%的燃料組分，攝動分析第一階段和主燃燒階段的滯燃期變化，敏感性系數(shù)分別定義為：

式中，S1表示第一階段著火滯燃期敏感度，S表示主燃燒階段著火滯燃期敏感度。

如圖11所示，不同的組分對于兩個階段滯燃期的影響有很大區(qū)別，正庚烷顯著地提高了第一階段和主燃燒階段的反應(yīng)活性，異辛烷、甲苯和2-戊烯則分別減慢了兩階段的反應(yīng)活性。需要指出的是，2-戊烯對于第一階段的著火滯燃期影響最大，即烯烴組分的添加，有利于靈活調(diào)控兩階段著火的滯燃期。它啟示我們，通過合理優(yōu)化表征燃料的組分比例可實現(xiàn)發(fā)動機(jī)燃燒相位的可控和優(yōu)化。

4　結(jié)論

本文提出了一套新的系統(tǒng)化多級機(jī)理簡化方法，通過分析各個方法的特點(diǎn)及其適用性，確定了機(jī)理簡化的順序，進(jìn)而形成了有效的簡化策略，并將其應(yīng)用于大規(guī)模汽油多組分表征燃料詳細(xì)反應(yīng)機(jī)理的簡化，構(gòu)建了適用于HCCI發(fā)動機(jī)燃燒邊界下，可直接用于CFD耦合的簡化機(jī)理模型。主要結(jié)論如下：

(1)在充分調(diào)研國內(nèi)外汽油表征燃料機(jī)理的基礎(chǔ)上，選擇以美國RD387(research grade)汽油為研究對象，其詳細(xì)機(jī)理共包含1389物種、5935反應(yīng)，在較寬廣的實驗條件下和各種基礎(chǔ)反應(yīng)器中，對單組分、兩組分及混合物進(jìn)行系統(tǒng)的驗證，表現(xiàn)出較高的預(yù)測精度，作為機(jī)理簡化的起點(diǎn)。

(2)DRGEP方法的快速高效性，可用于大規(guī)模機(jī)理簡化的第一步，簡化效果顯著；峰值濃度分析法可以直接衡量單一物種及其反應(yīng)的對于簡化目標(biāo)的貢獻(xiàn)，方便實現(xiàn)自動機(jī)理簡化；線性同分異構(gòu)化集總法針對低溫反應(yīng)的同分異構(gòu)體路徑，可以有效合并處理，而不損失機(jī)理的精度；主組分分析法用于移除冗余基元反應(yīng)；溫度敏感性分析和產(chǎn)率分析的耦合使用，作為機(jī)理簡化的最后一步，盡可能達(dá)到簡化機(jī)理的極限，效果讓人滿意，最終機(jī)理包含149個組分、414個反應(yīng)。

(3)與激波管和快速壓縮機(jī)的實驗比較表明，簡化機(jī)理可以滿意地預(yù)測高壓下著火滯燃期，同時，在EGR工況下，表現(xiàn)出良好的預(yù)測能力。

(4)與增壓HCCI發(fā)動機(jī)的實驗比較發(fā)現(xiàn)，簡化機(jī)理在不同進(jìn)氣壓力下計算所得到的缸壓曲線及CO、UHC變化曲線與實驗趨勢完全相同，成功捕捉到燃燒相位先推遲后提前的變化。通過與詳細(xì)機(jī)理的比較，簡化機(jī)理和詳細(xì)機(jī)理完全重合，很好地再現(xiàn)了詳細(xì)機(jī)理的燃燒特征。

(5)在與放熱率的實驗比較發(fā)現(xiàn)，簡化機(jī)理在典型的壓力和溫度下，均表現(xiàn)出良好的預(yù)測能力，通過進(jìn)一步的放熱率化學(xué)動力學(xué)分析，發(fā)現(xiàn)了R+O2反應(yīng)是控制中間溫度區(qū)放熱的關(guān)鍵基元反應(yīng)，在高壓低溫下IC8H18放熱起到?jīng)Q定性的作用。

(6)通過快速壓縮機(jī)中的壓力變化曲線對比研究，表明汽油四組分表征燃料模型要優(yōu)于三組分燃料模型，2-戊烯的作用明顯，尤其是對于第一階段著火滯燃期的有顯著影響，為進(jìn)一步探索調(diào)和組分比例控制燃燒提供了一條新思路。

References

(1)Kalghatgi,G.T.Proc.Combust.Inst.2015,35,101.doi: 10.1016/j.proci.2014.10.002

(2)Reitz,R.D.Combust.Flame 2013,160,1.doi:org/10.1016/j. combustflame.2012.11.002

(3)Battin-Leclerc,F.;Blurock,E.;Bounaceur,R.;Fournet,R.; Glaude,P.;Herbinet,O.;Sirjean,B.;Warth,V.Chem.Soc.Rev. 2011,40,4762.doi:10.1039/C0CS00207K

(4)Dryer,F.L.Proc.Combust.Inst.2015,35,117.doi:org/ 10.1016/j.proci.2014.09.008

(5)Pitz,W.J.;Cernansky,N.P.;Dryer,F.L.;Egolfopoulos,F.N.; Farrell,J.T.;Friend,D.G.;Pitsch,H.SAE Tech.Pap.Ser. 2007,2007-01-0175.doi:10.4271/2007-01-0175

(6)Pera,C.;Knop,V.Fuel 2012,96,59.doi:10.1016/j. fuel.2012.01.008

(7)Badra,J.A.;Bokhumseen,N.;Mulla,N.;Sarathy,S.M.; Farooq,A.;Kalghatgi,G.T.;Gaillard,P.Fuel 2015,160,458. doi:org/10.1016/j.fuel.2015.08.007

(8)Gauthier,B.M.;Davidson,D.F.;Hanson,R.K.Combust. Flame 2004,139,300.doi:10.1016/j. combustflame.2004.08.015

(9)Kukkadapu,G.;Kumar,K.;Sung,C.J.;Mehl,M.;Pitz,W.J. Combust.Flame 2012,159,3066.doi:org/10.1016/j. combustflame.2012.05.008

(10)Dec,J.E.;Yang,Y.SAE Tech.Pap.Ser.2010,2010-01-1086. doi:10.4271/2010-01-1086

(11)Yang,Y.;Dec,J.E.;Dronniou,N.;Sjoberg,M.;Cannella,W. SAE Tech.Pap.Ser.2011,2011-01-1359.doi:10.4271/2011-01-1359

(12)Mehl,M.;Chen,J.Y.;Pitz,W.J.;Sarathy,S.M.;Westbrook, C.K.Energy Fuels 2011,25,5215.doi:org/10.1021/ ef201099y

(13)Perez,P.L.;Boehman,A.L.Energy Fuels 2012,26,6106.doi: org/10.1021/ef300503b

(14)Naik,C.V.;Pitz,W.J.;Westbrook,C.K.;Sjoberg,M.;Dec,J. E.;Orme,J.;Curran,H.J.;Simmie,J.M.SAE Tech.Pap.Ser. 2005,2005-01-3741.doi:10.4271/2005-01-3741

(15)Fikri,M.;Herzler,J.;Starke,R.;Schulz,C.;Roth,P.; Kalghatgi,G.T.Combust.Flame 2008,152,276.doi:10.1016/ j.combustflame.2007.07.010

(16)Andrae,J.C.G.Fuel 2008,87,2013.doi:10.1016/j. fuel.2007.09.010

(17)Yahyaoui,M.;Djeba?li-Chaumeix,N.;Dagaut,P.;Paillard,C. E.;Gail,S.Proc.Combust.Inst.2007,31,385.doi:10.1016/j. proci.2006.07.179

(18)Sarathy,S.M.;Kukkadapu,G.;Mehl,M.;Wang,W.J.;Javed, T.;Park,S.;Oehlschlaeger,M.A.;Farooq,A.;Pitz,W.J.; Sung C.J.Proc.Combust.Inst.2015,35,249.doi:org/ 10.1016/j.proci.2014.05.122

(19)Ahmed,A.;Goteng,G.;Shankar,V.S.B.;Qurashi,K.A.; Roberts,W.L.;Sarathy,S.M.Fuel 2015,143,290.doi:org/ 10.1016/j.fuel.2014.11.022

(20)Lu,T.F.;Law,C.K.Prog.Energy Combust.Sci.2009,35, 192.doi:10.1016/j.pecs.2008.10.002

(21)Lu,T.F.;Law,C.K.Proc.Combust.Inst.2005,30,1333.doi: 10.1016/j.proci.2004.08.145

(22)Pepiot-Desjardins,P.;Pitsch,H.Combust.Flame 2008,154, 67.doi:10.1016/j.combustflame.2007.10.020

(23)Sun,W.T.;Chen,Z.;Gou,X.L.;Ju,Y.G.Combust.Flame 2010,157,1298.doi:10.1016/j.combustflame.2010.03.006

(24)Luo,Z.Y.;Lu,T.F.;Maciaszek,M.J.;Som,S.;Longman,D. E.Energy Fuels 2010,24,6283.doi:10.1021/ef1012227

(25)Fang,Y.M.;Wang,Q.D.;Wang,F.;Li,X.Y.Acta Phys.-Chim.Sin.2012,28,2536.

[方亞梅,王全德,王繁,李象遠(yuǎn).物理化學(xué)學(xué)報,2012,28,2536.]doi:10.3866/PKU. WHXB201208201

(26)Hua,X.X.;Wang,J.B.;Wang,Q.D.;Tan,N.X.;Li,X.Y. Acta Phys.-Chim.Sin.2011,27,2755.

[華曉筱,王靜波,王全德,談寧馨,李象遠(yuǎn).物理化學(xué)學(xué)報,2011,27,2755.]doi: 10.3866/PKU.WHXB20112755

(27)Liu,Y.D.;Jia,M.;Xie,M.Z.;Pang,B.Energy Fuels 2013, 27,4899.doi:org/10.1021/ef4009955

(28)Zhang,Q.F.;Zheng,Z.L.;He,Z.W.;Wang,Y.Acta Phys.-Chim.Sin.2011,27,530.

[張慶峰,鄭朝蕾,何祖威,王迎.物理化學(xué)學(xué)報,2011,27,530.]doi:10.3866/PKU. WHXB20110334

(29)Mehl,M.;Pitz,W.J.;Westbrook,C.K.;Curran,H.J.Proc. Combust.Inst.2011,33,193.doi:10.1016/j.proci.2010.05.027

(30)Lu,T.F.;Ju,Y.G.;Law,C.K.Combust.Flame 2001,126, 1445.doi:10.1016/S0010-2180(01)00252-8

(31)Lu,T.F.;Law,C.K.Combust.Flame 2008,154,153.doi: 10.1016/j.combustflame.2007.11.013

(32)Turanyi,T.J.Math.Chem.1990,5,203.doi:10.1007/ BF01166355

(33)Maroteaux,F.;Noel,L.Combust.Flame 2006,146,246.doi: 10.1016/j.combustflame.2006.03.006

(34)Lutz,A.E.;Kee,R.J.;Miller,J.A.SENKIN:a Fortran Program for Predicting Homogeneous Gas Phase Chemical Kinetics with SensitivityAnalysis.Report SAND87-8248. Sandia,1987.

(35)Kee,R.J.;Grear,J.F.;Smooke,M.D.;Miller,J,A.Chemkin-II:a Fortran Chemical Kinetics Package for theAnalysis of Gas-Phase Chemical Kinetics.Report SAND89-8009.Sandia, 1989.

(36)Luo,Z.Y.;Plomer,M.;Lu,T.F.;Som,S.;Longman,D.E. Combust.Theor.Model.2012,16,369.doi:org/10.1080/ 13647830.2011.631034

(37)Xiao,G..;Zhang,Y.S.;Lang,J.Chinese Internal Combustion Engine Engineering 2013,34,20.

[肖干,張煜盛,郎靜.內(nèi)燃機(jī)工程,2013,34,20.]

(38)Sjoberg,M.;Dec,J.;Hwang,J.Y.SAE Tech.Pap.Ser.2007, 2007-01-0207.doi:10.4271/2007-01-0207.

(39)Vuilleumier,D.;Kozarac,D.;Mehl,M.;Saxena,S.;Pitz,W.; Dibble,R.;Chen,J.Y.;Sarathy,M.;Combust.Flame 2014, 161,680.doi:10.1016/j.combustflame.2013.10.008

(40)Yang,Y.;Dec,J.;Sjoberg,M.;Ji,C.S.Combust.Flame 2015, 162,4008.doi:10.1016/j.combustflame.2015.07.040

(41)Mehl,M.;Pitz,W.;Westbrook,C.K.;Yasunag,K.;Conroy, C.;Curran,J.Proc.Combust.Inst.2011,33,201.doi:10.1016/ j.proci.2010.05.040

Systematic Construction and Validation of the Reduced Chemical Kinetic Model of Gasoline Multi-Component Surrogate Fuel

XIAO Gan1ZHANG Yu-Sheng1,*JIANG Guang-Jun1,2
(1College of Energy and Power Engineering,Huazhong University of Science&Technology,Wuhan 430074,P.R.China;2College of Machine and Electron Engineering,Nanchang University,Nanchang 330031,P.R.China)

Asystematic multi-stage mechanism reduction strategy for performing skeletal reductions of gasoline four-component surrogate fuel is presented.The approach includes the directed relation graph with error propagation,peak concentration analysis,linear isomer lumping,principal component analysis,temperature sensitivity analysis and rate of production analysis.The final reduced mechanism comprises 149 species and 414 reactions with embedded cross-reactions,which is suitable for homogeneous charge compression ignition (HCCI)engine application.Comparisons between computational and experimental data including the shock tube and rapid compression machine,indicate that the new reduced mechanism can provide good predictability of the ignition delay over extensive parameter space.Applying the reduced mechanism to the HCCI single zone model also shows satisfactory combustion and emission characteristics of the boosted HCCI combustion. Further heat release analysis demonstrates that R+O2are the key reactions controlling the intermediate temperature heat release and under high pressure and low temperature conditions,iso-octane is the most important species resulting in a large portion of heat release.After the addition of 2-pentene,the new four component model displays better predictability than the three component model,especially relative to the firststage ignition delay.Based on these new findings,we can use different composition ratios to arbitrarily control the combustion phasing of HCCI combustion.

Gasoline surrogate fuel;Multi-stage mechanism reduction;HCCI engine; Intermediate temperature heat release;Chemical kinetic modeling

November 17,2015;Revised:January 25,2016;Published on Web:January 26,2016.*Corresponding author.Email:yszhang@hust.edu.cn;Tel:+86-18986143063. The project was supported by the National Natural Science Foundation of China(51176057).

O643

10.3866/PKU.WHXB201601261

國家自然科學(xué)基金(51176057)資助項目