劉含笑，袁建國，郭 鏈，姚宇平，酈建國

(浙江菲達(dá)環(huán)?？萍脊煞萦邢薰?，浙江 諸暨 311800)

?

濕法脫硫裝置除霧器流場及除霧特性的數(shù)值研究

劉含笑，袁建國，郭 鏈，姚宇平，酈建國

(浙江菲達(dá)環(huán)保科技股份有限公司，浙江 諸暨 311800)

基于商業(yè)CFD軟件對折線形除霧器流場及除霧特性進(jìn)行了數(shù)值模擬，計(jì)算中，采用離散相模型模擬兩相流，湍流場采用k-ε兩方程湍流模型計(jì)算，顆粒運(yùn)動(dòng)軌跡采用顆粒云模型計(jì)算，通過調(diào)整入口條件及結(jié)構(gòu)參數(shù)，揭示了除霧器除霧效率和壓力損失的變化規(guī)律，可為用于結(jié)構(gòu)設(shè)計(jì)和優(yōu)化。進(jìn)一步對模型進(jìn)行優(yōu)化和補(bǔ)充，探討了不同時(shí)刻除霧器流場及顆粒運(yùn)動(dòng)的瞬時(shí)特性，對于深入了解其除霧過程具有一定的參考價(jià)值。

濕法脫硫；除霧器；流場；顆粒運(yùn)動(dòng)軌跡；數(shù)值模擬

Foundation items：National High-tech R&D Program of China (2013AA065002)

《火電廠大氣污染物排放標(biāo)準(zhǔn)》(GB13223-2011)進(jìn)一步降低了燃煤電廠大氣污染物的排放限值，其中重點(diǎn)控制地區(qū)，要求煙塵排放限值20 mg/m3，2014年7月1日起現(xiàn)役機(jī)組均執(zhí)行此標(biāo)準(zhǔn)。由于環(huán)境容量有限等原因，長三角、珠三角等地部分燃煤電廠已參考燃?xì)廨啓C(jī)組標(biāo)準(zhǔn)限值，要求在基準(zhǔn)氧含量6%條件下，煙塵、SO2、NOx排放濃度分別不高于5、35、50 mg/m3，即燃煤電廠大氣污染物需達(dá)到超低排放的要求。國家發(fā)改委、環(huán)保部和能源局三部委聯(lián)合印發(fā)的《煤電節(jié)能減排升級(jí)與改造行動(dòng)計(jì)劃(2014—2020年)》中要求東部等地區(qū)新建燃煤機(jī)組排放達(dá)到或接近燃機(jī)標(biāo)準(zhǔn)，進(jìn)一步促進(jìn)煤炭的清潔高效利用。

隨著行業(yè)內(nèi)外超低排放呼聲的愈演愈烈，業(yè)內(nèi)對于濕法脫硫裝置(以下簡稱WFGD)的協(xié)同除塵效果的要求也越來越高。據(jù)報(bào)道，國外煙氣協(xié)同治理技術(shù)路線中，WFGD的協(xié)同除塵效果均可達(dá)到70%以上[1-3]。增強(qiáng)WFGD協(xié)同除塵效果可從提高WFGD煙塵脫除效率和減少石膏漿液攜帶兩個(gè)方面來考慮，而增強(qiáng)除霧器除霧效果是減少WFGD石膏漿液攜帶的關(guān)鍵。除霧器內(nèi)流場的流動(dòng)狀態(tài)復(fù)雜，傳統(tǒng)的實(shí)驗(yàn)方法雖然可以粗略測到除霧器的霧滴脫除效率，但難以確定最佳除霧器結(jié)構(gòu)尺寸，即同時(shí)具有較低系統(tǒng)壓力和較高霧滴脫除效率的除霧器葉片幾何尺寸，且極難觀測除霧器內(nèi)液滴的運(yùn)動(dòng)及捕集過程。因此，很有必要發(fā)展一種有效的數(shù)值計(jì)算方法[4-7]，較為準(zhǔn)確預(yù)測除霧器最佳結(jié)構(gòu)形式及參數(shù)配置。

