朱慶霞馮 青王化能付海丹

（1.景德鎮(zhèn)陶瓷學(xué)院材料科學(xué)與工程學(xué)院，江西省先進(jìn)陶瓷材料重點(diǎn)實(shí)驗(yàn)室，江西景德鎮(zhèn) 333001；2.杭州諾貝爾集團(tuán)有限公司，浙江杭州 311122）

梭式窯富氧燃燒高溫階段換熱特性的模擬研究

朱慶霞1馮 青1王化能2付海丹1

（1.景德鎮(zhèn)陶瓷學(xué)院材料科學(xué)與工程學(xué)院，江西省先進(jìn)陶瓷材料重點(diǎn)實(shí)驗(yàn)室，江西景德鎮(zhèn) 333001；2.杭州諾貝爾集團(tuán)有限公司，浙江杭州 311122）

采用計(jì)算流體力學(xué)（CFD）軟件——FLUENT研究了富氧濃度對(duì)高溫階段梭式窯內(nèi)換熱特性的影響。結(jié)果表明采用富氧空氣助燃技術(shù)，提高窯內(nèi)溫度的同時(shí)也增大了窯內(nèi)溫差。分析了產(chǎn)生截面溫差的原因，為梭式窯富氧燃燒的穩(wěn)定運(yùn)行奠定理論基礎(chǔ)。

梭式窯；富氧燃燒；富氧濃度；換熱特性；溫差

0 引言

富氧燃燒具有很多優(yōu)點(diǎn),如提高熱效率、提高處理速率、減少煙塵排放、減小裝置尺寸等，是一種高效節(jié)能的燃燒技術(shù)，被發(fā)達(dá)國(guó)家稱之為“資源的創(chuàng)造性技術(shù)”[1～2]。富氧燃燒在回轉(zhuǎn)窯和玻璃熔窯中的應(yīng)用研究較多[3～4]，在梭式窯方面的應(yīng)用研究相對(duì)較少。梭式窯高溫段的流動(dòng)與換熱狀況是影響產(chǎn)品質(zhì)量和節(jié)能效率等的重要因素，因此掌握梭式窯預(yù)混富氧燃燒高溫階段換熱特性對(duì)于窯爐的節(jié)能降耗和結(jié)構(gòu)優(yōu)化等起著至關(guān)重要的作用。由于陶瓷窯爐內(nèi)部結(jié)構(gòu)十分復(fù)雜，對(duì)流場(chǎng)、氣氛?qǐng)龊蜏囟葓?chǎng)的測(cè)量顯得十分困難，實(shí)驗(yàn)研究仍然不足以全面揭示流動(dòng)和傳熱規(guī)律。計(jì)算流體力學(xué)FLUENT軟件是一個(gè)很好的解決辦法。FLUENT軟件設(shè)計(jì)針對(duì)每一種流動(dòng)的物理問(wèn)題的特點(diǎn)，采用適合它的數(shù)值解法在計(jì)算速度、穩(wěn)定性和精度等各方面達(dá)到最佳，省卻了科研工作者在計(jì)算方法、編程、前后處理等方面投入的重復(fù)的、低效的勞動(dòng)，而將主要的精力、智慧用于物理本身的探索上[5]。目前，F(xiàn)LUENT軟件在陶瓷窯爐模擬研究中的應(yīng)用日益廣泛[6-7]，宮小龍[8]等驗(yàn)證了利用計(jì)算流體力學(xué)(CFD)數(shù)值模擬技術(shù)模擬研究梭式窯的可行性。因此，本文在此研究基礎(chǔ)上，利用FLUENT軟件對(duì)采用不同氧濃度的富氧空氣整體助燃時(shí)梭式窯內(nèi)的換熱特性進(jìn)行數(shù)值模擬研究，以期為梭式窯富氧燃燒的穩(wěn)定運(yùn)行奠定理論基礎(chǔ)。

1 模型的建立

