曾小軍，汪小憨，趙黛青?，蔣利橋，楊衛(wèi)斌

（1. 中國(guó)科學(xué)院廣州能源研究所，廣州 510640；2. 中國(guó)科學(xué)院可再生能源重點(diǎn)實(shí)驗(yàn)室，廣州 510640；3. 廣東省新能源和可再生能源研究開(kāi)發(fā)與應(yīng)用重點(diǎn)實(shí)驗(yàn)室，廣州 510640）

兩種噴管中心間距下微陣列火焰場(chǎng)燃燒特性數(shù)值研究*

曾小軍1，2，3，汪小憨1，2，3，趙黛青1，2，3?，蔣利橋1，2，3，楊衛(wèi)斌1，2，3

（1. 中國(guó)科學(xué)院廣州能源研究所，廣州 510640；2. 中國(guó)科學(xué)院可再生能源重點(diǎn)實(shí)驗(yàn)室，廣州 510640；3. 廣東省新能源和可再生能源研究開(kāi)發(fā)與應(yīng)用重點(diǎn)實(shí)驗(yàn)室，廣州 510640）

構(gòu)建了處于無(wú)限陣列中管徑為1.5 mm的3 × 3噴管微陣列火焰模型，采用甲烷簡(jiǎn)化機(jī)理，分別對(duì)噴管中心間距為2 mm和4 mm的微陣列火焰場(chǎng)燃燒特性進(jìn)行數(shù)值研究。結(jié)果表明，在相同輸入工況下，隨著管中心間距的減小，噴管火焰間的相互影響增強(qiáng)，燃燒反應(yīng)區(qū)域變長(zhǎng)，最高溫度增大；同時(shí)，計(jì)算區(qū)域出口CO排放增加、燃燒效率下降。噴管中心間距2 mm和4 mm兩種工況下，距離管口5 mm后的區(qū)域溫度都開(kāi)始均勻分布，為研究溫度均勻的微陣列火焰加熱器提供了設(shè)計(jì)依據(jù)。

微小火焰；陣列火焰；溫度均勻性；火焰加熱；甲烷

0　 引 言

火焰直接加熱是常用的加熱方式，對(duì)厘米級(jí)、毫米級(jí)加熱對(duì)象及特殊均溫要求的加熱體，傳統(tǒng)的常規(guī)大尺度火焰存在火焰長(zhǎng)、加熱對(duì)象受熱不均勻、NOx和碳煙排放高等缺點(diǎn)。如何將長(zhǎng)火焰改善成為溫度均勻的清潔火焰，并開(kāi)發(fā)溫度可控的清潔火焰加熱系統(tǒng)，對(duì)厘米級(jí)以下有均溫要求的加熱對(duì)象有著重要的實(shí)際意義。研究表明，微尺度火焰藍(lán)色透明、清潔，隨著尺度變小，其最高發(fā)熱率增大［1，3］、最高溫度絕對(duì)值降低、溫度梯度減小［2，4］，同時(shí)透明的微火焰也不易產(chǎn)生碳煙。對(duì)于小型和微小型加熱對(duì)象，微尺度火焰可以很好地克服常規(guī)擴(kuò)散火焰的缺點(diǎn)，形成清潔、溫度均勻的加熱場(chǎng)。單個(gè)微火焰釋放的熱量絕對(duì)值很小，很難滿足高熱負(fù)荷的加熱要求，而采用多組微噴管來(lái)實(shí)現(xiàn)微火焰的疊加，最終形成滿足高負(fù)荷要求的加熱系統(tǒng)，是實(shí)現(xiàn)微尺度火焰替代大尺度火焰用于均勻加熱的有效途徑。采用多組微噴管構(gòu)建微陣列火焰加熱系統(tǒng)時(shí)，噴管中心間距對(duì)微陣列火焰場(chǎng)燃燒特性有著關(guān)鍵的影響作用，是微噴管陣列火焰加熱器設(shè)計(jì)優(yōu)化的重要參數(shù)，因此需要進(jìn)行深入的研究。

