曾海翔， 王 平， SHROTRIYA Prshnt， 姜霖松 MURUGESAN Meentchidevi

(江蘇大學(xué) a. 能源與動(dòng)力工程學(xué)院; b. 能源研究院，江蘇 鎮(zhèn)江 212013)

部分預(yù)混燃燒的特點(diǎn)是混合氣在空間和時(shí)間上具有可變的當(dāng)量比，且反應(yīng)發(fā)生在若干子區(qū)域，因而同時(shí)兼具擴(kuò)散火焰及預(yù)混火焰的特性[1].通過(guò)部分預(yù)混燃燒技術(shù)既可以降低燃燒過(guò)程中污染物的排放，又避免了預(yù)混火焰容易發(fā)生的回火問(wèn)題[2-3].

早期工作中，對(duì)部分預(yù)混火焰主要通過(guò)實(shí)驗(yàn)測(cè)量的方法加以研究.文獻(xiàn)[4]在高低不同雷諾數(shù)下的甲烷部分預(yù)混對(duì)沖火焰實(shí)驗(yàn)中，揭示了火焰熄滅主要受大尺度的渦流控制.文獻(xiàn)[5]研究了正常和微重力條件對(duì)部分預(yù)混層流火焰的舉升特性及火焰結(jié)構(gòu)的影響.文獻(xiàn)[6]在同心錐形燃燒器中產(chǎn)生丙烷射流部分預(yù)混湍流火焰，研究了伴流對(duì)火焰穩(wěn)定性和火焰結(jié)構(gòu)的影響.文獻(xiàn)[7]就進(jìn)口速度和當(dāng)量比振蕩對(duì)部分預(yù)混旋流火焰的擾動(dòng)進(jìn)行了實(shí)驗(yàn)探究，結(jié)果表明通過(guò)改變噴油位置、噴油器阻抗和平均噴嘴速度，可以有效地控制燃燒室進(jìn)口處的當(dāng)量比擾動(dòng)和進(jìn)口速度之間的相位差，進(jìn)而抑制或放大自激壓力振蕩的強(qiáng)度.文獻(xiàn)[8]采用粒子圖像測(cè)速法和平面激光誘導(dǎo)熒光(PLIF)技術(shù)，研究了稀薄部分預(yù)混湍流旋流火焰中流動(dòng)與反應(yīng)之間的非穩(wěn)態(tài)相互作用，結(jié)果表明非定常流場(chǎng)的主導(dǎo)結(jié)構(gòu)受剪切層中的進(jìn)動(dòng)渦核控制，并對(duì)火焰的穩(wěn)定起著重要作用.

得益于計(jì)算機(jī)性能的提高和數(shù)值模擬技術(shù)的發(fā)展，對(duì)湍流火焰的數(shù)值模擬逐漸成為一種常用的手段.由于燃燒過(guò)程中的火焰鋒面區(qū)域發(fā)生強(qiáng)烈化學(xué)反應(yīng)[9]，同時(shí)伴隨著流場(chǎng)的輸運(yùn)和混合，這對(duì)火焰面的建模提出了極大的挑戰(zhàn).文獻(xiàn)[10]利用進(jìn)度變量、混合分?jǐn)?shù)和標(biāo)量耗散率的聯(lián)合概率密度函數(shù)(PDF)統(tǒng)計(jì)將非預(yù)混湍流燃燒模擬中層流火焰的概念推廣到部分預(yù)混區(qū)域，并通過(guò)湍流擴(kuò)散舉升火焰和預(yù)混本生火焰的數(shù)值模擬證明了這一方法的廣泛適用性.文獻(xiàn)[11]將大渦模擬與條件矩封閉(CMC)方法進(jìn)行耦合，通過(guò)計(jì)算湍流部分預(yù)混二甲醚射流火焰，證明了CMC方法在復(fù)雜化學(xué)條件下仍可以精確地預(yù)測(cè)湍流-化學(xué)之間的相互作用.文獻(xiàn)[12]在小火焰控制方程中通過(guò)引入火焰指數(shù)這一概念來(lái)識(shí)別不同區(qū)域的火焰狀態(tài)，建立了適用于湍流噴霧燃燒大渦模擬的部分預(yù)混小火焰模型.

