李景冠，張林進，陳川輝，葉旭初

（南京工業(yè)大學材料化學工程國家重點實驗室，江蘇南京210009）

鉻鐵礦回轉窯氧體積分數(shù)分布的模擬仿真

李景冠，張林進，陳川輝，葉旭初

（南京工業(yè)大學材料化學工程國家重點實驗室，江蘇南京210009）

摘 要：針對鉻鐵礦氧化焙燒過程中氧對窯內(nèi)溫度分布的影響和對提高鉻氧化率的重要作用，利用數(shù)值仿真（CFD）方法研究了回轉窯內(nèi)氧體積分數(shù)的分布情況。主要分析了空氣焙燒過程中回轉窯內(nèi)氧體積分數(shù)分布的一般規(guī)律，討論了空氣過剩系數(shù)和一、二次風中氧體積分數(shù)的變化對焙燒溫度、氧體積分數(shù)的影響。結果表明，空氣過剩系數(shù)為1.7時，窯內(nèi)壁附近氧體積分數(shù)沿窯長方向不斷減小，在高溫氧化帶之后窯內(nèi)氧分布趨于均勻，窯尾煙氣中氧體積分數(shù)約為7.5%；空氣過剩系數(shù)為1.7左右時較為合理；隨著一、二次風中氧體積分數(shù)的增加，窯內(nèi)溫度先是明顯升高，后略有降低。

關鍵詞：鉻鐵礦；回轉窯；模擬仿真

回轉窯是鉻鐵礦氧化焙燒工藝的關鍵設備。窯內(nèi)氧體積分數(shù)是影響鉻鐵礦氧化的重要因素［1］。氧既可保證燃料充分燃燒，又有利于鉻鐵礦氧化。盡管鉻鐵礦氧化需要的氧含量不多，但氧體積分數(shù)在整個焙燒過程中對鉻氧化率的影響較大［2］。通過對回轉窯內(nèi)氧體積分數(shù)分布規(guī)律的研究，可以探索氧對焙燒及鉻鐵礦氧化的影響，進一步提高對鉻鐵礦焙燒工藝的認識。目前國內(nèi)外有許多關于回轉窯的研究［3-5］，但針對鉻鐵礦焙燒回轉窯的研究較少。筆者通過對鉻鐵礦回轉窯氧體積分數(shù)分布的模擬仿真，旨在為鉻鐵礦焙燒工藝的發(fā)展提供參考。

1 數(shù)學模型的建立

實驗中模擬的回轉窯以天然氣作為燃料。窯長為32m，內(nèi)半徑為0.95m，燃燒器的主要操作參數(shù)如表1所示。

表1 燃燒器操作參數(shù)

2 數(shù)學模型

選擇標準k-ε模型對回轉窯內(nèi)的氣體流動、傳熱傳質(zhì)和燃燒過程進行模擬。燃燒反應采用非預混燃燒模型［6］。采用離散坐標P1輻射模型來描述窯內(nèi)的氣體與窯壁之間的輻射傳熱過程［7］。因為研究的主要對象為窯內(nèi)煙氣，在對結果影響不大的情況下，數(shù)值模擬過程中忽略了回轉窯的轉動過程，以及生料在回轉窯中的運動和存在方式。文獻［8］證明，鉻鐵礦生料在窯內(nèi)化學反應過程中，放熱和吸熱反應的總熱效應近似為零。將生料因溫度升高所需的吸熱、生料反應吸收的O2和放出的CO2以面積平均的方式在窯內(nèi)壁的源項中計算。窯內(nèi)煙氣中O2及CO2等組分的質(zhì)量守恒方程如式1所示［11］：

式中：cs為組分S的體積濃度，ρ為質(zhì)量濃度，Ds為擴散系數(shù)，Ss為系統(tǒng)內(nèi)單位時間內(nèi)單位體積通過化學反應產(chǎn)生的該組分的質(zhì)量。

3 結果及討論

3.1 溫度場

空氣過剩系數(shù)為1.7時，中心截面（X=0）上回轉窯內(nèi)溫度場分布如圖1所示。根據(jù)窯內(nèi)各區(qū)域?qū)ξ锪系牟煌饔?，按溫度可大致分為冷卻帶、高溫氧化帶、碳酸鹽分解帶和干燥預熱帶［10］，其中高溫氧化帶是鉻鐵礦氧化的主要區(qū)域。從圖1可見，窯內(nèi)煙氣溫度分布比較規(guī)整，從左到右呈逐級遞減的趨勢。

圖1 縱截面X=0處空氣焙燒時窯內(nèi)溫度分布

3.2 氧體積分數(shù)