本文以0.1m3液化石油氣梭式窯為模擬對(duì)象，在建立模型之前，先做如下簡(jiǎn)化：（1）幾何體尺寸以實(shí)際梭式窯的尺寸為準(zhǔn)；（2）窯門與窯墻之間不存在任何縫隙，即燒成過(guò)程中不會(huì)漏入冷空氣；（3）梭式窯燒嘴入口僅僅為固定空燃比的一個(gè)入口（模擬預(yù)混燃燒），不考慮二次空氣；（4）料垛以空心圓柱體表示，在窯內(nèi)按模擬要求分布。

由于窯體幾何結(jié)構(gòu)比較復(fù)雜，僅用結(jié)構(gòu)化網(wǎng)格來(lái)離散計(jì)算區(qū)域是非常困難的；倘若完全使用非結(jié)構(gòu)化網(wǎng)格，由于窯內(nèi)料垛排列緊密造成網(wǎng)格數(shù)目的龐大而大大增加計(jì)算機(jī)的計(jì)算工作量，從而增加了計(jì)算時(shí)間，同時(shí)也難于控制網(wǎng)格疏密造成局部網(wǎng)格過(guò)細(xì)，加大數(shù)值計(jì)算的舍入誤差。為了克服以上困難，本文采用結(jié)構(gòu)化和非結(jié)構(gòu)化網(wǎng)格相結(jié)合的方法，將梭式窯計(jì)算區(qū)域分為燒嘴區(qū)域、料垛區(qū)域和周邊區(qū)域三個(gè)部分，分別生成三對(duì)彼此獨(dú)立的網(wǎng)格（即燒嘴區(qū)域和料垛區(qū)域采用結(jié)構(gòu)化網(wǎng)格離散，周邊區(qū)域采用非結(jié)構(gòu)化網(wǎng)格離散），燒嘴出口、煙氣出口、料垛的周邊等速度梯度較大的區(qū)域進(jìn)行局部加密。具體網(wǎng)格劃分圖見(jiàn)圖1。

表1 不同氧濃度條件下的燒嘴出口速度值及模擬的邊界條件Tab.1 Exit velocities at different oxygen concentrations and boundary conditions for simulation

1.1 邊界條件及物性參數(shù)的確定

為了便于邊界條件的確定，現(xiàn)對(duì)模型進(jìn)行如下假設(shè)：（1）窯車（或窯底）、窯門與窯墻的材料近似為全陶瓷纖維材料或耐火粘土磚；（2）環(huán)境溫度為300K（27℃）；（3）在模擬過(guò)程中忽略坯體和釉料的物理化學(xué)反應(yīng)熱效應(yīng)，即僅僅包括燃料的燃燒，窯墻，坯體本身及煙氣的換熱過(guò)程；（4）梭式窯內(nèi)某一離散溫度點(diǎn)的燃燒為三維穩(wěn)態(tài)預(yù)混湍流燃燒；（5）氣體為不可壓理想氣體；（6）僅考慮煙囪底部的壓力，不計(jì)算煙囪內(nèi)部煙氣的流動(dòng)情況。

1.2 速度入口及壓力出口邊界條件

據(jù)相關(guān)文獻(xiàn)表明,富氧燃燒技術(shù)節(jié)能效果最好是在富氧空氣中氧濃度在21%～30%范圍內(nèi)[9]，因此，本文將分別模擬研究氧濃度為21%、24%、27%時(shí)梭式窯內(nèi)的溫度場(chǎng)。為了貼近實(shí)驗(yàn)工況，離散成3個(gè)溫度點(diǎn)（1050℃、1200℃、1300℃），保持各溫度點(diǎn)燃料入口速度不變。

經(jīng)過(guò)計(jì)算[10]，梭式窯熱態(tài)模擬所需的邊界條件如表1所示，單位均采用國(guó)際標(biāo)準(zhǔn)單位。針對(duì)不同氧濃度，根據(jù)離散溫度點(diǎn)的燃料入口速度，通過(guò)燃料燃燒學(xué)和流體力學(xué)的相關(guān)知識(shí)可計(jì)算出相應(yīng)狀態(tài)下的空氣入口速度和燒嘴出口速度。煙氣出口的壓力即煙囪底部的壓力，可根據(jù)窯內(nèi)不同氣氛取適當(dāng)?shù)膲毫χ?。在FLUENT里對(duì)坯體的邊界條件設(shè)置時(shí)，隨著窯內(nèi)溫度的升高，由于輻射所傳遞的熱量急劇上升，換熱快且均勻，此時(shí)制品表面取等熱流密度會(huì)對(duì)整個(gè)模擬計(jì)算造成很大的誤差，所以坯體表面在高溫時(shí)取等溫度邊界條件或耦合邊界條件[11]。應(yīng)用傳熱學(xué)的相關(guān)知識(shí)，可估算窯墻的熱流密度。