文獻(xiàn)[13]首次提出“人工增厚火焰”(ATF)模型，多年以來(lái)，世界各地的專(zhuān)家和學(xué)者廣泛對(duì)其進(jìn)行了論證和改進(jìn).文獻(xiàn)[14]通過(guò)引入效率函數(shù)解決了增厚后的火焰褶皺問(wèn)題.文獻(xiàn)[15]考慮到需要增厚的火焰位置只發(fā)生在火焰的鋒面處，對(duì)于遠(yuǎn)離火焰區(qū)的位置無(wú)需增厚處理，于是引進(jìn)了動(dòng)態(tài)增厚因子，這極大地拓展了增厚火焰模型的適用性.文獻(xiàn)[16]使用動(dòng)態(tài)增厚火焰(DTF)模型結(jié)合小火焰生成流形(FGM)化學(xué)制表方法對(duì)旋流預(yù)混燃燒器進(jìn)行模擬驗(yàn)證，結(jié)果與實(shí)驗(yàn)值吻合良好.文獻(xiàn)[17]為研究燃?xì)廨啓C(jī)內(nèi)部的燃燒不穩(wěn)定性與燃料混合狀態(tài)的關(guān)聯(lián)，分別采用完全預(yù)混和部分預(yù)混兩種假設(shè)，并利用DTF模型對(duì)PRECCINSTA燃燒器當(dāng)量比為0.7和0.83兩種工況進(jìn)行驗(yàn)證；在部分預(yù)混方案中，針對(duì)部分預(yù)混火焰的特性對(duì)化學(xué)反應(yīng)指前因子進(jìn)行修正，進(jìn)而調(diào)整火焰的傳播速度和厚度.結(jié)果顯示在當(dāng)量比為0.7(接近貧燃極限)的工況下，采用部分預(yù)混假設(shè)方案捕捉到了自激不穩(wěn)定現(xiàn)象和頻譜，而完全預(yù)混假設(shè)方案與實(shí)驗(yàn)結(jié)果差距較大.

DTF模型在保留火焰層流燃燒速度的同時(shí)，將火焰面加厚到了網(wǎng)格可解的尺度，通過(guò)直接求解化學(xué)反應(yīng)方程即可獲得火焰的詳細(xì)信息，避免了模型封閉的困難.理論上，DTF模型廣泛適用于各種燃燒模式，然而，DTF模型對(duì)部分預(yù)混燃火焰的預(yù)測(cè)還少有報(bào)道.這是由于傳統(tǒng)DTF模型在使用過(guò)程中需要預(yù)設(shè)增厚因子、層流火焰?zhèn)鞑ニ俣燃皩恿骰鹧婧穸?個(gè)參數(shù).在處理擴(kuò)散火焰時(shí)，可以認(rèn)為著火只發(fā)生在化學(xué)恰當(dāng)比的對(duì)沖區(qū)域，故上述3個(gè)參數(shù)可以按照化學(xué)恰當(dāng)比的層流火焰進(jìn)行設(shè)置；當(dāng)處理著火界限內(nèi)的某個(gè)預(yù)混火焰時(shí)，同樣可以根據(jù)其具體的化學(xué)當(dāng)量比來(lái)設(shè)置對(duì)應(yīng)的參數(shù).然而在部分預(yù)混火焰中，可燃物的化學(xué)當(dāng)量比隨時(shí)間和流場(chǎng)位置不斷地演化，依靠單一化學(xué)當(dāng)量比的建模方式已不太適合，此時(shí)需要對(duì)著火界限內(nèi)的層流火焰燃燒特性進(jìn)行動(dòng)態(tài)修正.

