楊成禹

(中冶賽迪電氣技術(shù)有限公司，重慶 400010)

高爐渣是高爐煉鐵過(guò)程中的主要副產(chǎn)品，其排放溫度在1 450～1 650 ℃[1]，是一種高品位的余熱資源。近年來(lái)，高爐渣干法?；坝酂釙?huì)水技術(shù)得到了廣泛關(guān)注[2]。針對(duì)粒化后的高溫渣粒，自流床余熱鍋爐直接回收高溫渣粒余熱的技術(shù)具有顯著優(yōu)勢(shì)。與其他方法相比，該技術(shù)降低了空氣鼓風(fēng)的能量損失，減少了余熱回收過(guò)程造成的環(huán)境污染，使余熱回收效率得到明顯提高。

在自流床余熱鍋爐內(nèi)，高溫顆粒在重力的作用下向下流動(dòng)，顆粒流繞流鍋爐管束進(jìn)行換熱。在這個(gè)過(guò)程中，顆粒的流動(dòng)行為將直接影響著高溫顆粒與鍋爐管束的換熱特性。Kurochkin[3]通過(guò)實(shí)驗(yàn)對(duì)顆粒環(huán)繞圓管流動(dòng)特性進(jìn)行研究，發(fā)現(xiàn)管束頂部存在停滯區(qū)，且顆粒直徑越大，停滯區(qū)范圍越小。Niegsch[4]等人通過(guò)實(shí)驗(yàn)研究了小直徑顆粒繞流叉排管束時(shí)的流動(dòng)行為，對(duì)管束停滯區(qū)與空隙區(qū)范圍的計(jì)算方法進(jìn)行分析，并搭建了漏斗流模型用以描述顆粒繞流管束時(shí)的速度場(chǎng)分布。于慶波等人[5]在進(jìn)行高溫渣粒繞流管束流動(dòng)與傳熱特性的數(shù)值模擬研究中，將流動(dòng)的高爐渣顆粒當(dāng)作連續(xù)的黏性流體，發(fā)現(xiàn)增大顆粒的入口速度，換熱效果提高。Bartsch[6]研究了埋有水平換熱管的移動(dòng)床內(nèi)顆粒的流動(dòng)換熱過(guò)程，并分別建立了CFD模型和DEM模型進(jìn)行了對(duì)比分析，發(fā)現(xiàn)兩種模型都能較好描述顆粒繞流圓管的過(guò)程。

目前，針對(duì)高溫顆粒繞管流動(dòng)與傳熱耦合特性的研究尚需進(jìn)一步進(jìn)行完善與補(bǔ)充，而顆粒的流動(dòng)特性直接決定了其與管束間的傳熱效果，因此有必要深入研究高溫顆粒繞流換熱管的流動(dòng)與傳熱特性。本文采用MFiX軟件，分析高爐渣顆粒繞流換熱管束的流動(dòng)與傳熱耦合特性，研究換熱管型及其排列方式對(duì)顆粒流動(dòng)特性與傳熱特性的影響，進(jìn)而指導(dǎo)自流床余熱鍋爐的結(jié)構(gòu)設(shè)計(jì)。

1 模型建立

1.1 幾何模型

選取自流床余熱鍋爐結(jié)構(gòu)的特征單元，并將其簡(jiǎn)化為一個(gè)帶管束的矩形模型[7]，如圖1所示。其中換熱管的外徑為38 mm，從上至下共設(shè)置4排，換熱管橫向與縱向中心間距均為76 mm，即2倍管外徑。模型頂端距離首排管束頂端57 mm，底端距離最后一排管束底端19 mm。在本研究中，換熱管選用3種管型，分別為圓管、橢圓管與菱形管，且其中心定位位置保持一致。

圖1 幾何模型圖

1.2 數(shù)學(xué)模型

顆粒繞流換熱管的流動(dòng)過(guò)程中，顆粒的運(yùn)動(dòng)狀態(tài)可以分為平動(dòng)與轉(zhuǎn)動(dòng)[8]，故對(duì)于單個(gè)顆粒i來(lái)說(shuō)，其運(yùn)動(dòng)的控制方程如式(1)～式(2)所示。

平動(dòng)方程：

(1)

轉(zhuǎn)動(dòng)方程

(2)