1 除霧器主要類型

WFGD布置的除霧器一般為折流板型除霧器，其利用慣性，通過液滴與除霧器壁面碰撞而將其捕集。通常, 折流板除霧器的板間距為20～30 mm，且煙氣流速不宜過高，氣速過高會(huì)引起二次攜帶，當(dāng)垂直安置時(shí)，氣體的平均流速一般為2～3 m/s，而水平放置時(shí)，氣體的流速可以高些，一般為6～10 m/s。

除霧器葉片是除霧器中最核心的部件，其材料一般由高分子材料(如聚丙稀、FRP等)或不銹鋼(如317L)兩大類材料制作而成。形狀大致可分為折線型和流線型兩類。

2 理論模型

2.1 假設(shè)條件

除霧器在實(shí)際運(yùn)行過程中，液滴及氣流在彎曲的除霧器通道內(nèi)的運(yùn)動(dòng)過程較為復(fù)雜，屬于非定常、黏性可壓縮流體的流動(dòng)過程?？紤]計(jì)算機(jī)內(nèi)存量及計(jì)算時(shí)間的復(fù)雜度，在誤差允許的范圍內(nèi)，根據(jù)實(shí)際情況對模型進(jìn)行適當(dāng)簡化和假設(shè)。

(1) 以空氣代替煙氣進(jìn)行計(jì)算，且流速較低，可假設(shè)為不可壓縮流體。

(2) 計(jì)算為定常工況，并忽略溫度對流場的影響。

(3) 以水滴代替石膏漿液進(jìn)行計(jì)算，并假設(shè)為球形顆粒。

(4) 忽略顆粒的凝聚、破碎、沉降、變形及相變等。

(5) 液滴碰到除霧器壁面即認(rèn)為被捕集，到達(dá)出口時(shí)，認(rèn)為是逃逸，忽略石膏漿液的二次攜帶。

2.2 連續(xù)相計(jì)算模型

氣體相的控制方程可用雷諾平均后的動(dòng)量及連續(xù)性方程，考慮到模型入口為射流,因此選擇處理射流問題較好的可實(shí)現(xiàn)k-ε兩方程湍流模型求解本模型。

連續(xù)方程：

(1)

動(dòng)量守恒方程：

(2)

2.3 顆粒相(液滴)計(jì)算模型

除霧器中液滴相的體積分?jǐn)?shù)小于10%，采用拉格朗日法來描述液滴顆粒的運(yùn)動(dòng)較為準(zhǔn)確，將其作為惰性離散相顆粒處理，即選用DPM模型來跟蹤顆粒運(yùn)動(dòng)，利用顆粒云模型計(jì)算湍流擴(kuò)散對射流顆粒的影響，顆粒相的作用力平衡方程在直角坐標(biāo)系下的形式為：

(3)

式中FD為連續(xù)相阻力；G為重力；Fsaffman為saffman力；Fother為其它相間作用力，在本計(jì)算中不作考慮。

3 物理模型

WFGD內(nèi)除霧器各工作區(qū)域具有一定的對稱性和空間重復(fù)性，忽略來流非均勻性、非對稱結(jié)構(gòu)對除霧器各工作區(qū)域的影響，可將幾何模型簡化為二維結(jié)構(gòu)，以折線形除霧器為例進(jìn)行試算，其幾何模型及主要計(jì)算參數(shù)如表1所示。網(wǎng)格劃分如圖1所示，均為結(jié)構(gòu)化網(wǎng)格。

表1 幾何模型及主要計(jì)算參數(shù)

圖1 網(wǎng)格劃分

4 數(shù)值計(jì)算

4.1 流場分析

采用上述計(jì)算方法，通過CFD計(jì)算軟件對圖1網(wǎng)格文件進(jìn)行計(jì)算，并將計(jì)算結(jié)果進(jìn)行后處理，分別以a=30 mm，b=30 mm，c=20 mm，α=90°，入口流速為2 m/s，液滴直徑為10 μm和a=30 mm，b=30 mm，c=20 mm，α=90°，入口流速為3 m/s，液滴直徑為30 μm為例，計(jì)算除霧器內(nèi)顆粒運(yùn)動(dòng)軌跡及速度、壓力云圖如圖2、3所示。

