翟亞軍，張懷宇，周 密，陳 寧

（1. 中遠(yuǎn)海運(yùn)重工有限公司，上海 200120；2. 江蘇科技大學(xué)，江蘇鎮(zhèn)江 212003）

0 引言

為了有效地控制船舶硫氧化物的排放，1997年 9月在海上環(huán)境保護(hù)委員會(huì)（MEPC）第40屆會(huì)議上通過(guò)了國(guó)際防止船舶造成污染公約（MARPOL）73/781997年議定書（包括MARPOL附則Ⅵ—防止船舶造成空氣污染規(guī)則）及8個(gè)決議案。MARPOL附則Ⅵ第14條嚴(yán)格規(guī)定了船舶燃油含硫量標(biāo)準(zhǔn)和強(qiáng)制執(zhí)行的時(shí)間節(jié)點(diǎn)[1]。自2003年MEPC 50屆會(huì)議至2013年MEPC 65屆會(huì)議，10年共16屆環(huán)保會(huì)議上制定和修改了國(guó)際海事組織（IMO）的MARPOL附則Ⅵ，在全球范圍內(nèi)確定了一些特殊區(qū)域和特別敏感海區(qū)，設(shè)為船舶廢氣排放控制區(qū)（ECA）[2]。因此，脫硫塔作為船用煙氣脫硫系統(tǒng)中的關(guān)鍵設(shè)備正廣泛應(yīng)用于船舶上。但由于船舶上空間有限，傳統(tǒng)噴淋脫硫塔體積過(guò)于龐大，使得脫硫塔在船舶上的安裝運(yùn)行受到很大限制[3-5]。因此，本文提出一種利用旋風(fēng)分離器原理的新型脫硫塔，塔內(nèi)的旋轉(zhuǎn)氣流不僅增大了湍流強(qiáng)度使得氣液混合更充分、反應(yīng)更快速?gòu)氐?，而且大大增加了煙氣在塔?nèi)的行程，使得煙氣在有限尺寸的脫硫塔內(nèi)依然有足夠的反應(yīng)時(shí)間。

1 旋流脫硫塔結(jié)構(gòu)設(shè)計(jì)

旋流脫硫塔的結(jié)構(gòu)參數(shù)包括：塔徑D、煙氣進(jìn)口尺寸W×L、煙氣出口直徑D0、底流口直徑Du、柱段長(zhǎng)度H、煙氣出口插入深度H0、錐角θ，其結(jié)構(gòu)參數(shù)如圖1所示，具體數(shù)值如表1所示。

圖1 旋流脫硫塔結(jié)構(gòu)參數(shù)

表1 旋流脫硫塔結(jié)構(gòu)尺寸

2 氣-液兩相流的數(shù)值模擬

2.1 湍流模型

旋流脫硫塔內(nèi)的流場(chǎng)是一種高雷諾數(shù)、強(qiáng)旋流的流動(dòng)，Bernardo研究發(fā)現(xiàn)各相同性標(biāo)準(zhǔn)的k-ε模型已不能準(zhǔn)確描述強(qiáng)旋流流場(chǎng)特征。Erdal采用雷諾應(yīng)力模型準(zhǔn)確預(yù)測(cè)了旋流器內(nèi)的流場(chǎng)分布。本文采用雷諾應(yīng)力模型（RSM）模擬旋流脫硫塔內(nèi)的強(qiáng)旋湍流流場(chǎng)。雷諾應(yīng)力模型所對(duì)應(yīng)的方程包括：

1）連續(xù)性方程

2）運(yùn)動(dòng)方程

式中：Ui為時(shí)均速度；gi為重力加速度分量；P為介質(zhì)壓力；u′為脈動(dòng)速度；xi為笛卡爾坐標(biāo)分量；ρ為介質(zhì)密度；μ為介質(zhì)黏度。

3）雷諾應(yīng)力運(yùn)輸方程

4）湍動(dòng)能方程

5）湍動(dòng)能耗散方程

2.2 組分傳輸模型