另外一個(gè)問(wèn)題是對(duì)部分預(yù)混火焰的火焰位置的捕捉.在之前的研究中，通過(guò)監(jiān)測(cè)CO2的生成反應(yīng)[18]或是基于CO2的質(zhì)量分?jǐn)?shù)構(gòu)建進(jìn)度變量[19-20]所捕捉的火焰面都不甚理想.而由文獻(xiàn)[21]提出的基于熱放率捕捉火焰面的方法表現(xiàn)良好，故本文所使用的DTF模型同樣利用反應(yīng)放熱率方法來(lái)識(shí)別火焰位置.

由于常溫常壓下，火焰的層流燃燒速度隨當(dāng)量比波動(dòng)較為明顯，鑒于此，所使用的DTF模型通過(guò)判斷當(dāng)?shù)乜扇嘉锏幕旌戏謹(jǐn)?shù)，進(jìn)而自動(dòng)調(diào)整層流火焰的燃燒速度和厚度.最后，利用動(dòng)態(tài)k方程亞格子模型，結(jié)合DTF模型和甲烷多步化學(xué)反應(yīng)機(jī)理，對(duì)悉尼部分預(yù)混燃燒裝置下FJ200-5GP-Lr75-103接近吹熄極限的工況進(jìn)行了模擬驗(yàn)證，下文中此工況根據(jù)其速度大小簡(jiǎn)記為U103，證明了DTF模型在計(jì)算部分預(yù)混火焰方面具有足夠的適應(yīng)性和準(zhǔn)確性.

1 湍流燃燒模型

1.1 湍流模型

對(duì)三維納維-斯托克斯方程進(jìn)行Favre濾波得到的大渦模擬控制方程參考文獻(xiàn)[21]，由濾波帶來(lái)的亞網(wǎng)格應(yīng)力項(xiàng)及亞網(wǎng)格標(biāo)量通量項(xiàng)需要模型進(jìn)行處理，其中本文計(jì)算中采用文獻(xiàn)[22]提到的動(dòng)態(tài)k方程模型對(duì)亞格子應(yīng)力項(xiàng)加以封閉：

(1)

1.2 動(dòng)態(tài)增厚火焰燃燒模型

文獻(xiàn)[13]考察了層流火焰的燃燒速度和火焰面厚度，發(fā)現(xiàn)其均勻與組分的擴(kuò)散系數(shù)和化學(xué)反應(yīng)有某種關(guān)聯(lián)：

(2)

(3)

文獻(xiàn)[14]認(rèn)為增厚火焰模型降低了火焰對(duì)流場(chǎng)的敏感性，于是引入褶皺因子E的概念.針對(duì)部分預(yù)混火焰的特性，首先基于GRII-Mech3.0機(jī)理得到了可燃界限內(nèi)甲烷常溫常壓下，層流燃燒速度sL和層流火焰面厚度δL隨當(dāng)量比φ的變化情況，如圖1所示.隨后，通過(guò)MATLAB提出二者隨當(dāng)量比變化規(guī)律的擬合函數(shù)，并植入DTF模型中.在所使用的DTF模型中，火焰面的位置通過(guò)多步化學(xué)反應(yīng)的放熱量dq來(lái)捕捉，使得火焰的增厚系數(shù)只在燃燒火焰區(qū)內(nèi)有效，修正后的火焰增厚因子Fdyn和褶皺因子Edyn分別由下式表示：

(4)

式中：ndtf為預(yù)設(shè)的增厚因子；Δe為網(wǎng)格尺度；a1、a2為調(diào)整流場(chǎng)中動(dòng)態(tài)增厚曲線形態(tài)的縮放因子，此處a1=10,a2=25.褶皺因子可由下式表示：

(5)

(6)

圖1 層流火焰的燃燒特性Fig.1 Combustion characteristics of laminar flame

2 算例設(shè)置