1.3 數(shù)值計(jì)算方法

在數(shù)值方法上，選用FLUENT提供的分離求解器。為了提高計(jì)算精度同時(shí)加快收斂速度，在初始計(jì)算中先計(jì)算流場(chǎng)，達(dá)到收斂以后引入能量方程和動(dòng)量方程。并且能量方程和動(dòng)量方程選用一階迎風(fēng)格式，收斂以后選用二階迎風(fēng)格式，再迭代至最終收斂。在算法上選用SIMPLEC算法以提高收斂性。

2 模擬結(jié)果分析與討論

在坯體燒制過(guò)程中，600℃～1050℃屬氧化階段，1050℃～1200℃屬還原階段，1200℃～1300℃以上為燒結(jié)階段，此時(shí)窯內(nèi)換熱由預(yù)熱階段的對(duì)流換熱為主逐漸隨溫度的升高轉(zhuǎn)為以輻射換熱為主。一般在這一階段，燒嘴全開(kāi)。為了貼近實(shí)驗(yàn)工況要求，高溫階段模擬研究也在梭式窯四角點(diǎn)“點(diǎn)”四只燒嘴。

本文以1300℃離散溫度點(diǎn)為例，對(duì)采用不同氧濃度的富氧空氣進(jìn)行整體增氧助燃時(shí)梭式窯內(nèi)換熱特性進(jìn)行分析。從圖2、圖3的窯爐中部（坯體頂部）和窯爐底部（坯體底面）兩截面的溫度等值線圖可以看出，當(dāng)采用普通空氣助燃時(shí)，整個(gè)窯爐中部及底部的截面溫度均勻，截面溫差較小。當(dāng)采用富氧空氣助燃時(shí)，隨著富氧濃度的增加，截面溫度逐漸提高，這主要是因?yàn)樵诟邷仉A段，煙氣與制品的換熱主要以輻射換熱為主，由氣體輻射特點(diǎn)可知，只有三原子和多原子氣體具有輻射能力，而雙原子氣體幾乎無(wú)輻射能力，富氧燃燒減少了無(wú)輻射能力的氮?dú)饬?，也就相?yīng)增加了煙氣對(duì)制品的換熱，從而使窯內(nèi)的溫度普遍提高；但當(dāng)采用富氧空氣進(jìn)行助燃時(shí)，窯內(nèi)截面的溫差也會(huì)逐漸變大，這主要是因?yàn)殡S著氧濃度的增加，火焰的長(zhǎng)度變短，火焰燃燒變急，造成窯內(nèi)的輻射換熱不均勻，進(jìn)而造成窯內(nèi)各個(gè)截面溫差變大。

從圖2還可以看出，在氧濃度為24%時(shí)，y=65mm處窯爐中部截面噴火口處的溫度最高。這是因?yàn)椴捎酶谎蹩諝庵紩r(shí),隨著富氧濃度的增大，噴火口處火焰的溫度場(chǎng)溫度梯度逐漸增大[12]。而富氧濃度在21～24%變化范圍內(nèi)，火焰增加幅度又比較大（窯內(nèi)平均溫度增加90℃），加之其火焰燃燒的溫度場(chǎng)溫度梯度較氧濃度為27%時(shí)的溫度梯度小，因而出現(xiàn)了在窯爐中部截面上氧濃度為24%時(shí)，噴火口處溫度最高的反常現(xiàn)象。而氧濃度為27%時(shí)雖然平均溫度最高，但由于其火焰溫度梯度大，窯爐中部截面（y=65mm）可能偏離其火焰最高溫度面較遠(yuǎn)。

