顏浩 陳強(qiáng) 安澤文 張龍龍 王翠蘋

摘要： 針對(duì)脫硫塔后煙氣中攜帶著大量水汽問題，本文采用數(shù)值模擬的方法，對(duì)除霧器內(nèi)部流場(chǎng)流動(dòng)及液體分布與團(tuán)聚情況進(jìn)行研究。采用Solidworks軟件，建立葉輪除霧器單元管的三維模型，對(duì)新型的被動(dòng)式葉輪除霧器一個(gè)單元管的三維模型進(jìn)行兩相流數(shù)值模擬。同時(shí)，為了更好地模擬實(shí)際流動(dòng)，選用準(zhǔn)確的湍流模型kε模型、群體平衡模型（population banlance model，PBM）和6DOF模型，并選擇SAMPLE算法。模擬結(jié)果表明，煙氣經(jīng)過除霧器時(shí)，葉片區(qū)域存在強(qiáng)烈湍流耗散，液滴被旋轉(zhuǎn)氣流拋向壁面，從而實(shí)現(xiàn)液滴團(tuán)聚和氣液分離;除霧器內(nèi)對(duì)流體切向速度、液滴顆粒的團(tuán)聚等作用都優(yōu)于固定式葉片，且液滴粒徑較小時(shí)更易團(tuán)聚，除霧器單元管內(nèi)除霧效率隨流速增加而提高。該研究為優(yōu)化除霧器結(jié)構(gòu)提供了理論依據(jù)。

關(guān)鍵詞： 被動(dòng)旋轉(zhuǎn); 葉片除霧器; 液滴團(tuán)聚; 數(shù)值模擬; 除霧效率

中圖分類號(hào)： X701; TB126? 文獻(xiàn)標(biāo)識(shí)碼： A

近年來，在電力行業(yè)中濕法煙氣脫硫系統(tǒng)得到廣泛應(yīng)用，濕法脫硫可以有效的控制SO2等氣態(tài)污染物排放，但是在脫除SO2的同時(shí)，煙氣中會(huì)夾帶大量飽和水汽，除少量從煙囪底部凝結(jié)排出，大部分液滴隨煙氣從煙囪頂部排出。液滴的排放會(huì)造成其溶解攜帶的PM2.5等粒子污染物增加，同時(shí)液滴夾帶是造成石膏雨的重要原因，所以脫除煙氣中攜帶的細(xì)小液滴非常必要[17]。目前，國內(nèi)常用的電廠除霧設(shè)備主要有折流板型、離心式型及正在推廣使用的濕式電除塵。根據(jù)霧滴的碰撞和粘附等原理，攜帶液滴的煙氣經(jīng)過密集的折流板時(shí)，煙氣在彎曲管道內(nèi)流動(dòng)產(chǎn)生的離心力使氣液分離，液滴黏附在板面上形成很薄的液膜，經(jīng)重力作用緩慢流到集液槽中，從而實(shí)現(xiàn)氣液分離[810]，但是這種板式除霧器對(duì)小于5 μm，甚至更小的細(xì)霧無法有效地去除。為保證除霧效率，對(duì)攜帶液滴煙氣的速度有一定限制，當(dāng)超過臨界氣速時(shí)，因?yàn)槎螉A帶使脫霧效率降低[1113]。離心式除霧器在洗滌器內(nèi)部設(shè)置葉片使氣流旋轉(zhuǎn)，一種是葉片固定，氣流從葉片縫隙穿過，形成旋流;另一種是葉片不固定，通過電機(jī)驅(qū)動(dòng)使葉片旋轉(zhuǎn)。由于電驅(qū)動(dòng)式阻力小，可提供比較大的圓周速度，其除霧效率比固定式高。但由于其提升除霧效率主要是提高能耗獲得，且這兩種除霧器不能根據(jù)風(fēng)速同步調(diào)節(jié)葉輪轉(zhuǎn)速，所以當(dāng)風(fēng)速波動(dòng)較大時(shí)，除霧效果會(huì)明顯降低[1417]?；诖耍疚奶岢鲆环N新型自旋轉(zhuǎn)葉輪除霧器，葉輪由氣流的流動(dòng)帶動(dòng)旋轉(zhuǎn)，從而增強(qiáng)除霧效果，且不需要其他額外動(dòng)力。通過Solidworks建立新型除霧器的單元管三維模型，選取適當(dāng)?shù)臄?shù)學(xué)模型及恰當(dāng)?shù)倪\(yùn)行參數(shù)，對(duì)除霧器單元管內(nèi)氣液兩相流動(dòng)、液顆粒的碰撞、粘附以及液滴的團(tuán)聚等進(jìn)行數(shù)值模擬，并對(duì)其氣液分離性能進(jìn)行初步預(yù)測(cè)。該研究為優(yōu)化除霧器結(jié)構(gòu)提供了理論依據(jù)。