2.1 悉尼燃燒器

為研究入口處燃料的混合情況對(duì)火焰的影響，文獻(xiàn)[23]精心設(shè)計(jì)了一套燃燒器，結(jié)構(gòu)如圖2(a)所示，其中：Lr為回縮距離.燃燒器的入流管道是由兩個(gè)壁厚為0.25 mm的同心管組成，內(nèi)管通入甲烷氣體，內(nèi)管直徑為4 mm，外管通入空氣，外管直徑為7.5 mm.混合氣體通過(guò)外圍導(dǎo)管中的高溫燃燒廢氣點(diǎn)燃，廢氣導(dǎo)管的直徑為18 mm，整個(gè)裝置放置于25 mm × 25 mm的伴流風(fēng)洞中.燃燒器具有凹入的內(nèi)管，該內(nèi)管可以從外管出口平面向上游回縮.當(dāng)Lr=0時(shí)，內(nèi)管與外管平齊(內(nèi)管位于外管出口平面)，此時(shí)混合氣燃燒模式為完全非預(yù)混，即擴(kuò)散燃燒模式；當(dāng)Lr≥300 mm(40倍的外管直徑)時(shí)，此時(shí)燃?xì)馀c空氣充分混合，燃燒模式為完全預(yù)混燃燒模式.當(dāng)Lr處于兩個(gè)極限距離之間時(shí)，燃燒模式為部分預(yù)混燃燒模式.為驗(yàn)證DTF模型對(duì)部分預(yù)混火焰的預(yù)測(cè)能力，計(jì)算內(nèi)管回縮距離選取 75 mm，對(duì)應(yīng)混合氣體處于部分預(yù)混狀態(tài)，速度入流參數(shù)如表1所示，其中：U為內(nèi)外管混合段的混合氣平均速度；UA為外管入流段空氣的平均速度；UF為內(nèi)管入流段燃料的平均速度.為了點(diǎn)燃射流火焰，在射流火焰的外圍制造了一個(gè)環(huán)形的值班火焰區(qū)域，其中注入的是由C2H2、H2、CO2、9N2和空氣這5種組分按照甲烷-空氣火焰在化學(xué)恰當(dāng)比時(shí)的元素配比來(lái)配置的混合氣.計(jì)算過(guò)程中值班火焰入口取為實(shí)驗(yàn)測(cè)量結(jié)果，具體信息如表2所示.其中：w(N2)、w(CO2)、w(H2O)、w(O2)和w(CO)分別為對(duì)應(yīng)組分的質(zhì)量分?jǐn)?shù)；Up為值班火焰環(huán)空內(nèi)混合氣的平均速度；T為值班火焰峰值溫度.

圖2 悉尼部分預(yù)混燃燒器及計(jì)算網(wǎng)格示意圖[23]Fig.2 Schematic diagram of Sydney partially premixed burner and calculation grid[23]

表1 U103工況下入流速度信息Tab.1 Inflow speed information under U103 condition

表2 U103工況下值班火焰尾氣中各組分質(zhì)量分?jǐn)?shù)、流速和溫度信息

2.2 計(jì)算設(shè)置

本文計(jì)算基于開(kāi)源的計(jì)算流體力學(xué)軟件OpenFOAM，其向來(lái)以其代碼的高效率和優(yōu)良的可拓展性著稱(chēng).在原有的rhoReactingFoam 求解器基礎(chǔ)上，本課題組植入了DTF模型.計(jì)算采用一套分塊的包含3.82×106個(gè)單元的結(jié)構(gòu)網(wǎng)格，如圖2(b)所示.其中，徑向、周向、軸向分別部署了76×80×628個(gè)網(wǎng)格單元，火焰區(qū)域內(nèi)的網(wǎng)格最大寬度為0.2 mm.

