梁毅， 張玉銀，， 李世琰， 徐斌

(1.上海交通大學(xué)機(jī)械與動(dòng)力工程學(xué)院， 上海 200240；2.河南科技大學(xué)交通與車(chē)輛工程學(xué)院， 河南 洛陽(yáng) 471000)

?

·性能研究·

基于化學(xué)反應(yīng)動(dòng)力學(xué)耦合G方程的定容彈爆震燃燒研究

梁毅1， 張玉銀1，2， 李世琰1， 徐斌2

(1.上海交通大學(xué)機(jī)械與動(dòng)力工程學(xué)院， 上海 200240；2.河南科技大學(xué)交通與車(chē)輛工程學(xué)院， 河南 洛陽(yáng) 471000)

由于爆震受多方面的因素共同影響，且這些因素往往是相互耦合在一起，直接在內(nèi)燃機(jī)上對(duì)各因素進(jìn)行解耦進(jìn)而研究單因素對(duì)爆震的影響幾乎不可能。針對(duì)上述問(wèn)題，首先基于化學(xué)反應(yīng)動(dòng)力學(xué)和G方程火焰面模型建立了長(zhǎng)方體燃燒彈內(nèi)的爆震燃燒三維CFD模型；然后在該模型基礎(chǔ)上把初始?jí)毫?、初始?xì)鉁睾统跏急跍?末端)設(shè)為獨(dú)立變量，系統(tǒng)地研究了各變量對(duì)正庚烷與空氣混合氣的爆震界限及爆震強(qiáng)度的影響。結(jié)果表明：當(dāng)初始?jí)毫Σ桓哂?.25 MPa時(shí)，初始?xì)鉁睾湍┒吮跍卦?00～700 K之間無(wú)論如何變化都不會(huì)發(fā)生爆震；當(dāng)初始?xì)鉁夭桓哂?50 K時(shí)，初始?jí)毫υ?.2～0.5 MPa之間、末端壁溫在300～700 K之間無(wú)論如何變化也都不會(huì)發(fā)生爆震。同時(shí)發(fā)現(xiàn)最大壓力振幅可作為發(fā)生爆震與否的判據(jù)，本研究中只要最大壓力振幅不高于0.02 MPa即不會(huì)發(fā)生爆震。

爆震; 定容燃燒彈; 化學(xué)反應(yīng)動(dòng)力學(xué)

對(duì)于內(nèi)燃機(jī)來(lái)說(shuō)爆震(Knocking,又稱(chēng)Deflagration)是一種非正常燃燒現(xiàn)象[1]。實(shí)踐表明，嚴(yán)重的爆震會(huì)影響整個(gè)內(nèi)燃機(jī)的性能，甚至對(duì)內(nèi)燃機(jī)造成不可挽回的機(jī)械破壞。

鑒于爆震現(xiàn)象對(duì)內(nèi)燃機(jī)的危害，學(xué)者對(duì)爆震現(xiàn)象做了很多的研究工作[2-6]。1968年，Affeleck等指出任何影響末端氣體傳熱的因素都對(duì)爆震產(chǎn)生影響[2]；2006年何鑫等在快速壓縮機(jī)上進(jìn)行了高速攝影試驗(yàn)，結(jié)果顯示燃料的摩爾分?jǐn)?shù)是決定自燃模式的關(guān)鍵因素[3]；2007年M. Poschl等在快速壓縮機(jī)上進(jìn)行研究，發(fā)現(xiàn)引起爆震的主要因素是溫度分布的不均勻性，且不同的分布情況會(huì)導(dǎo)致不同程度的爆震[4]；2011年Zahdeh等直接在內(nèi)燃機(jī)上研究進(jìn)氣溫度、壓縮比和增壓比等多個(gè)因素對(duì)爆震的影響，結(jié)果表明進(jìn)氣溫度的提高會(huì)導(dǎo)致末端混合氣更容易發(fā)生自燃[5]；2013年王志等在直噴汽油機(jī)上針對(duì)高負(fù)荷下的爆震進(jìn)行試驗(yàn)研究，證實(shí)了降低進(jìn)氣門(mén)關(guān)閉時(shí)的缸內(nèi)氣體溫度可抑制爆震[6]。然而，上述對(duì)爆震現(xiàn)象的研究均是直接在內(nèi)燃機(jī)或者快速壓縮機(jī)上實(shí)現(xiàn)的，導(dǎo)致直接在內(nèi)燃機(jī)上對(duì)影響爆震且相互耦合的各因素進(jìn)行解耦分析，研究單一因素對(duì)爆震的影響幾乎是不可能的。但是，單一因素對(duì)爆震界限及爆震強(qiáng)度的影響對(duì)于探究爆震燃燒機(jī)理和抑制內(nèi)燃機(jī)爆震又是非常重要的。

