閆 申，高 峰

(1.山西經(jīng)濟(jì)管理干部學(xué)院，山西 太原 030024；2.太原理工大學(xué) 材料科學(xué)與工程學(xué)院，山西 太原 030024)

高溫空氣燃燒技術(shù)(high temperature air combustion(HTAC))是20世紀(jì)90年代開發(fā)成功的一項(xiàng)燃料燃燒領(lǐng)域中的新技術(shù)。HTAC包括兩項(xiàng)基本技術(shù)手段：一是通過性能優(yōu)良的蓄熱回收裝置使燃燒產(chǎn)物顯熱極限回收；二是燃料在低氧氣氛下燃燒，達(dá)到溫室氣體排放量減少的目的。HTAC技術(shù)可以降低30%的二氧化碳排放和50%的污染物排放，同時(shí)減小25%的設(shè)備尺寸[1]。

在我國，大部分倒焰窯爐的熱利用效率較低，同時(shí)排放大量的二氧化硫、氮氧化物等污染物。目前高溫空氣燃燒技術(shù)已經(jīng)在很多行業(yè)的工業(yè)窯爐中有所應(yīng)用，例如冶金和玻璃窯爐等[5]?，F(xiàn)在經(jīng)濟(jì)發(fā)展新常態(tài)下，節(jié)能減排及環(huán)保檢查的常態(tài)化，高溫空氣燃燒技術(shù)作為一項(xiàng)節(jié)能、環(huán)保技術(shù)也一定會(huì)引起關(guān)注。

通過查閱工業(yè)窯爐數(shù)值模擬、HTAC技術(shù)的大量相關(guān)資料，本文利用FLUENT軟件對(duì)采用高溫空氣燃燒技術(shù)的新型倒焰窯爐進(jìn)行數(shù)值模擬，通過改變?nèi)剂蠂娙胨俣?，研究窯爐內(nèi)部的溫度變化、組分分布的情況。

1 模型的建立

1.1 模擬高溫窯爐的構(gòu)造

模擬對(duì)象為一個(gè)小型窯爐，有效容積為2m3，其尺寸是4403mm×3800mm×4104mm，該窯爐兩側(cè)分別設(shè)置有一個(gè)燃料燒嘴；底部分成兩部分，空氣通過一側(cè)底部煙道進(jìn)入窯爐內(nèi)與燃料反應(yīng)，煙氣從另一側(cè)煙道排出。

燃料燒嘴的尺寸大小為Φ30mm；燃料燒嘴中心和底部煙道的距離為570mm；底部分布煙道24個(gè)，長×寬=115×30mm。窯爐的結(jié)構(gòu)如圖1。

圖1 窯爐模型 Fig.1 Kiln model

1.2 窯爐模型的網(wǎng)格劃分

網(wǎng)格的劃分質(zhì)量對(duì)數(shù)值模擬的計(jì)算精度和結(jié)果的收斂性有很大影響，一般情況下網(wǎng)格質(zhì)量越高計(jì)算出的結(jié)果精度也越高。在對(duì)計(jì)算區(qū)域進(jìn)行網(wǎng)格劃分時(shí)，主要考慮的內(nèi)容有：初始化的時(shí)間、計(jì)算花費(fèi)和數(shù)值耗散[2]。

在進(jìn)行網(wǎng)格劃分時(shí)，根據(jù)實(shí)際情況，首先將模型劃分為不同的多個(gè)區(qū)域，然后分別對(duì)每一個(gè)區(qū)域進(jìn)行網(wǎng)格劃分[3]。計(jì)算模型對(duì)燒嘴局部區(qū)域進(jìn)行網(wǎng)格加密，保證計(jì)算的精度。計(jì)算區(qū)域的網(wǎng)格劃分見圖2。

圖2 計(jì)算區(qū)域的網(wǎng)格劃分 Fig.2 Grid division of regional

2 模型條件設(shè)置

2.1 基本條件設(shè)置

湍流模型采用Reynolds Stress模型和κ-ε兩方程紊流模型。因?yàn)镽eynolds Stress模型可以精確的預(yù)測復(fù)雜流動(dòng)模型，而κ-ε兩方程模型計(jì)算簡單、經(jīng)濟(jì)、方便和成熟[3-4]。燃燒模型選用甲烷兩步完成反應(yīng)機(jī)制，燃料燃燒系統(tǒng)用finite-rate化學(xué)模型分析[3]。

模型的幾何結(jié)構(gòu)簡單，故熱輻射模型采用計(jì)算精度適中，占用計(jì)算機(jī)內(nèi)存不大的P-1輻射模型[3]。

方程求解器采用SIMPLEC方法[3]?？刂品匠痰碾x散格式選用精度較高的二階迎風(fēng)格式，各方程的松弛因子采用默認(rèn)值[3]。所有參數(shù)設(shè)定完畢后，最后對(duì)方程迭代求解。

2.2 計(jì)算條件

燃?xì)膺M(jìn)口采用速度邊界條件，空氣入口和出口都采用壓力邊界條件，燃燒室墻壁使用壁面邊界條件，溫度為固定值300K。

