基于計(jì)算流體動(dòng)力學(xué)的礦熱爐水冷型煙道模擬研究

摘 """""要： 為了研究鎳鐵冶煉礦熱爐水冷型煙道中的煙氣流動(dòng)特性，利用FLUENT對(duì)礦熱爐水冷型煙道內(nèi)煙氣沉積行為進(jìn)行模擬。采用RNG k-ε湍流模型計(jì)算了流場(chǎng)和傳熱。同時(shí)，采用DPM模型，并考慮了Grant和Tabakoff顆粒反彈模型及JKR理論，建立了顆粒沉積的判斷準(zhǔn)則，分析了爐氣顆粒在碰撞過(guò)程中的受力情況及其運(yùn)動(dòng)軌跡。

關(guān) "鍵 "詞：礦熱爐；水冷型煙道；顆粒沉積；數(shù)值模擬

中圖分類號(hào)：TF806.7 """""文獻(xiàn)標(biāo)志碼：A """"文章編號(hào)：1004-0935（2024）10-1564-03

回轉(zhuǎn)窯-礦熱爐工藝因其工藝成熟的特點(diǎn)廣泛應(yīng)用于鎳鐵的實(shí)際生產(chǎn)^[1]。礦熱爐在實(shí)際生產(chǎn)中會(huì)產(chǎn)生高溫爐氣，爐氣中含有的礦料顆粒會(huì)隨爐氣一同進(jìn)入煙道，并在流動(dòng)中與煙道壁面產(chǎn)生碰撞與摩擦，這易于造成煙道的沉積、磨蝕、腐蝕和堵塞等問(wèn)題，進(jìn)而影響煙道的使用壽命。為了確保高溫爐氣在煙道中的高效輸運(yùn)，對(duì)礦熱爐水冷型煙道內(nèi)氣固兩相流中顆粒沉積行為和去除特性進(jìn)行了模擬研究。

在目前的研究中，學(xué)者發(fā)現(xiàn)顆粒的沉積行為與氣體、顆粒及碰撞壁面的特性有較大的關(guān)系。CROSBY等^[2]研究發(fā)現(xiàn)顆粒的沉積受到氣體溫度影響較大，并且發(fā)現(xiàn)顆粒的沉積速率隨著氣體溫度的降低而降低。HAN等^[3]則利用改進(jìn)的模型研究了流速和顆粒直徑對(duì)顆粒沉積速率的影響，研究發(fā)現(xiàn)流速越大顆粒沉積速率越小，粒徑越大顆粒的沉積速率越大。BRAVO等^[4]研究了粒徑分布和運(yùn)動(dòng)軌跡對(duì)顆粒沉積的影響，發(fā)現(xiàn)粒徑分布和顆粒的初始軌跡對(duì)預(yù)測(cè)顆粒沉積起著至關(guān)重要的作用。HONG等^[5]研究了顆粒尺寸、粗糙結(jié)構(gòu)及表面能等因素對(duì)顆粒沉積特性的影響，發(fā)現(xiàn)不同粗糙結(jié)構(gòu)對(duì)顆粒沉積有較為顯著的影響，并且表面能的增大會(huì)導(dǎo)致顆粒沉積速率的增大。

以48 MW礦熱爐煙道為研究對(duì)象，用FLUENT對(duì)礦熱爐水冷型煙道內(nèi)煙氣沉積行為進(jìn)行模擬，預(yù)測(cè)了煙道內(nèi)煙氣的流動(dòng)特性。采用RNG k-ε湍流模型計(jì)算了水冷型煙道內(nèi)的流場(chǎng)和傳熱。同時(shí)，采用DPM模型考慮了Grant和Tabakoff顆粒反彈模型及JKR理論，建立了顆粒沉積的判斷準(zhǔn)則，分析了爐氣顆粒在碰撞過(guò)程中的受力情況及其運(yùn)動(dòng)軌跡。

1 "模型建立與網(wǎng)格劃分

1.1 "模型建立

礦熱爐水冷型煙道模型示意圖如圖1所示。煙道除了爐口段為漸縮段，其余各管段的直徑均為1 000 mm，Ⅰ段煙道長(zhǎng)為 5"300"mm、Ⅱ段煙道長(zhǎng)為3"700"mm、Ⅲ段煙道長(zhǎng)為10"000"mm，盲管段長(zhǎng)為1"500"mm，出口段長(zhǎng)為 8"000"mm，出口段與Ⅲ段煙道夾角為定值，且α=80°。水冷方向與煙氣流動(dòng)方向相反，各段水冷管數(shù)均為10。

1.2 "網(wǎng)格無(wú)關(guān)性驗(yàn)證

網(wǎng)格數(shù)量會(huì)影響沉積速率的大小，為了計(jì)算效率最大化，物理模型整體網(wǎng)格數(shù)量取1"429"133個(gè)。將模擬數(shù)據(jù)與HAN等的研究結(jié)果作對(duì)比，模型驗(yàn)證結(jié)果如圖2所示。

2 "數(shù)學(xué)模型

為了研究礦熱爐煙道中煙氣的流動(dòng)特性，利用歐拉-拉格朗日法來(lái)求解含有顆粒爐氣的流動(dòng)問(wèn)題。其中利用歐拉模型描述氣相流動(dòng)，利用拉格朗日模型預(yù)測(cè)顆粒的運(yùn)動(dòng)。爐氣是含有顆粒、不可壓縮的三維湍流流體，冷卻水為常溫常壓、不可壓縮流體。在模擬爐氣流動(dòng)的過(guò)程中，采用RNG k-ε湍流模型計(jì)算了流場(chǎng)和傳熱。由于連續(xù)相的質(zhì)量、動(dòng)量和能量的控制方程以及RNG k-ε湍流模型方程在大量文獻(xiàn)中存在，因此不再贅述^[6-7]。

采用DPM模型求解固相顆粒的運(yùn)動(dòng)，基于顆粒粒徑分布特性，顆粒輸運(yùn)過(guò)程中受力分析主要考慮布朗力、虛擬質(zhì)量力、Saffman力、重力、曳力和慣性力等，采用DEFINE_DPM_DRAG宏拓展曳力計(jì)算，準(zhǔn)確模擬爐氣與顆粒發(fā)生相對(duì)運(yùn)動(dòng)時(shí)產(chǎn)生的相互作用力。

顆粒碰撞示意圖如圖3所示。

為了能準(zhǔn)確預(yù)測(cè)顆粒的運(yùn)動(dòng)軌跡，提出了顆粒沉積模型。該模型考慮了顆粒與壁面和顆粒與顆粒之間的相互作用。此外還考慮了顆粒反彈模型，通過(guò)顆粒碰撞過(guò)程中能量的損耗來(lái)預(yù)測(cè)顆粒反彈速度^[8]。

3 "結(jié)果與討論

3.1 "煙氣流速對(duì)沉積的影響

煙氣的流速是影響顆粒沉積速率的一個(gè)主要因素。煙氣流速對(duì)沉積速率的影響如圖4所示。水冷型煙道隨著煙氣流速的增加，沉積速率逐漸減小^[9]。煙氣流速越低，顆粒的跟隨性就越差，顆粒的慣性力和擴(kuò)散力都小于阻力。隨著煙氣流速的增加，顆粒的沉積會(huì)使煙氣的尾跡渦面積減小，使得煙道內(nèi)的流場(chǎng)發(fā)生變化，顆粒難以進(jìn)一步發(fā)生沉積行為。

煙道沉積云圖如圖5所示。隨著煙氣速度的增加，煙道的入口段、Ⅰ段、Ⅱ段、Ⅲ段和出口段沉積速率逐漸減小。隨著煙氣流速的增加，顆粒受煙氣的夾帶作用，慣性作用更加明顯，進(jìn)而減弱了顆粒的湍流渦和布朗擴(kuò)散作用。

不同煙氣流速對(duì)各管段沉積情況的影響如圖6所示。顆粒在入口段、Ⅰ段、Ⅱ段、Ⅲ段和出口段不易產(chǎn)生沉積行為。隨著煙氣流速的增加，盲管段和三通段的沉積速率逐漸增大。當(dāng)煙氣流速較大時(shí)，顆粒動(dòng)能增大，使得顆粒所受慣性作用更大，因此顆粒不易跟隨煙氣流動(dòng)。煙氣流速的增大會(huì)使顆粒在三通處流到盲管段的趨勢(shì)更加明顯，因此在較大的煙氣流速下顆粒在盲管和三通處更易產(chǎn)生沉積行為。

3.2 "顆粒粒徑對(duì)沉積速率的影響

顆粒粒徑是影響顆粒沉積速率的另一個(gè)主要因素?？紤]了顆粒直徑為5～140 μm時(shí)對(duì)顆粒沉積速率的影響。不同顆粒粒徑對(duì)沉積速率的影響如圖7所示。

由圖7可知，隨著顆粒粒徑的增大，沉積速率逐漸增大。顆粒沉積主要受到渦流擴(kuò)散碰撞、湍流擴(kuò)散和慣性碰撞的影響，并且不同粒徑受到的影響各不相同。對(duì)于小粒徑顆粒，容易被煙道中的渦流夾帶，增加了顆粒與壁面碰撞的概率。大粒徑顆粒自身質(zhì)量較大，所以在煙氣中所具有的慣性力較大，因此大粒徑顆粒受慣性碰撞的影響產(chǎn)生沉積行為。

4 "結(jié) 論

本文利用FLUENT對(duì)礦熱爐煙道內(nèi)煙氣顆粒沉積行為進(jìn)行模擬，預(yù)測(cè)了礦熱爐水冷型煙道煙氣流動(dòng)特性。研究了不同煙氣流速和不同顆粒粒徑對(duì)水冷型煙道顆粒沉積速率的影響，研究結(jié)論如下：影響顆粒沉積速率的因素有煙氣流速和顆粒粒徑，顆粒沉積速率隨著煙氣流速的增加而減少。盲管段和三通段的顆粒沉積速率隨著煙氣流速的增加而增加，其余各管段的顆粒沉積速率隨著煙氣流速的增加而減少。顆粒的沉積速率隨著顆粒粒徑的增加而增加。隨著顆粒粒徑的增大，盲管段和三通段的沉積速率增大，其余各管段的顆粒沉積速率逐漸減少。

參考文獻(xiàn)：

[1] LIU P， LI B K， CHEUNG S C P， et al. Material and energy flows in rotary kiln-electric furnace smelting of ferronickel alloy with energy saving[J]. Applied Thermal Engineering， 2016， 109： 542-559.

[2] CROSBY J M， LEWIS S， BONS J P， et al. Effects of temperature and particle size on deposition in land based turbines[J]. Journal of Engineering for Gas Turbines and Power， 2008， 130（5）： 051503.

[3] HAN Z M， XU Z M， SUN A D， et al. The deposition characteristics of micron particles in heat exchange pipelines[J]. Applied Thermal Engineering， 2019， 158： 113732.

[4] BRAVO L G， JAIN N， KHARE P， et al. Physical aspects of CMAS particle dynamics and deposition in turboshaft engines[J]. Journal of Materials Research， 2020， 35（17）： 2249-2259.

[5] HONG W P， WANG B H， ZHENG J X. Numerical study on the influence of fine particle deposition characteristics on wall roughness[J]. Powder Technology， 2020， 360： 120-128.

[6] ABD-ELHADY M S， RINDT C C M， WIJERS J G， et al. Particulate fouling in waste incinerators as influenced by the critical sticking velocity and layer porosity[J]. Energy， 2005， 30（8）： 1469-1479.

[7] LU H， LU L， JIANG Y. Numerical study of monodispersed particle deposition rates in variable-section ducts with different expanding or contracting ratios[J]. Applied Thermal Engineering， 2017， 110： 150-161.

[8] BRACH R M， DUNN P F. A mathematical model of the impact and adhesion of microsphers[J]. Aerosol Science and Technology， 1992， 16（1）： 51-64.

[9] TANG S Z， WANG F L， REN Q L， et al. Fouling characteristics analysis and morphology prediction of heat exchangers with a particulate fouling model considering deposition and removal mechanisms[J]. Fuel， 2017， 203： 725-738.

Simulation Study on Water-Cooled Flue of Submerged Arc Furnace"Based on Computational Fluid Dynamics

GONG"Ke-xin¹， LIU"Wei¹， LIU Peng¹， WANG Cui-hua¹， GUAN Qing-chen¹， SUN Chao¹

（1.School of Mechanical and Power Engineering， Shenyang University of Chemical Technology，"Shenyang Liaoning 113142，China）

Abstract：""In order to study the flow characteristics of flue gas in water-cooled flue in submerged arc furnace，"FLUENT was used to simulate the behavior of flue gas"particlenbsp;deposition in the water-cooled flue of submerged arc furnace."The RNG k-ε turbulence model was used to calculate the flow field and heat transfer."At the same time， using DPM model， considering Grant and Tabakoff particle rebound model and JKR theory， the judgment criterion of particle deposition was"established， and the stress and motion trajectory of furnace gas particles in the collision process were analyzed.

Key words：""Electric furnace; Water cooled flue; Particle deposition; Numerical simulation