化學(xué)反應(yīng)機(jī)理對(duì)模擬的精度有著明顯的影響，文獻(xiàn)[24]考量了5種總包機(jī)理與文獻(xiàn)[25]提出的骨架機(jī)理在貧燃部分預(yù)混旋流火焰中的表現(xiàn)，所獲得的結(jié)論為，總包機(jī)理可以大大降低計(jì)算成本，但同時(shí)也影響結(jié)果的質(zhì)量；而包含73步化學(xué)反應(yīng)的Lu骨架機(jī)理表現(xiàn)與實(shí)驗(yàn)值吻合的最好，但計(jì)算成本也遠(yuǎn)高于總包機(jī)理.本文中化學(xué)反應(yīng)機(jī)理選用文獻(xiàn)[26]提出的甲烷19步骨架反應(yīng)機(jī)理，該機(jī)理同樣是由GRI-Mech1.2機(jī)理簡(jiǎn)化而來(lái).計(jì)算中燃料和空氣入口施加強(qiáng)度為10%的隨機(jī)擾動(dòng)，時(shí)間離散為二階Backward格式，對(duì)流項(xiàng)采用二階有界線性格式，拉普拉斯項(xiàng)選用高斯線性正交.

3 結(jié)果與討論

流場(chǎng)在不同軸向位置處的流場(chǎng)信息如圖3所示.其中：r為軸向距離；x為徑向距離；d為鈍體直徑；Ux為徑向平均速度；Ux_rms為徑向均方根速度.從圖3中可以看出，可燃混合氣剛進(jìn)入燃燒室時(shí)，中心軸向處的速度較高，并向外逐漸衰減，混合氣與周?chē)榱鳉怏w接觸的剪切層中的速度梯度很大.隨著混合氣沿著下游方向不斷的發(fā)展，主流的平均速度不斷降低，而周?chē)黧w的速度卻在不斷增大.這是由于剪切層的作用，使得速度較高的可燃?xì)獠粩嗟貙⒅車(chē)目諝饩砣胫髁髦?，以降低剪切層中的速度梯度，而周?chē)黧w由此受益使其平均速度逐漸增大.觀察流場(chǎng)中的均方根脈動(dòng)，可以看到在整個(gè)流場(chǎng)中，均方根脈動(dòng)的峰值幾乎沒(méi)有發(fā)生變化，但脈動(dòng)的頻譜卻越來(lái)越寬，整個(gè)流場(chǎng)越來(lái)越趨于一致.

圖3 不同軸向位置處的速度場(chǎng)信息Fig.3 Velocity field information at different axial positions

為了辨別流場(chǎng)中的渦旋結(jié)構(gòu)與位置，利用無(wú)量綱化的Q準(zhǔn)則和渦量ω提取流場(chǎng)中的擬序結(jié)構(gòu)，如圖4所示.Q準(zhǔn)則云圖中可以看到在火焰的根部r/d約為0.5的位置處達(dá)到峰值，然后向下游緩緩的擴(kuò)張開(kāi)來(lái)，在x/d=5到x/d=20的區(qū)間內(nèi)渦旋較為密集.渦量圖顯示了流體微團(tuán)的旋轉(zhuǎn)規(guī)律.從圖4中可以看到，沿軸向的渦旋幾乎總是與反向渦旋成對(duì)出現(xiàn)，在x/d>20的下游位置渦量逐漸減弱.對(duì)于不同軸向位置，其渦旋峰值所在的位置正好與圖4流場(chǎng)中速度沿x方向的均方根波動(dòng)峰位相吻合，這進(jìn)一步表明由剪切層產(chǎn)生的渦旋加劇了流場(chǎng)的摻混和軸向速度的傳遞.