為了克服以往工作的不足，本研究設(shè)計(jì)了一種可以對(duì)影響內(nèi)燃機(jī)爆震的因素進(jìn)行解耦分析的長(zhǎng)方體爆震燃燒彈(以下稱(chēng)為長(zhǎng)方體燃燒彈)，該長(zhǎng)方體燃燒彈可以把初始混合氣溫度、壓力、初始末端壁面溫度解耦?；诨瘜W(xué)反應(yīng)動(dòng)力學(xué)耦合G方程建立了該長(zhǎng)方體燃燒彈內(nèi)的三維CFD模型，研究了初始混合氣溫度、壓力、初始末端壁面溫度分別作為單一因素變化時(shí)對(duì)爆震燃燒的影響規(guī)律，并獲得發(fā)生爆震與否的界限。本研究將為該類(lèi)燃燒裝置的設(shè)計(jì)與優(yōu)化和內(nèi)燃機(jī)抑制爆震燃燒提供理論依據(jù)。

1 模型建立與驗(yàn)證

1.1 長(zhǎng)方體燃燒彈仿真模型建立

本研究通過(guò)Converge軟件建立三維CFD的長(zhǎng)方體燃燒彈仿真模型。采用燃燒室為20 mm×20 mm×100 mm的小型臥式長(zhǎng)方體燃燒彈實(shí)現(xiàn)各因素分離，在固定其他因素不變的情況下，研究單一因素對(duì)爆震的影響。長(zhǎng)方體燃燒彈外形見(jiàn)圖1a，燃燒室見(jiàn)圖1b，其中A點(diǎn)為末端混合氣壓力(以下稱(chēng)為末端壓力)數(shù)據(jù)采樣位置，距離末端壁面1 mm。

為了捕捉到爆震時(shí)的快速燃燒過(guò)程，末端壓力采樣時(shí)間間隔為5×10-6s?；鸹ㄈ贾迷陂L(zhǎng)方體燃燒彈的B點(diǎn)位置(見(jiàn)圖1b)，距離右端壁面5 mm，點(diǎn)火時(shí)燃燒室內(nèi)為化學(xué)當(dāng)量比的正庚烷(C7H16)與空氣均質(zhì)混合氣，且混合氣溫度均勻、無(wú)壓力間斷。本研究分別采用兩個(gè)熱點(diǎn)模擬實(shí)際火花塞點(diǎn)火時(shí)的電容放電和電感放電過(guò)程，總能量為500 mJ的高能量點(diǎn)火?？紤]到既要節(jié)約計(jì)算成本又要保證計(jì)算的精度，采用笛卡爾網(wǎng)格，初始網(wǎng)格設(shè)置為2 mm ×2 mm×2 mm，計(jì)算過(guò)程中程序根據(jù)溫度梯度和速度梯度自適應(yīng)加密到最小網(wǎng)格邊長(zhǎng)0.5 mm，點(diǎn)火附近則固定加密到0.125 mm。有關(guān)網(wǎng)格設(shè)置參考了Yang[11]等建立的模型。

