燒結礦余熱回收豎罐內氣固傳熱特性

馮軍勝，董輝，劉靖宇，梁凱

（東北大學材料與冶金學院，國家環境保護生態工業重點實驗室，遼寧 沈陽 110819）

燒結礦余熱回收豎罐內氣固傳熱特性

馮軍勝，董輝，劉靖宇，梁凱

（東北大學材料與冶金學院，國家環境保護生態工業重點實驗室，遼寧 沈陽 110819）

以自制氣固傳熱實驗裝置為操作平臺，實驗研究了燒結礦顆粒填充床內的氣固傳熱特性。結果表明：氣體表觀流速和燒結礦顆粒直徑是影響顆粒床層內氣固傳熱過程的主要因素。氣體表觀流速越大，顆粒直徑越小，床層內氣固傳熱系數就越大。當燒結礦顆粒直徑和氣體表觀流速一定時，床層內燒結礦顆粒溫度的升高導致床層氣固傳熱系數的增加和傳熱Nusselt數的減小。由于計算誤差較大，現有的經驗關聯式不適用于求解燒結礦顆粒床層內的氣固傳熱過程。基于量綱分析法，并結合實驗測量數據擬合得出了能夠描述燒結礦顆粒床層內氣固傳熱特性的實驗關聯式，平均計算誤差為4.22%，顯示了良好的預測性能。

燒結；填充床；豎罐；多孔介質；傳熱系數

引　言

燒結過程余熱資源的高效回收利用是目前降低燒結工序能耗的主要途徑之一[1]。燒結礦余熱罐式回收是針對現有燒結余熱回收系統存在難以克服的弊端[2]，借鑒干熄焦爐的結構和工藝提出的一種燒結余熱高效回收方式[3-6]。就床層特點而言，燒結礦余熱回收豎罐是一種隨機非結構化大顆粒填充床，床層內燒結礦的填充結構近似于多孔介質。豎罐床層內的氣固傳熱特性是決定燒結礦豎罐式余熱回收可行性的關鍵問題。其中，床層氣固傳熱系數是影響豎罐內燒結礦和冷卻介質熱量傳遞，進而影響燒結礦豎罐結構和操作參數的主要因素之一。

由于顆粒填充床在冶金、化工和農業等領域應用較廣，其床層內氣固傳熱問題也被廣泛地研究。迄今，國內外學者已經對顆粒填充床內氣固傳熱特性進行了大量的研究。Wakao等[7]對前人的實驗結果進行了總結，得出了Reynolds數在20～10000范圍內適用于空氣和水的多孔介質換熱關系式。Whitaker[8]在前人實驗的基礎上分別推導出Reynolds數在 22～8000范圍內球體填充床內的換熱準則關系式。Kays等[9]給出了Reynolds數在10～5000范圍內，以空氣為換熱介質的球體多孔介質換熱關系式。Nie等[10]研究了填充顆粒與熱空氣之間的換熱，并給出Reynolds數在5～280范圍內以空氣為介質的傳熱關系式。Kuwahara等[11]采用實驗方法得收了雙能量方程中的對流換熱Nusselt數，這一對流換熱公式是由空隙率、Prandtl數和 Reynolds數組成。Abdulmohsin等[12]實驗測量了球體填充床內的對流傳熱系數及其徑向分布，分析了氣體速度對對流傳熱系數的影響，并得收了傳熱Nusselt數沿著床層徑向和軸向分布的預測公式。Yang等[13]實驗研究了非球體顆粒填充床內的流動和傳熱特性，并根據實驗數據收正了床層壓降公式和換熱 Nusselt數。Thakur等[14]采用實驗方法，收正了低空隙率填充床內傳熱系數和顆粒摩擦因子，認為二者均是床層幾何參數的函數。Pe?i?等[15]實驗研究了大的冷球體放入熱球體顆粒填充床內的傳熱規律，并根據實驗測量數據得收了顆粒床層內的傳熱系數。Cong 等[16]實驗研究了恒壁面熱流填充床內氣固兩相流之間的換熱特性，并采用對數平均溫差法獲得了床層內的綜合傳熱系數。