微噴管火焰間的影響程度是決定微火焰疊加的主要因素，NAKAMURA等［5］探索了兩個(gè)微噴嘴擴(kuò)散火焰的相互影響和合并條件，LEE等［6］對(duì)不同形狀單排陣列擴(kuò)散火焰的火焰特征進(jìn)行了詳細(xì)研究，這些研究多集中于擴(kuò)散燃燒方式，而針對(duì)微尺度陣列預(yù)混火焰的研究報(bào)道較少。在微尺度陣列火焰燃燒場(chǎng)中，由于噴管間火焰的相互影響，其燃燒特性和自由空間內(nèi)單噴管火焰有著明顯的差異，故本文構(gòu)建了處于無(wú)限陣列中管徑為1.5 mm的3 × 3噴管方形陣列火焰模型，采用甲烷空氣簡(jiǎn)化機(jī)理，在相同輸入工況下，針對(duì)不同噴管中心間距對(duì)微陣列預(yù)混火焰的影響規(guī)律進(jìn)行研究，對(duì)燃燒場(chǎng)的溫度、速度以及組分濃度進(jìn)行分析，并對(duì)不同間距下噴管間火焰相互影響程度和控制因素進(jìn)行討論。

1　 物理及數(shù)學(xué)模型

本文構(gòu)建的微陣列火焰物理模型如圖1所示，采用9組管徑Ф1.5 mm的微噴管等間距正方形布置，噴管高度為5 mm，管壁面無(wú)限薄，噴管進(jìn)口通入常溫常壓的預(yù)混甲烷-空氣；計(jì)算區(qū)域高度為40 mm，區(qū)域側(cè)面均設(shè)置為symmetry面對(duì)稱(chēng)邊界條件， 用于模擬無(wú)限擴(kuò)展的微噴管陣列結(jié)構(gòu)。預(yù)混氣當(dāng)量比設(shè)定為1，流速為0.5 m/s，底部噴管外區(qū)域有同向0.02 m/s的伴流空氣，模擬實(shí)際噴管底部存在的自然卷吸條件；噴管中心間距S設(shè)定為2 mm和4 mm兩個(gè)條件，其中C1為中心噴管圓心點(diǎn)，C2為中心噴管與水平方向相鄰噴管間距的中點(diǎn)，C3為中心噴管與對(duì)角線方向相鄰噴管間距的中點(diǎn)；后續(xù)分析中所取的中心軸面為C1、C2、C3點(diǎn)所在的中心切面。

圖1　物理模型示意圖Fig. 1 Schematic of the physical model

計(jì)算控制方程由質(zhì)量、動(dòng)量、能量及組分守恒方程構(gòu)成；模型采用組分輸運(yùn)和層流有限速率模型進(jìn)行燃燒反應(yīng)模擬；化學(xué)反應(yīng)機(jī)理為CHEMKIN導(dǎo)入、包含16種組分和41步反應(yīng)的甲烷空氣簡(jiǎn)化機(jī)理［7］（16種組分包括：H2O、CO2、O2、CH4、CO、H、H2、OH、O、CH3、HCO、HO2、H2O2、CH2O、CH3O、N2）；離散格式采用二階中心差分；采用壓力速度的耦合算法進(jìn)行層流穩(wěn)態(tài)火焰的計(jì)算求解。模型采用結(jié)構(gòu)化網(wǎng)格進(jìn)行劃分，對(duì)主要燃燒區(qū)域的噴管出口進(jìn)行加密處理，總的計(jì)算網(wǎng)格數(shù)分別為24萬(wàn)個(gè)和41.4萬(wàn)個(gè)。

2　 結(jié)果與討論

2.1中心間距對(duì)微陣列火焰溫度場(chǎng)的影響