長(zhǎng)方體燃燒彈燃燒室壁面溫度的設(shè)置參考了實(shí)際汽油機(jī)的燃燒室、活塞頂壁面溫度值。Nakagawa等通過(guò)試驗(yàn)發(fā)現(xiàn)自燃始于燃燒室壁面，因此，初始末端壁面溫度(以下稱(chēng)為末端壁溫，Tw)可能是爆震重要的影響參數(shù)之一，本研究把末端壁溫設(shè)為變量[12]。燃燒室靠近A點(diǎn)的端面為離末端混合氣最近的壁面(類(lèi)似排氣門(mén)頂面，相當(dāng)于實(shí)際內(nèi)燃機(jī)中的局部熱點(diǎn))。除末端壁面外，其他壁面溫度與汽油機(jī)正常運(yùn)行時(shí)氣缸壁面溫度相當(dāng)，設(shè)置固定值440 K(見(jiàn)圖1b)。

本研究將系統(tǒng)地研究初始混合氣溫度(以下稱(chēng)為初始?xì)鉁兀琓g)、初始混合氣壓力(以下稱(chēng)為初始?jí)毫Γ琾i)和末端壁溫作為單一變量變化對(duì)爆震界限及爆震強(qiáng)度的影響。

圖1 長(zhǎng)方體燃燒彈仿真模型

1.2 長(zhǎng)方體燃燒彈仿真模型驗(yàn)證

為了保證該長(zhǎng)方體燃燒彈仿真模型模擬的可靠性，對(duì)模型中的化學(xué)反應(yīng)機(jī)理、子模型(如燃燒模型、湍流模型等)及其參數(shù)設(shè)置的正確性進(jìn)行了驗(yàn)證。

參考Tinaut等的試驗(yàn)結(jié)果，利用與長(zhǎng)方體燃燒器模型一致的化學(xué)反應(yīng)機(jī)理、子模型及其參數(shù)設(shè)置建立了球形燃燒器三維CFD模型[13]。Tinaut的試驗(yàn)是在直徑為200 mm的球形燃燒器上進(jìn)行，基本參數(shù)見(jiàn)表1[13]。在球形燃燒器中壓力隨時(shí)間變化的實(shí)測(cè)值與本研究計(jì)算值的比較情況見(jiàn)圖2，CFD仿真模型模擬壓力和試驗(yàn)結(jié)果基本吻合，說(shuō)明球形燃燒器的三維CFD仿真模型基本上能正確反映出試驗(yàn)時(shí)的燃燒特性，從而驗(yàn)證了長(zhǎng)方體燃燒彈仿真模型的子模型及其參數(shù)設(shè)置具有充足的可靠性，可用于進(jìn)一步的燃燒過(guò)程模擬研究。

表1 球形燃燒器初始參數(shù)

圖2 燃燒器壓力仿真值與試驗(yàn)值的對(duì)比

2 結(jié)果分析

參考實(shí)際汽油機(jī)發(fā)生爆震時(shí)的缸內(nèi)環(huán)境，考慮到長(zhǎng)方體燃燒彈火花點(diǎn)火后形成的火焰沿長(zhǎng)軸方向傳播時(shí)壓力升高，未燃混合氣不斷受到壓縮，這個(gè)過(guò)程可以看作汽油機(jī)內(nèi)活塞的壓縮過(guò)程。長(zhǎng)方體燃燒彈內(nèi)爆震發(fā)生前的壓力與汽油機(jī)壓縮終了的壓力相當(dāng)，如果把初始?jí)毫υO(shè)置在0.4 MPa，爆震前的壓力大約為2 MPa。本研究的初始?jí)毫υO(shè)置為0.2～0.5 MPa，最小間隔為0.05 MPa。自然吸氣汽油機(jī)進(jìn)氣溫度較低，而增壓并引入EGR的汽油機(jī)進(jìn)氣溫度較高，本研究中把初始?xì)鉁卦O(shè)置在300～700 K之間，間隔為50 K。排氣門(mén)頂面溫度在冷起動(dòng)時(shí)較低，而在正常運(yùn)行一段時(shí)間后溫度又很高，因此將末端壁溫設(shè)在300～700 K之間，間隔為50 K。

2.1 爆震評(píng)價(jià)指標(biāo)