綜上所述，目前關于顆粒填充床內氣固傳熱特性的研究較多，但是這些研究大部分集中在球形顆粒或者均勻顆粒床層內，涉及大顆粒非均勻顆粒床層內氣固傳熱特性的研究較少，而且這些有限的研究中所提供的關聯式不適用于求解燒結礦豎罐內的氣固傳熱過程。基于此，本文在自制氣固傳熱實驗裝置上，實驗研究了影響豎罐內氣固傳熱系數的主要因素及其影響規律，并在此基礎上，通過量綱分析法并結合實驗測量數據擬合得出描述燒結礦顆粒床層內氣固傳熱特性的實驗關聯式，旨在為燒結礦余熱豎罐的設計和優化提供理論依據。

1　實驗方法及裝置

1.1實驗方法

在 dτ時間內，氣體在固定顆粒填充床微元段dV內的氣固傳熱量可表示為

其中

在dτ時間內，微元段內燒結礦和冷卻空氣的溫度變化分別為：

對于燒結礦

對于冷卻空氣

式（3）與式（4）相減可得

將式（1）代入式（4）整理得

兩邊積分得

由此得出氣固綜合傳熱系數求解公式如下

Nusselt數的求解公式如下

式（9）中Ts為豎罐床層內燒結礦的平均溫度；Tg為豎罐進出口空氣溫度的平均值。

1.2實驗裝置

氣固傳熱實驗裝置如圖1所示。實驗豎罐的橫截面為圓形，其內徑和高度分別為450 mm和1000mm。冷卻空氣在鼓風機的驅動下，經過調節閥和孔板流量計后進入豎罐料層內，最后從豎罐頂部流出。實驗所有工況的料層高度為 800 mm，在豎罐本體外壁高度200、400、600 mm和空氣出口處各有一個測溫孔，通過K-type (Ni/Al)型熱電偶對不同高度處的燒結礦溫度進行檢測，通過抽氣式熱電偶對空氣出口的溫度進行檢測。實驗通過調節閥來控制進入豎罐內的空氣流量，通過孔板流量計來測量空氣流量。所有實驗工況均在常溫下進行，即豎罐的進口風溫為20℃。實驗所用篩分后的燒結礦顆粒平均直徑、顆粒形狀因子[17]、床層空隙率[18]和顆粒當量直徑的變化范圍如表1所示。其中，顆粒形狀因子，即顆粒球形度，是采用等體積的測算方法求出的。因此，計算所得的顆粒當量直徑為等體積當量直徑。

圖1　氣固傳熱實驗裝置Fig. 1　Schematic diagram of gas-solid heat transfer experimental setup

表1　燒結礦顆粒相關參數Table 1　Related parameters of sinter particles

影響燒結礦豎罐內氣固傳熱過程的因素可分為2個方面：床層因素，如燒結礦顆粒直徑dp、床層孔隙率ε（通過顆粒直徑dp反映）等；流動介質因素，如氣體表觀流速u、氣體黏度μ等。通過對式（2）進行分析可知，顆粒直徑dp是影響床層空隙率ε和顆粒比表面積av的主要因素，而氣體表觀流速是影響燒結礦平均溫度和空氣出口溫度的主要因素。基于此，本文擬討論氣體表觀流速和燒結礦顆粒直徑對罐體內氣固傳熱過程的影響。

實驗中，分別對表1中3種不同顆粒直徑燒結礦在不同氣體表觀流速條件下的實驗工況進行測試，測量出不同實驗工況下的燒結礦溫度和出口氣體溫度，國后取其平均值，采用式（8）計算出不同實驗工況下床層內的氣固傳熱系數。本文以顆粒當量直徑dp作為特征長度，以空氣進出口溫度的平均值Tg作為定性溫度，文中Rep定義為

2　實驗結果與分析

通過測量顆粒填充床層內燒結礦的平均溫度和空氣出口溫度，得出了顆粒直徑d=18 mm條件下床層內氣固傳熱系數隨氣體表觀流速變化關系，如圖2所示。

圖2　氣體表觀流速對氣固傳熱系數的影響Fig.2　Effect of gas superficial velocity on gas-solid heat transfer coefficient

由圖2可知，燒結礦顆粒溫度越高，氣體表觀流速越大，床層內燒結礦與冷卻空氣之間的綜合傳熱系數就越大。這是由于，隨著燒結礦溫度的升高，燒結礦與冷卻空氣間的輻射傳熱會大幅增加，導致與燒結礦顆粒進行傳熱的冷卻空氣溫度的升高。空氣溫度的升高導致空氣熱導率的增加，從而導致燒結礦顆粒與空氣之間傳導傳熱的增加，同時空氣溫度的升高也會導致輻射傳熱的增加。因此，床層內燒結礦顆粒與冷卻空氣間的傳熱系數會隨著顆粒溫度的升高而增大。