窯爐底部（坯體底面）截面(y=－310mm)溫度比窯爐中部（坯體頂部）截面(y=65mm)平均溫度要低，局部區(qū)域兩截面溫差最高達(dá)到100～200℃左右，這主要是因?yàn)楦G爐底面氣流流速低，換熱狀況較窯爐中部弱，加之窯爐底部截面又處于燒嘴燃燒火焰的根部，從而產(chǎn)生窯爐底面溫度比窯爐中部溫度低的現(xiàn)象；采用富氧空氣助燃時(shí)，兩截面溫差更大，這同樣是由于富氧燃燒會(huì)增大火焰燃燒的溫度梯度。因此，在采用富氧空氣助燃技術(shù)，有必要加設(shè)氣流擾動(dòng)裝置，在整體提高窯內(nèi)溫度的同時(shí)，均勻溫度場(chǎng)，避免局部高溫惡化產(chǎn)品質(zhì)量以及破壞砌體材料。

從圖2～3還可以看出，當(dāng)氧濃度為21%時(shí)，窯尾溫度要比窯門方向溫度高出40～85℃；當(dāng)氧濃度為24%、27%時(shí)，窯尾溫度比窯門方向高出50～120℃。這主要是因?yàn)榕艧熆谠O(shè)在窯尾，窯頭煙氣比窯尾煙氣在窯內(nèi)滯留時(shí)間長(zhǎng)，導(dǎo)致梭式窯內(nèi)局部對(duì)流換熱速度不同以及窯內(nèi)換熱不均勻；氧濃度越大，火焰燃燒越急，窯內(nèi)熱輻射就越不均勻，所以截面溫差越大。因此建議無(wú)論采用普通空氣助燃方法還是富氧空氣助燃方法，窯內(nèi)煙氣出口的位置都不應(yīng)全部集中在窯尾底部排煙，應(yīng)該在窯頭設(shè)計(jì)適當(dāng)?shù)呐艧熆?，使氣流排出口均勻分布在窯底，這樣能使煙氣射流在窯內(nèi)分布均勻，減少窯內(nèi)截面溫差，有利于對(duì)坯體的加熱。

在圖3中，當(dāng)采用普通空氣助燃時(shí)，窯爐底面換熱較弱，溫度較低，在窯爐底面同一水平截面上，左右兩側(cè)中間部位出現(xiàn)了溫度偏高的區(qū)域帶，造成上述現(xiàn)象的主要原因是研究的梭式窯有效容積較?。?.1m3），料垛排列相對(duì)較密，煙氣不能充分沖刷燒嘴附近區(qū)域料垛，加之燒嘴煙氣射流頂部回流的出現(xiàn)，造成部分料垛間隙區(qū)域的氣流流動(dòng)緩慢，湍流度不高；而兩側(cè)燒嘴是呈對(duì)稱布置，在窯爐左右兩側(cè)對(duì)稱面處區(qū)域煙氣射流匯集，氣體流速較大，湍流度大，換熱較強(qiáng)而導(dǎo)致間隙溫度偏高。因此建議無(wú)論采用普通空氣助燃方法還是富氧空氣助燃方法，料垛中部應(yīng)適當(dāng)碼密些，可在一定程度上克服料垛區(qū)域左右方向中間部位溫度偏高的現(xiàn)象，緩解局部換熱不均勻的問(wèn)題。

3 結(jié)論

（1）運(yùn)用FLUENT軟件建立了0.1m3梭式窯的幾何模型，并模擬了不同氧濃度的富氧空氣助燃時(shí)高溫階段窯內(nèi)的換熱特性，為梭式窯的富氧燃燒運(yùn)行提供了有益的參考。

（2）采用富氧空氣助燃技術(shù)，能使窯內(nèi)的截面溫度普遍得到提高，但是同時(shí)也使窯內(nèi)同一水平截面溫差增加，這主要是因?yàn)殡S著氧濃度的增加，火焰的長(zhǎng)度將會(huì)變短，火焰燃燒變急，造成窯內(nèi)的輻射換熱不均勻，進(jìn)而造成窯內(nèi)溫差變大。

