衣淑立,周勝名

（金川集團股份有限公司，甘肅金昌737100）

氯氣吸收塔氣液兩相流的數(shù)值模擬研究

衣淑立,周勝名

（金川集團股份有限公司，甘肅金昌737100）

利用Fluent軟件對鎳鈷冶煉系統(tǒng)產(chǎn)生的廢氯氣堿液吸收工藝中的吸收塔進行了模擬計算?？疾炝艘欢ú僮鳁l件下吸收塔內(nèi)氣液兩相流的速度、壓力及氯氣濃度分布。模擬計算的數(shù)據(jù)表明吸收液噴射量的變化對吸收塔出口氯氣濃度影響顯著，而對塔內(nèi)壓降影響不大，與實驗測試值吻合較好，驗證了氯氣吸收塔氣液兩相流數(shù)學模型的正確性。

氯氣吸收塔；數(shù)值模擬；Fluent;吸收液噴射量;

氯氣是一種強氧化劑,主要用于鹽酸、農(nóng)藥、炸藥、有機染料、有機溶劑及化學試劑的制備，有色及稀有金屬冶煉的電解、氯化焙燒、沸騰氯化或氯化揮發(fā)及其它氯化冶金過程等。氯氣又是一種有毒氣體，排入環(huán)境后會對人類、動植物和其他器物造成危害。隨著計算機技術(shù)的迅速發(fā)展，計算流體力學(computational fluid dynamic,CFD)模擬高復雜度的多相流動成為可能。近年來,國內(nèi)外學者應(yīng)用CFD軟件對氣體吸收塔內(nèi)的多相流進行了模擬計算[1-6],研究結(jié)果用于改進吸收塔設(shè)計,取得了較好的效果。某公司鎳、鈷冶煉車間電解陽極有大量的廢氯氣產(chǎn)生,這些廢氯氣通常采用一級或多級吸收塔以堿液吸收[7-9]。對吸收塔內(nèi)氯氣與堿液兩相流進行數(shù)值模擬研究,對改進氯氣吸收塔具有指導意義。本文作者擬采用目前比較流行的商用CFD軟件Fluent對氯氣吸收塔內(nèi)氣液兩相流的流動特性進行數(shù)值模擬研究。

1 數(shù)學模型

1.1 物理模型

為便于建模及計算,對氯氣吸收塔內(nèi)煙氣與吸收液兩相流動狀況作假設(shè)和簡化：將煙氣視為不可壓縮牛頓流體,不考慮塔內(nèi)噴淋管等組件對氣液流場的影響,假設(shè)液滴為球形,不考慮液滴的碰撞、破碎及聚并,暫不考慮煙氣與液滴之間的傳質(zhì)、傳熱和反應(yīng)。采用歐拉法將煙氣處理為連續(xù)相對其進行描述。

1.2 控制方程

基于上述假設(shè)與簡化,對于不可壓煙氣流場的連續(xù)控制方程可表示為：

動量守恒控制方程為：

式中：p為靜壓，F(xiàn)x、Fy、Fz為液滴對煙氣流場的反作用力。根據(jù)塔內(nèi)煙氣流動湍流實際情況,本文采用標準～.湍流模型來模擬塔內(nèi)煙氣湍流運動。

由于噴淋液滴的分散體積率較小(遠＜10%),故采用拉格朗日法來描述液滴的運動較為精確,即選用離散相模型(discrete phase models,DPM)來跟蹤液滴運動。吸收塔內(nèi)的液滴受到諸如重力、曳力、浮力等力的作用,本文主要考慮重力和曳力對液滴的作用。根據(jù)液滴在煙氣流場中受力分析,液滴顆粒的運動方程可表示為：

式中：FD(ug-up)為顆粒受煙氣的單位質(zhì)量曳力，其表達式為：

式（6）中：Re為相對雷諾數(shù)（或顆粒雷諾數(shù)），其定義為：

CD為液滴受煙氣的阻力系數(shù)與Re相關(guān)。出于計算收斂的考慮,當Re≤1 000時采用了工程上應(yīng)用較多的CD取值。Re的取值如下：

由于吸收塔內(nèi)煙氣流動為湍流,故需考慮煙氣湍流對塔內(nèi)液滴運動的影響。本文采用隨機漫步模型(discrete random walk，DRW)來考察湍流對液滴軌跡的影響。

液滴的粒徑分布對氣液兩相流計算影響很大。在塔內(nèi)液滴的直徑大約分布在0.2 mm～5 mm之間。為了較精確模擬液滴在噴淋塔內(nèi)的運動,本文采用Rosin-Rammler模型來描述液滴的粒徑分布：

根據(jù)噴嘴試驗數(shù)據(jù)得出：d=0.002 5，n=3。