當燒結礦溫度一定時，氣體表觀流速的增加，導致床層內空氣的熱運動加劇，冷卻空氣與燒結礦顆粒之間的對流傳熱增加，從而導致燒結礦顆粒與空氣之間傳熱系數的增加。

圖3為氣體表觀流速u=1.04 m·s-1條件下床層內氣固傳熱系數隨燒結礦顆粒直徑的變化。從圖3中可以看出，燒結礦顆粒溫度越高，顆粒直徑越小，床層內燒結礦與冷卻空氣之間的綜合傳熱系數就越大。這是由于顆粒直徑越小，床層內燒結礦顆粒的比表面積就越大，導致燒結礦顆粒與空氣之間的對流傳熱增加。同時，顆粒直徑越小，燒結礦顆粒內部的導熱熱阻也就越小，顆粒與空氣間的傳導傳熱就會增加。因此，床層內燒結礦顆粒與冷卻空氣間的傳熱系數會隨著顆粒直徑的減小而增大。

圖3　燒結礦顆粒直徑對氣固傳熱系數的影響Fig.3　Effect of sinter particle diameter on gas-solid heat transfer coefficient

在以往球體顆粒填充床氣固傳熱研究中，Gupta等[19]、Ranz等[20]和Chaube等[21]分別提出了描述顆粒床層內氣固傳熱特性的預測關聯式，如下所示。

預測關聯式的計算平均誤差η采用式（14）進行計算

由式（9）與式（11）～式（13）分別得收了氣體表觀流速u=1.04 m·s-1條件下不同顆粒直徑填充床內傳熱Nusselt數隨顆粒Reynolds數Rep的變化關系，實驗值與經驗關聯式計算值的對比如圖4所示。

圖4　不同顆粒床層內Nu-Rep關系實驗值與計算值的對比Fig.4　Comparison between experimental values and calculation values with Nu vs Repin different particle bed layer

從圖 4可以得出，通過式（9）所求得的實驗Nusselt數隨著顆粒Reynolds數Rep的增加而增加。Gupta關聯式[19]求解不同顆粒直徑填充床內氣固傳熱系數的計算值與實驗值相差較大，η值達收了234.2%，而Ranz關聯式[20]和Chaube關聯式[21]求解不同顆粒直徑填充床內氣固傳熱系數的計算值與實驗值相差較小，尤其在圖 4(c)中，Ranz關聯式[20]求解顆粒直徑d =36 mm填充床內氣固傳熱系數的計算值與實驗值之間的平均誤差僅為 6.33%，但基于整個實驗工況，Ranz關聯式[20]和Chaube關聯式[21]的η值達收了21.68%和40.53%。因此，以上3個氣固傳熱經驗關聯式不適用于求解燒結礦豎罐內的氣固傳熱過程。

基于以上分析，采用球體顆粒填充床研究所得收的經驗關聯式不適用于求解燒結礦顆粒床層內的氣固傳熱過程。因此，本文通過量綱分析法并結合實驗測量數據擬合得出描述燒結礦顆粒床層內氣固傳熱特性的實驗關聯式。

由文獻[7-16]可知，影響顆粒填充床內氣固傳熱特性的因素有：床層內氣固傳熱系數 h、床層空隙率ε、流體密度ρ、流體比熱容cp、流體熱導率λ、動力黏度μ、氣體表觀流速u和顆粒直徑dp。基于此，描述顆粒床層內流動阻力特性的一般關系式如下

采用量綱分析法，并結合π定理計算可得如下量綱1方程

式（16）是描述燒結礦顆粒床層內氣固傳熱特性的實驗關聯式，其中，k、α、β和γ為實驗常數，需要通過實驗數據擬合得收。

式（16）是多元非線性方程，將其兩邊取對數轉變為線性方程為

利用最小二乘法，將實驗數據代入式（17）中，采用Excel軟件進行線性回歸計算，擬合所得的實驗關聯式如下所示。

其中，540≤Rep≤3724，0.674≤Pr≤0.703。

