高爐鼓風脫濕系統(tǒng)除霧器擋板結(jié)構(gòu)兩相流場數(shù)值模擬

羅雨慧　謝建中　左可

摘要：

對高爐鼓風脫濕系統(tǒng)除霧器的除霧機理進行了詳細分析，在此基礎(chǔ)上參考實際運行參數(shù)，采用計算流體力學(xué)（CFD）方法對除霧器擋板結(jié)構(gòu)內(nèi)的氣液兩相流動進行了數(shù)值模擬，得到了液滴的運動軌跡以及液滴質(zhì)量濃度、壓力、速度和旋渦分布情況，并對除霧器擋板結(jié)構(gòu)內(nèi)部兩相流場進行了深入分析.結(jié)果表明：鼓風脫濕系統(tǒng)除霧器能夠?qū)y霧氣流中的液滴進行有效的分離，從而保證鼓風脫濕系統(tǒng)連續(xù)可靠地運行.研究對除霧器的優(yōu)化設(shè)計和運行具有指導(dǎo)意義.

關(guān)鍵詞：

除霧器； 鼓風脫濕； 兩相流； 壓降； 流場

中圖分類號： TQ 051.8文獻標志碼： A

Abstract：

Based on the analysis of the demisting mechanism of the demister in the blower dehumidification system of a blast furnace，the computational fluid dynamics（CFD） method was used to simulate numerically the twophase flow of the gas and liquid in the baffle structure of a demister with the actual operation conditions as a reference.The detailed distributions of the pressure，velocity，vortex，and droplet concentration were provided.The simulation also showed the tracks of the droplets.In addition，an indepth analysis and research on the twophase flow were conducted.The results showed that the demister could separate the droplets from the airflow effectively，so as to ensure that the blower dehumidification system could run continuously and reliably.This study is helpful to the optimization and operation of the demister.

Keywords：

demister； blast dehumidification； two phase flow； pressure drop； flow field

鋼鐵企業(yè)高爐鼓風站的電能消耗十分巨大，其中，鼓風脫濕系統(tǒng)中除霧器阻損對鼓風機的電能消耗影響較大[1].高爐鼓風系統(tǒng)設(shè)置的除霧器主要用于清除由鼓風冷凝時產(chǎn)生的游離霧滴，提高鼓風質(zhì)量.

本文擬對鼓風脫濕系統(tǒng)中絲網(wǎng)與擋板組合式除霧器的分離機理進行詳細分析，在此基礎(chǔ)上對除霧器擋板結(jié)構(gòu)內(nèi)部氣液兩相流動進行數(shù)值模擬，得到流道內(nèi)液滴的運動軌跡及壓力、速度、旋渦分布，為除霧器的進一步優(yōu)化設(shè)計提供指導(dǎo).

1除霧器的結(jié)構(gòu)及工作原理

在脫濕除霧器中，擋板主要起導(dǎo)流作用，究其原因在于多層絲網(wǎng)結(jié)構(gòu)的除霧效率很高（根據(jù)現(xiàn)場測試數(shù)據(jù)，對于粒徑大于1 μm的粒子，分離效率高達99%；對于粒徑大于2 μm的粒子，分離效率高達99.6%）.攜霧氣流流過絲網(wǎng)結(jié)構(gòu)后，其中所夾帶的液滴幾乎全部被捕集除去.液滴在重力作用下沿絲網(wǎng)的絲徑向向下運動，同時繼續(xù)吸附氣體中夾帶的霧滴，凝聚變大的霧滴滴落在擋板上，沿擋板流入水槽中，進而排出除霧器外，實現(xiàn)除霧.

本文所研究的鼓風脫濕系統(tǒng)中除霧器布置在脫濕器內(nèi)二級冷卻器之后，為壓板固定的框架式8層×8塊鋼絲網(wǎng)結(jié)構(gòu)（7塊斜擋板）.除霧器除霧構(gòu)件如圖1所示.在除霧器中，攜霧氣流流過鋼絲網(wǎng)層實現(xiàn)氣液分離，而被絲網(wǎng)捕集分離和與擋板撞擊分離的液滴則經(jīng)擋板導(dǎo)流引入水槽中，進而被排出除霧器外.