1 除霧器物理模型的建立

采用Solidworks軟件，建立葉輪除霧器單元管的三維模型，尺寸設(shè)計(jì)參照離心式除霧器的類比設(shè)計(jì)。葉片為被動(dòng)式旋轉(zhuǎn)。葉片的扇形與中心線成一定角度，該夾角即為葉片夾角，新型除霧器單元管的物理模型如圖1所示。

該模型直徑200 mm，長300 mm，內(nèi)部進(jìn)水管直徑20 mm。在除霧器的每個(gè)單元管中，葉片數(shù)量m=18，葉片與水平夾角，即仰角為20°。以上幾何參數(shù)的選擇主要考慮如下幾點(diǎn)：

1） 葉片形狀為平板扇形，葉片的扭曲角為0°，這是考慮在實(shí)際過程中，制造、加工及建模的方便[18]。

2） 葉片仰角的經(jīng)驗(yàn)值為20°～ 30°，由于該除霧器為被動(dòng)旋轉(zhuǎn)，所以取仰角為最小值20°，以提高葉輪轉(zhuǎn)速。

3） 葉片數(shù)量，在旋流板除霧器設(shè)計(jì)中，根據(jù)氣液負(fù)荷與開孔面積，通常取24片或18片，考慮攜帶液滴量不太大，所以取18片，以減小氣液穿過時(shí)的流動(dòng)阻力。

4） 葉片的徑向角度取零，有利于旋流強(qiáng)度的增加及氣液分離。

2 除霧器單元管內(nèi)液滴團(tuán)聚與氣液分離的模擬

2.1 除霧器單元管內(nèi)流場(chǎng)的簡化與假設(shè)

除霧器性能受內(nèi)部流動(dòng)規(guī)律和液滴分離規(guī)律的制約。在除霧器中，由于攜帶液滴的煙氣葉片間的流動(dòng)較為復(fù)雜，為了提高計(jì)算效率，可在允許的誤差范圍內(nèi)對(duì)模型進(jìn)行簡化。

1） 本文的模擬以空氣代替煙氣，進(jìn)入除霧器的煙氣流速較小，一般為3～8 m/s，可把氣體視為理想氣體。

2） 在理想情況下，煙氣的流動(dòng)參數(shù)不隨時(shí)間變化，故流動(dòng)可視為定常流動(dòng)。

3） 模擬計(jì)算以水滴代替漿液滴，除霧器內(nèi)液滴占煙氣體積分?jǐn)?shù)較小，粒徑較?。ㄒ旱蔚闹睆綖槲⒚准?jí)），但由于要考慮細(xì)小液滴因氣流攜帶，通過葉片時(shí)發(fā)生旋轉(zhuǎn)而團(tuán)聚成較大液滴，使液滴的直徑發(fā)生變化，因此液滴間的相互作用不可忽略。

4） 在分析過程中，考慮氣體流動(dòng)對(duì)顆粒的粘性作用，忽略其反作用。

2.2 數(shù)學(xué)模型的選取