為了定量分析爆震，需要確定爆震的評(píng)價(jià)指標(biāo)，判定是否發(fā)生爆震及爆震強(qiáng)度。Oppenheim[14]給出了較全面的爆震定義：爆震是一種特殊的燃燒方式，它是自發(fā)和隨機(jī)地產(chǎn)生，導(dǎo)致銳利的壓力脈沖。本研究提出以長(zhǎng)方體燃燒彈內(nèi)末端混合氣(圖1b中的A點(diǎn))的最大壓力振幅作為爆震評(píng)價(jià)指標(biāo)。最大壓力振幅定義為燃燒室末端壓力與平均壓力之差的絕對(duì)值的最大值，平均壓力為燃燒室內(nèi)腔壓力平均值。

圖3和圖4分別示出了長(zhǎng)方體燃燒彈內(nèi)爆震與非爆震工況下的壓力變化情況。由圖可知，發(fā)生爆震時(shí)末端壓力出現(xiàn)劇烈的壓力振蕩，壓力振幅曲線(xiàn)存在多個(gè)壓力脈沖；而非爆震工況下則不出現(xiàn)劇烈的壓力振蕩，壓力振幅曲線(xiàn)中不存在壓力脈沖。

圖3 爆震工況壓力曲線(xiàn)

圖4 非爆震工況壓力曲線(xiàn)

圖5示出了爆震與非爆震工況下長(zhǎng)方體燃燒彈燃燒室內(nèi)火焰的溫度分布對(duì)比。由圖5可看出，火花點(diǎn)火形成的火焰面往前推進(jìn)一段時(shí)間后，火焰前鋒面呈郁金花狀向前傳播，將燃燒室劃分為己燃區(qū)和末燃區(qū)；初始?xì)鉁?00 K工況中的末端混合氣從7.82 ms開(kāi)始發(fā)生自燃，同時(shí)，圖3中長(zhǎng)方體燃燒彈的末端壓力開(kāi)始出現(xiàn)劇烈振蕩，這是由于爆震壓力波在燃燒室內(nèi)來(lái)回反射的結(jié)果。顯然，明顯的壓力振幅標(biāo)志著爆震的開(kāi)始，因此采用最大壓力振幅來(lái)評(píng)價(jià)爆震開(kāi)始是合理的，后面還將進(jìn)一步說(shuō)明用最大壓力振幅作為判定是否發(fā)生爆震和爆震強(qiáng)度的依據(jù)。

圖5 爆震與非爆震工況長(zhǎng)方體燃燒彈燃燒室內(nèi)溫度分布

2.2 初始?jí)毫?duì)爆震的影響

為了研究初始?jí)毫?duì)爆震燃燒特性的影響，在同一末端壁溫、初始?xì)鉁貤l件下選取不同的初始?jí)毫M(jìn)行計(jì)算。長(zhǎng)方體燃燒彈的初始參數(shù)見(jiàn)表2。

表2 初始?jí)毫Φ淖兓秶?/p>

初始?jí)毫?duì)爆震壓力波的影響見(jiàn)圖6。圖6表明：在初始?jí)毫Σ怀^(guò)0.25 MPa的工況，末端壓力不出現(xiàn)明顯振蕩；在初始?jí)毫?.3 MPa及更高的工況時(shí)，隨初始?jí)毫μ岣撸饡r(shí)刻明顯提前，壓力振幅明顯增大。

圖6 初始?jí)毫?duì)爆震壓力波的影響

初始?jí)毫?duì)爆震強(qiáng)度的影響見(jiàn)圖7。從圖7可以看出：初始?jí)毫υ?.3～0.5 MPa之間變化時(shí)，隨初始?jí)毫μ岣咦畲髩毫φ穹焖僭龃蟆＝Y(jié)合圖6和圖7可知：初始?jí)毫?.25 MPa，0.2 MPa時(shí)末端壓力無(wú)明顯的振蕩，不發(fā)生爆震，最大壓力振幅均不大于0.004 MPa；初始?jí)毫Σ坏陀?.3 MPa時(shí)末端壓力有明顯的振蕩，發(fā)生爆震，最大壓力振幅均不小于0.24 MPa，且初始?jí)毫υ礁?，末端壓力振蕩幅度越大，爆震越?qiáng)烈，對(duì)應(yīng)的最大壓力振幅越大。顯然，以初始?jí)毫ψ鳛閱我蛔兞?，爆震與非爆震的工況點(diǎn)的最大壓力振幅分別在兩個(gè)數(shù)量級(jí)上。這里取最大壓力振幅0.02 MPa為閾值，認(rèn)為最大壓力振幅超過(guò)0.02 MPa時(shí)為爆震，小于0.02 MPa時(shí)為非爆

