曾海翔，王 平，余 倩，李偉超，陳明敏

(1.江蘇大學(xué)能源與動(dòng)力工程學(xué)院，鎮(zhèn)江 212013；2.江蘇大學(xué)能源研究院，鎮(zhèn)江 212013；3.上海電氣燃?xì)廨啓C(jī)有限公司，上海 201199)

分層燃燒是一種當(dāng)量比分布不均勻的燃燒方式，其優(yōu)點(diǎn)之一是可擴(kuò)大可燃混合物的燃燒極限[1].這一概念已經(jīng)在實(shí)際燃燒器中實(shí)現(xiàn)[2].按照火焰的傳播方式可將分層火焰分為“焰前支持”(front-supported)火焰和“焰后支持”(back-supported)火焰.“焰前支持”火焰，即火焰的傳播方向由稀預(yù)混分支向濃預(yù)混分支蔓延；“焰后支持”火焰的傳播規(guī)律則恰好相反，即火焰由濃預(yù)混分支向稀預(yù)混分支發(fā)展.而按照當(dāng)量比變化進(jìn)行分類(lèi)，當(dāng)量比介于富燃極限和貧燃極限之間時(shí)，僅能發(fā)生預(yù)混燃燒；一旦超過(guò)可燃極限，火焰中會(huì)同時(shí)存在預(yù)混和擴(kuò)散兩種燃燒方式，此時(shí)的分層燃燒又稱(chēng)為部分預(yù)混燃燒[3].

在燃料/空氣混合層中傳播的火焰，當(dāng)燃料濃度梯度不均勻分布時(shí)，便會(huì)出現(xiàn)三岔火焰結(jié)構(gòu)[4].三岔火焰作為一種典型的分層燃燒火焰，兼具稀預(yù)混燃燒、濃預(yù)混燃燒和擴(kuò)散燃燒3 種燃燒模式.濃預(yù)混分支剩余的燃料和稀預(yù)混分支剩余的氧化劑擴(kuò)散到中間，形成擴(kuò)散火焰，擴(kuò)散火焰分支中混合物接近化學(xué)平衡[5].三岔火焰廣泛存在于多種流場(chǎng)條件下，例如協(xié)同流、自由射流和對(duì)沖流等.不同流場(chǎng)和燃燒裝置下形成的三岔火焰具有兩種形式，對(duì)沖流下平行三岔火焰和協(xié)同流下分叉型三岔火焰.分叉型三岔火焰3個(gè)火焰分支匯聚于化學(xué)當(dāng)量比為1 的三岔點(diǎn)上.

研究三岔火焰對(duì)了解非穩(wěn)態(tài)火焰的燃燒特性如著火、發(fā)展以及熄滅機(jī)制有很大幫助.自從1965 年P(guān)hillips 發(fā)現(xiàn)三岔火焰現(xiàn)象以來(lái)[6]，國(guó)內(nèi)外學(xué)者對(duì)三岔火焰產(chǎn)生了濃厚的興趣.由于部分預(yù)混火焰具有多反應(yīng)區(qū)的特征，數(shù)值模擬中經(jīng)常采用的總包反應(yīng)機(jī)理難以描述其火焰結(jié)構(gòu)，因而大多數(shù)學(xué)者在計(jì)算過(guò)程中會(huì)選擇較為詳細(xì)的化學(xué)反應(yīng)機(jī)理，Qin 等[7]針對(duì)重力對(duì)甲烷/空氣部分預(yù)混層流火焰的影響開(kāi)展了研究，計(jì)算過(guò)程中采用81 步甲烷反應(yīng)機(jī)理；Yamamoto等[8]使用25 步甲烷化學(xué)反應(yīng)機(jī)理來(lái)研究甲烷/空氣共流射流下三岔火焰形態(tài)；Yuan 等[9]分析流場(chǎng)的壁面對(duì)甲烷/空氣三岔火焰?zhèn)鞑ヅc幾何結(jié)構(gòu)特性的影響，其中的數(shù)值計(jì)算部分采用了甲烷184 步化學(xué)反應(yīng)機(jī)理，對(duì)應(yīng)的計(jì)算精度有一定提升，但其計(jì)算量也大大增加.

