王建平, 田萬一
( 中冶長天國際工程有限責(zé)任公司, 湖南 長沙 410007)
摘 要: 通過對(duì)環(huán)冷機(jī)上燒結(jié)礦冷卻過程進(jìn)行數(shù)值仿真分析, 以提高余熱綜合利用率為目標(biāo), 研究了冷卻過程運(yùn)行參數(shù)對(duì)燒結(jié)礦余熱利用的影響, 并通過正交實(shí)驗(yàn)找到了燒結(jié)礦余熱利用率最優(yōu)的冷卻過程運(yùn)行參數(shù)組合。
關(guān)鍵詞: 余熱利用; 數(shù)值仿真; 環(huán)冷機(jī)
1 前 言
鋼鐵工業(yè)是高能耗工業(yè), 鋼鐵總能耗占全國工業(yè)總能耗的 15% 左右, 而鋼鐵生產(chǎn)過程中的能源有效率僅為 30%[1, 2]。當(dāng)前我國燒結(jié)生產(chǎn)中余熱利用效率與國際先進(jìn)水平相比約有20% 以上的差距。因此, 我國鋼鐵行業(yè)燒結(jié)系統(tǒng)節(jié)能的潛力非常大。由于燒結(jié)礦冷卻廢氣的回收利用相較于燒結(jié)煙氣的回收利用在技術(shù)和流程上較為方便和可行, 因此, 燒結(jié)礦冷卻過程中的余熱回收利用成為了當(dāng)前降低燒結(jié)工序能耗、提高余熱利用率的熱點(diǎn)。本文通過建立燒結(jié)礦冷卻過程的傳熱計(jì)算模型, 對(duì)環(huán)冷機(jī)上冷卻過程進(jìn)行數(shù)值模擬, 并與現(xiàn)場測試相結(jié)合, 對(duì)影響燒結(jié)礦冷卻過程以及余熱利用的因素進(jìn)行了系統(tǒng)研究, 得到單一運(yùn)行參數(shù)對(duì)環(huán)冷機(jī)余熱回收的影響趨勢, 并通過正交實(shí)驗(yàn)的方法, 得到了環(huán)冷機(jī)余熱利用最優(yōu)運(yùn)行參數(shù)組合。
2 研究方法
采用數(shù)值模擬與現(xiàn)場實(shí)測相結(jié)合的方法進(jìn)行, 通過現(xiàn)場實(shí)測的結(jié)果來驗(yàn)證數(shù)學(xué)計(jì)算模型的合理性。
2、1 數(shù)值模擬方法
根據(jù)環(huán)冷機(jī)單個(gè)臺(tái)車的實(shí)際幾何尺寸建立物理模型, 其尺寸為實(shí)測值。本文對(duì)計(jì)算模型進(jìn)行網(wǎng)格劃分時(shí), 首先將環(huán)冷機(jī)臺(tái)車整體進(jìn)行結(jié)構(gòu)化網(wǎng)格劃分, 然后在環(huán)冷機(jī)臺(tái)車入口等重要位置局部加密, 這樣可以同時(shí)保證計(jì)算精度和計(jì)算量的有機(jī)結(jié)合。環(huán)冷機(jī)臺(tái)車網(wǎng)格系統(tǒng)如圖 1。
采用多孔介質(zhì)模型來模擬燒結(jié)礦的冷卻過程, 以 Ergun 方程計(jì)算氣流所受到的料層阻力。由于 FLUENT 中的能量方程為熱力學(xué)平衡狀態(tài)下的單能量方程, 燒結(jié)礦在冷卻過程中與氣流存在著換熱, 溫度是不等的, 即燒結(jié)礦與氣體之間是局部非熱力學(xué)平衡狀態(tài), 因此需要對(duì)能量方程進(jìn)行一定的修正。本文采用局部非熱平衡雙能量方程模型[ 3, 4], 即分別對(duì)燒結(jié)礦和冷卻空氣建立能量方程[ 5]如下:
以上兩式中, hv 為燒結(jié)礦與冷卻空氣之間的傳熱系數(shù)。其余參數(shù)的選取通過測試及在測試基礎(chǔ)上的計(jì)算得到, 分別是: Qs= 1 600 kg/m3, <= 014, hv= 01199 w/ ( m3#K) , cs= 920 J/( kg# e ) 。
根據(jù)燒結(jié)礦的冷卻過程, 研究對(duì)象的邊界條件為: 環(huán)冷機(jī)臺(tái)車?yán)鋮s空氣入口邊界、環(huán)冷機(jī)臺(tái)車壁面、環(huán)冷機(jī)臺(tái)車頂部出口。對(duì)于環(huán)冷機(jī)臺(tái)車?yán)鋮s空氣入口邊界, 由于單位時(shí)間進(jìn)入環(huán)冷機(jī)臺(tái)車的冷卻空氣的速度是固定的, 所以采用速度入口邊界條件, 數(shù)值以測試值為準(zhǔn)。環(huán)冷機(jī)臺(tái)車壁面采用定熱流壁面。對(duì)于環(huán)冷機(jī)臺(tái)車頂部出口, 采用壓力出口邊界條件。
2、2 現(xiàn)場測試方法
結(jié)合現(xiàn)場生產(chǎn)條件, 利用熱電偶測量環(huán)冷機(jī)余熱利用區(qū)出口煙氣溫度, 通過對(duì)測量位置選取不同的測點(diǎn), 并對(duì)不同的測點(diǎn)進(jìn)行多次測量, 取其平均值。
3 各因素對(duì)余熱利用的影響
在前述建立的物理數(shù)學(xué)模型的基礎(chǔ)上, 利用流體動(dòng)力學(xué)計(jì)算軟件 FLUENT 對(duì)燒結(jié)礦冷卻過程進(jìn)行數(shù)值模擬計(jì)算, 在保持其他參數(shù)不變的條件下, 考察單一因素對(duì)環(huán)冷機(jī)余熱利用的影響。以下分別對(duì)影響環(huán)冷機(jī)余熱回收利用量的五個(gè)主要因素進(jìn)行仿真計(jì)算: 余熱利用區(qū)內(nèi)環(huán)冷機(jī)臺(tái)車入口風(fēng)速、燒結(jié)礦料層高度、燒結(jié)礦粒徑、孔隙率、入口冷卻空氣溫度。
3、1 環(huán)冷機(jī)臺(tái)車入口風(fēng)速對(duì)余熱利用的影響
不同環(huán)冷機(jī)臺(tái)車入口風(fēng)速條件下, 環(huán)冷機(jī)出口煙氣平均溫度隨時(shí)間變化曲線如圖2 所示。
由圖 2 可知, 隨入口風(fēng)速的增加, 出口煙氣平均溫度逐漸降低。說明入口風(fēng)速越大, 單位時(shí)間內(nèi)氣體帶走的熱量越多, 有效余熱利用量增加。
3、2 燒結(jié)礦料層高度
不同燒結(jié)礦料層高度下, 環(huán)冷機(jī)出口煙氣平均溫度隨時(shí)間的變化曲線如圖 3 所示。
由圖 3 可知, 環(huán)冷機(jī)出口煙氣的平均溫度隨著燒結(jié)礦層高度的增加而升高。由此可知,環(huán)冷機(jī)出口煙氣的平均溫度和有效余熱利用量隨著燒結(jié)礦層高度的增加而增大。因此, 在考慮鼓風(fēng)機(jī)動(dòng)力消耗允許的情況下, 應(yīng)增加料層高度, 以提高余熱回收率。
3、3 燒結(jié)礦粒徑
在燒結(jié)礦粒徑不同的條件下, 環(huán)冷機(jī)出口煙氣平均溫度隨時(shí)間的變化曲線如圖 4。
由圖 4 可知, 燒結(jié)礦粒徑的增大將會(huì)導(dǎo)致環(huán)冷機(jī)臺(tái)車出口煙氣的平均溫度降低。由此可知, 環(huán)冷機(jī)出口煙氣平均溫度和有效余熱利用量隨著燒結(jié)礦粒徑的增大而減小。
3、4 孔隙率
不同孔隙率的情況下, 環(huán)冷機(jī)出口煙氣平均溫度隨時(shí)間變化曲線如圖 5。
由圖 5 可知, 在冷卻過程初期階段(200 s 之前), 隨著孔隙率的增加, 出口煙氣的平均溫度增高; 在 200 s 之后, 則趨勢相反。
3、5 入口冷卻空氣溫度
在入口冷卻空氣溫度不同的條件下, 環(huán)冷機(jī)出口煙氣平均溫度隨時(shí)間的變化曲線如圖 6。
由圖 6 可知, 入口冷卻空氣溫度升高, 使得出口煙氣的平均溫度升高。說明隨著入口冷卻空氣溫度的增加, 出口煙氣平均溫度和有效余熱利用量也隨之增加。
4 數(shù)值仿真結(jié)果的驗(yàn)證
在常規(guī)工況下, 利用 origin 軟件對(duì)數(shù)值模擬計(jì)算所得到的結(jié)果進(jìn)行曲線擬合, 圖 7 為環(huán)冷機(jī)出口煙氣平均溫度隨燒結(jié)礦冷卻時(shí)間的變化曲線。如圖所示, 擬合曲線與實(shí)測曲線能夠理想重合。由此說明, 在工程實(shí)際中, 可以采用擬合曲線的數(shù)學(xué)公式對(duì)不同冷卻時(shí)間的出口煙氣溫度進(jìn)行預(yù)測。
5 正交優(yōu)化仿真實(shí)驗(yàn)
利用正交實(shí)驗(yàn)的方法, 對(duì)影響環(huán)冷機(jī)上燒結(jié)礦余熱利用的七個(gè)主要運(yùn)行參數(shù)進(jìn)行優(yōu)化分析。七個(gè)運(yùn)行參數(shù)分別是: 環(huán)冷機(jī)臺(tái)車上下層物料高度、上中下層的燒結(jié)礦粒徑、環(huán)冷機(jī)臺(tái)車入口風(fēng)速、冷卻空氣入口平均溫度。七個(gè)因素之間相互獨(dú)立, 本實(shí)驗(yàn)采用標(biāo)準(zhǔn)正交表 L18( 37) 。優(yōu)化指標(biāo)為燒結(jié)礦每小時(shí)的余熱利用量;根據(jù)燒結(jié)礦進(jìn)入環(huán)冷機(jī)的實(shí)際溫度, 將燒結(jié)礦的初始溫度設(shè)定為 1 023115 K( 750 e ) 。表 1列出了正交實(shí)驗(yàn)各工況仿真計(jì)算結(jié)果。
通過正交實(shí)驗(yàn)結(jié)果可得到, 七個(gè)不同運(yùn)行參數(shù)按其對(duì)燒結(jié)礦余熱利用的影響, 從強(qiáng)到弱的排序依次為: 環(huán)冷機(jī)臺(tái)車入口風(fēng)速、環(huán)冷機(jī)臺(tái)車上層燒結(jié)礦粒徑、環(huán)冷機(jī)臺(tái)車中層燒結(jié)礦粒徑、環(huán)冷機(jī)臺(tái)車下層燒結(jié)礦粒徑、冷卻空氣入口平均溫度、環(huán)冷機(jī)臺(tái)車下層物料高度、環(huán)冷機(jī)臺(tái)車上層物料高度。
本實(shí)驗(yàn)結(jié)果表明, 實(shí)驗(yàn)中的七個(gè)不同參數(shù)的最優(yōu)組合為: 環(huán)冷機(jī)臺(tái)車下層物料高度為 014m, 環(huán)冷機(jī)臺(tái)車下層物料粒徑大小 01035 m; 環(huán)冷機(jī)臺(tái)車中層物料高度為 015 m, 環(huán)冷機(jī)臺(tái)車中層物料粒徑大小 01025 m; 環(huán)冷機(jī)臺(tái)車上層物料高度為 015 m, 環(huán)冷機(jī)臺(tái)車上層物料粒徑大小0103 m; 環(huán)冷機(jī)臺(tái)車入口風(fēng)速為 7165 m/ s; 冷卻空氣入口平均溫度為 404 K。在此最優(yōu)參數(shù)組合的工況下, 模擬仿真計(jì)算得到燒結(jié)礦余熱利用量為 176 783 265 kJ/ h, 與正常工況條件下余熱回收量相比, 提高了 2516%, 環(huán)冷機(jī)出口煙氣平均溫度為 660186 K, 比正常工況條件下環(huán)冷機(jī)出口煙氣平均溫度提高了 49185 K。
6 結(jié) 論
1) 利用 FLU ENT 軟件作為計(jì)算平臺(tái), 采用局部非熱力學(xué)平衡的雙能量方程模型, 建立了燒結(jié)礦冷卻仿真模型。
2) 利用數(shù)學(xué)物理模型進(jìn)行模擬仿真, 在保持其它因素不變的情況下, 研究單一因素對(duì)環(huán)冷機(jī)出口煙氣溫度和環(huán)冷機(jī)余熱利用量的影響。
3) 采用 origin 軟件將利用 FLU ENT 軟件得到的數(shù)值仿真計(jì)算結(jié)果和在現(xiàn)場實(shí)測數(shù)據(jù)分別進(jìn)行曲線擬合, 得到了環(huán)冷機(jī)出口煙氣平均溫度隨冷卻時(shí)間變化的曲線。經(jīng)過對(duì)比, 兩條曲線重合理想, 所以在工程實(shí)際中可以采用擬合曲線的數(shù)學(xué)公式對(duì)不同冷卻時(shí)間的環(huán)冷機(jī)出口煙氣溫度進(jìn)行預(yù)測。
4) 利用正交試驗(yàn)方法, 對(duì)影響燒結(jié)礦冷卻過程的七個(gè)運(yùn)行參數(shù)進(jìn)行正交試驗(yàn)分析, 找出了燒結(jié)礦余熱利用率最大的參數(shù)組合, 從而達(dá)到提高燒結(jié)礦余熱利用率的目標(biāo)。