為了研究燃料流速度對(duì)對(duì)窯爐內(nèi)的溫度場的影響。對(duì)模型進(jìn)行了燃料射流速度為10m/s、20m/s、30m/s、35m/s、40m/s等五個(gè)工況下爐內(nèi)溫度場，組分場數(shù)值模擬，計(jì)算工況見表1。

表1 計(jì)算工況Tab.1 Calculation conditions

3 計(jì)算結(jié)果

3.1 燃料射流速度對(duì)窯爐內(nèi)部溫度場的影響

燃料射流速度不同時(shí)，燃料燒嘴所在截面的溫度分布見圖3～圖7。

圖3 A的溫度分布圖 Fig.3 Temperature distribution of A

圖4 B的溫度分布圖 Fig.4 Temperature distribution of B

圖5 C的溫度分布圖 Fig.5 Temperature distribution of C

圖6 D的溫度分布圖 Fig.6 Temperature distribution of D

圖7 E的溫度分布圖 Fig.7 Temperature distribution of E

從圖3～圖7可以看出，提高燃料噴入速度，燒嘴截面上的最高溫度升高，同時(shí)高溫區(qū)域逐步擴(kuò)大，溫差逐步減小。當(dāng)燃料以較低速度噴入時(shí)，整個(gè)窯爐內(nèi)部溫度較低，只有在燒嘴的附近有較高溫度存在。這主要是因?yàn)槿剂蠂娙胨俣容^低，燃燒發(fā)生在燒嘴周圍，爐內(nèi)氣體流動(dòng)性差，擾動(dòng)少，傳熱較差，故此，高溫區(qū)域較集中，大部分區(qū)域溫度較低，溫度分部不均勻。燃料噴入速度的進(jìn)一步增大，燃燒放出的熱能量增加，最高溫度在升高，區(qū)域范圍在擴(kuò)大，這樣達(dá)到燃燒條件的燃?xì)饩驮蕉?。同時(shí)，由于燃料速度的增大，卷吸大量熱煙氣，提高氣體的湍流強(qiáng)度，帶動(dòng)周圍氣體的流動(dòng)，有利于傳熱，從而使溫度的差值減小。下表2列出了不同工況下模擬得到的爐內(nèi)溫度。其最高溫度、平均溫度為模擬計(jì)算得到的窯爐內(nèi)的煙氣溫度的最大值、平均值。

表2 最高溫度與平均溫度比較Tab2 Comparison of the maximum temperature and average temperature

表2可以看出，燃料噴出速度從10m/s～40m/s，窯爐內(nèi)部的溫差在減少，窯爐內(nèi)部的溫度逐漸均勻化。

3.2 燃料射流速度對(duì)窯爐內(nèi)部組分場的影響

燃料的射流速度會(huì)對(duì)組分場產(chǎn)生影響。對(duì)CH4、O2和CO2組分場的影響分布圖見圖8～圖26。

圖8 A的甲烷濃度分布 Fig.8 Methane concentration of A

圖9 B的甲烷濃度分布 Fig.9 Methane concentration of B

圖10 C的甲烷濃度分布 Fig.10 Methane concentration of C

圖11 D的甲烷濃度分布 Fig.11 Methane concentration of D

圖12 E的甲烷濃度分布 Fig.12 Methane concentration of E

圖13 A的氧氣濃度分布 Fig.13 Oxygen concentration of A

圖14 B的氧氣濃度分布 Fig.14 Oxygen concentration of B

圖15 C的氧氣濃度分布 Fig.15 Oxygen concentration of C

圖16 D的氧氣濃度分布 Fig.16 Oxygen concentration of D

圖17 E的氧氣濃度分布 Fig.17 Oxygen concentration of E

圖18 A的二氧化碳濃度分布 Fig.18 Carbon dioxide concentration of A

圖19 B的二氧化碳濃度分布 Fig.19 Carbon dioxide concentration of B

圖20 C的二氧化碳濃度分布 Fig.20 Carbon dioxide concentration of C

圖21 D的二氧化碳濃度分布 Fig.21 Carbon dioxide concentration of D

圖22 E的二氧化碳濃度分布 Fig.22 Carbon dioxide concentration of E

圖23 燃料燒嘴軸線上CH4的濃度圖 Fig.23 Methane concentration of central axis

圖8～圖12可知，燃料噴口速度的加快，帶動(dòng)周圍氣體，使燃燒的區(qū)域遠(yuǎn)離窯爐的內(nèi)壁，這樣就保護(hù)了窯爐內(nèi)壁。不同工況下，燃料都可以很快發(fā)生燃燒反應(yīng)，且燃燒完全。當(dāng)燃料以較高速度噴入時(shí)，加大對(duì)爐內(nèi)煙氣和空氣的擾動(dòng)，有利于熱量的散發(fā)，從而降低窯爐內(nèi)部的溫差。