近些年來(lái)國(guó)內(nèi)外同時(shí)針對(duì)甲烷發(fā)展了大量的化學(xué)反應(yīng)機(jī)理.Hughes 等[10]通過(guò)對(duì)甲烷反應(yīng)的動(dòng)力學(xué)特性的測(cè)量，提出一套包含37 種組分351 個(gè)正向反應(yīng)的甲烷詳細(xì)化學(xué)反應(yīng)機(jī)理，反應(yīng)的詳細(xì)參數(shù)以CHEMKIN 格式置于網(wǎng)上.Slavinskaya 等[11]基于敏感性分析的方法將烴類(lèi)化學(xué)動(dòng)力學(xué)數(shù)據(jù)庫(kù)中C0～C1氧化的一般動(dòng)力學(xué)方案進(jìn)行推廣，得到一套低壓甲烷燃燒的骨架反應(yīng)機(jī)理，并通過(guò)著火延遲和層流火焰速度對(duì)比對(duì)其進(jìn)行了評(píng)估.Novosselov 等[12]為預(yù)測(cè)工業(yè)燃?xì)廨啓C(jī)的污染排放特性開(kāi)發(fā)一套八步反應(yīng)的甲烷總包反應(yīng)機(jī)理，相較之前提出的五步反應(yīng)機(jī)理進(jìn)一步增加了壓力梯度的適應(yīng)范圍.劉合等[13]在研究甲烷/空氣預(yù)混超聲速燃燒時(shí)提出一套包含14 種組分18 個(gè)基元反應(yīng)的簡(jiǎn)化機(jī)理，該機(jī)理在預(yù)測(cè)點(diǎn)火延遲方面具有較高精度且在高、低壓力范圍內(nèi)均表現(xiàn)良好.胡賢忠等[14]基于GRI 3.0 詳細(xì)機(jī)理構(gòu)建了甲烷在O2/CO2氣氛下燃燒的框架機(jī)理和總包機(jī)理，為甲烷在O2/CO2氣氛下的燃燒過(guò)程的數(shù)值模擬提供高效的簡(jiǎn)化模型.

為了降低計(jì)算成本的同時(shí)有效提高計(jì)算精度，在過(guò)去幾十年中已經(jīng)開(kāi)發(fā)了多種簡(jiǎn)化動(dòng)力學(xué)方法，主要有計(jì)算奇異(CSP)方法[15]、自適應(yīng)建表(ISAT)方法[16]和固有低維流形(ILDM)方法[17].其中ILDM 方法是基于對(duì)化學(xué)系統(tǒng)非線性響應(yīng)動(dòng)態(tài)行為的直接數(shù)值分析，通過(guò)對(duì)化學(xué)反應(yīng)項(xiàng)的Jacobian 矩陣進(jìn)行特征值分析，把反應(yīng)空間分解為快反應(yīng)子空間和慢反應(yīng)子空間.在慢反應(yīng)子空間中構(gòu)造流形，從而實(shí)現(xiàn)用部分參數(shù)來(lái)表征整個(gè)化學(xué)系統(tǒng)的組成.ILDM 方法在簡(jiǎn)化機(jī)理的過(guò)程中假設(shè)化學(xué)反應(yīng)占主導(dǎo)地位，忽略了低溫區(qū)對(duì)流擴(kuò)散對(duì)化學(xué)反應(yīng)造成的影響.因此，只在高溫區(qū)有很好的模擬精度.增加子空間維度是提高低溫區(qū)的計(jì)算精度的方法之一，但同時(shí)極大地增加了計(jì)算量.由van Oijen 等[18]提出的火焰面生成流形(FGM)和Gicquel 等[19]提出的火焰固有流形延伸(PFI)方法，正是基于ILDM 方法并考慮了擴(kuò)散效應(yīng)，得到了廣泛的應(yīng)用.這兩種方法成功應(yīng)用于預(yù)混和擴(kuò)散燃燒的數(shù)值模擬后，越來(lái)越多的學(xué)者嘗試將其用于部分預(yù)混火焰.Fiorina 等[20]將二維PFI 表格應(yīng)用于部分預(yù)混火焰，發(fā)現(xiàn)當(dāng)混合物當(dāng)量比位于可燃界限內(nèi)時(shí)，PFI 模型獲得優(yōu)良結(jié)果；反之，出現(xiàn)FPI 結(jié)果與詳細(xì)化學(xué)反應(yīng)機(jī)理結(jié)果的偏離.Wu 等[5]在計(jì)算二維三岔火焰時(shí)，發(fā)現(xiàn)FPI 制表模型預(yù)測(cè)的CO 質(zhì)量分?jǐn)?shù)遠(yuǎn)大于詳細(xì)反應(yīng)機(jī)理下的結(jié)果.宋帥[21]通過(guò)引入火焰指數(shù)將FGM 預(yù)混表和FGM 擴(kuò)散表應(yīng)用于部分預(yù)混火焰，得到了較好的結(jié)果，但計(jì)算量有所增加.Steinhilber 和Maas[22]發(fā)展的“Reaction-Diffusion Mainfold”方法則較好地克服了計(jì)算量上升這個(gè)問(wèn)題.王蔡軍等[23]利用二維REDIM 表結(jié)合假定概率密度函數(shù)(PFDF)對(duì)Darmstadt 大學(xué)分層燃燒裝置進(jìn)行了計(jì)算，與實(shí)驗(yàn)結(jié)果吻合很好，但計(jì)算所設(shè)置的混合物當(dāng)量比均在可燃界限內(nèi).