攜霧氣流在除霧構(gòu)件內(nèi)流動，由于流線偏折，在慣性作用下，液滴不能隨氣流偏轉(zhuǎn)而撞擊到擋板上，其中：動量較小的液滴粘附在擋板表面被捕獲；動量較大的大液滴撞擊擋板表面發(fā)生濺射，產(chǎn)生多個小液滴.聚集在擋板表面的水滴受三種力的作用，即拉力、重力和表面張力.當重力占主導(dǎo)地位時，液滴在重力作用下沿擋板表面流入集液槽排出.氣流沖刷擋板表面的液膜，將其卷起、帶走.撞擊在擋板表面的液滴由于自身的動量過大而破裂、飛濺，均可能導(dǎo)致霧沫的二次夾帶.

2模型建立及網(wǎng)格劃分

攜霧氣流在除霧器內(nèi)的流動實際上是一種可壓縮的黏性流體的三維、非定常的復(fù)雜流動.對這一實際流動情況無法采用數(shù)學(xué)形式精確描述.本文從既能較好反映實際情況又力求模型構(gòu)造簡單角度出發(fā)，在合理誤差范圍內(nèi)對流場作簡化：

（1） 由于氣流速度較小，馬赫數(shù)遠小于0.1，故可把氣體視為不可壓縮氣體處理[2]；

（2） 考慮在實際的穩(wěn)定工作條件下，流動參數(shù)與時間的關(guān)系及氣流的振蕩對流場的影響可忽略，故將流動視為定常流動；

（3） 采用冷態(tài)條件，過程中不考慮溫度的影響，不考慮氣液傳熱和液體蒸發(fā)等現(xiàn)象；

（4） 由于液滴粒徑很小，故可作球形處理，考慮流動中其重力和氣相對液滴的曳力；同時假定其在運動中直徑不變，不考慮液滴之間的碰撞、聚合等現(xiàn)象，且忽略蒸發(fā)、摩擦、撕裂及熱效應(yīng)的影響，不考慮氣液兩相之間的任何能量交換；

（5） 液滴無濺射，不考慮壁面反彈、液膜形成與撕裂，忽略二次帶水的影響；

（6） 液滴接觸擋板，即認為被捕集；液滴到達除霧器出口時，即認為液滴逃逸.

本文采用計算流體力學(xué)（CFD）應(yīng)用軟件Fluent進行數(shù)值計算，應(yīng)用其前處理軟件Gambit生成網(wǎng)格，并參考以往的數(shù)值計算結(jié)果和經(jīng)驗[3]，整個計算區(qū)域采用非均勻的網(wǎng)格布置方式.為了提高整體網(wǎng)格質(zhì)量，網(wǎng)格劃分時先對除霧器部分進行網(wǎng)格加密，采用內(nèi)部面將除霧器所在體與其他體分離.

將模型分割為四個計算區(qū)域：方管段、除霧器、漸縮段和圓管段.先劃分除霧器的相關(guān)網(wǎng)格，在擋板結(jié)構(gòu)上生成線網(wǎng)格，采用相似邊界軟連接，再采用混合結(jié)構(gòu)網(wǎng)格從已有邊界網(wǎng)格生成擋板結(jié)構(gòu)體網(wǎng)格，最后生成除霧器的體網(wǎng)格.除霧器簡化三維模型如圖2所示.其四個計算區(qū)域的網(wǎng)格劃分均采用混合結(jié)構(gòu)網(wǎng)格，共生成75 965個節(jié)點，353 251個混合形式網(wǎng)格.除霧器網(wǎng)格劃分如圖3所示.

3數(shù)值計算方法及邊界條件

在除霧器內(nèi)流場的數(shù)值模擬中，通道內(nèi)包括互相之間交換質(zhì)與能的三相——氣相、氣流夾帶的液滴相和擋板表面的液膜相.本文采用兩相流模型，只考慮氣相和液滴相.

對于連續(xù)相（氣相），由雷諾數(shù)可知流動為湍流.本文采用SIMPLE算法進行壓力-速度耦合，采用有限體積法對算例進行離散處理.壓力采用Standard離散格式，動量、湍流動能、湍流耗散率均采用二階迎風離散格式，以獲得較準確的解.