由于煙氣粘性和其它性質(zhì)存在，故在吸收塔內(nèi)由噴嘴向下噴出液滴會受到向上流動煙氣的阻力和其它力的作用,因而需要考慮煙氣與液滴之間的作用力。在DPM模型中煙氣(連續(xù)相)與液滴(離散相)之間的作用可表示為：

式中：Fo為其它作用力。

1.3 數(shù)值求解

本文以某冶煉廠廢氯氣吸收塔進行模擬：吸收塔直徑為1.5 m,塔高6 m;煙氣進口距吸收塔底部2 m,管徑400 mm;吸收液噴嘴距吸收塔底部4.5 m,采用動力波噴頭;煙氣出口位于吸收塔縱向軸線方向,管徑400 mm。

模擬計算所用計算機配置為Intel Pentium 2.66 G,DDR RAM 512MB，80 GB,操作系統(tǒng)為Windows XP professional sp2,CFD軟件為Fluent 6.3.13。采用Fluent專用的前置處理器Gambit建立幾何形狀并生成混合網(wǎng)格(圖1),網(wǎng)格數(shù)量約1×105萬個,然后由Fluent進行求解。

圖1 吸收塔網(wǎng)格

模擬計算所用的進口煙氣組分見表1。模擬中改變3個變量：煙氣流量、氯氣深度和吸收液流量,共得到27組模擬結(jié)果。

表1 模擬計算所用參數(shù)

2 結(jié)果與討論

2.1 塔內(nèi)速度場分布

取煙氣流量5 000 Nm3/h,氯氣深度8×104mg/Nm3,吸收液噴淋流量60 m3/h進行模擬計算,得到吸收塔內(nèi)流場的跡線圖、速度矢量圖及不同界面的流體速度分布,分別見圖2(a)、圖2(b)、圖3。

圖2 塔內(nèi)煙氣速度場—跡線圖(a)和速度矢量圖(b)

圖3 沿吸收塔縱向不同截面的速度分布/(m/s)

從圖2(a)可以看出,吸收塔內(nèi)存在兩個較大的漩渦區(qū)域。位于吸收塔底部的漩渦區(qū)域流體呈順時針方向。不同吸收液噴淋速度條件下的模擬都表明，此漩渦區(qū)域不會隨著吸收液噴淋速度的改變而發(fā)生改變。這說明這個漩渦區(qū)域的產(chǎn)生主要是由于煙氣進口對面壁面的阻擋作用引起。位于吸收液噴嘴下方的漩渦區(qū)域呈逆時針方向，這個漩渦的產(chǎn)生主要是由于噴嘴噴出的液滴對煙氣的曳力作用而產(chǎn)生。從傳質(zhì)角度來說,渦旋的存在有利于強化氣液兩相間的傳質(zhì)，從而有利于吸收反應(yīng)的進行。而從能量損失的角度來看,渦旋的產(chǎn)生必然會造成壓損的增加。實際工程當中的理想目標就是尋找兩者最佳結(jié)合點的流場。另外,由于煙氣出口位于吸收塔軸線位置,使得一部分煙氣沒有與吸收液充分混合就直接從出口排出，降低了氯氣的吸收效率。因此，這種出口布局應(yīng)進行改進,如將煙氣進口與出口布置在吸收塔同一側(cè),可改善氣液傳質(zhì)過程。從圖2(b)可以看出,煙氣流速在塔內(nèi)發(fā)生了較大的變化,靠近噴嘴上方煙氣流速有一個極小區(qū)域,速度接近0。塔內(nèi)組件如吸收液噴嘴等對流場的均勻性產(chǎn)生影響,在吸收塔設(shè)計時須充分考慮此因素的影響,通過合理調(diào)整噴嘴的位置可以將此影響降低。從圖3可以看出,隨著高度增加,塔內(nèi)的煙氣流場漸趨均勻,在煙氣進口方向截面上的流場變化最大。

2.2 塔內(nèi)壓力分布

吸收塔內(nèi)壓降變化情況一直是實際工程設(shè)計中比較關(guān)注的方面。圖4是煙氣流量為5 000 Nm3/h、氯氣深度100 000 mg/Nm3條件下塔內(nèi)沿縱向界面的壓力分布圖。

圖4 不同吸收液流量時塔內(nèi)沿縱向截面的壓力分布/Pa

對比圖4(a)、圖4(b)、圖4(c)可以看出,塔內(nèi)的壓力分布變化不大,吸收液噴嘴對壓力分布不產(chǎn)生顯著影響，吸收液流量變化對塔內(nèi)壓降影響很小,不同吸收液噴射量下塔內(nèi)的壓降大約都在130 Pa左右,而構(gòu)成塔內(nèi)壓降的主要區(qū)域是煙氣進口和出口區(qū)域。這說明在煙氣流量不變的情況下,改變吸收液噴射量,對塔內(nèi)壓降的影響不明顯。因此,為了提高氣液傳質(zhì)效率和氯氣吸收率，可以盡量增大吸收液的循環(huán)噴射量。在實際工程設(shè)計與運行過程中，可通過優(yōu)化進出口結(jié)構(gòu)設(shè)計，有效地降低壓損,節(jié)約能耗。