FLUENT軟件提供了通用有限速度模型、非預(yù)混和燃燒模型、預(yù)混和燃燒模型、部分預(yù)混和燃燒模型4種模擬反應(yīng)的方法。本文采用的通用有限速度模型是基于組分質(zhì)量分?jǐn)?shù)的輸運(yùn)方程解，可以根據(jù)定義的化學(xué)反應(yīng)機(jī)制，對(duì)化學(xué)反應(yīng)進(jìn)行模擬。這種模型的應(yīng)用十分廣泛，可包含預(yù)混、部分預(yù)混、非預(yù)混以及燃燒等，主要用于模擬化學(xué)組分混合、輸運(yùn)和反應(yīng)的問(wèn)題。

2.3 三維模型及網(wǎng)格劃分

根據(jù)旋流脫硫塔的結(jié)構(gòu)尺寸，應(yīng)用PROE軟件建立旋流脫硫塔三維結(jié)構(gòu)模型，如圖2所示，再將模型導(dǎo)入ICEM軟件中劃分結(jié)構(gòu)網(wǎng)格，如圖3所示。

圖2 幾何模型

圖3 網(wǎng)格圖

2.4 初始條件和邊界條件設(shè)置

1）模擬的初始條件

采用RSM湍流模型求解瞬態(tài)的雷諾應(yīng)力方程，計(jì)算塔內(nèi)強(qiáng)旋流場(chǎng)；采用能量模型和組分傳輸模型模擬強(qiáng)旋流場(chǎng)內(nèi)的化學(xué)反應(yīng)。煙氣溫度為373 K，其中SO2的體積分?jǐn)?shù)為0.12%。CaCO3溶液在進(jìn)口處被霧化成微小液滴，按照與SO2氣體摩爾比4∶1的比例和煙氣混合進(jìn)入脫硫塔內(nèi)。

2）模擬的邊界條件

入口邊界條件：入口邊界類型選擇速度進(jìn)口（velocity inlet），煙氣流量Q=28.49 kg/s，計(jì)算入口流速為25 m/s。

出口邊界條件：出口邊界類型選擇流動(dòng)出口（outflow）。

壁面邊界條件：計(jì)算壁面的邊界條件采用無(wú)滑移邊界。

3 結(jié)果與分析

3.1 旋流脫硫塔內(nèi)的速度場(chǎng)模擬

在旋流反應(yīng)器的速度場(chǎng)中，切向速度占據(jù)著一個(gè)重要的地位，這不僅因?yàn)榍邢蛩俣仍跀?shù)值上遠(yuǎn)大于其他兩個(gè)方向的速度，而且因?yàn)榍邢蛩俣仁请x心力產(chǎn)生的前提。

為了利于觀察不同截面的切向速度分布，在旋流脫硫塔柱段和錐段分別取2個(gè)截面進(jìn)行分析，即y1=2 m，y2=0.6 m，y3=?1.4 m，y4=?3.4 m，其分布示意圖如圖4所示。

圖4 旋流脫硫塔截面分布示意圖

圖5 中，a）、b）、c）、d）分別為y1、y2、y3、y4旋流脫硫橫截面的切向速度分布圖。由圖5可知：塔內(nèi)流場(chǎng)的切向速度分布基本呈軸對(duì)稱，速度由中心處沿徑向方向隨著半徑的增大而不斷增大，到達(dá)最大值以后又逐漸減小，在塔壁附近降至較低，由此可以看出塔壁處和中心軸線處為切向速度低速區(qū)，旋渦在這兩處形成，塔壁處為自由渦，中心軸線處為強(qiáng)制渦，雙渦交匯處切向速度最大，約在r=0.5 m處取得。

3.2 旋流脫硫塔內(nèi)的壓力場(chǎng)模擬

圖6為z=0 m、y=3 m處的截面靜壓分布云圖。從圖6中z=0 m處的剖面可以看出：靜壓沿軸向由塔頂向塔底，壓力沒(méi)有明顯的變化，等壓的位置幾乎與中心軸線平行，所以可以得出軸向位置對(duì)壓力的影響很?。粡膱D6中y=3 m處的截面可以看出，壓力沿著徑向方向隨半徑的減小而減小，塔壁處壓力最高，中心軸線處壓力最低，這一現(xiàn)象與理論上旋流器中心軸線處易形成低壓氣柱相吻合。低壓區(qū)使得煙氣能夠沿著中心軸線處上升直至塔頂處的煙氣出口排出。