（3）為了減少窯頭窯尾的溫差，應(yīng)均勻布置排煙口。料垛中部適當(dāng)碼密些，在一定程度上可克服料垛左右方向中部區(qū)域溫度偏高現(xiàn)象，緩解料垛局部換熱不均勻的問(wèn)題。

1 劉毅,李光界,孫鶴等.富氧助燃技術(shù)及其應(yīng)用.節(jié)能與環(huán)保,2005,2:28～29

2 鄧偉強(qiáng),曾令可,稅安澤等.富氧燃燒技術(shù)在陶瓷燒成工藝中的應(yīng)用.工業(yè)爐,2006,28(6):12～15

3 陳文仲,王春華,田遠(yuǎn)航等.回轉(zhuǎn)窯熱工狀況影響參數(shù)的數(shù)值模擬.化工學(xué)報(bào),2011,62(1):47～52

4 韓韜,劉宗明,李賢松等.單元玻璃窯爐富氧燃燒空間的數(shù)值模擬研究.硅酸鹽通報(bào),2010,29(1):77～82

5 李勇,劉志友,安亦然.介紹通用流體力學(xué)計(jì)算軟件——FLUENT.水動(dòng)力研究與進(jìn)展,2001,16(2):255～259

6 曾令可,鄧偉強(qiáng)等.NOx在陶瓷窯爐內(nèi)生成的數(shù)值模擬.華南理工大學(xué)學(xué)報(bào)(自然科學(xué)版),2007,35(2):102～106

7 童劍輝,馮青,汪和平.輥道窯燒成帶氣體流場(chǎng)與溫度場(chǎng)的三維數(shù)值模擬.中國(guó)陶瓷,2006,42(5):27～31

8 宮小龍,湯文菊,馮青等.液化氣梭式窯高溫段流場(chǎng)與換熱的數(shù)值模擬研究.中國(guó)陶瓷,2008,44(4):35～37

9 潘小勇,宮小龍,馮青等.基于輥道窯富氧燃燒計(jì)算與分析.陶瓷學(xué)報(bào),2010,31(4):627～631

10 付海丹.梭式窯富氧預(yù)混燃燒火焰特性及窯內(nèi)換熱特性研究.景德鎮(zhèn):景德鎮(zhèn)陶瓷學(xué)院,2011

11 孫晉濤.硅酸鹽工業(yè)熱工基礎(chǔ).湖北:武漢工業(yè)大學(xué)出版社,1999

12 朱慶霞,馮青,付海丹.預(yù)混富氧燃燒火焰特性的實(shí)驗(yàn)研究.陶瓷學(xué)報(bào),2012,33(2):627～631

Numerical Research on Heat-transfer Characteristics of Oxygenenriched Combustion in the High-temperature Stage of Shuttle Kiln

ZHU Qingxia1FENG Qing1WANG Huaneng2FU Haidan1
(1.Key Laboratory of Advanced Ceramics of Jiangxi Province,School of Materials Science and Engineering,Jingdezhen Ceramic Institute,Jingdezhen Jiangxi 333001,China;2.Hangzhou Nabel Group.co.,Ltd,Hangzhou Zhejiang 311122,China)

The effect of oxygen concentration on the heat-transfer characteristics in the high-temperature stage of shuttle kiln was investigated by computational fluid dynamics (CFD)software-FLUENT.The results show that the oxygen-enriched combustion technology can raise the temperature,and increase the temperature difference.Causes of cross-section temperature difference were analyzed to find the theory basis for the stable operation of oxygen-enriched combustion in the shuttle kiln.

shuttlekiln;oxygen-enrichedcombustion;oxygenconcentration;heat-transfercharacteristics;temperaturedifference

on Mar.26,2012

T Q 1 7 4.4+3

A

1000-2278（2012）03-0315-05

2012-03-26

江西省教育廳青年科學(xué)基金項(xiàng)目（編號(hào)：GJJ11201）

朱慶霞，E-mail:qingxia-zhu@163.com

ZHU Qingxia,E-mail:qingxia-zhu@163.com