為了研究不同當(dāng)量比條件下三岔火焰及分層火焰的燃燒特性，同時(shí)驗(yàn)證二維REDIM 表格對(duì)部分超出可燃極限燃燒狀態(tài)的描述能力，本文考慮了4 種不同當(dāng)量比布局的層流火焰，并采用詳細(xì)甲烷化學(xué)反應(yīng)對(duì)其進(jìn)行計(jì)算.計(jì)算穩(wěn)定后，對(duì)比分析了這4 種火焰的火焰結(jié)構(gòu)及火焰面上的燃燒狀態(tài)，并將計(jì)算域內(nèi)每個(gè)網(wǎng)格點(diǎn)的N2和CO2值投影到二維REDIM 表格中，觀察REDIM 表格是否完全容納投影的結(jié)果.

1 反應(yīng)-擴(kuò)散流形方法

REDIM 方法在生成低維流形時(shí)同樣考慮了化學(xué)反應(yīng)對(duì)分子輸運(yùn)過(guò)程的影響，并引入不變流形理論.最終得到的低維流形方程[24]如式(1)所示：

式中：ρ為密度；χ＝grad(θ)為θ的空間梯度矢量為Hessian 矩陣；d 為分子擴(kuò)散系數(shù)矩陣，符號(hào)“?”為具有三階張量的矢量運(yùn)算.REDIM 表格的生成方法可以分為以下幾個(gè)步驟[22]：

(1) 指定方程(1)初始條件和邊界條件；

(2) 將梯度grad(θ)指定為θ的函數(shù)；

(3) 對(duì)方程(1)進(jìn)行時(shí)間積分，直到達(dá)到穩(wěn)定狀態(tài)；

本文中采用Dirichlet 邊界條件，通過(guò)INSFLA[25]程序使用詳細(xì)化學(xué)反應(yīng)機(jī)理求解不同一維層流小火焰單元，繼而生成REDIM 表格的計(jì)算邊界和初始流形.這些小火焰單元包括Air-CH4冷態(tài)的純混合流、Air-pilot 對(duì)沖火焰和fuel-pilot 對(duì)沖火焰等.如圖1(a)所示，紅色線條表示fuel-pilot 對(duì)沖火焰，綠色線條表示pilot-Air 對(duì)沖火焰，底部藍(lán)色線條表示Air-CH4對(duì)沖流中無(wú)燃燒的純混合過(guò)程，這3 條線構(gòu)成二維REDIM 表格的邊界.邊界內(nèi)部的虛線表示不同拉伸率下的fuel-Air 對(duì)沖火焰，用于指定REDIM 表的初始流型.曲線從上到下表示拉伸率從小到大，當(dāng)拉伸率過(guò)大時(shí)，火焰會(huì)呈現(xiàn)熄火狀態(tài)，初始流型只要在計(jì)算邊界內(nèi)均可用于REDIM 表的計(jì)算.圖1(b)表示的是由N2-CO2-OH 所表示的三維骨架圖，圖中綠色的虛線表示瞬間熄火的小火焰單元.由此可見(jiàn)二維REDIM 表格能夠描述熄火狀態(tài)，比較適用于計(jì)算部分預(yù)混火焰.將所計(jì)算的小火焰單元 Profilecounterdiff 文件中最后一個(gè)DATA 數(shù)據(jù)導(dǎo)入REDIM表格生成程序中，指定方程的初始條件和邊界條件，得到如圖1(c)所示的初始流形.