已有的對(duì)旋流的計(jì)算，多是在固定容器內(nèi)發(fā)生簡單的旋風(fēng)分離。本文模擬的是混合氣流沿管道進(jìn)入除霧器，受到葉片的阻力而推動(dòng)葉片轉(zhuǎn)動(dòng)，同時(shí)葉片帶動(dòng)氣流由平行流動(dòng)變?yōu)樾?，使液滴在垂直于軸線的方向上做離心運(yùn)動(dòng)，從而產(chǎn)生離心力將液滴甩向壁面，在重力作用下沿壁面緩慢流下來;并且強(qiáng)烈的湍流運(yùn)動(dòng)使液滴與液滴、葉片之間發(fā)生劇烈的撞擊，由于分子間的相互作用力，使液滴團(tuán)聚在一起，變?yōu)檩^大顆粒，更容易甩到壁面，從而脫除霧滴。為了更好地模擬實(shí)際流動(dòng)，需要選擇準(zhǔn)確的湍流模型。

2.2.1 標(biāo)準(zhǔn)kε模型

目前，湍流模型主要包括標(biāo)準(zhǔn)kε模型、RNG模型、標(biāo)準(zhǔn)kω模型、雷諾應(yīng)力模型、SpalartAlmaras模型等。為了適應(yīng)復(fù)雜多變的模擬環(huán)境，本文選用標(biāo)準(zhǔn)kε模型，模擬液滴顆粒因劇烈旋轉(zhuǎn)而產(chǎn)生的團(tuán)聚作用。

2.2.2 群體平衡模型（PBM）

在考慮粒徑分布的兩相體系中，氣體和液滴的粒徑分布會(huì)隨著兩相體系的反應(yīng)和傳遞現(xiàn)象的發(fā)生而變化。因此，不僅要考慮能量、動(dòng)量和質(zhì)量守恒，還需要添加一個(gè)可以描述粒子平衡的方程，即群體平衡模型。離散法中的非均一離散法可以應(yīng)用于多個(gè)離散粒子群，調(diào)用不同的相速度，將顆粒群的粒徑范圍離散為有限的粒徑間隔。本文采用非均一離散法，模擬計(jì)算各相的粒徑分布。

液滴的團(tuán)聚主要包括庫倫團(tuán)聚、熱團(tuán)聚、布朗團(tuán)聚和湍流團(tuán)聚。由于本研究中的顆粒為電中性，庫倫團(tuán)聚不加考慮;同時(shí)，熱團(tuán)聚對(duì)整個(gè)過程中團(tuán)聚影響較小，故忽略，因此只考慮湍流團(tuán)聚和布朗團(tuán)聚，自行編譯液滴團(tuán)聚的用戶自定義函數(shù)（user definefunction，UDF）進(jìn)行模擬計(jì)算。研究的液滴顆粒粒徑為微米級(jí)別，可視為自由分子區(qū)，此時(shí)布朗團(tuán)聚核為

2.2.3 6DOF模型

本文研究被動(dòng)式葉輪的旋轉(zhuǎn)，邊界運(yùn)動(dòng)規(guī)律是未知的，在動(dòng)網(wǎng)格計(jì)算中，網(wǎng)格變化需要多次調(diào)整，并且這類運(yùn)動(dòng)是由內(nèi)部流動(dòng)造成涉及到力和力矩的變化。對(duì)于此問題，在Fluent中采用6DOF模型進(jìn)行計(jì)算，6DOF即六自由度模型，一般用于解決被動(dòng)網(wǎng)格問題，涉及力和力矩的計(jì)算。利用模型時(shí)需要確定運(yùn)動(dòng)部件的重量、3個(gè)方向的轉(zhuǎn)動(dòng)慣量及慣性矩、重心坐標(biāo)等，以上參數(shù)的獲取可通過Solidworks軟件獲得，在模擬時(shí)加入自行編譯的被動(dòng)型動(dòng)網(wǎng)格UDF進(jìn)行計(jì)算。