圖5是煙氣流量為5 000Nm3/h、氯氣深度1×105mg/Nm3條件下,不同吸收液流量時壓力沿煙氣進口管軸線方向上的分布,x=0為吸收塔中心位置?？梢钥闯?隨著吸收液流量增加,這些特定方向的壓力有所增加。在煙氣進口方向,煙氣離開進口管進入吸收塔后壓力逐漸上升,在距吸收塔中心處約一半距離時達到最大，此后下降，在吸收塔中心位置達到一個極小值，隨后由于碰到壁面，壓力又迅速上升。

圖5 不同吸收液流量時沿煙氣進口軸線方向的壓力分布

圖6是同樣條件下不同吸收液流量時沿吸收塔縱向軸線方向的壓力分布，x=0為吸收塔底部位置。在沿吸收塔軸線方向,壓力總體上不斷下降,但在進氣管附近,壓力達到一個極大值。在高度4.5 m的位置,由于存在噴頭,產(chǎn)生了一個壓力盲區(qū),使壓力分布產(chǎn)生不連續(xù)。

圖6 不同吸收液流量時沿吸收塔軸線方向的壓力分布

圖7是氯氣深度為1×105mg/Nm3、吸收液流量為60 m3/h時不同煙氣流量下塔內(nèi)沿縱向截面的壓力分布圖。隨著煙氣量增加,塔內(nèi)平均壓力逐漸升高,但壓力分布趨勢變化不大。

圖7 不同煙氣量時塔內(nèi)沿縱向截面的壓力分布

圖8是氯氣深度為1×105mg/Nm3、吸收液流量為60 m3/h時,塔內(nèi)壓降與煙氣量的關(guān)系。隨著煙氣量增加,塔內(nèi)壓降相應(yīng)升高。但在實際的煙氣量范圍內(nèi),這個壓降的增加是可以接受的。另外，氯氣深度對塔內(nèi)壓降的影響不大,這是因為煙氣中氯氣含量較低,因氯氣與吸收液的反應(yīng)而造成的壓力降低較小。

圖8 煙氣量與塔內(nèi)壓降關(guān)系

2.3 氯氣濃度分布

圖9是當煙氣流量為5 000 Nm3/h、氯氣深度為1×105mg/Nm3時,不同吸收液流量條件下塔內(nèi)沿縱截面的Cl2濃度分布。很明顯20 m3/h的吸收液流量是偏低的，煙氣出口處的Cl2濃度仍然較高。隨著吸收液流量增大,Cl2與吸收液的傳質(zhì)和吸收反應(yīng)得到強化,吸收液噴嘴上方的低Cl2濃度區(qū)域相應(yīng)增大,煙氣出口處的Cl2濃度逐漸降低;同時,整個吸收塔內(nèi)Cl2的平均濃度也在下降,Cl2的吸收效率提高。

圖9 不同吸收液流量時塔內(nèi)沿縱截面的Cl2濃度分布

圖10是吸收液流量變化時塔內(nèi)沿著吸收塔軸線方向Cl2濃度的變化曲線?？傮w上來看,增大吸收液噴射量,可使Cl2濃度明顯降低。由于吸收塔軸線穿過吸收液噴嘴，噴嘴區(qū)域無法進行模擬計算,因此曲線圖上靠近噴嘴區(qū)域存在一個濃度分布間斷區(qū)。對比3條曲線可以看出,吸收液流量為20～40 m3/h時,噴嘴以上Cl2濃度有一個較小的增加,這可能是由于在較大的煙氣流量下,由于液體噴射流量較小,在噴嘴上方形成了一個較小的霧化區(qū),這一區(qū)域氣液傳質(zhì)過程受到抑制,氯氣濃度稍有上升。當吸收液噴射量達到60 m3/h時,靠近噴嘴附近,Cl2濃度顯著降低。由于吸收液噴射量的增加對塔內(nèi)壓降不產(chǎn)生顯著影響,為了提高Cl2吸收率,應(yīng)盡可能提高吸收液的循環(huán)噴射量。

圖10 不同吸收液流量時沿吸收塔軸線方向Cl2的濃度分布

圖11是不同煙氣量時塔內(nèi)沿縱截面的Cl2濃度分布。當氯氣深度為1×105mg/Nm3、吸收液流量為20 m3/h時,煙氣流量對Cl2的縱向分布影響不大。但仍然可以看出,煙氣流量的增加對Cl2的濃度還是產(chǎn)生了影響,隨著煙氣流量增加,氯氣濃度亦稍有增加。在這樣的煙氣流量和Cl2深度條件下,20 m3/h的吸收液噴射量不夠的,煙氣出口處Cl2的濃度依然較高。