圖7 初始?jí)毫?duì)爆震強(qiáng)度的影響

震, 并認(rèn)為壓力振幅首次達(dá)到0.02 MPa時(shí)刻為爆震時(shí)刻，且最大壓力振幅越大爆震越強(qiáng)烈。因此，使用最大壓力振幅來(lái)判定爆震強(qiáng)度是合理的，后續(xù)結(jié)果也將證實(shí)這一判據(jù)的合理性。

2.3 初始?xì)鉁貙?duì)爆震的影響

為了研究初始?xì)鉁貙?duì)爆震燃燒特性的影響，在同一末端壁溫、初始?jí)毫l件下選取不同的初始?xì)鉁胤植键c(diǎn)進(jìn)行研究。初始?xì)鉁刈兓瘯r(shí)長(zhǎng)方體燃燒彈的初始參數(shù)見(jiàn)表3。

表3 初始?xì)鉁氐淖兓秶?/p>

初始?xì)鉁貙?duì)爆震壓力波和末端壓力的影響見(jiàn)圖8。圖8表明：當(dāng)初始?xì)鉁夭怀^(guò)500 K時(shí)，末端壓力不出現(xiàn)明顯振蕩，上升平緩；當(dāng)初始?xì)鉁剡_(dá)到550 K時(shí)，隨初始?xì)鉁靥岣撸饡r(shí)刻明顯提前，壓力振幅明顯增大。

圖8 初始?xì)鉁貙?duì)爆震壓力波的影響

初始?xì)鉁貙?duì)爆震強(qiáng)度的影響見(jiàn)圖9。從圖9可以看出：初始?xì)鉁卦?00～700 K之間變化時(shí)，隨初始?xì)鉁靥岣咦畲髩毫φ穹焖僭龃?。結(jié)合圖8和圖9可知：初始?xì)鉁氐陀诨虻扔?00 K時(shí)末端壓力無(wú)明顯的振蕩，不發(fā)生爆震，最大壓力振幅均不大于0.016 MPa；當(dāng)初始?xì)鉁馗哂?00 K時(shí)末端壓力有明顯的振蕩，發(fā)生爆震，最大壓力振幅均不小于0.18 MPa。顯然，以初始?xì)鉁刈鳛閱我蛔兞繒r(shí)，最大壓力振幅閾值0.02 MPa可以用于判斷爆震與否，且最大壓力振幅越大爆震越強(qiáng)烈。

圖9 初始?xì)鉁貙?duì)爆震強(qiáng)度的影響

2.4 末端壁面溫度對(duì)爆震的影響

為了研究末端壁溫對(duì)爆震燃燒特性的影響，在同一初始?xì)鉁?、初始?jí)毫l件下選取不同的末端壁溫分布點(diǎn)進(jìn)行研究。末端壁溫變化時(shí)長(zhǎng)方體燃燒彈的初始參數(shù)見(jiàn)表4。

表4 末端壁溫的變化范圍

末端壁溫對(duì)爆震壓力波的影響見(jiàn)圖10。圖10表明：末端壁溫不超過(guò)450 K時(shí)，各工況下末端混合氣的壓力曲線(xiàn)變化不大；末端壁溫升高到500 K及以上時(shí)，末端壓力出現(xiàn)明顯的振蕩。

圖10 末端壁溫對(duì)爆震壓力波的影響

末端壁溫對(duì)爆震強(qiáng)度的影響見(jiàn)圖11。從圖11可以看出：末端壁溫在450 K及以下時(shí)，最大壓力振幅均較小(不超過(guò)0.01 MPa)；末端壁溫為500 K及以上時(shí)，最大壓力振幅均較大(不低于0.59 MPa)。