2.3 初始條件與邊界條件的設(shè)定

初始條件與邊界條件設(shè)定如下：

1） 進(jìn)、出口相對(duì)壓力為一個(gè)標(biāo)準(zhǔn)大氣壓。

2） 進(jìn)口速度為5 m/s。

3） 進(jìn)口氣體中液滴所占的體積分?jǐn)?shù)為5%。

4） 液滴粒徑按大小分為6組，分別為Bin1～Bin6等級(jí)，其中Bin5設(shè)置為1，即煙氣進(jìn)入除霧器時(shí)的初始液滴粒徑為均勻的1 μm，進(jìn)入后經(jīng)過團(tuán)聚，可能粒徑增長為Bin1（10 μm），Bin2（6.6 μm），Bin3（4.7 μm），Bin4（2.6 μm），也有可能破碎為Bin6（0.3 μm）。而設(shè)為初始液滴粒徑為5 μm時(shí)，則Bin1（50 μm），Bin2（34 μm），Bin3（21 μm），Bin4（12 μm），也有可能破碎為Bin6（2.7 μm）。

5） 選擇示蹤粒子液滴數(shù)量為1 000。

3 模擬結(jié)果及分析

3.1 被動(dòng)式葉輪除霧器內(nèi)部流場(chǎng)分布及固定式除霧器的對(duì)比

除霧器內(nèi)，被動(dòng)式除霧器切向的速度分布如圖2所示?；旌蠚庖簝上噙M(jìn)入除霧器時(shí)平行向前流動(dòng)，可以看做管內(nèi)層流流動(dòng)，當(dāng)氣流向前流動(dòng)時(shí)，由于遇到葉片的阻力而轉(zhuǎn)變?yōu)閺?qiáng)烈的湍流流動(dòng)，并改變流動(dòng)方向旋轉(zhuǎn)前進(jìn)。由圖2可以看出，軸線平面上的切向速度梯度分布，而且由于葉片的作用，使靠近壁面處的混合氣流速度增大，更利于攜帶液滴對(duì)壁面的碰撞和粘連，在壁面形成水膜匯集沿壁面流下，則實(shí)現(xiàn)了除霧的目的。攜帶液滴示蹤粒子軌跡如圖3所示，由圖3可以看出，混合氣流在經(jīng)過葉片后呈旋流前進(jìn)，在垂直軸線的平面上，混合氣流繞軸心作離心運(yùn)動(dòng)，并產(chǎn)生離心力，有利于氣液分離。固定式除霧器切向速度分布如圖4所示。

由圖4可以看出，大部分液滴在穿過葉輪后，被動(dòng)式葉輪下的切向速度比固定式葉輪大。切向速度越大，則液滴受到的離心力越大，因而有利于液滴粒子向容器壁面轉(zhuǎn)移，對(duì)強(qiáng)化液滴與液滴、壁面間的碰撞非常有利，這是液滴在相鄰葉片之間的縫隙繞流時(shí)獲得的切向速度，而自由式葉輪又受到氣流沖撞，旋轉(zhuǎn)時(shí)對(duì)氣流施加反作用，在葉片的帶動(dòng)下獲得更大的切向速度，利于氣液分離。同時(shí)，被動(dòng)式葉輪對(duì)氣流的阻力比固定式葉片小。

由圖3中顆粒的運(yùn)動(dòng)軌跡可以看出，前進(jìn)過程中粒徑較大的顆粒易于被甩向壁面，同時(shí)較小的顆粒之間發(fā)生相互碰撞，在各種力的作用下進(jìn)行布朗團(tuán)聚與湍流團(tuán)聚，更易于提高除霧效率。

3.2 固定式及被動(dòng)式葉片液滴粒徑分布