本文研究的高爐渣顆粒溫度高于1 000 ℃，所以輻射換熱不可忽略。顆粒在自流床余熱鍋爐中的移動(dòng)速度相對(duì)較小，顆粒間隙內(nèi)的空氣被顆粒攜帶移動(dòng)，故對(duì)流換熱作用小，所以換熱過(guò)程中忽略氣固兩相間的對(duì)流換熱。因此換熱過(guò)程包括顆粒之間的導(dǎo)熱和輻射兩大部分，其中導(dǎo)熱又可以分為：通過(guò)顆粒與顆粒之間直接接觸面積的導(dǎo)熱、通過(guò)顆粒與顆粒之間氣膜的導(dǎo)熱、通過(guò)顆粒與換熱管壁面接觸面積的導(dǎo)熱以及通過(guò)顆粒與換熱管壁面之間氣膜的導(dǎo)熱；輻射換熱部分同樣可以分為顆粒與顆粒和顆粒與換熱管壁面之間的輻射換熱。故單個(gè)顆粒i傳熱的控制方程為

(3)

式中：cp,i為顆粒i的比熱容，J/(kg·℃)；Ti為顆粒i在t時(shí)刻的溫度，℃；Qi,j,cond、Qi,j,rad分別為顆粒i與其他顆粒的導(dǎo)熱量和輻射換熱量，J；Qi,tube,cond、Qi,tube,rad分別為顆粒i與換熱管壁面的導(dǎo)熱量和輻射換熱量，J。

在邊界條件的設(shè)置中，模型上方為高溫顆粒的入口，設(shè)置為質(zhì)量流量邊界，給定顆粒的入口溫度與顆粒的質(zhì)量流量；模型下方為顆粒的出口，同樣設(shè)置為質(zhì)量流量邊界，給定顆粒出口流速；換熱管束外壁為固定溫度的無(wú)滑移壁面；其余表面均為絕熱的無(wú)滑移壁面。

1.3 網(wǎng)格劃分及參數(shù)設(shè)置

計(jì)算網(wǎng)格為均勻矩形網(wǎng)格，網(wǎng)格尺寸為2倍顆粒直徑。根據(jù)目前的檢測(cè)水平，不同溫度條件下高爐渣顆粒的相關(guān)物性參數(shù)難于確定。所以，在本研究中高爐渣顆粒的物性參數(shù)保持恒定，具體參數(shù)如表1所示。

表1 高爐渣顆粒的參數(shù)

2 結(jié)果與分析

2.1 換熱管管型及排列方式對(duì)顆粒流動(dòng)特性的影響

圖2為不同管型和管排列方式條件下顆粒的運(yùn)動(dòng)軌跡。在遠(yuǎn)離換熱管束的上部區(qū)域，無(wú)論何種管排列方式，顆粒的流動(dòng)過(guò)程為平推流[9]，其運(yùn)動(dòng)軌跡為直線(xiàn)。這是由于運(yùn)動(dòng)過(guò)程中，顆粒間的相互作用較小，且不受管束等其他物體的阻礙作用，因此顆粒的運(yùn)動(dòng)不會(huì)發(fā)生橫向偏移。當(dāng)顆?？拷着艙Q熱管時(shí)，顆粒的運(yùn)動(dòng)軌跡開(kāi)始發(fā)生變化，這是由于在換熱管上部區(qū)域處顆粒的運(yùn)動(dòng)發(fā)生了較大變化，然后通過(guò)力鏈的傳遞作用，進(jìn)而影響靠近換熱管的顆粒。顆粒與首排管束接觸碰撞后，便產(chǎn)生了明顯的水平速度，故發(fā)生橫向偏移，此時(shí)貼壁顆粒沿著換熱管的外輪廓進(jìn)行移動(dòng)，其運(yùn)動(dòng)軌跡與管型相似。

圖2 顆粒流線(xiàn)圖

當(dāng)換熱管順排排列時(shí)，顆粒流過(guò)首排管束后，其運(yùn)動(dòng)軌跡基本為一條直線(xiàn)。繞流區(qū)的顆粒在下落過(guò)程中與管壁接觸，并發(fā)生繞流，其運(yùn)動(dòng)軌跡近似為半“S”形，非繞流區(qū)顆粒大多未與管壁接觸，且受到其他顆粒的擾動(dòng)較小，其運(yùn)動(dòng)軌跡近似為一條直線(xiàn)。從圖2(a)～圖2(c)可以看出，當(dāng)換熱管為圓管或者橢圓管時(shí)，由于換熱管的阻礙作用，顆粒在換熱管正上方形成顆粒“停滯區(qū)”；在重力的作用下，顆粒繞流換熱管移動(dòng)至換熱管下方，然后顆粒與換熱管分離，在換熱管正下方形成顆?！翱障秴^(qū)”；當(dāng)換熱管為菱形管時(shí)，顆粒沿著菱形管的垂直壁面直接向下流動(dòng)，然后沿著斜面向菱形管的正下方移動(dòng)。