圖2 顆粒軌跡及流場計(jì)算結(jié)果(a=30 mm，b=30 mm，c=20 mm，α=90°，入口流速為2 m/s，液滴直徑為10 μm)

圖3 顆粒軌跡及流場計(jì)算結(jié)果(a=30 mm，b=30 mm，c=20 mm，α=90°，入口流速為3 m/s，液滴直徑為30 μm)

由圖2(a)所示，除霧器中部分顆粒跟隨流體流出除霧器通道，而一部分顆粒碰到除霧器壁面而被捕集，且大部分捕集區(qū)域發(fā)生在折角處(a圖的框處)，此處連續(xù)相的的流速也較高(b圖的框處)，容易發(fā)生氣液分離；除霧器的壓力損失主要發(fā)生在除霧器折角處(c圖的框處)；圖3(a)的顆粒已被完全捕集，表明除霧器對于大粒徑、高流速的液滴具有較好的捕集效果。

通過流場及顆粒運(yùn)動(dòng)軌跡的分析，可較直觀的了解除霧器的除霧原理。

4.2 除霧效率的影響參數(shù)

4.2.1 除霧效率計(jì)算

本計(jì)算均為單一粒徑工況計(jì)算，為計(jì)算方便，此處除霧效率暫以液滴數(shù)濃度脫除率代替，即進(jìn)口顆粒個(gè)數(shù)為n個(gè)，除霧器出口液滴個(gè)數(shù)為m個(gè)，則被除霧器壁面捕獲的顆粒個(gè)數(shù)為(n-m)個(gè)，則此時(shí)除霧器的除霧效率計(jì)算公式為：

(4)

4.2.2 液滴直徑對除霧效率的影響

以結(jié)構(gòu)參數(shù)a=30 mm，b=30 mm，c=30 mm，α=90°除霧器為例，計(jì)算不同入口流速條件下液滴直徑與除霧效率的關(guān)系如圖4所示。

圖4 除霧效率—液滴直徑關(guān)系曲線

由圖4可知，在相同入口流速時(shí)，液滴粒徑越大，除霧器對其脫出效率越高，最高可達(dá)100%，這是因?yàn)橐旱沃睆皆酱?，其慣性越大，在煙氣中運(yùn)動(dòng)的響應(yīng)時(shí)間也就越長，即液滴運(yùn)動(dòng)的跟隨性變差，因此更容易碰到除霧器壁面而被捕獲。當(dāng)流速為2～4 m/s時(shí)，對于粒徑在10～50 μm范圍內(nèi)的液滴脫除率均達(dá)不到100%，即此時(shí)除霧器的極限粒徑大于50 μm。當(dāng)流速大于5 m/s時(shí)，除霧器極限粒徑大約為30 μm。

4.2.3 煙氣流速對除霧效率的影響

以結(jié)構(gòu)參數(shù)a=30 mm，b=30 mm，c=30 mm，α=90°除霧器為例，不同液滴直徑時(shí)除霧器入口流速與除霧效率的關(guān)系如圖5所示。

由圖5可知，對于50 μm以上的液滴，幾乎已被完全捕獲。粒徑為20、30 μm的液滴，在低流速情況下除霧效果就已經(jīng)較高。隨著流速升高，除霧效果也迅速提升，很快便達(dá)到或接近100%，這是因?yàn)榱魉僭礁撸旱坞x心慣性力越大，越容易撞擊到除霧器壁面。對于10 μm顆粒，隨著流速增加，除霧效果時(shí)增時(shí)減，推測是因?yàn)榱魉僭黾?，流場湍流?qiáng)度增強(qiáng)，亞微觀粒子(1～10 μm)所受的Saffman升力作用更加明顯[8-9]，顆粒運(yùn)動(dòng)軌跡變得不穩(wěn)定，從而導(dǎo)致除霧效率變化不定。

圖5 除霧效率—入口流速關(guān)系曲線