固定式葉片液滴粒徑分布如圖5所示，被動(dòng)式葉片液滴粒徑分布如圖6所示。由圖5和圖6可以看出，相同粒徑的液滴進(jìn)入除霧器，在通過葉輪后，粒徑逐漸增大。但被動(dòng)自由式葉輪的液滴粒徑比固定式葉片的粒徑更大，且大顆粒所占比例更高，這是因?yàn)闅饬髟谕ㄟ^被動(dòng)自由式葉輪時(shí)的擾動(dòng)更加強(qiáng)烈，而形成更強(qiáng)的旋流，加劇液滴間的碰撞，使液滴間的湍流團(tuán)聚與布朗團(tuán)聚增強(qiáng)，更易形成較大粒徑的液滴顆粒，從而被甩向容器壁面完成氣液分離，說明新型被動(dòng)自由式除霧器在增強(qiáng)旋流程度及液滴團(tuán)聚方面更具優(yōu)勢(shì)。

3.3 新型除霧器的模擬計(jì)算

對(duì)不同粒徑大小及不同流速下的新型除霧器進(jìn)行模擬計(jì)算，了解不同工況下的液滴分離效果及影響氣液分離的因素。當(dāng)液滴粒徑為1 μm和5 μm時(shí)，進(jìn)入除霧器的初始粒徑液滴團(tuán)聚后顆粒粒徑分布如圖7所示。

由圖7可以看出，當(dāng)初始液滴為1 μm時(shí)，流過葉片后Bin6 和Bin5 等所占比例比5 μm粒徑時(shí)要小，并且Bin2和Bin1等較大粒徑液滴顆粒所占比例比粒徑為5 μm時(shí)要大，說明初始小粒徑液滴的氣體在經(jīng)過旋流運(yùn)動(dòng)時(shí)，容易發(fā)生湍流團(tuán)聚與布朗團(tuán)聚，形成更大粒徑的液滴，從而被壁面收集，這也反映了被動(dòng)葉輪除霧器可以更好的脫除較小液滴，具有優(yōu)異的除霧性能。

當(dāng)初始粒徑為1 μm，流速分別為5 m/s和8 m/s時(shí)，流經(jīng)葉片后液滴體積分布如圖8所示。

由圖8 b可以看出，在8 m/s流速下，Bin6和Bin5等小粒徑所占體積分?jǐn)?shù)更小，Bin3，Bin2和Bin1等較大粒徑所占比例加大，說明液滴粒徑整體呈現(xiàn)向較大粒徑轉(zhuǎn)變。這是因?yàn)槿肟跉饬鞯脑龃螅谷~輪被動(dòng)旋轉(zhuǎn)更快，同時(shí)葉輪對(duì)氣流旋流的反作用更強(qiáng)，容器內(nèi)部呈現(xiàn)劇烈的湍流耗散，使細(xì)小液滴之間的相互碰撞更加激烈，因而大大增強(qiáng)了湍流團(tuán)聚效果，把小液滴團(tuán)聚成較大液滴，更利于液滴脫離氣流在壁面匯集，實(shí)現(xiàn)除霧作用。

4 結(jié)束語

本文通過對(duì)固定式和被動(dòng)自由式葉輪除霧器單元管內(nèi)流場(chǎng)進(jìn)行模擬，重點(diǎn)分析了本課題組提出的被動(dòng)自由式葉輪除霧器除霧的特性。由除霧器內(nèi)部流場(chǎng)可以看出，液滴群的最大切向速度在管內(nèi)壁面附近，氣流穿過葉輪發(fā)生明顯的旋流運(yùn)動(dòng)，有利于液滴群間的碰撞和團(tuán)聚;通過對(duì)比兩種除霧器除霧計(jì)算結(jié)果，被動(dòng)自由式葉輪在液滴的切向速度和團(tuán)聚作用上都優(yōu)于固定式葉輪除霧器，驗(yàn)證了新型除霧器的除霧效率更高;當(dāng)初始液滴粒徑分別為1 μm和5 μm時(shí)，粒徑較小的液滴更容易發(fā)生團(tuán)聚，形成較大顆粒而分離脫除;將不同初始流速5 m/s和8 m/s進(jìn)行比較，流速越大，流經(jīng)葉輪后液滴的團(tuán)聚效果越好，使氣液分離的效率更高。該研究為除霧器設(shè)計(jì)及優(yōu)化提供了理論依據(jù)。