當(dāng)換熱管叉排排列時(shí)，顆粒繞管換熱管的運(yùn)動(dòng)軌跡近似為“S”型曲線(xiàn)[9]。在這種布置方式下，顆粒的整體流動(dòng)狀態(tài)則是由入口處的平推流過(guò)渡為管束外的“S”形流動(dòng)。在y方向上，整個(gè)流動(dòng)區(qū)域由非繞流區(qū)和繞流區(qū)交替組合形成，致使管束外的每個(gè)顆粒都會(huì)繞管流動(dòng)，具有較大的x方向偏移量，故顆粒間的摻混效果更好。此外，在相同管間距的情況下，由于在叉排排列時(shí)具有“S”形運(yùn)動(dòng)軌跡，顆粒的實(shí)際流動(dòng)距離要比順排排列時(shí)顆粒的流動(dòng)距離更長(zhǎng)，這說(shuō)明叉排時(shí)顆粒與管壁的接觸時(shí)間長(zhǎng)，利于顆粒與換熱管壁的充分換熱。

圖3為不同管束排列方式條件下首排換熱管處顆粒在y方向的速度分布圖。從圖中可以看出，順排時(shí)顆粒的y方向速度分布由正反的“U”形峰交替組合而成，而叉排時(shí)顆粒的y方向速度分布則是由正反的“V”形峰交替組合而成。二者在非繞流區(qū)的流速相差無(wú)幾，故使顆粒的豎直流動(dòng)距離相當(dāng)。但是在“停滯區(qū)”附近，叉排時(shí)顆粒具有更大的流速，這也表明叉排時(shí)顆粒的整體流動(dòng)性更好，“停滯區(qū)”內(nèi)顆粒也在參與流動(dòng)，進(jìn)而能強(qiáng)化“停滯區(qū)”內(nèi)顆粒的擾動(dòng)，強(qiáng)化該區(qū)域顆粒與流動(dòng)區(qū)域內(nèi)顆粒的混合、傳熱。

圖3 水平方向上顆粒的縱向速度分布

圖4為不同管型和管排列方式條件下顆粒與管壁的接觸面積。顆粒與換熱管壁的接觸面積越大，則顆粒與換熱管之間的導(dǎo)熱量則增多。由于顆粒在換熱管周?chē)淖饔昧Ψ较虿煌趽Q熱管的圓周方向，顆粒與管壁的接觸面積發(fā)生著顯著變化。在換熱管的正上方，顆粒與管壁的接觸面積最大，而在換熱管的正下方，顆粒與換熱管壁的接觸面積為0，說(shuō)明該區(qū)域?yàn)轭w?？障秴^(qū)。當(dāng)換熱管順排排列時(shí)，顆粒與管周向壁面的接觸面積基本相似，其中六邊形管的空隙區(qū)最大，橢圓管的空隙區(qū)最小。當(dāng)換熱管叉排排列時(shí)，顆粒與管上部的接觸面積較大，其中同樣是六邊形管的空隙區(qū)最大，橢圓管的空隙區(qū)最小。因此，顆粒與橢圓換熱管壁間的接觸更緊密，從而更有利于顆粒與管壁間的熱交換。

2.2 換熱管型及排列方式對(duì)顆粒傳熱特性的影響

圖5為不同管型和管排列方式條件下顆粒的溫度分布圖。通過(guò)前面的分析可知，換熱管排列方式的不同使顆粒在管間的運(yùn)動(dòng)軌跡不同，從而影響管間的顆粒溫度分布。換熱管順排排列時(shí)，在換熱管的上下相鄰管間存在明顯的顆粒流動(dòng)“停滯區(qū)”，此區(qū)域內(nèi)的顆粒流動(dòng)緩慢，進(jìn)而與換熱管接觸換熱的時(shí)間較長(zhǎng)，所以在溫度分布圖中可以觀察到顆粒的低溫區(qū)域出現(xiàn)在下層換熱管的正上方。換熱管叉排排列時(shí)，管間不存在明顯的流動(dòng)死區(qū)，故在同一水平高度處顆粒的溫度分布相對(duì)均勻。從圖5中可以看出，換熱管管型的變化對(duì)顆粒溫度分布的影響并不明顯。