圖5為擬合實驗關聯式的計算值Nucal與實驗值Nuexp之間的對比。從圖5中可以看出，采用擬合實驗關聯式計算所得床層內氣固傳熱系數與實驗值能較好地吻合，對于整個實驗工況，η值也僅為4.22%。因此，式（18）可以用來描述燒結礦顆粒床層內的氣固傳熱特性。

圖5　擬合實驗關聯式計算值與實驗值的對比Fig.5　Comparison between calculation values and experimental values

3　結　論

本文在自制氣固傳熱裝置上，實驗研究了燒結礦顆粒床層內氣固傳熱特性，所得結論如下。

（1）影響燒結礦顆粒床層內氣固傳熱特性的主要因素為氣體表觀流速和燒結礦顆粒直徑。氣體表觀流速越大，顆粒直徑越小，床層內氣固傳熱系數就越大。另外，床層內燒結礦溫度對氣固傳熱系數的影響也很大，床層內氣固傳熱系數隨燒結礦溫度的升高而增大。

（2）對于不同顆粒直徑填充床內的氣固傳熱過程，當氣體表觀流速一定時，床層內氣固傳熱Nusselt數隨顆粒Reynolds數Rep的增加而增加，隨床層內燒結礦顆粒溫度的升高而減小。

（3）由于計算誤差較大，基于球體顆粒填充床研究所得的經驗關聯式不適用于求解燒結礦顆粒床層內的氣固傳熱過程。采用量綱分析法，并結合實驗測量數據擬合得出的實驗關聯式，平均計算誤差僅為4.22%，能夠很好地描述燒結礦顆粒床層內的氣固傳熱特性。

符號說明

av——床層內顆粒比表面積，m2·m-3

c——比熱容，J·kg-1·K-1

d——顆粒平均直徑，mm

dp——顆粒當量直徑，mm

h——氣固傳熱系數，W·m-2·K-1

Nu——氣固傳熱Nusselt數

Pr——Prandtl數

Rep——顆粒Reynolds數

T——溫度，K

u——氣體表觀流速，m·s-1

ε——空隙率

η——計算誤差

λf——空氣熱導率，W·m-1·K-1

μ——運動黏度，Pa·s

ρ——密度，kg·m-3

Ф——顆粒形狀因子

下角標

cal——計算值

exp ——實驗值

g ——空氣

s ——燒結礦

References

[1] Cai Jiuju (蔡九菊), Wang Jianjun (王建馮), Chen Chunxia (陳春霞), Lu Zhongwu (陸鐘武). Waste heat recovery and utilization in iron and steel industry [J]. Iron and Steel (鋼鐵), 2007, 42 (6): 1-7.

[2] Dong Hui (董董), Zhao Yong (趙勇), Cai Jiuju (蔡九菊), Zhou Jiewang (周節旺), Ma Guangyu (馬光宇). On the air leakage problem in sintering cooling system [J]. Iron and Steel (鋼鐵), 2012, 47 (1):95-99.

[3] Marcelo R E, Luiz F M. Thermodynamic analysis of a coke dry quenching unit [J]. Energ. Convers. Manage., 2000, 41 (2): 109-127.

[4] Dong Hui, Jia Fengrui, Zhao Yong, Wang Meng, Cai Jiuju. Experimental investigation on the drying process of the sinter mixture [J]. Powder Technol., 2012, 218: 1-4.

[5] Cai Jiuju (蔡九菊), Dong Hui (董董). The method and device of sintering waste heat recovery and utilization with vertical tank [P]: CN, 200910187381.8. 2009-09-15.

[6] Dong Hui (董董), Li Lei (李磊), Liu Wenjun (劉文馮), Wang Bo (王博), Suo Yanshuai (索延帥), Cai Jiuju (蔡九菊). Process of waste heat recovery and utilization for sinter in vertical tank [J]. China Metallurgy (中國冶金), 2012, 22 (1): 6-11.

[7] Wakao N, Kaguei S. Heat and Mass Transfer in Packed Beds [M]. New York: McGraw-Hill, 1982.

[8] Whitaker S. Forced convection heat transfer correlations for flow in pipes, past flat plates, single cylinders, single spheres, and for flow in packed beds and tube bundles [J]. AIChE J., 1972, 18 (2): 361-371.

[9] Kays W M, London A L. Compact Heat Exchangers [M]. New York: McGraw-Hill, 1984.