4.2.4 葉片間距對除霧效率的影響

分別以直徑10 μm和20 μm的液滴顆粒為例，除霧器葉片間距對除霧效率的影響如圖6所示。

圖6 除霧效率—葉片間距關(guān)系曲線

由圖6可知，除霧器葉片間距越大，其除霧效果越差，這是因?yàn)槿~片間距增大，氣流在通道中的流通面積變大，在煙氣流方向速度變化趨于平緩，顆粒對氣流的跟隨性變好，不容被除霧器壁面捕獲。

4.2.5 折彎角度對除霧效率的影響

分別以直徑10 μm和20 μm的液滴顆粒為例，除霧器葉片折彎角度對除霧效率的影響如圖7所示。由圖7可知，除霧器葉片折彎角度越大，其除霧效果越差，這是因?yàn)檎蹚澖嵌茸兇?，氣流在通道?nèi)的沿來流方向發(fā)生偏轉(zhuǎn)的角度變小，顆粒的跟隨性變好，不容易被除霧器壁面捕獲。

圖7 除霧效率—折彎角度關(guān)系曲線

4.3 除霧器壓力損失的影響參數(shù)

壓力損失是除霧器的重要性能參數(shù)，壓力損失越大，系統(tǒng)能耗越高，經(jīng)濟(jì)性也就越差。除霧器的壓力損失主要表現(xiàn)在除霧器的沿程阻力損失和折角處的局部壓力損失，如圖2、3中圖(c)所示。

4.3.1 流場特性(流速、粒徑)

以結(jié)構(gòu)參數(shù)a=30 mm，b=30 mm，c=30 mm，α=90°除霧器為例，計(jì)算不同入口流速、液滴直徑條件下壓力損失如圖8所示。液滴直徑對除霧器的壓力損失幾何沒有影響，這是因?yàn)橛?jì)算中未考慮液滴對連續(xù)相的反作用。隨著流速的增加，除霧器的壓力損失明顯增加，這是因?yàn)閴毫p失與速度的平方呈正比。

圖8 不同流速、液滴粒徑時(shí)除霧器壓力損失

4.3.2 結(jié)構(gòu)特性(葉片間距、折彎角度)

以結(jié)構(gòu)參數(shù)a=30 mm，b=30 mm，α=90°除霧器為例，計(jì)算不同葉片間距下壓力損失如圖9所示。由圖9可知，葉片間距越大，除霧器壓力損失越小，且流速越高，不同間距除霧器的壓力損失差別越大。

以結(jié)構(gòu)參數(shù)a=30 mm，b=30 mm，c=30 mm除霧器為例，計(jì)算不同葉片折彎角度下壓力損失如圖10所示。由圖10可知，折彎角度越大，除霧器壓力損失越小，且流速越高，不同間距除霧器的壓力損失差別越大。

圖9 不同葉片間距時(shí)除霧器壓力損失

圖10 不同折彎角度時(shí)除霧器壓力損失

5 計(jì)算模型的進(jìn)一步優(yōu)化及補(bǔ)充

上述計(jì)算結(jié)果可較好反應(yīng)除霧器內(nèi)流場及顆粒運(yùn)動(dòng)的特性，且計(jì)算規(guī)律與相關(guān)文獻(xiàn)的結(jié)果具有較好的一致性[4，7，10]，可為除霧器結(jié)構(gòu)及參數(shù)優(yōu)化提供參考。但隨著理論研究的不斷深入，僅時(shí)均、定性的規(guī)律已經(jīng)不能滿足需求，人們越來越關(guān)注流場、顆粒運(yùn)動(dòng)的瞬時(shí)特性，本文在原有研究基礎(chǔ)上進(jìn)步一對計(jì)算模型進(jìn)行適當(dāng)優(yōu)化及補(bǔ)充，探討不同時(shí)刻除霧器內(nèi)流場及顆粒運(yùn)動(dòng)特性。

5.1 理論模型優(yōu)化及補(bǔ)充