參考文獻(xiàn)：

[1] 楊柳， 王世和， 王小明. 濕式脫硫塔除霧器流場(chǎng)的數(shù)值模擬[J]. 華東電力， 2004， 32（10）： 46.

[2] 石振晶， 陶明， 何育東， 等. 噴淋脫硫塔內(nèi)除霧器性能數(shù)值模擬[J]. 熱力發(fā)電， 2016， 45（3）： 9297， 104.

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

[4] 郝雅潔. 濕法煙氣脫硫除霧效率和流場(chǎng)優(yōu)化的數(shù)值模擬研究[D]. 南京： 東南大學(xué)， 2015.

[5] 沈浩. 濕式電除塵器用于控制燃煤煙氣污染物的測(cè)評(píng)[J]. 電力與能源， 2014， 35（1）： 5458.

[6] 張龍龍， 安澤文， 姜旭， 等. 團(tuán)聚液霧化氣泡粒徑分布特性的可視化研究[J]. 化工進(jìn)展， 2016， 35（10）： 30863091.

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

[8] 魏宏鴿， 葉偉平， 柴磊， 等. 濕法脫硫系統(tǒng)除塵效果分析與提效措施[J]. 中國電力， 2015， 48（8）： 3336.

[9] 魏書洲， 張國新， 郝劍， 等. 脫硫吸收塔流場(chǎng)優(yōu)化及研究[J]. 華北電力技術(shù)， 2016（2）： 6570.

[10] 徐文勝， 馬志剛. 濕法脫硫塔入口煙道流場(chǎng)的優(yōu)化分析[J]. 能源工程， 2011， 3（6）： 4851， 57.

[11] 熊偉. 葉輪除霧器氣液分離過程的數(shù)值模擬[D]. 武漢： 武漢科技大學(xué)， 2007.

[12] 郝雅潔. 濕法煙氣脫硫除霧效率和流場(chǎng)優(yōu)化的數(shù)值模擬研究[D]. 南京： 東南大學(xué)， 2016.

[13] 陳明功， 付永強(qiáng)， 袁細(xì)寧. 濕法噴霧脫硫塔內(nèi)霧化器的應(yīng)用[J]. 環(huán)境工程， 2006， 24（4）： 4041， 44.

[14] 孫志春， 郭永紅， 肖海平， 等. 鼓泡脫硫塔除霧器除霧特性數(shù)值研究及實(shí)驗(yàn)驗(yàn)證[J]. 中國電機(jī)工程學(xué)報(bào)， 2010， 30（8）： 6875.

[15] Babai L， Nisant N， Szegedy M. Multiparty protocols， pseudorandom generators for logspace， and timespace tradeoffs[J]. Journal of Computer & System Sciences， 1992， 45（2）： 204232.

[16] 姚杰， 仲兆平， 周山明. 濕法煙氣脫硫帶鉤波紋板除霧器結(jié)構(gòu)優(yōu)化數(shù)值模擬[J]. 中國電機(jī)工程學(xué)報(bào)， 2010， 30（14）： 6167.

[17] 皮威. 折流板除霧器流場(chǎng)數(shù)值模擬[D]. 武漢： 武漢工程大學(xué)， 2011.

[18] Charinpanitkul T， Tanthapanichakoon W. Deterministic model of openspace dust removal system using water spray nozzle： Effects of polydispersity of water droplet and dust particle[J]. Separation & Purification Technology， 2011， 77（3）： 382388.

[19] 趙海波. 離散系統(tǒng)動(dòng)力學(xué)演變過程中的顆粒群平衡模擬（精）[M]. 北京： 科學(xué)出版社， 2008.