[10] Nie Xiaodong, Evitts R, Besant R, Bolster J. A new technique to determine convection coefficients with flow through particle beds [J]. J. Heat Transfer, 2011, 133 (4): 41-60.

[11] Kuwahara F, Shirota M, Nakayama A. A numerical study of interfacial convective heat transfer coefficient in two-energy equation model for convection in porous media [J]. Int. J. Heat Mass Tran., 2001, 44 (6): 1153-1159.

[12] Abdulmohsin R S, Al-Dahha M H. Characteristics of convective heat transport in a packed pebble-bed reactor [J]. Nucl. Eng. Des., 2015, 284: 143-152.

[13] Yang Jian, Wang Jing, Bu Shanshan, Zeng Min, Wang Qiuwang, Nakayama A. Experimental analysis of forced convective heat transfer in novel structured packed beds of particles [J]. Chem. Eng. Sci., 2012, 71 (26): 126-137.

[14] Thakur N S, Saini J S, Solanki S C. Heat transfer and friction factor correlations for packed bed solar air heater for a low porosity system [J]. Sol. Energy, 2003, 74 (4): 319-329.

[15] Pe?i? R, Radoi?i? T K, Bo?kovi?-Vragolovi? N, Arsenijevi? Z, Grbav?i? ?. Heat transfer between a packed bed and a larger immersed spherical particle [J]. Int. J. Heat Mass Tran., 2014, 78: 130-136.

[16] Cong T N, He Yurong, Chen Haisheng, Ding Yulong, Wen Dongsheng. Heat transfer of gas-solid two-phase mixtures flowing through a packed bed under constant wall heat flux conditions [J]. Chem. Eng. J., 2007, 130 (1): 1-10.

[17] Geldart D. Estimation of basic particle properties for use in fluid-particle process calculations [J]. Powder Technol., 1990, 60 (1): 1-13.

[18] Klerk A D. Voidage variation in packed beds at small column to particle diameter ratio [J]. AIChE J., 2003, 49 (8): 2022-2029.

[19] Gupta A S, Thodos G. Mass and heat transfer in the flow of fluids through fixed and fluidized beds of spherical particles [J]. AIChE J., 1962, 8 (5): 608-610.

[20] Ranz W E, Marshall W R. Evaporation from drops [J]. Chem. Eng. Prog., 1952, 48: 141-146.

[21] Gupta S N, Chaube R B, Upadhyay S N. Fluid-particle heat transfer in fixed and fluidized beds [J]. Chem. Eng. Sci., 1974, 29 (3): 839-843.

Gas-solid heat transfer characteristics in vertical tank for sinter waste heat recovery

FENG Junsheng, DONG Hui, LIU Jingyu, LIANG Kai
(SEPA Key Laboratory on Eco-industry, School of Materials and Metallurgy, Northeastern University, Shenyang 110819, Liaoning, China)

With a homemade gas-solid heat transfer experimental setup, an experimental study was conducted on gas-solid heat transfer characteristics in a packed bed with sinter particles. The results show that the main factors affecting the gas-solid heat transfer process in sinter bed layer are gas superficial velocity and sinter particle diameter. With the increase of gas superficial velocity and the decrease of sinter particle diameter, the gas-solid heat transfer coefficient increases in the sinter bed. For given gas superficial velocity and sinter particle diameter, higher sinter particle temperature in the bed gives higher gas-solid heat transfer coefficient and lower Nusselt number of heat transfer. Due to the large calculation error, the existing prediction correlations are not suitable for the gas-solid heat transfer process in sinter bed. Based on the method of dimensional analysis, an experimental correlation for describing the gas-solid heat transfer characteristics in sinter bed is obtained by fitting the experimental data, and its mean deviation from experimental data is 4.22%, giving good prediction.

sintering; packed bed; vertical tank; porous media; heat transfer coefficient

date: 2015-05-04.

Prof. DONG Hui, Dongh@smm.neu.edu.cn

supported by the National Natural Science Foundation of China (51274065).

10.11949/j.issn.0438-1157.20150548

TK 11+5

A

0438—1157（2015）11—4418—06

2015-05-04收到初稿，2015-06-04收到修改稿。

聯系人：董輝。第一作者：馮軍勝(1988—)，男，博士研究生。

國家自然科學基金項目(51274065)。