圖13～圖17可以發(fā)現(xiàn)，燃料噴口速度的加快，更多的氧氣參與燃燒反應(yīng)，并且使氧氣的流動(dòng)更加的快速，促進(jìn)熱量的散發(fā)。各工況下，由于燃燒反應(yīng)的徹底進(jìn)行，在燒嘴附近的區(qū)域形成第一個(gè)低氧區(qū)，但是在當(dāng)燃料速度大于30m/s時(shí)，在窯爐的出口處，形成第二個(gè)低氧區(qū)。隨著工況燃料速度的加大，第二個(gè)低氧區(qū)的范圍逐漸的加大，主要原因是流速加大，促進(jìn)煙氣的流動(dòng)，從而形成一個(gè)小的漩渦，進(jìn)一步稀釋空氣而造成的。

圖18～圖22可以發(fā)現(xiàn)，圖形外觀上看，二氧化碳的濃度分布和氧氣的濃度分布大小正好互補(bǔ)。

綜合以上，燃料噴口速度加快，燃料量增加，燃料在窯爐內(nèi)的濃度提高，有利于燃燒反應(yīng)和稀釋氧氣，放出大量的熱；同時(shí)加強(qiáng)了窯爐內(nèi)部煙氣的擾動(dòng)，有利于進(jìn)一步減少窯爐內(nèi)部的溫差，可以減少污染物的排放。

燃料的射流速度不同時(shí)，燃料燒嘴軸線上的CH4、O2和CO2的分布分別見圖23～圖26。其中圖24是圖23中A區(qū)域的局部放大圖。

圖24 燃料燒嘴軸線上CH4的濃度局部放大圖 Fig.24 Methane concentration

圖25 燃料燒嘴軸線上O2的濃度圖 Fig.25 Oxygen concentration of central axis enlarged drawing of central axis

從圖23、圖24可以看出，燃料噴出燒嘴后，燒嘴軸線上的燃料的濃度迅速降低。這主要是由于來自底部的空氣對(duì)燃料的吹動(dòng)作用，使燃料向上偏移。燃料的速度越快，越穩(wěn)定，空氣越難吹動(dòng)燃料，使其脫離軸線。

圖26 燃料燒嘴軸線上CO2的濃度圖 Fig.26 Carbon dioxide concentration of central axis

圖26可以看出，第一階段二氧化碳的濃度隨著燃燒反應(yīng)的進(jìn)行而提高；在后期，由于出口壓力的吸力作用，二氧化碳的濃度降低。整個(gè)階段，燃料噴口的速度越高，反應(yīng)產(chǎn)生的二氧化碳就越多，變化也越大，可見，提高燃料的速度，能夠加強(qiáng)窯爐內(nèi)部煙氣的流動(dòng)。

4 結(jié)論

(1)燃料噴口速度的加快，可以強(qiáng)化窯爐內(nèi)部擾動(dòng)，燒嘴所在面上的最高溫度升高，高溫區(qū)進(jìn)一步擴(kuò)大，窯爐內(nèi)部溫差逐步減小；

(2)隨著燃料噴口速度的提高，爐內(nèi)煙氣流動(dòng)加強(qiáng)，氧氣在燃料燒嘴所在面上高濃度的區(qū)域進(jìn)一步減少，低濃度氧氣的區(qū)域增加；在低濃度氧氣區(qū)域二氧化碳的濃度較高，高氧濃度區(qū)域二氧化碳的濃度低。與傳統(tǒng)的倒焰窯相比，二氧化碳的排放量整體上是降低了。

[1] 劉曉東.常溫空氣MILD的實(shí)驗(yàn)研究與火焰面模型的模擬[D].武漢:華中科技大學(xué),2011.

[2] 陶文銓.數(shù)值傳熱學(xué)[M].2版.西安:西安交通大學(xué)出版社,2001.

[3] 韓占忠,王 敬,蘭小平,等.FLUENT流體工程仿真計(jì)算實(shí)例與應(yīng)用[M].北京.北京理工大學(xué)出版社,2004.

[4] 于 蕾,李 敏,林 豹,等.FLUENT軟件在射流模擬中的應(yīng)用[J].科技情報(bào)開發(fā)與經(jīng)濟(jì),2005,15(13)：154-155.

[5] 程小蘇,柯善軍,曾令可,等.高溫空氣燃燒技術(shù)在陶瓷窯爐的應(yīng)用分析[J].工業(yè)爐,2009，30(3):1-4.

[6] 呂 軍,唐一科,歐陽奇.高溫空氣對(duì)燃燒回流區(qū)流動(dòng)影響的數(shù)值模擬[J].物理測試,2008,26(1):22-25.

[7] 羅國民,郭漢杰,溫志紅.高溫空氣燃燒和富氧燃燒在加熱爐生產(chǎn)應(yīng)用中的對(duì)比研究[J].工業(yè)爐2008,37(5):44-48.

[8] 高力明.中國陶瓷熱工技術(shù)和窯爐的發(fā)展水平與前景[J].中國陶瓷,2005,41(1):1-6.

[9] Tong J H,Feng Q,Wang H P.Numerical simulation of the gas flow sand temperature field in the firing zone of roller heath kiln [J]. China Ceramics,2006,42(5):27-31 .

[10] 宋新南.倒焰窯熱工特性與改進(jìn)方向研究[J].工業(yè)加熱,2002(2):39-44.