圖5 不同旋流脫硫橫截面的切向速度分布圖

圖6 z=0 m、y=3 m處?kù)o壓分布云圖

3.3 旋流脫硫塔內(nèi)的濃度場(chǎng)模擬

3.3.1 旋流脫硫塔內(nèi)SO2濃度分布

圖7為z=0 m、y=3 m處SO2濃度分布云圖。從圖7中z=0 m處的剖面可以看出：塔內(nèi)柱段SO2濃度變化明顯，從上到下、從塔壁到中心處，SO2濃度急劇減小，說(shuō)明堿液與SO2氣體充分接觸并反應(yīng)，錐段處SO2濃度穩(wěn)定在一個(gè)極低的范圍內(nèi)；從圖7中y=3 m處的截面可以看出：煙氣由進(jìn)口切向進(jìn)入沿塔壁環(huán)流后濃度急劇減少，其變化的趨勢(shì)呈現(xiàn)旋轉(zhuǎn)渦流的形態(tài)，其間，SO2幾乎全部被霧化的堿液吸收，脫除效率達(dá)到了99.5%以上。

圖7 z=0 m、y=3 m處SO2濃度分布云圖

3.3.2 旋流脫硫塔內(nèi)堿液濃度分布

圖8為z=0 m和y=3 m處堿液濃度分布云圖。從圖8中可以看出：塔內(nèi)堿液濃度與SO2濃度變化趨勢(shì)一致。由進(jìn)口切向進(jìn)入塔內(nèi)后沿塔壁環(huán)流后濃度快速減少，變化趨勢(shì)為旋轉(zhuǎn)渦流的形態(tài)。

圖8 z=0 m、y=3 m處堿液濃度分布云圖

4 結(jié)論

為了更好地滿足IMO關(guān)于船舶硫氧化物的排放標(biāo)準(zhǔn)，本文提出一種利用旋風(fēng)分離器原理的新型脫硫塔，并利用FLUENT軟件，采用雷諾應(yīng)力模型和組分傳輸模型對(duì)塔內(nèi)的強(qiáng)旋流場(chǎng)和化學(xué)反應(yīng)進(jìn)行三維數(shù)值模擬，基于模擬結(jié)果可以得到以下結(jié)論：

1）塔內(nèi)流場(chǎng)的切向速度分布基本符合蘭金組合渦特征。塔中塔壁處和中心軸線處切向速度較小，在這兩處形成了旋渦，塔壁處為自由渦，中心軸線處為強(qiáng)制渦，雙渦交匯處切向速度最大，大約是在r=0.5 m處。

2）軸向位置對(duì)壓力的影響很小，塔壁處壓力最高，中心軸線處壓力最低，這一現(xiàn)象與理論上旋流器中心軸線處易形成低壓氣柱相吻合。

3）SO2和堿液經(jīng)進(jìn)口切向進(jìn)入塔內(nèi)以后，濃度沿塔壁進(jìn)行環(huán)流的過(guò)程中快速減少，其變化的趨勢(shì)呈渦流狀；SO2在塔內(nèi)與堿液充分接觸并幾乎被全部吸收，其脫出效率達(dá)到了99.5%以上。

[1] 劉新華. 控制船舶柴油機(jī)排放的最佳方法[J]. 航海技術(shù), 1995(1): 42-46.

[2] 馮春林, 盧明. 柴油機(jī)廢氣排放控制技術(shù)開發(fā)[J].柳州職業(yè)技術(shù)學(xué)院學(xué)報(bào), 2007(3): 68-72.

[3] 馬明志. 船舶柴油機(jī)減排新形勢(shì)與技術(shù)發(fā)展[J].上海造船, 2009(4): 34-37, 41.

[4] 龔為佳, 沈衛(wèi)東, 邵錦平, 等. 水洗消煙系統(tǒng)在柴油機(jī)后處理中的應(yīng)用[J]. 柴油機(jī)設(shè)計(jì)與制造,2010(3): 22-25.

[5] 加拿大優(yōu)加公司. 加拿大優(yōu)加公司陶瓷膜廢氣凈化系統(tǒng)介紹說(shuō)明書[Z]. 2010.