圖2是中心距為4 mm和2 mm工況下，噴管出口10 mm范圍內(nèi)800 K以上溫度區(qū)域分布。從圖中可以看出，中心距4 mm時(shí)各噴管火焰間獨(dú)立性較強(qiáng)，而中心距2 mm工況下火焰出現(xiàn)相互影響和合并現(xiàn)象；距離噴管出口5 mm位置，中心距2 mm工況溫度明顯高于中心距4 mm工況。在中心距4 mm工況時(shí)，中心噴管出口軸線1.1 mm左右到3.5 mm左右高度，溫度由800 K（1）上升到1 832 K（15）；中心距2 mm工況時(shí)，中心噴管出口軸線0.8 mm左右到約2.7 mm高度，溫度由800 K（1）升到1 979 K（17），噴口附近的溫度梯度明顯大于中心距4 mm工況。其主要原因?yàn)橹行木嘤? mm減小到2 mm工況時(shí)，噴管整體布置空間更加緊湊，使得火焰下游相同位置的火焰溫度升高，同時(shí)由于反應(yīng)區(qū)間在空間上的縮小，促使燃燒反應(yīng)區(qū)間在軸向的長(zhǎng)度增加，從而使噴管火焰的長(zhǎng)度變長(zhǎng)；另一方面噴管中心間距的減小使得噴管間火焰出現(xiàn)較強(qiáng)的相互影響作用，火焰出現(xiàn)疊加合并，燃燒反應(yīng)區(qū)域的疊加干涉導(dǎo)致了噴管出口附近溫度梯度的增大。

圖2　噴管出口中心軸面溫度分布Fig. 2 Temperature distributions in center section

圖3為兩種噴管間距工況下以中心點(diǎn)C1為起點(diǎn)的中心軸線溫度分布。中心軸線上的溫度分布更為清晰地表現(xiàn)了兩種結(jié)構(gòu)下火焰及燃燒下游的溫度變化，從圖中可以看出，噴管出口附近的溫度梯度隨著噴管中心間距的縮小而增大，在燃燒區(qū)域下游兩種工況的溫度都趨于平穩(wěn)，但中心距2 mm工況最終溫度比4 mm工況高約200 K。噴口附近溫度梯度的變化主要是因?yàn)楣荛g距的縮小使得中心噴管火焰受到四周?chē)姽芑鹧娴寞B加影響，多噴管火焰向一個(gè)合并火焰發(fā)展，使噴管出口附近整體溫度升高，導(dǎo)致離噴口相同距離內(nèi)的溫度梯度增大；而燃燒最終均溫隨著中心間距的縮小而增大則是噴管整體布置空間更加緊湊所致。

圖3　沿C1點(diǎn)軸線溫度分布Fig. 3 Temperature distributions along the axis line from C1

2.2中心間距對(duì)微陣列火焰速度場(chǎng)影響

圖4為噴口附近的速度分布情況。從圖中可以看出工況4 mm條件下，噴管出口附近的燃燒氣流表現(xiàn)出明顯的獨(dú)立性，相互間不存在明顯的疊加影響區(qū)域；而工況2 mm條件下，噴管出口的燃燒氣流存在明顯的疊加區(qū)域，燃燒區(qū)域下游的速率也明顯大于4 mm工況。其主要原因，一方面是噴管間距的縮小使得燃燒區(qū)間更加緊湊，使得噴管間速度干涉影響作用明顯增強(qiáng)；另一方面是因?yàn)閲姽芑鹧嬖? mm工況下出現(xiàn)了明顯的疊加合并區(qū)域，促進(jìn)了燃燒氣流之間的混合和擴(kuò)散，從而表現(xiàn)出明顯的速度疊加合并區(qū)。