以二氧化碳的瞬時(shí)生成速率云圖RCO2及火焰指數(shù)ξ來(lái)分析火焰面的信息，兩種方案的結(jié)果如圖5所示.其中：虛線為RCO2=3 kg/(m3·s) 的等值線，用來(lái)勾勒火焰面的位置.由圖5可知，由CO2生成的速率云圖中不同時(shí)刻的火焰面的底部x/d<10的部分均沒(méi)有產(chǎn)生中斷，這意味著火焰底部附近的燃燒較為穩(wěn)定.這是由于外部的值班火焰(根部左、右兩支角狀火焰分支)持續(xù)作用的結(jié)果.由于過(guò)量的甲烷和空氣在導(dǎo)管內(nèi)來(lái)不及充分預(yù)混便進(jìn)入到燃燒室，火焰根部的混合分?jǐn)?shù)較高，只有在外層較為稀薄的地方被伴流火焰所點(diǎn)燃，所以根部附近的火焰面較窄，通過(guò)分析火焰指數(shù)云圖可以看到，兩種方案在此區(qū)域均為預(yù)混燃燒模式.然而在火焰的內(nèi)部區(qū)域，混合氣體此時(shí)由擴(kuò)散模式主導(dǎo).在火焰中部1020)火焰繼續(xù)向下游發(fā)展，孔洞出現(xiàn)的越來(lái)越多，其范圍也逐漸增大，局部熄火現(xiàn)象進(jìn)一步加劇，這一區(qū)域的火焰由擴(kuò)散模式所主導(dǎo).

圖4 燃燒室內(nèi)無(wú)量綱化的Q準(zhǔn)則及渦量ω圖Fig.4 Dimensionless Q criterion and vorticity ω in combustion chamber

圖5 CO2瞬時(shí)生成速率云圖和火焰指數(shù)云圖Fig.5 Cloud maps of CO2 instantaneous generation rate and flame index

圖6 不同軸向位置處混合分?jǐn)?shù)均值及均方根波動(dòng)值對(duì)比Fig.6 Comparison of mean mixture fraction and root mean square fluctuation at different axial positions

可燃?xì)獾幕旌蠣顟B(tài)決定了火焰的溫度及燃燒產(chǎn)物，故而計(jì)算中必須對(duì)混合分?jǐn)?shù)做出精準(zhǔn)的預(yù)測(cè).圖6顯示了不同軸向位置處混合物的質(zhì)量分?jǐn)?shù)對(duì)比，其中：Zmean為混合分?jǐn)?shù)均值；ZRMS為混合分?jǐn)?shù)均方根波動(dòng)值.由圖6可以看到，無(wú)論是火焰的根部亦或是下游區(qū)域，采用動(dòng)態(tài)k方程湍流模型及結(jié)構(gòu)化的382萬(wàn)個(gè)網(wǎng)格對(duì)可燃?xì)獾幕旌戏謹(jǐn)?shù)的計(jì)算結(jié)果與實(shí)驗(yàn)值吻合良好.在整個(gè)流場(chǎng)中，混合氣的質(zhì)量分?jǐn)?shù)呈現(xiàn)出了較好的對(duì)稱(chēng)性，在r/d=0的中心線上，混合氣的質(zhì)量分?jǐn)?shù)達(dá)到最大，然后向兩邊逐漸遞減.然而混合分?jǐn)?shù)的均方根脈動(dòng)峰值卻沒(méi)有出現(xiàn)在中心線上，而是在中心線附近并沿中心線對(duì)稱(chēng)，呈現(xiàn)出兩個(gè)峰狀.隨著高度的增加，混合氣的質(zhì)量分?jǐn)?shù)峰值逐漸降低，火焰面隨之也變得越來(lái)越寬.這是由于邊緣較為稀薄的可燃?xì)獠粩嗟谋幌?，混合氣從濃度高的中心區(qū)域不斷向兩邊擴(kuò)散，整個(gè)混合氣的濃度逐漸被稀釋的結(jié)果.