結(jié)合圖10和圖11可知，末端壁溫為450 K及以下時(shí)末端壓力無(wú)明顯的振蕩，不發(fā)生爆震，最大壓力振幅均不大于0.008 MPa；末端壁溫在500 K及以上時(shí)末端壓力出現(xiàn)明顯的振蕩，發(fā)生爆震，最大壓力振幅均不小于0.59 MPa。顯然，以末端壁溫作為單一變量時(shí)，最大壓力振幅閾值0.02 MPa可以用于判斷爆震與否，且最大壓力振幅越大爆震越強(qiáng)烈。

圖11 末端壁溫對(duì)爆震強(qiáng)度的影響

2.5 爆震閾值與極限

上述仿真計(jì)算結(jié)果表明：在所有非爆震工況中，最大壓力振幅均小于0.02 MPa；在所有爆震工況中，最大壓力振幅均大于0.02 MPa。因此，判斷爆震與否的依據(jù)可以采用最大壓力振幅，認(rèn)為最大壓力振幅超過(guò)0.02 MPa為爆震，最大壓力振幅越大爆震越強(qiáng)烈，同時(shí)，爆震開(kāi)始時(shí)刻可以從壓力振幅首次達(dá)到0.02 MPa時(shí)判定。

通過(guò)上述閾值獲得了爆震界限(見(jiàn)圖12)，本研究計(jì)算點(diǎn)均在黑色線(xiàn)圍成的四方體框圖內(nèi)。從圖12可以看出：當(dāng)初始?jí)毫Σ桓哂?.25 MPa時(shí)，初始?xì)鉁睾湍┒吮跍卦?00～700 K范圍內(nèi)不管怎樣變化都不會(huì)發(fā)生爆震；在初始?xì)鉁夭桓哂?50 K時(shí)，初始?jí)毫υ?.2～0.5 MPa之間、末端壁溫在300～700 K范圍內(nèi)無(wú)論如何變動(dòng)也都不會(huì)發(fā)生爆震。因此，初始?jí)毫统跏細(xì)鉁厥菦Q定爆震與否的主要因素，而末端壁溫對(duì)爆震界限的影響相對(duì)較小。

圖12 爆震界限

3 結(jié)論

a) 初始?jí)毫统跏細(xì)鉁厥菦Q定爆震與否的關(guān)鍵因素，末端壁溫對(duì)爆震界限的影響相對(duì)較小;當(dāng)初始?jí)毫Σ桓哂?.25 MPa時(shí)，初始?xì)鉁睾湍┒吮跍卦?00～700 K范圍內(nèi)無(wú)論如何變化都不會(huì)發(fā)生爆震；當(dāng)初始?xì)鉁夭桓哂?50 K時(shí)，初始?jí)毫υ?.2～0.5 MPa范圍、末端壁溫在300～700 K范圍內(nèi)無(wú)論如何變化時(shí)都不會(huì)發(fā)生爆震;

b) 隨著初始?jí)毫?、初始?xì)鉁靥岣撸饡r(shí)刻提前，爆震強(qiáng)度(最大壓力振幅)增大;

c) 判斷爆震與否的依據(jù)可以采用最大壓力振幅，本研究中最大壓力振幅在0.02 MPa以上時(shí)發(fā)生爆震；爆震與非爆震時(shí)的最大壓力振幅相差一個(gè)數(shù)量級(jí)以上。

[1] 武得鈺, 顧笑映, 李繼軍, 等. 火花點(diǎn)火發(fā)動(dòng)機(jī)爆震燃燒研究的發(fā)展與現(xiàn)狀 (Ⅰ)[J].車(chē)用發(fā)動(dòng)機(jī),1995(5):1-6.

[2] Affleck W S, Fish A. Knock: Flame Acceleration or Spontaneous Ignition[J].Combustion and Flame,1968,12(3):243-252.

[3] Walton S M, He X, Zigler B T, et al. An Experimental Investigation of Iso-octane Ignition Phenomena[J].Combustion and Flame,2007,150(3):246-262.

[4] Poschl M, Sattelmayer T. Influence of Temperature Inhomogeneities on Knocking Combustion[J].Combustion and Flame,2008,153(4):562-573.