圖1 二維REDIM表格生成過(guò)程示意Fig.1 Schematic of generation of 2D REDIM table

通過(guò)執(zhí)行步驟(2)、(3)和(4)，最終得到以N2和CO2質(zhì)量分?jǐn)?shù)為自變量的二維REDIM 查詢(xún)表，如圖1(d)所示表征反應(yīng)進(jìn)程，wN2表征混合過(guò)程，表示CO2生成速率.其他熱力學(xué)變量信息可以通過(guò)二維REDIM 表查詢(xún)得到.無(wú)燃燒的混合過(guò)程中CO2的質(zhì)量分?jǐn)?shù)為零；N2的質(zhì)量分?jǐn)?shù)可以反映混合氣的混合狀態(tài)，純空氣流對(duì)應(yīng)N2質(zhì)量分?jǐn)?shù)取最大，而當(dāng)量比達(dá)到最大值時(shí)，此時(shí)的 N2質(zhì)量分?jǐn)?shù)最小.根據(jù)其物理意義，這兩個(gè)組分的質(zhì)量分?jǐn)?shù)主要集中在圖1(d)所示的三角形范圍內(nèi).

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

2.1 邊界條件設(shè)定

本文一共計(jì)算了4 種狀態(tài)下的火焰，如圖2 所示，其中狀態(tài)A 和B 兩種分層火焰，混合氣當(dāng)量比分別超出可燃上限和可燃下限.由當(dāng)量比等于1 的甲烷/空氣混合物燃燒后的高溫廢氣分別點(diǎn)燃當(dāng)量比在富燃區(qū)間(狀態(tài)A)和貧燃區(qū)間(狀態(tài)B)線性變化的新鮮混合氣體，廢氣的入口速度為0.5 m/s，新鮮混合氣入口速度為1 m/s.其中，狀態(tài)A 為“焰前支持”火焰，預(yù)混氣體位于5 mm＜y≤15 mm 處，混合氣當(dāng)量比從1 到2 線性增加；狀態(tài)B 為“焰后支持”火焰，預(yù)混氣體位于5 mm＜y≤15 mm 處，混合氣當(dāng)量比從1～0 線性減小.狀態(tài)C 和D 是兩種特殊的分層火焰，具有典型的三岔火焰結(jié)構(gòu)，混合氣入口溫度為300 K，計(jì)算域的內(nèi)部為高溫燃燒廢氣.狀態(tài)C 的入口混合氣的速度為0.5 m/s，當(dāng)量比從0～2 線性增加，y＝7.5 mm 處當(dāng)量比為1；狀態(tài)D 的入口混合氣的速度為0.4 m/s，當(dāng)量比從0.5～1.5 線性增加，y＝7.5 mm 處當(dāng)量比為1.除了入口條件的設(shè)置不同，這4 個(gè)狀態(tài)的火焰具有相同的計(jì)算條件.其中壁面采用滑移邊界條件；壓力出口設(shè)置為常壓，p＝101 325 Pa.所有計(jì)算采用15 mm×25 mm 的二維矩形計(jì)算域，在x 和y 方向分別均勻設(shè)置600 和1 000 個(gè)網(wǎng)格單元，網(wǎng)格尺寸為25 μm，足夠求解火焰結(jié)構(gòu).

圖2 4個(gè)層流火焰的計(jì)算域及邊界條件示意Fig.2 Schematic of calculation domain and boundary conditions of four laminar flames

2.2 計(jì)算方法