圖5為沿軸向高度方向C1、C2和C3點(diǎn)的軸線速度分布，由于圖1中噴管是對(duì)稱(chēng)分布，同一高度三點(diǎn)的速度分布可以代表所在平面的速度分布。從圖中可以看出，兩種中心距工況下，三點(diǎn)軸線速度最終在下游趨于一致；不同的是中心距2 mm工況下，在軸向高度約10 mm處三點(diǎn)軸線速度趨于1.67 m/s，而中心距4 mm工況下在軸向高度約15 mm處三點(diǎn)軸線速度趨于0.47 m/s，可見(jiàn)火焰速度趨于一致的軸向位置隨著噴管中心間距的縮小而降低；從 C2、C3點(diǎn)的速度分布可以看出，中心距2 mm工況下噴管出口的C2、C3點(diǎn)速度梯度急劇增大，而中心距4 mm工況下速度梯度則較為平緩，其主要原因是噴管間間距的縮小使得噴口附近的火焰間影響強(qiáng)度增加，進(jìn)而溫度梯度變大，高溫區(qū)向噴口附近移動(dòng)，燃燒氣氛被加熱從而使速度增加更快，最終在離噴口更近的地方速度達(dá)到一致。

圖4　噴口附近速度分布Fig. 4 Velocity distributions of nozzle outlet

圖5　沿C1、C2、C3點(diǎn)軸線的速度分布Fig. 5 Velocity distributions along the axis lines from C1， C2 and C3

從速度分布的分析結(jié)果可知，噴管中心間距的縮小使得噴管間氣流在噴口附近出現(xiàn)明顯的疊加干涉區(qū)域，從而促使沿軸向高度方向水平面上的速率更快趨于一致，且速率大小增加，這些因素的改變將增加燃燒加熱場(chǎng)與加熱對(duì)象間的對(duì)流換熱系數(shù)，從而進(jìn)一步提升系統(tǒng)的加熱效率。

2.3微陣列火焰場(chǎng)CO組分濃度分布

為了進(jìn)一步說(shuō)明微陣列火焰環(huán)境中噴管間火焰的相互影響，下面對(duì)兩種結(jié)構(gòu)下的燃燒組分CO的濃度分布進(jìn)行分析。圖6是噴口附近CO的摩爾濃度分布。從圖中可以看出，管間距2 mm工況下火焰燃燒氣氛表現(xiàn)出明顯的擴(kuò)散和相互疊加現(xiàn)象，這進(jìn)一步說(shuō)明隨著噴管中心間距的縮小，噴管火焰之間的疊加干涉作用增強(qiáng)，使得燃燒產(chǎn)物和氣氛相互擴(kuò)散，增加了燃燒場(chǎng)溫度和速度的混合和擴(kuò)散程度。同時(shí)可以從圖中看出，噴管間距的減小導(dǎo)致燃燒反應(yīng)區(qū)間變長(zhǎng)，CO被消耗完的位置明顯向燃燒下游移動(dòng)。

圖6 噴口附近CO摩爾濃度分布Fig. 6 CO mole fraction distributions of nozzle outlet

圖7是沿軸向高度方向C1、C2和C3點(diǎn)的軸線CO濃度分布。從圖中可以看出，在中心間距2 mm工況下，以C2、C3為出發(fā)點(diǎn)的軸線上都有CO分布，距離噴口約1 mm范圍內(nèi)，CO濃度達(dá)到峰值，其值與C1點(diǎn)CO濃度峰值相差較小，說(shuō)明在此位置火焰處于強(qiáng)疊加區(qū)域；在距離噴管出口約2 mm位置以后，三點(diǎn)間的CO濃度趨于一致且數(shù)值逐漸減小，并且在計(jì)算區(qū)域高度的40 mm范圍內(nèi)，仍然有CO濃度的分布；而在中心間距4 mm工況下，因?yàn)閲姽芑鹧孑^為獨(dú)立，以 C2、C3為出發(fā)點(diǎn)的軸線上CO分布除噴管出口區(qū)域外幾乎為零，而噴管出口附近C1點(diǎn)CO濃度最大值遠(yuǎn)高于C2、C3點(diǎn)，在距離噴管出口位置約12 mm后，C1點(diǎn)CO濃度開(kāi)始趨于零。出現(xiàn)這種情況的主要原因是噴管中心間距增大時(shí)，火焰為獨(dú)立的單火焰，在本實(shí)驗(yàn)化學(xué)當(dāng)量比和卷吸空氣的條件下容易燃燒完全；而在噴管中心間距減小時(shí)，一方面因?yàn)榻Y(jié)構(gòu)布置緊湊使得燃燒反應(yīng)區(qū)域變長(zhǎng)，另一方面是由于噴管火焰間相互干涉作用增強(qiáng)，火焰疊加使得噴管間燃燒組分相互擴(kuò)散，火焰合并向大火焰發(fā)展，加上噴管間速度疊加的影響，導(dǎo)致CO被消耗完的距離增加，即計(jì)算高度范圍內(nèi)都有CO濃度分布。