鉻鐵礦的氧化反應主要發(fā)生在窯壁處，窯壁處的氧體積分數(shù)的分布更為關鍵?？諝膺^剩系數(shù)為1.7時，中心截面上各特征線上的氧體積分數(shù)曲線見圖2。由圖2可知，在窯頭氣體入口處（Z=0～4m），因為燃料還沒有完全燃燒，窯壁處（Y=0.94m）氧體積分數(shù)較大；在高溫氧化帶（Z=8～20m），窯壁處氧體積分數(shù)減小比較明顯，主要是受到燃料燃燒的影響，同時生料反應也會產(chǎn)生大量CO2并消耗少量O2，使得窯壁處氧體積分數(shù)從窯頭至窯尾不斷減??；在碳酸鹽分解帶和干燥預熱帶處（Z=20～32m），窯壁處氧體積分數(shù)變化較?。桓G內(nèi)部其他3條特征線上的氧體積分數(shù)，因受到燃料燃燒、氧擴散、煙氣速度等因素的影響不同而變化各異，在24m處氧體積分數(shù)趨于穩(wěn)定，最終窯內(nèi)部煙氣中的氧體積分數(shù)約為7.5%，與空氣過剩系數(shù)為1.7時理論計算結論一致［10］。

圖2 各特征線上的氧體積分數(shù)曲線

3.3 空氣過剩系數(shù)

模擬計算過程中，通過調(diào)整二次風的流量使空氣過剩系數(shù)（a）分別為1.3、1.7、2.3、3。圖3為空氣過剩系數(shù)對窯壁處特征線上的氧體積分數(shù)的影響。從圖3可以看出，空氣過剩系數(shù)的提高使窯壁處的氧體積分數(shù)顯著增大。圖4為空氣過剩系數(shù)對窯壁處溫度分布的影響。由圖4可以看出，隨著空氣過剩系數(shù)的提高，窯壁處的溫度依次明顯降低。

圖3 空氣過剩系數(shù)對窯壁處特征線上的氧體積分數(shù)的影響

圖4 空氣過剩系數(shù)對 窯壁處溫度分布的影響

結合圖3、4可知，空氣過剩系數(shù)為3時，窯壁處氧體積分數(shù)較大，但溫度較低，達不到理想的焙燒溫度，而且較大空氣過剩系數(shù)也會導致燃料熱效率低；空氣過剩系數(shù)為1.3時，窯內(nèi)氧體積分數(shù)低，且溫度過高會使液相增多，將阻礙氧向物料的擴散，不利于氧化反應。因而合理的空氣過剩系數(shù)應為1.7左右，與現(xiàn)行的鉻鐵礦焙燒工藝要求相符［10］。

3.4 風中氧含量

控制一、二次風中的氧體積分數(shù)分別為21%、30%、40%、74%。窯內(nèi)壁處特征線上的氧體積分數(shù)分布如圖5所示。從圖5可以看出，風中的氧體積分數(shù)越大，窯內(nèi)氧體積分數(shù)越高。風中的氧體積分數(shù)不同時，窯內(nèi)各氧體積分數(shù)曲線的形狀非常相似。這是因為燃料及鉻鐵礦耗氧量不變的情況下，窯內(nèi)氧體積分數(shù)的變化相同。

圖5 風中氧體積分數(shù)對窯壁處特征線上氧體積分數(shù)的影響

圖6為風中氧體積分數(shù)對窯壁處特征線上的溫度分布的影響。從圖6中可以看出，風中氧體積分數(shù)為21%、30%、40%時，煙氣高溫點隨氧體積分數(shù)的提高而升高，且氧體積分數(shù)在21%～30%時對窯內(nèi)煙氣溫度的影響較大；風中氧體積分數(shù)為30%、40%、74%時，煙氣高溫點較高，高溫點靠前，溫度分布相似，這是因為氧體積分數(shù)較大時，火焰燃燒較為集中，而且足夠的氧可以使火焰向四周拓展；風中氧體積分數(shù)大于40%時，因為高溫時（1 000～1 200℃）O2比N2的定容比熱容高約5%［11］，因而此時煙氣吸熱較多，所以φ（風氧）=74%時煙氣高溫點和溫度比φ（風氧）=40%時低。

圖6 風中氧體積分數(shù)對窯壁處特征線上的溫度分布的影響

4 結論

1）窯內(nèi)煙氣溫度分布比較規(guī)整，在窯長方向呈逐級遞減的趨勢，其中高于1 100℃的區(qū)域較長，高溫氧化帶位于窯內(nèi)8～20m處，窯尾出口溫度約為700℃。