為了更細(xì)致的模擬不同時(shí)刻流場的湍流特性，引入2d-LES模型計(jì)算湍流流場[11]；兩相流模型依然采用PBM模型，為了考慮瞬時(shí)湍流強(qiáng)度對顆粒運(yùn)動(dòng)的影響，改用隨機(jī)軌道模型計(jì)算顆粒運(yùn)動(dòng)軌跡；液滴顆粒粒徑從1～100 μm分為10個(gè)區(qū)間，符合R-R分布：

R(dk)=exp[-(dk/d)n]

(5)

R(dk)表示尺寸大于dk的那部分顆粒占所有顆粒群的質(zhì)量百分?jǐn)?shù)，為R(dk)=e-1時(shí)的顆粒直徑。顆粒粒徑分布如圖11所示。

圖11 顆粒質(zhì)量百分?jǐn)?shù)隨粒徑變化

考慮液滴除霧器內(nèi)運(yùn)動(dòng)的過程中顆粒的凝并及二次破碎，液滴碰撞即認(rèn)為凝并，液滴破碎模型選用泰勒類比破碎模型。

5.2 計(jì)算結(jié)果及分析

采用非定常計(jì)算方法，時(shí)間步長設(shè)為0.001 s，以結(jié)構(gòu)參數(shù)a=30 mm，b=30 mm，c=30 mm，α=120°除霧器為例，經(jīng)迭代計(jì)算，不同時(shí)刻除霧器內(nèi)速度場、顆粒運(yùn)動(dòng)軌跡的計(jì)算結(jié)果如圖12～圖14所示。

圖12 速度云圖(t=0.01～0.10 s)

圖13 顆粒運(yùn)動(dòng)軌跡(t=0.01～0.10 s)

圖14 速度云圖及顆粒運(yùn)動(dòng)軌跡(t=0.15 s)

圖12、圖13分別為t=0.01、0.03、0.06、0.10 s時(shí)刻除霧器內(nèi)速度云圖和顆粒運(yùn)動(dòng)軌跡圖。由圖12可知，除霧器葉片折角處有明顯的速度擾流，而且該擾流隨著時(shí)間不斷變化。圖14為t=0.15 s時(shí)刻除霧器內(nèi)速度云圖和顆粒運(yùn)動(dòng)軌跡，其中(c)圖為對(b)圖線框內(nèi)部分的局部放大。由圖14可知，液滴的主要捕集區(qū)域?yàn)檎劢巧媳诿嫣帲掖笠旱屋^容易被捕集，紅色液滴出口已經(jīng)看不到，即已經(jīng)被完全捕集，而藍(lán)色液滴(小液滴)仍有較多逃逸出去。

6 結(jié)論

運(yùn)用離散相模型(DPM)對不同工況條件下折線形除霧器的流場及顆粒運(yùn)動(dòng)軌跡進(jìn)行了數(shù)值模擬研究。結(jié)果表明：煙氣流速、液滴直徑與除霧效率正相關(guān)，即流速越高、液滴執(zhí)行越大，除霧效果越好；除霧器葉片間距和折彎角度與除霧效率成反比；液滴直徑與除霧器壓力損失無關(guān)，流速越大，葉片間距、折彎角度越小，除霧器壓力損失越大，能耗越高，經(jīng)濟(jì)性也就越差。計(jì)算結(jié)果與文獻(xiàn)結(jié)果吻合較好，除霧效率和壓力損失時(shí)除霧器設(shè)計(jì)和應(yīng)用中的重要參數(shù)，本文展示的計(jì)算方法和結(jié)果可為除霧器的選型設(shè)計(jì)提供有益參考。

另外，引入2d-LES模型模擬湍流流場，隨機(jī)軌道模型計(jì)算顆粒運(yùn)動(dòng)軌跡，入口液滴粒徑分布定義為R-R分布，并考慮液滴的凝并、破碎，優(yōu)化、補(bǔ)充后的計(jì)算模型可得到實(shí)驗(yàn)極難觀測到的不同時(shí)刻流場及顆粒運(yùn)動(dòng)的瞬時(shí)特性，對于研究人員深入了解除霧器除霧過程具有一定的參考價(jià)值。