本文計(jì)算基于 OpenFOAM 開(kāi)源軟件，選用rhoReactingBuoyantFoam 作為求解器.計(jì)算過(guò)程中采用甲烷-空氣詳細(xì)反應(yīng)機(jī)理[26]，對(duì)控制方程中的時(shí)間項(xiàng)采用了默認(rèn)的一階隱式Euler 差分格式，對(duì)流項(xiàng)采用限制差分型TVD 格式，擴(kuò)散項(xiàng)采用二階線性修正高斯差分格式.在OpenFOAM 中可以通過(guò)庫(kù)朗數(shù)自動(dòng)調(diào)整時(shí)間步長(zhǎng)，這里取為0.4，計(jì)算過(guò)程中迭代誤差統(tǒng)一設(shè)置為1×10-9.

3 結(jié)果分析

3.1 火焰結(jié)構(gòu)分析

為了區(qū)分燃燒過(guò)程中的預(yù)混和擴(kuò)散區(qū)域，以進(jìn)一步判斷當(dāng)?shù)氐娜紵龣C(jī)制，Yamashita 等[27]提出了火焰指數(shù)的概念，火焰指數(shù)定義為

圖3 狀態(tài)A溫度、當(dāng)量比、CO2生成速率以及火焰指數(shù)示意Fig.3 Skematic of temperature，equivalence ratio，formation rate of CO2，and flame index in state A

圖4 狀態(tài)B溫度、當(dāng)量比、CO2生成速率以及火焰指數(shù)示意Fig.4 Skematic of temperature，equivalence ratio，formation rate of CO2，and flame index in state B

圖3 和圖4 分別展示了兩種分層火焰下的溫度(T)、CO2生成速率(rCO2)、當(dāng)量比(φ)和火焰指數(shù)云圖.對(duì)比圖3、圖4 發(fā)現(xiàn)，兩種分層火焰狀態(tài)只有一個(gè)火焰分支，且火焰面的位置基本一致.通過(guò)CO2的生成來(lái)反映火焰面的位置，由火焰面處的火焰指數(shù)發(fā)現(xiàn)：狀態(tài)A 和狀態(tài)B 火焰面處火焰指數(shù)都接近1，因此這兩種燃燒模式下的火焰均為預(yù)混火焰主導(dǎo).雖然入口反應(yīng)物濃度范圍分別為[0，1]和[1，2]，均部分超出可燃極限，但并未產(chǎn)生擴(kuò)散火焰分支.由當(dāng)量比分布云圖知，狀態(tài)A 和狀態(tài)B 中計(jì)算域內(nèi)混合物當(dāng)量比呈現(xiàn)分層分布現(xiàn)象，符合入口混合物的當(dāng)量比設(shè)置.

圖5 和圖6 中OH 質(zhì)量分?jǐn)?shù)分布云圖清晰地展現(xiàn)了狀態(tài)C 和狀態(tài)D 中的3 個(gè)火焰分支.結(jié)合圖5、6 中的CO2生成速率rCO2云圖可知，在這兩個(gè)狀態(tài)下兩側(cè)的火焰分支均為預(yù)混火焰，中間的火焰分支為擴(kuò)散火焰.濃預(yù)混火焰分支位于擴(kuò)散火焰分支的上方，稀預(yù)混火焰分支位于擴(kuò)散火焰分支的下方.同時(shí)發(fā)現(xiàn)，在三岔點(diǎn)上和擴(kuò)散火焰分支處OH 質(zhì)量分?jǐn)?shù)遠(yuǎn)高于兩個(gè)預(yù)混火焰分支，稀預(yù)混火焰分支處OH 質(zhì)量分?jǐn)?shù)稍高于濃預(yù)混火焰分支.狀態(tài)C 混合物入口濃度梯度大于狀態(tài)D，由CO2生成速率圖可知：狀態(tài)C 三岔火焰擴(kuò)散火焰分支相較于狀態(tài)D 更加清晰，與文獻(xiàn)[28]中所提及濃度梯度增加，擴(kuò)散火焰分支逐漸增強(qiáng)這一結(jié)論相符.狀態(tài)C 混合氣入口速度為0.5 m/s，在計(jì)算中火焰能夠穩(wěn)定在某一位置；對(duì)狀態(tài)D 設(shè)置相同的入口速度時(shí)，發(fā)現(xiàn)火焰不停地向下游傳播，直至離開(kāi)計(jì)算域；將預(yù)混氣的入流速度設(shè)置為0.4 m/s時(shí)，火焰駐定在計(jì)算域中，但穩(wěn)定高度距離入口約8 mm 處，高于狀態(tài)C 的3.7 mm 舉升高度.由此認(rèn)為狀態(tài)C 三岔火焰的傳播速度大于狀態(tài)D.狀態(tài)C 和狀態(tài)D 的火焰?zhèn)鞑ニ俣染哂诰鶆蝾A(yù)混條件下的層流火焰的傳播速度(約為0.3 m/s)，考慮到狀態(tài)C 的預(yù)混氣體的濃度梯度高于狀態(tài)D，可以判斷在均勻濃度梯度分布的條件下，增大預(yù)混氣的濃度梯度有助于火焰的傳播，同時(shí)火焰的舉升高度也發(fā)生顯著的降低.