圖11 不同煙氣流量時塔內(nèi)沿縱截面的Cl2濃度分布

圖12是同樣條件下煙氣流量變化時塔內(nèi)沿著吸收塔軸線方向Cl2濃度的變化曲線?？梢钥闯?，在2.5 m以下,煙氣流量對這個方向的Cl2濃度影響不大,3條曲線非常接近；在2.5 m以上,Cl2濃度稍有差別,煙氣流量增加,Cl2濃度也有所增加。

圖12 不同煙氣流量時沿吸收塔軸線方向Cl2的濃度分布

圖13是不同Cl2深度時塔內(nèi)沿縱截面的Cl2濃度分布。當煙氣流量為5 000 Nm3/h、吸收液流量為40 m3/h時，隨著Cl2深度增加，煙氣進口處的Cl2濃度逐漸增加，整個塔內(nèi)的Cl2濃度也相應(yīng)上升。

圖13 不同Cl2深度時塔內(nèi)沿縱截面的Cl2濃度分布

圖14是同樣條件下Cl2深度變化時塔內(nèi)沿著吸收塔軸線方向Cl2濃度的變化曲線。Cl2深度的提高顯然使塔內(nèi)的Cl2濃度升高。在煙氣進口以下,3條曲線差別較大,而在進口以上,沿吸收塔軸線方向Cl2的濃度差別減小,特別是在出口處,含氯深度8× 104mg/Nm3與1×105mg/Nm3條件下Cl2的濃度幾乎相同。

圖14 不同Cl2深度時沿吸收塔軸線方向Cl2的濃度分布

3 結(jié)論

本文對比分析了氯氣吸收塔內(nèi)不同工況下的流體速度場分布、壓力變化和氯氣濃度分布。塔內(nèi)壓降、出口氯氣濃度等指標的模擬計算結(jié)果與實驗測試值吻合較好,驗證了氯氣吸收塔氣液兩相流數(shù)學模型的正確性。由上述分析可以看出,吸收液噴射量的變化對吸收塔出口氯氣濃度影響顯著,而對塔內(nèi)壓降影響不大,因此可通過增大吸收液噴射量有效提高氣液傳質(zhì)效率,提高氯氣吸收效果。數(shù)值模擬計算結(jié)果為優(yōu)化氯氣吸收塔設(shè)計與工況條件提供了依據(jù)。[J].山東化工,2014(1)：28-30.

[2] 郭欣.低溫甲醇洗吸收塔的計算機模擬[J].煤炭轉(zhuǎn)化,2013(1)：89-92.

[3] 李斌,劉強,吳其榮.180度雙入口脫硫吸收塔的優(yōu)化模擬[J].四川環(huán)境,2012(3)：122-125.

[4] 李艷.高壓水洗沼氣凈化吸收塔的模擬研究[J].可再生能源,2011 (5)：71-74.

[5] 謝書圣.天然氣三甘醇脫水系統(tǒng)吸收塔模擬計算研究[J].計算機與應(yīng)用化學,2011(3)：343-346.

[6] 劉佳,黃偉.合成氣一步法制二甲醚吸收塔的模擬研究[J].煤化工, 2009(3)：7-10.

[7] 鄭紀東.廢氯氣吸收裝置運行總結(jié)[J].中國氯堿,2013(8)：30-32.

[8] 吳培德.氯氣吸收塔的腐蝕與防護[J].中國有色冶金,2011(1)：42-43,48.

[9] 周穎.雙塔處理吸收廢氯氣的探討[J].氯堿工業(yè),2012(7)：18-19.

[1] 湯民波,周新明,王桂清.有風情況下氯氣泄漏擴散數(shù)值模擬研究

Numerical Simulation Research on Gas-liquid Two-phase Flows in a Chlorine Absorption Tower

YI Shuli,ZHOU Shengming

(Jinchuan Group Co.,Ltd.,Jinchang,Gansu 737100,China)

Make numerical simulation on chlorine absorption column for nickelamp;cobalt metallurgical process by using Fluent software.Under the certain operating conditions,the inside tower gas-liquid flow rate,pressure,chlorine distribution were investigated. The data from analog computation shows that the change of absorption liquid emitted dose has great influence on chlorine density of absorption outlet,while it cause a minor impact on tower inside pressure drop.The test result is finely matched with the experimental measurement,which proves the accuracy of mathematical model of gas-liquid two-phase flows in a chlorine absorption tower.

chlorine absorption tower;numeral simulation;Fluent;emitted dose of absorption liquid

TF815；TF816

B

1004-4345(2014)05-0030-06

2014-05-22

衣淑立(1966—),男,高級工程師,主要從事有色冶金管理與技術(shù)研究工作。