對于離散相（液滴相），由于除霧器中氣流內(nèi)液滴的體積分數(shù)小于10%，故可忽略顆粒之間的相互作用、顆粒對氣相的影響、顆粒的運動軌跡[3-5].本文選用基于LagraianEulerian法的DPM（discrete phase model）模型，按拉格朗日方法對各個顆粒方程進行積分求解.

交替求解離散相與連續(xù)相的控制方程，直到兩者均收斂從而實現(xiàn)離散相與連續(xù)相雙向耦合.

攜霧氣流中液滴流動軌跡計算時只需考慮阻力和重力的影響[6].

氣相作為連續(xù)相在歐拉坐標系中描述；液滴相作為離散相在拉格朗日坐標系中描述[7].初始條件及邊界條件分別為：

（1） 連續(xù)相（氣相）

介質(zhì)：空氣，密度為1.205 kg·m-3，動力黏度（20℃）為18.1 μPa·s.

進口條件：給定氣流速度uy=1.2 m·s-1，ux=0，假設(shè)其在進口截面為均勻分布，湍流度為0.05[8].

出口條件：自由出流.

壁面條件：無滑移，絕熱.

（2） 離散相（液滴相）：

介質(zhì)：水，呈細小液滴狀，密度為1.0×103 kg·m-3，平均液滴直徑為5 μm，鼓風含濕量為10 g·Nm-3.

進口條件：設(shè)定液滴速度與氣流速度相同[9]，噴射類型選為表面，使液滴在進口截面均勻分布.

壁面條件：選擇捕集類型，即不考慮反彈，液滴觸及壁面即認為被捕集，不考慮二次夾帶效應(yīng).

當計算連續(xù)性殘差、速度殘差、湍流動能殘差、湍流耗散率殘差均降至10-3，且入口與出口氣相流量相差小于3%時認為計算收斂[10].

4流場分布及分析

4.1液滴質(zhì)量濃度分布

圖4為計算得到的液滴質(zhì)量濃度c的分布.由圖可知：在捕集液滴的除霧構(gòu)件區(qū)域及起導(dǎo)流作用的擋板附近液滴質(zhì)量濃度明顯較高，而在整個除霧構(gòu)件中第四至第六塊斜擋板之間區(qū)域的液滴質(zhì)量離效果已足以保證系統(tǒng)連續(xù)、可靠地運行.

4.2液滴顆粒的運動軌跡

采用Fluent軟件對除霧器進行模擬，得到的連續(xù)相和離散相的運動軌跡如圖5所示，其中圖5（b）為計算得到的除霧器內(nèi)液滴的隨機運動軌跡.由圖5（b）可知，液滴運動軌道主要集中在除霧器除霧構(gòu)件中的中間五塊擋板部分區(qū)域，這對于除霧器的工程優(yōu)化設(shè)計有著重要的指導(dǎo)意義.此外，液滴主要在除霧器中被捕集，小部分進入漸縮段與壁面發(fā)生碰撞而被分離.被捕集的液滴主要源自攜霧氣流的主流兩側(cè)，由于受到氣流速度梯度影響，其在流線偏折時所受離心慣性作用較大，容易發(fā)生偏轉(zhuǎn)，撞擊在壁面上而被捕集.故除霧器入口處液滴參數(shù)的設(shè)計對除霧效率的提高具有較大影響.

由圖5（b）亦可知，鼓風脫濕系統(tǒng)中的絲網(wǎng)與擋板組合式除霧器可有效地實現(xiàn)氣液分離.

4.3壓力及速度場分析

除霧器中壓力和速度分布云圖如圖6所示.由圖可知，方管段和圓管段壓力、速度分布較為均勻.攜霧氣流流過除霧器擋板結(jié)構(gòu)的過程中未出現(xiàn)顯著的壓降，而流過除霧器出口處的漸縮段時發(fā)生了明顯的壓降，降幅達142 Pa左右.其產(chǎn)生的主要原因是流體的流通面積減小，并且截面形狀突然改變，流體流線被迫發(fā)生改變[11].由于旋流作用較強，在漸縮段方圓突變截面后出現(xiàn)了明顯的低壓區(qū)，在漸縮段后方圓突變截面處及其后一段距離內(nèi)亦出現(xiàn)了顯著的壓降.