圖7　沿C1、C2、C3點(diǎn)軸線CO摩爾濃度分布Fig. 7 CO mole fraction distributions along the axis lines from C1， C2 and C3

2.4兩種工況微陣列火焰場(chǎng)特性參數(shù)對(duì)比

不同管間距條件下燃燒區(qū)域體積變化較大，為更加準(zhǔn)確地比較兩種工況下微陣列火焰燃燒場(chǎng)的特性變化，以管間距2 mm為計(jì)算基準(zhǔn)，將具有可比性的火焰場(chǎng)特性參數(shù)列在表1中并計(jì)算其變化率，其中出口 CO質(zhì)量流量為計(jì)算區(qū)域出口平面的 CO絕對(duì)質(zhì)量流量積分所得，出口燃燒效率為出口平面CO2質(zhì)量流量除以理論燃燒后CO2質(zhì)量流量所得。

從表1中各參數(shù)的變化率來(lái)看，在相同輸入工況下，管中心間距由2 mm增加為4 mm時(shí)，燃燒區(qū)域體積增加300%，但是出口溫度只下降11.64%；所生成的CO質(zhì)量流量下降92.77%，燃燒效率增加5.59%。即在管中心間距變化時(shí)，同一工況下的燃燒空間體積、單位時(shí)間輸入熱負(fù)荷、污染物CO的排放質(zhì)量將有非常大的變化，而計(jì)算區(qū)域的出口溫度下降率和燃燒效率變化較小。因此從燃料的綜合利用、提升加熱效率和加熱負(fù)荷、減少污染物排放等角度考慮，在一定的化學(xué)當(dāng)量比條件下，微陣列火焰燃燒器的噴管間距、燃燒空間熱負(fù)荷分布、CO污染物排放和出口燃燒效率間應(yīng)有個(gè)匹配的最佳范圍，這是影響微陣列火焰燃燒場(chǎng)特性的重要因素，也是設(shè)計(jì)微陣列火焰加熱系統(tǒng)時(shí)必須考慮的重要因素。

表1　兩種中心間距下的甲烷微陣列燃燒場(chǎng)特性參數(shù)Table 1 Characteristic parameters of combustion field under two nozzle center distance conditions

3　 結(jié) 論

本文對(duì)相同輸入工況下的微陣列火焰燃燒特性進(jìn)行了計(jì)算研究，計(jì)算結(jié)果表明在微陣列火焰中：

（1）在相同管徑和輸入條件下，隨著噴管中心距的減小，噴管出口附近的溫度梯度增大；沿軸向方向水平面的平均溫度增加；噴管分布更加緊湊，使得燃燒下游整體溫度得到提高，有利于提升整體加熱效率和實(shí)現(xiàn)更高的加熱負(fù)荷；

（2）隨著管間距的縮小，噴管間氣流在噴口附近出現(xiàn)明顯的疊加干涉區(qū)域，從而促使沿軸向高度方向水平面上的速率更快趨于一致，且速率大小增加，促進(jìn)了燃燒反應(yīng)的相互影響和火焰疊加合并的實(shí)現(xiàn)；

（3）微小噴管中心間距逐漸減小時(shí)，微陣列火焰由獨(dú)立向合并疊加變化，燃燒反應(yīng)區(qū)間也變長(zhǎng)，CO被消耗完的位置明顯向燃燒下游移動(dòng)，導(dǎo)致燃燒效率下降和CO排放增加；