圖5 狀態(tài)C溫度、CO2生成速率、OH質(zhì)量分?jǐn)?shù)以及火焰指數(shù)示意Fig.5 Skematic of temperature，formation rate of CO2，OH mass fraction，and flame index in state C

圖6 狀態(tài)D溫度、CO2生成速率、OH質(zhì)量分?jǐn)?shù)以及火焰指數(shù)示意Fig.6 Skematic of temperature，formation rate of CO2，OH mass fraction，and flame index in state D

結(jié)合4種狀態(tài)下火焰指數(shù)分布圖發(fā)現(xiàn)，對(duì)于簡(jiǎn)單的預(yù)混和擴(kuò)散火焰，火焰指數(shù)能夠很好地用來(lái)判斷當(dāng)?shù)氐娜紵Ｊ?，如狀態(tài)A 和狀態(tài)B 兩種火焰的火焰面處均只存在一個(gè)火焰指數(shù)0.9，這表明此時(shí)狀態(tài)A和狀態(tài)B 火焰面上的燃料和氧化劑的位移夾角非常小(約為3°以?xún)?nèi)).而對(duì)于本文中的三岔火焰狀態(tài)C和狀態(tài)D，結(jié)合OH 基生成云圖，可以清晰地看見(jiàn)此時(shí)上下兩端的預(yù)混火焰分支和中間的擴(kuò)散火焰分支于y＝7.5 mm 附近的三岔點(diǎn)處相遇，對(duì)應(yīng)化學(xué)當(dāng)量比等于1.此時(shí)，該點(diǎn)處的燃燒模式既不屬于預(yù)混火焰，也不屬于擴(kuò)散火焰.然而在圖5 和圖6 火焰指數(shù)分布示意圖中，位于三岔點(diǎn)的火焰指數(shù)顯示為0.99.根據(jù)火焰指數(shù)定義，此處的燃燒模式應(yīng)該為預(yù)混燃燒主導(dǎo).這與該點(diǎn)的實(shí)際燃燒模式不符，故此火焰指數(shù)的表達(dá)式在三岔點(diǎn)處不盡合理，有必要進(jìn)一步完善.

3.2 計(jì)算域內(nèi)N2、CO2質(zhì)量分?jǐn)?shù)散點(diǎn)圖分析