[5] Zahdeh A,Rothenberger P,Nguyen W,et al. Fun-damental Approach to Investigate Pre-Ignition in Boosted SI Engines[J].SAE Paper 2011-01-0340.

[6] 王志，徐雅齊，王建昕，等.增壓直噴汽油機(jī)掃氣抑制爆震試驗(yàn)研究及模擬解析[J].內(nèi)燃機(jī)工程,2013,34(4):13-17.

[7] Liang L,Reitz R D, Iyer C O, et al.Modeling knock in spark-ignition engines using a G-equation combustion model incorporating detailed chemical kinetics[C].SAE Paper 2007-01-0165.

[8] LIU Yaodong,JIA Ming,XIE Maozhao,et al. En-hancement on a skeletal kinetic model for primary reference fuel oxidation by using a semidecoupling methodology[J].Energy & Fuels,2012,26(12):7069-7083.

[9] Han Z, Reitz R D.Turbulence Modeling of Internal Combustion Engines using RNG κ-ε Models[J].Combustion Science and Technology,1995,106(4-6):267-295.

[10] Han Z,Reitz R D.A Temperature Wall Function Formulation for Variable-density Turbulent Flows with Application to Engine Convective Heat Transfer Modeling[J].International Journal of Heat and Mass Transfer,1997,40(3):613-625.

[11] Yang X,Solomon A,Kuo T W. Ignition and combustion simulations of spray-guided SIDI engine using Arrhenius combustion with spark-energy deposition model[C].SAE Paper 2012-01-0147.

[12] Nakagawa Y,Takagi Y,Itoh T.Laser shadow graphic analysis of knocking in SI engine[C].SAE Paper 845001,1984.

[13] Tinaut F V,Reyes M,Giménez B, et al. Measurements of OH and CH Chemiluminescence in Premixed Flames in a Constant Volume Combustion Bomb under Autoignition Conditions[J].Energy & Fuels,2010,25(1):119-129.

[14] Oppenheim A K. The knock syndrome-its cures and its victims[C].SAE Paper 841339,1984.

[編輯: 李建新]

Knock Combustion of Constant Volume Bomb Based on Chemical Kinetics Coupling G-equation

LIANG Yi1, ZHANG Yuyin1,2, LI Shiyan1, XU Bin2

(1.School of Mechanical Engineering, Shanghai Jiaotong University, Shanghai 200240, China; 2. School of Vehicle and Transportation, Henan University of Science and Technology, Luoyang 471000, China)

It was extremely difficult to analyze the effect of single factor on knock due to the collective action and mutual coupling of knock influencing factors. For the above problem, a 3D CFD model of knock combustion in a rectangular constant volume bomb was built based on the chemical kinetics and G-equation flame model. Then the influence of single variable such as initial pressure, initial mixture and wall temperature on knock limit and intensity of n-heptane and air mixture was researched systematically. The results show that the knock will not occur when the initial temperature of mixture and wall is between 300 K and 700 K and the initial pressure is lower than 0.25 MPa; and when the initial temperature of mixture is lower than 450 K, the initial temperature of wall is between 300 K and 700 K and the initial pressure is between 0.2 MPa and 0.5 MPa. Besides, the maximum fluctuation amplitude of pressure can be used as knock criterion and its upper limit is 0.02 MPa.

knock combustion; constant volume bomb; chemical kinetics

2016-01-11;

2016-01-31

內(nèi)燃機(jī)燃燒學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室開(kāi)放基金資助項(xiàng)目(K2014-2)；上海市教育委員會(huì)科研創(chuàng)新項(xiàng)目(14ZZ022)

梁毅(1988—)，男，碩士，主要研究方向?yàn)楸鹑紵^(guò)程；liangyisjtu@126.com。

通訊簡(jiǎn)介： 張玉銀(1962—)，男，教授，博士，主要研究方向?yàn)閲婌F燃燒及發(fā)動(dòng)機(jī)燃燒過(guò)程等；yuyinzhang@sjtu.edu.cn。

10.3969/j.issn.1001-2222.2016.02.004

TK411.22

B

1001-2222(2016)02-0017-06