圖4 顆粒與管周向壁面的接觸面積

圖5 顆粒溫度分布隨換熱管管型和排列方式的變化

圖6為不同顆粒流速時(shí)各換熱管管型的熱回收率。從圖6中可以看出，在所研究顆粒流速范圍內(nèi)，系統(tǒng)熱回收率隨顆粒流速的增大而不斷下降。雖然顆粒流速的增大會(huì)使換熱系數(shù)增大，但同時(shí)會(huì)有更多比例的顆粒未能與換熱管進(jìn)行充分的換熱，造成綜合熱回收率的下降，所以在較大顆粒流速工況下，換熱裝置需要裝備更多的換熱管。從圖6中還可以看出，在低顆粒流速下，采用圓形換熱管時(shí)系統(tǒng)的余熱回收率相對(duì)較高，這是由于在低流速范圍時(shí)，顆粒在換熱管周?chē)亩逊e相對(duì)密實(shí)，換熱管下方的顆?？障秴^(qū)均較小。當(dāng)換熱管的當(dāng)量直徑相等時(shí)，圓形換熱管的表面積要大于橢圓形換熱管的表面積，因而圓形換熱管與渣粒的換熱量相對(duì)較大。隨著顆粒流速的增加，顆粒運(yùn)動(dòng)到換熱管下半周的時(shí)候會(huì)很快脫離換熱管表面，換熱管下方的顆?？障秴^(qū)面積逐漸增大。而對(duì)于橢圓形換熱管來(lái)說(shuō)，由于顆粒的流動(dòng)方向與橢圓的長(zhǎng)軸方向相同，顆粒與橢圓換熱管表面的接觸時(shí)間相對(duì)較長(zhǎng)，即實(shí)際的換熱面積相對(duì)大一些，所以，在高顆粒流速下會(huì)出現(xiàn)橢圓形換熱管的熱回收率偏高。但是在實(shí)際應(yīng)用當(dāng)中，橢圓形換熱管的加工難度相對(duì)較大，該類(lèi)型的換熱裝置的制造成本會(huì)相對(duì)較高，所以，在確定換熱管類(lèi)型的時(shí)候要結(jié)合顆粒的流速及熱回收裝置的制造成本進(jìn)行綜合考量。

圖6 換熱管管型對(duì)熱回收率的影響

3 結(jié) 論

本文以開(kāi)源代碼MFiX為平臺(tái)，基于CFD-DEM方法，對(duì)自流床內(nèi)顆粒物料繞流換熱管束的流動(dòng)與傳熱過(guò)程進(jìn)行了數(shù)值模擬研究，獲得了顆粒繞流管束流動(dòng)與傳熱的耦合特性。

(1)在換熱管順排排列時(shí)，顆粒在流經(jīng)首排換熱管時(shí)產(chǎn)生橫向偏移，而后運(yùn)動(dòng)軌跡近似為直線(xiàn)；而在換熱管順排排列時(shí)叉排時(shí)，顆粒繞管呈“S”型流動(dòng)，顆粒在通道內(nèi)的流動(dòng)時(shí)間相對(duì)較長(zhǎng)，且加強(qiáng)了顆粒在流動(dòng)過(guò)程中的擾動(dòng)作用，強(qiáng)化了顆粒與管壁的換熱作用。

(2)在換熱管順排排列時(shí)，由于在換熱管的下方存在明顯的顆粒“停滯區(qū)”，換熱管的上下相鄰管間存在低溫顆粒區(qū)域；在換熱管順排排列時(shí)叉排時(shí)，由于顆粒的流型強(qiáng)化了顆粒的擾動(dòng)作用，顆粒在管間同一水平位置的溫度分布更加均勻。

(3)由于換熱管結(jié)構(gòu)形式的不同，菱形管下方的空隙區(qū)最大，橢圓管下方的空隙區(qū)最小，因此顆粒與橢圓換熱管壁間的接觸更緊密，且在高顆粒流速下，橢圓形換熱管的熱回收率相對(duì)較高。但考慮橢圓形換熱管的加工難度大等問(wèn)題，換熱管的選型要根據(jù)熱回收系統(tǒng)的經(jīng)濟(jì)效益進(jìn)行綜合考量。