圖7為計(jì)算結(jié)果與實(shí)驗(yàn)值的溫度對(duì)比，其中：Tmean為溫度均值；TRMS為溫度均方根波動(dòng)值.整體上看，改進(jìn)后的DTF模型計(jì)算結(jié)果較為良好的展現(xiàn)出了火焰中的溫度分布及波動(dòng)趨勢(shì).當(dāng)x/d<10時(shí)，平均溫度峰值出現(xiàn)在r/d=1附近，而均方根波動(dòng)的峰值出現(xiàn)在r/d=0.5附近；隨著下游熄火位點(diǎn)的增多，溫度場(chǎng)的均值不斷下降，而其波動(dòng)也逐漸蔓延到外側(cè)區(qū)域.

圖7 不同軸向位置處的平均溫度及均方根波動(dòng)值對(duì)比Fig.7 Comparison of mean value and root mean square fluctuation of temperature at different axial positions

文獻(xiàn)[2]在層流部分預(yù)混射流火焰的直接數(shù)值模擬中指出，部分預(yù)混燃燒火焰中包含兩個(gè)反應(yīng)層，內(nèi)層中的甲烷和氧氣不斷消耗產(chǎn)生CO、H2、OH等中間產(chǎn)物，然后在外層的氧化層中被反應(yīng)為CO2和H2O.這里選擇CO作為中間產(chǎn)物的代表，并與實(shí)驗(yàn)值進(jìn)行對(duì)比，結(jié)果如圖8所示，其中:w(CO)mean為一氧化碳質(zhì)量分?jǐn)?shù)均值;w(CO)RMS為一氧化碳質(zhì)量分?jǐn)?shù)均方根波動(dòng)值.對(duì)于中上游中間產(chǎn)物的計(jì)算結(jié)果，DTF模型結(jié)果較好的展現(xiàn)了中間產(chǎn)物CO的含量及波動(dòng)特征.同時(shí)可以看到在火焰的下游區(qū)域，計(jì)算結(jié)果較實(shí)驗(yàn)值出現(xiàn)了一定的偏差，此區(qū)域的CO質(zhì)量分?jǐn)?shù)均方根波動(dòng)也變得較為劇烈，這可能是由于采用簡(jiǎn)化機(jī)理無(wú)法較為全面地描述CO的反應(yīng)過(guò)程.從圖8左圖可以看到，不同平面高度上CO的峰值均出現(xiàn)在r/d=0.75附近，然后向兩邊緩慢擴(kuò)散；到了下游區(qū)域，曲線逐漸變得平滑.與圖7中溫度曲線的峰值位置進(jìn)行對(duì)比，可以發(fā)現(xiàn)放熱峰值出現(xiàn)在r/d=0.9附近，稍遠(yuǎn)于CO峰值面，這表明Suresh的假設(shè)是比較合理的.由于避免了環(huán)境因素的影響，計(jì)算值比實(shí)驗(yàn)值呈現(xiàn)出了更好的對(duì)稱(chēng)性.