4.4旋渦分布

脫濕除霧器內(nèi)的湍流動能和湍流耗散率分布如圖7所示.在整個脫濕除霧器中，湍流耗散最強烈的區(qū)域為漸縮段后一段距離.在除霧器除霧構(gòu)件中，中間五塊擋板部分區(qū)域的湍流動能最大，湍流耗散強烈，是除霧器內(nèi)實現(xiàn)氣液分離的關(guān)鍵區(qū)域，攜霧氣流流過除霧器通道時，氣流主要流道在該區(qū)域.由于慣性力的作用，氣流中液滴的跟隨性變差，速度遲豫時間延長，所以易于碰到壁面而被捕集[12].在第三、四塊擋板下部區(qū)域的高壓低速區(qū)是除霧器擋板結(jié)構(gòu)中湍流動能最大、湍流耗散最強烈的區(qū)域.該區(qū)域除霧效率較高.氣體流過除霧器后在最上及最下處豎直擋板后易形成回流區(qū)，此處亦可能產(chǎn)生角渦，加劇了流場的擾動，增大了系統(tǒng)壓力損失.上部豎直擋板附近流域局部放大圖如圖8所示.

5結(jié)論

本文運用兩相流模型對鼓風脫濕系統(tǒng)除霧器擋板結(jié)構(gòu)的主要性能進行了數(shù)值模擬，得到了液滴質(zhì)量濃度、壓力、速度及旋渦分布情況，展示了其三維內(nèi)部流場，揭示了實驗手段難以獲得的數(shù)據(jù)和現(xiàn)象.計算結(jié)果表明：鼓風脫濕系統(tǒng)除霧器能夠有效實現(xiàn)攜霧氣流中液滴的分離.模擬結(jié)果對除霧器結(jié)構(gòu)設(shè)計、降低除霧器內(nèi)流動阻力具有一定的參考性意義.

參考文獻：

[1]張永紅，袁熙志，羅冬梅，等.我國鋼鐵行業(yè)節(jié)能降耗現(xiàn)狀與發(fā)展[J].工業(yè)爐，2013，35（3）：12-17.

[2]周光炯，嚴宗毅，許世雄，等.流體力學(xué)[M].2版.北京：高等教育出版社，2006.

[3]李鵬飛，徐敏義，王飛飛.精通CFD工程仿真與案例實戰(zhàn)[M].北京：人民郵電出版社，2011.

[4]孫秀君，孫海鷗，姜任秋.油氣分離器內(nèi)油滴軌跡的數(shù)值模擬[J].應(yīng)用科技，2006，33（10）：69-72.

[5]李文艷，徐妍，申林艷.數(shù)值模擬技術(shù)在濕法脫硫除霧器優(yōu)化設(shè)計上的應(yīng)用[J].熱力發(fā)電，2007，29（5）：10-15.

[6]韓占忠，王敬，蘭小平.FLUENT——流體工程仿真計算實例與應(yīng)用[M].2版.北京：北京理工大學(xué)出版社，2010.

[7]郭鵬宇，楊震.除霧器通道內(nèi)二維兩相流場的數(shù)值模擬[J].電力科學(xué)與工程，2005（2）：31-33.

[8]王福軍.計算流體動力學(xué)分析[M].北京：清華大學(xué)出版社，2004.

[9]趙毅，華偉，王亞君，等.濕式煙氣脫硫塔中折線型擋板除霧器分離效率的數(shù)值模擬[J].動力工程學(xué)報，2005，25（2）：293-297.

[10]陳凱華，宋存義，李強，等.濕法煙氣脫硫系統(tǒng)中折板式除霧器性能的數(shù)值模擬[J].環(huán)境工程學(xué)報，2007，1（7）：91-96.

[11]CRANE工程部編.流體流經(jīng)閥門、管件和管道的流體計算TP410[M].北京：化學(xué)工業(yè)出版社，2013.

[12]RUDINGER G.氣體-顆粒流基礎(chǔ)[M].張遠君譯.北京：國防工業(yè)出版社，1986.