[1]Yoshio Nakayama，Satoshi Nakamura，Yasuhiro Takeuchi，et al. MHI High Efficiency System-Proven technology for multi pollutant removal[R]. Hiroshima Research & Development Center, 2011: 1-11.

[2]藤島英勝，土屋喜重，西田定二，等. 低低溫EP適用による石炭火力用排煙処理システムの合理化設(shè)計(jì)[J]. 三菱重工技報(bào),1994,31(4):247-251.

[3]名嶋慎司. 石炭火力用低低溫電気集塵裝置[J].住友重機(jī)械技報(bào), 2001,146:35-38.

[4]王 霄，閔 健，高正明，等．脫硫吸收塔除霧器性能的實(shí)驗(yàn)研究和數(shù)值模擬[J]．環(huán)境工程學(xué)報(bào)，2008，11(3)：1529-1534．

WANG Xiao，MIN Jian，GAO Zheng-ming，et al．Experimental study an numerical simulation of characteristics of demister of absorbing tower in desulfurization system[J]．Chinese Journal of Environmental Engineering，2008，11(3)：1529-1534．

[5]楊 柳，王世和，王小明. 濕式脫硫塔除霧器流場的數(shù)值模擬[J]. 華東電力, 2004，32(10)：4-6．

YANG Liu ，WANG Shi-he ，WANG Xiao-ming. Numerical sim ulation of flow field for mist eliminator in wet desulphurization tower[J]. East China Electric Power, 2004，32(10)：4-6．

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

YAO Jie，ZHONG Zhao-ping，ZHOU Shan-ming. Numerical Simulation on Wave-plate Demister with Hamulus in Wet Flue Gas Desulfurization[J]. Proceedings of the CSEE, 2010， 30(14): 61-67.

[7]王政允. 濕法脫硫系統(tǒng)除霧器除霧特性的實(shí)驗(yàn)研究與數(shù)值模擬[D]. 北京：華北電力大學(xué)，2009.

[8]周力行. 多相湍流反應(yīng)流體力學(xué)[M]. 北京：國防工業(yè)出版社, 2002.

[9]孔 瓏. 兩相流體力學(xué)[M]. 北京：高等教育出版社, 2004.

[10]黃新長. 脫硫除霧器性能實(shí)驗(yàn)研究與優(yōu)化設(shè)計(jì)[D]. 北京：華北電力大學(xué)，2009.

[11]酈建國，劉含笑，姚宇平. 微顆粒捕集增效裝置二維單擾流柱流場計(jì)算[J]. 華東電力, 2013, 41(11): 2404-2407

LI Jian-guo, LIU Han-xiao, YAO Yu-ping. The Flow Field Calculation of 2D Single Turbulence Column in the PM2.5 Capture Efficiency Device[J]. East China Electric Power, 2013, 41(11): 2404-2407.

(本文編輯：趙艷粉)

Numerical Research on Flow Field and Demisting Features of Wet Desulphurization Device Demister

LIU Han-xiao, YUAN Jian-guo, GUO Lian, YAO Yu-ping, LI Jian-guo

(Zhejiang Feida Enviro-tech Co., Ltd., Zhuji 311800, China)

Based on the commercial CFD software, numerical simulation was conducted for foldline-shaped demister flow field and demisting characteristics. In the calculation, the discrete phase model simulated two-phase flow; turbulent flow field used k-ε two-equation turbulence model; particle trajectory was calculated by using particle cloud model; by adjusting the population condition and structure parameters, the demister demisting efficiency and pressure loss were revealed. This research can be used for structural design and optimization, and further optimize and complement the model. It also discusses demister instantaneous characteristics of flow field and particle movement at the different time, providing a certain reference valu for further research into the demisting processe.

wet desulphurization; mist eliminator; the flow field; particle trajectory; numerical simulation

10.11973/dlyny201506027

國家高技術(shù)研究發(fā)展計(jì)劃(2013AA065002)

劉含笑(1987)，男，碩士，從事PM2.5捕集增效技術(shù)研發(fā)工作。

X51；X773

B

2095-1256(2015)06-0863-06

2015-08-13