綜合以上結(jié)論，微陣列火焰燃燒加熱系統(tǒng)設(shè)計(jì)時(shí)，需綜合考慮燃燒效率、提升加熱效率和加熱負(fù)荷、減少污染物排放等因素進(jìn)行噴管中心間距的最優(yōu)設(shè)計(jì)。

［1］ MATTA L M， NEUMEIER Y， LEMON B， et al. Characteristics of microscale diffusion flames［J］. Proceedings of combustion institute， 2002， 29（1）: 993-939. DOI: 10.1016/S1540-7489（02）80118-3.

［2］ 趙黛青， 蔣利橋， 黃顯峰， 等. 微尺度預(yù)混合火焰結(jié)構(gòu)和熄火特性研究［J］. 工程熱物理學(xué)報(bào)， 2006， 27（4）: 711-713. DOI: 10.3321/j.issn:0253-231X.2006.04.055.

［3］ 蔣利橋， 趙黛青， 汪小憨. 微尺度甲烷擴(kuò)散火焰及其熄滅特性［J］. 燃燒科學(xué)與技術(shù)， 2007， 13（2）: 183-186. DOI: 10.3321/j.issn:1006-8740.2007.02.016.

［4］ MARUTA K. Micro and mesoscale combustion［J］. Proceedings of the combustion institute， 2011， 33（1）: 125-150. DOI: 10.1016/j.proci.2010.09.005.

［5］ NAKAMURA Y， HIRASAWA T. Study on microflame array combustor［C］//Proceedings of 2009 China-Japan Workshop on Combustion in Micro Power System. Guangzhou， China: ［s.n.］， 2009.

［6］ LEE B J， KIM J S， LEE S. Enhancement of blowout limit by the interaction of multiple nonpremixed jet flames［J］. Combustion science and technology， 2004，176（4）: 481-497. DOI: 10.1080/00102200490276700.

［7］ 李鵬飛， 徐敏義， 王飛飛. 精通 CFD工程仿真與案例實(shí)戰(zhàn)［M］. 北京: 人民郵電出版社， 2011: 493-505.

Combustion Characteristics of Array Combustor for Methane Premixed Microflame under Two Nozzle Center Distance Conditions

ZENG Xiao-jun1，2，3， WANG Xiao-han1，2，3， ZHAO Dai-qing1，2，3，JIANG Li-qiao1，2，3， YANG Wei-bin1，2，3
（1. Guangzhou Institute of Energy Conversion， Chinese Academy of Sciences， Guangzhou 510640， China；2. Key Laboratory of Renewable Energy， Chinese Academy of Sciences， Guangzhou 510640， China；3. Guangdong Provincial Key Laboratory of New and Renewable Energy Research and Development， Guangzhou 510640， China）

A numerical model of 3 × 3 array combustor consisting of Ф1.5 mm tubes for methane premixed microflame was established to learn the effects of nozzle center distance on the combustion characteristics by using the reduced mechanism for methane. The results show that the effect among micro flames is enhanced by decreasing the nozzle center distance， resulting in a lengthening region of combustion reaction and an increase of the highest combustion temperature as well as an increase in CO emissions and decrease in combustion efficiency at the outlet of calculating region. Temperature uniformity in the heating field was obtained after 5 mm from the nozzles under different nozzle center distances of 2 mm and 4 mm， which provides a design basis for the micro array combustion heater.

micro combustion； array combustor； temperature uniformity； flame heating； methane

趙黛青（1958-），女，博士，研究員，從事微尺度流動(dòng)和燃燒解析的研究。

TK16

A

10.3969/j.issn.2095-560X.2016.03.006

2095-560X（2016）03-0201-05

2016-03-02

2016-04-05

國(guó)家自然科學(xué)基金面上項(xiàng)目（51176193）；國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃（2014CB239600）

曾小軍（1980-），男，工學(xué)博士，助理研究員，主要從事微尺度燃燒及新型燃燒器開(kāi)發(fā)。