2）在回轉窯內(nèi)，氧由窯壁處向內(nèi)部進行擴散，在20m以后，窯內(nèi)氧體積分數(shù)趨于穩(wěn)定。當空氣過剩系數(shù)為1.7時，窯尾煙氣中的氧體積分數(shù)約為7.5%。

3）較大的空氣過剩系數(shù)可以提高窯內(nèi)氧體積分數(shù)，但降低了燃料熱效率?？諝獗簾龝r需要較大的空氣過剩系數(shù)來保證窯內(nèi)有足夠的氧進行氧化反應，合理的空氣剩系數(shù)應為1.7左右。

4）一、二次風中氧體積分數(shù)在21%～40%時，風中氧體積分數(shù)的提高使窯內(nèi)煙氣高溫點升高，且在φ（風氧）=21%～30%時對窯內(nèi)溫度分布影響較大，φ（風氧）＞40%后窯內(nèi)煙氣溫度略有降低。φ（風氧）= 30%～74%時，窯內(nèi)溫度分布相似。

參考文獻：

［1］ Antony M P，Jha A，Tathavadkar V D.Alkali roasting of Indian chromiteores：thermodynamicand kinetic considerations［J］.Mineral Processingand ExtractiveMetallurgy，2006，115（2）：71-79.

［2］ Tathavadkar V D，Antony M P，Jha A.The effect of Salt-phase composition on the rate of soda-ash roasting of chromite ores［J］. Metallurgical and Materials Transactions B，2003，34（5）：555-563.

［3］ 張志霄，池涌，李水清.回轉窯傳熱模型與數(shù)值模擬［J］.化學工程，2003，31（4）：27-31.

［4］ 王匯，羅申，溫治.回轉窯內(nèi)流動傳熱和燃燒過程數(shù)學模型的研究現(xiàn)狀及其發(fā)展趨勢［J］.礦冶，2006，15（2）：28-34.

［5］ 馬愛純，周孑民，李旺興.回轉窯內(nèi)煙氣溫度分布的數(shù)值研究［J］.金屬材料與冶金工程，2007，35（3）：19-22，31.

［6］ 周密，閆立峰，郭慶祥，等.非預混燃燒模型模擬流化床生物質(zhì)氣化器中富氫氣體的制備［J］.化學物理學報，2006，19（2）：131-136.

［7］ MujumdarKS，RanadeVV.CFDmodelingofrotary cementkilns［J］. Asia-pacific Journal of Chemical Engineering，2008，3（2）：106-118.

［8］ 利西安柯，紀柱.鉻礦氧化焙燒回轉窯熱工特性綜合研究［J］.鉻鹽工業(yè)，1997（2）：13-22.

［9］ 龍?zhí)煊?計算流體力學［M］.重慶：重慶大學出版社，2007：15-17.

［10］ 丁翼，紀柱.鉻化合物生產(chǎn)與應用［M］.北京：化學工業(yè)出版社，2003：44-88.

［11］ 劉光啟，馬連湘，劉杰.化學化工物性數(shù)據(jù)手冊無機卷［M］.北京：化學工業(yè)出版社，2002：110-112.

聯(lián)系方式：li0676@163.com

中圖分類號：TQ136.11

文獻標識碼：A

文章編號：1006-4990(2012)02-0033-03

收稿日期：2011-08-20

作者簡介：李景冠（1979— ），男，碩士，主要從事無機材料生產(chǎn)工藝的研究，已公開發(fā)表文章1篇。

Numericalsimulation ofoxygen volume fraction in chrom ite roasting process

Li Jingguan，Zhang Linjin，Chen Chuanhui，Ye Xuchu
（State Key Laboratory ofMaterials-Oriented ChemicalEngineering，NanjingUniversity of Technology，Nanjing210009，China）

Abstract：Oxygen volume fraction distribution in rotary kiln was studied by numerical simulation（CFD）technology for the importance ofoxygen on temperature distribution and increasing chromite oxidation rate during roasting process.The general law of oxygen volume fraction distribution in the air roasting processwas analyzed，then the effectof excess-air coefficient，different oxygen volume fractions of inlet wind on temperature distribution and oxygen volume fraction in rotary kiln was discussed.Result indicated that，when excess-air coefficientwas 1.7，oxygen volume fraction near wall continually reduced lengthwise along kiln；in this condition oxygen volume fraction of kiln tail gaswas about 7.5%and oxygen in rotary kiln tended to distribute evenly after high temperature oxidation area；the proper excess-air coefficient was about 1.7；with increasing oxygen volume fraction of inlet wind，kiln temperature would obviously rise in the beginning，then decreased slightly.

Keywords：chromite；rotary kiln；numericalsimulation