本文生成的二維REDIM 表格選取的參考變量為N2和CO2的質(zhì)量分?jǐn)?shù)，在燃燒數(shù)值模擬中可以通過(guò)查找上述參考變量來(lái)確定對(duì)應(yīng)狀態(tài)下其他變量的信息，有效規(guī)避了詳細(xì)求解復(fù)雜的化學(xué)反應(yīng).在制表過(guò)程中，當(dāng)量比選取的范圍在0～10 之間，涵蓋了本文中4 種狀態(tài)所設(shè)置的當(dāng)量比分布.如圖7～圖10所示，藍(lán)色線條表示由二維REDIM 表格生成的計(jì)算邊界，假如REDIM 表格可以很好地反映火焰的詳細(xì)信息，則實(shí)際的火焰信息通過(guò)查表法轉(zhuǎn)換的散點(diǎn)應(yīng)該全部被容納在計(jì)算邊界內(nèi)，即圖中的藍(lán)色線條封閉的區(qū)域內(nèi).為了加以驗(yàn)證，將上述4 種火焰計(jì)算域內(nèi)各點(diǎn)的N2和CO2質(zhì)量分?jǐn)?shù)值投影到二維REDIM 表格中，觀察這些點(diǎn)的分布情況.從圖7 和圖8 中不難看出，狀態(tài)A 和狀態(tài)B 散點(diǎn)圖中N2質(zhì)量分?jǐn)?shù)呈現(xiàn)出沿化學(xué)當(dāng)量比下N2質(zhì)量分?jǐn)?shù)wN2＝0.725 這條豎線的對(duì)稱(chēng)分布，這正好符合狀態(tài)A 和狀態(tài)B 火焰入口處的梯度設(shè)置(狀態(tài)A 的當(dāng)量比范圍1～2，狀態(tài)B 的當(dāng)量比范圍0～1).圖9 中狀態(tài)C(當(dāng)量比0～2)散點(diǎn)圖中N2質(zhì)量分?jǐn)?shù)分布最廣，說(shuō)明此狀態(tài)下火焰內(nèi)部的當(dāng)量比分布最廣.狀態(tài)D(當(dāng)量比0.5～1.5)散點(diǎn)圖中N2質(zhì)量分?jǐn)?shù)分布較狀態(tài)C 稍窄，但這兩種狀態(tài)當(dāng)量比分布都顯示出火焰從貧燃區(qū)跨度到富燃區(qū)，驗(yàn)證了濃度梯度不均勻分布的分層燃燒技術(shù)可以有效擴(kuò)大可燃物的燃燒極限.

4 種狀態(tài)N2和CO2質(zhì)量分?jǐn)?shù)對(duì)應(yīng)點(diǎn)都落在二維REDIM 表格中，說(shuō)明4 種狀態(tài)下的計(jì)算域內(nèi)各網(wǎng)格點(diǎn)對(duì)應(yīng)的熱化學(xué)參數(shù)都能被此二維REDIM 表格求解，圖7～10 可以比較直觀地說(shuō)明二維REDIM 表格適用于部分預(yù)混火焰.后期將采用二維REDIM 表格計(jì)算本研究中四種狀態(tài)火焰，并與GRI3.0 機(jī)理的結(jié)果進(jìn)行對(duì)比，詳細(xì)分析二維REDIM 表格的適用性.

圖7 狀態(tài)A火焰散點(diǎn)圖分布Fig.7 Scatter plot in state A

圖8 狀態(tài)B火焰散點(diǎn)圖分布Fig.8 Scatter plot in state B

圖9 狀態(tài)C火焰散點(diǎn)圖分布Fig.9 Scatter plot in state C

圖10 狀態(tài)D火焰散點(diǎn)圖分布Fig.10 Scatter plot in state D

4 結(jié)論

本文通過(guò)4 種層流火焰的直接數(shù)值模擬計(jì)算，對(duì)火焰結(jié)構(gòu)進(jìn)行分析，并將計(jì)算域內(nèi)各點(diǎn)N2和CO2質(zhì)量分?jǐn)?shù)投影到二維REDIM 表格中檢驗(yàn)其適用性.研究表明：

(1) 雖然兩種分層燃燒狀態(tài)A 和狀態(tài)B 混合物入口處的當(dāng)量比范圍均超出可燃界限，但均只存在一個(gè)預(yù)混火焰分支；狀態(tài)C 和狀態(tài)D 則呈現(xiàn)出三岔火焰結(jié)構(gòu).

(2) 通過(guò)計(jì)算火焰指數(shù)，發(fā)現(xiàn)該火焰指數(shù)能很好地表征預(yù)混或擴(kuò)散火焰的火焰狀態(tài)，對(duì)于復(fù)雜的部分預(yù)混火焰，如三岔火焰，則存在一定的缺陷.

(3) 將計(jì)算域內(nèi)各網(wǎng)格點(diǎn)的熱力學(xué)狀態(tài)投影到二維REDIM 表格中，發(fā)現(xiàn)4 種狀態(tài)計(jì)算域內(nèi)各點(diǎn)N2和CO2質(zhì)量分?jǐn)?shù)散點(diǎn)都落在二維REDIM 表格內(nèi)，這表明本文生成的REDIM 表格可擴(kuò)展到部分預(yù)混火焰的研究中，為此類(lèi)問(wèn)題的研究提供一個(gè)可能的解決方案.