將混合氣的溫度隨質(zhì)量分?jǐn)?shù)分布制成散點(diǎn)圖如圖9所示，其中：T為溫度；Z為混合分?jǐn)?shù)；藍(lán)色虛線為化學(xué)恰當(dāng)比所在線，將火焰分割為左、右兩支.左半支的混合分?jǐn)?shù)較低，而對(duì)應(yīng)的溫度較高，故其為高溫燃燒廢氣與伴流空氣的混合支，右側(cè)的混合分?jǐn)?shù)較高但對(duì)應(yīng)的溫度比較低，對(duì)應(yīng)可燃?xì)馀c高溫燃燒廢氣的混合支.由于入流速度過(guò)高(接近吹熄極限)，導(dǎo)致火焰根部x/d=1附近大量的可燃?xì)怏w來(lái)不及引燃便被帶到下游.進(jìn)一步在x/d=5的位置觀察火焰燃燒狀態(tài)，可以發(fā)現(xiàn)此時(shí)混合氣的燃燒狀態(tài)較為理想.從x/d>10，左、右兩個(gè)混合分支間開(kāi)始出現(xiàn)熄火點(diǎn)，越到下游熄火點(diǎn)出現(xiàn)的越為密集.由圖3中的火焰指數(shù)云圖可以看出，在流場(chǎng)的中下游，擴(kuò)散燃燒模式逐漸占據(jù)主導(dǎo)地位.對(duì)比實(shí)驗(yàn)和計(jì)算的散點(diǎn)圖可以發(fā)現(xiàn)，采用改進(jìn)后的動(dòng)態(tài)層流火焰燃燒速度及層流火焰厚度的DTF模型整體非常良好的捕捉到了這些變化趨勢(shì).散點(diǎn)圖中火焰的根部x/d=1處，無(wú)論是實(shí)驗(yàn)亦或是計(jì)算模擬，散點(diǎn)數(shù)據(jù)都比較少；且大部分散點(diǎn)都集中在化學(xué)當(dāng)量比分割線附近.當(dāng)混合分?jǐn)?shù)Z超過(guò)0.1，300 K附近出現(xiàn)一條水平直線.這是由于混合氣剛進(jìn)入燃燒室時(shí)，底部主流區(qū)域的可燃?xì)饣旌戏謹(jǐn)?shù)過(guò)高且流速過(guò)快，來(lái)不及與周?chē)娜紵龔U氣進(jìn)行熱量交換，這部分散點(diǎn)全部集中在300 K的帶狀分布附近.由于這部分的散點(diǎn)比例過(guò)大，導(dǎo)致整個(gè)散點(diǎn)圖中的點(diǎn)顯得比較稀少.

圖8 不同軸向位置處混合物的CO質(zhì)量分?jǐn)?shù)平均值及均方根波動(dòng)值對(duì)比Fig.8 Comparison of mean mass fraction and root mean square fluctuation of CO at different axial positions

圖9 不同軸向位置處混合物的溫度隨混合分?jǐn)?shù)分布散點(diǎn)圖Fig.9 Scatter plot of temperature with mixture fraction distribution at different axial positions

4 結(jié)論

為了應(yīng)用動(dòng)態(tài)增厚火焰模型計(jì)算當(dāng)量比不均勻分布的湍流部分預(yù)混火焰，需要考慮到當(dāng)量比變化對(duì)火焰褶皺函數(shù)中出現(xiàn)的層流火焰速度、厚度的影響.本文通過(guò)計(jì)算一系列的層流火焰構(gòu)造了兩個(gè)擬合函數(shù)，其可以在DTF模型中自適應(yīng)調(diào)整火焰參數(shù)，很好地解決了上述問(wèn)題.結(jié)合動(dòng)態(tài)k方程亞格子湍流模型，對(duì)接近吹熄極限的悉尼部分預(yù)混火焰進(jìn)行了大渦模擬研究.本文研究所獲得的結(jié)論如下.

(1) 通過(guò)流場(chǎng)中混合分?jǐn)?shù)均值及均方根脈動(dòng)對(duì)比，顯示了動(dòng)態(tài)k方程亞格子湍流模型可以很好地預(yù)測(cè)非均勻入流的射流流動(dòng)問(wèn)題.

(2) 觀察CO2生成速率云圖和OH基渲染的三維瞬時(shí)火焰面，不難發(fā)現(xiàn)采用改進(jìn)后的DTF模型對(duì)此類(lèi)帶有熄火現(xiàn)象的部分預(yù)混火焰具有較好的預(yù)測(cè)能力，溫度-混合分?jǐn)?shù)散點(diǎn)圖分布也很好的說(shuō)明了這一點(diǎn).

(3) 本文使用的甲烷多步簡(jiǎn)化機(jī)理對(duì)某些組分如CO的質(zhì)量分?jǐn)?shù)預(yù)測(cè)較實(shí)驗(yàn)值偏低，若采用更為詳細(xì)的化學(xué)反應(yīng)機(jī)理或結(jié)合化學(xué)制表方法，這種現(xiàn)象可能會(huì)有一些改善.