旋流式鼓泡塔內(nèi)部流場特性的實(shí)驗?zāi)M

曹喬喬，郝惠娣，王 瑾

工藝與裝備

旋流式鼓泡塔內(nèi)部流場特性的實(shí)驗?zāi)M

曹喬喬，郝惠娣，王 瑾

（西北大學(xué) 化工學(xué)院, 陜西 西安 710069）

利用CFD技術(shù)對新型旋流式鼓泡的流場進(jìn)行了模擬，考察了不同入口氣速為2,3,4,5 m/s對鼓泡塔內(nèi)流場特性的影響，并對其速度云圖，速度矢量圖及液體體積分?jǐn)?shù)圖進(jìn)行了分析。研究發(fā)現(xiàn)：當(dāng)入口氣速為4 m/s時，攪拌槽內(nèi)的混合效果較好。

新型旋流；鼓泡塔； 流體；CFD模擬

鼓泡塔廣泛應(yīng)用于石油化學(xué)工程、生物工程、環(huán)境以及能源工業(yè)領(lǐng)域中，如加氫、費(fèi)托合成、污水處理、煤液化等［1］。近年來，大量的科研工作者通過各種實(shí)驗?zāi)M研究了鼓泡塔反應(yīng)器內(nèi)的流體特性[2-8]。本文利用CFD技術(shù)中的FLUENT計算軟件對入口氣速分別為2,3,4,5 m/s時的旋流式鼓泡反應(yīng)器流場進(jìn)行了模擬，研究結(jié)果對鼓泡塔的優(yōu)化設(shè)計及工業(yè)放大有實(shí)際的意義。

1 旋流式鼓泡反應(yīng)器的結(jié)構(gòu)

新型旋流式鼓泡塔的總體結(jié)構(gòu)如圖1 所示。旋流式鼓泡反應(yīng)器基本由筒體，氣體分布器，氣液分離器三部分組成。

圖1 總體結(jié)構(gòu)圖Fig.1 The overall structure

2 數(shù)值模擬

2.1 建立模型

模擬對象為圓柱形旋流式鼓泡三相反應(yīng)器,它的內(nèi)徑為800 mm、高為1 400 mm。主要相為液相—水；次要相為氣相—空氣，噴射鼓氣管直徑為20 mm。噴氣徑與反應(yīng)器直徑的比值（環(huán)徑比）為0.4[9]。

2.2 網(wǎng)格劃分（圖2）

圖2 旋流式鼓泡塔的網(wǎng)格劃分Fig.2 The mash of the swirling bubbling reactor

2.3 邊界條件

在本算例中，邊界條件設(shè)定包括壁面（wall）、速 度 入 口 （ velocity-inlet）、 壓 力 出 口（pressure-outlet）。邊界條件的設(shè)定[10]，對于復(fù)雜流場的收斂性有著決定性的影響。模擬變量邊界條件為：入口為氣相，在模擬過程中，入口的湍流強(qiáng)度（Turbulent Intensity）為5%。水力直徑（Hydraulic Diameter）由模型中噴氣管的數(shù)量不同而不同；出口為壓力出口，壓力為表壓為0 Pa。出口的湍流強(qiáng)度（Turbulent Intensity）為1%，水力直徑（Hydraulic Diameter）為800 m；壁面定義為無滑移邊界條件。

3 模擬結(jié)果分析

本文通過CFD技術(shù)中的Fluent對旋流式鼓泡塔反應(yīng)器流場進(jìn)行數(shù)值模擬計算，再通過 CFD-Post對其結(jié)果進(jìn)行后處理[11]，得出不同氣速下反應(yīng)器中的速度矢量圖，速度云圖, 水的體積分?jǐn)?shù)云圖等并對其結(jié)果進(jìn)行分析比較，確定最佳的入口氣速。

圖3所示為不同入口氣速情況下的鼓泡塔中流體矢量圖分析可以看出：當(dāng)速度為2 m/s時，流體速度在反應(yīng)器底部向兩邊呈漩渦狀分布，漩渦附近流體速度分布密度較大，反應(yīng)器中部向上流體速度分布趨于均勻，旋流效果較差，流體速度為0.9 m/s左右;當(dāng)速度為3 m/s時，流體速度在反應(yīng)器底部氣體入口處向兩邊呈漩渦狀分布，氣體入口和出口處的流體速度較大且分布較為密集，反應(yīng)器整體的流體速度有所增大，為1.2～5.5 m/s;當(dāng)速度為4 m/s時，流體速度矢量圖類似速度為3 m/s時的速度矢量圖，只是流體速度增大了，為4.5～16 m/s；當(dāng)速度為5 m/s時，速度矢量圖類似氣速速度為4 m/s時的速度矢量圖，整體速度增大，為8～23 m/s。結(jié)果分析，可以看出鼓泡塔氣速為4 m/s時，流體速度分布較好。體速度有所增大，為1.2～5.5 m/s;當(dāng)速度為4 m/s時，流體速度矢量圖類似速度為3 m/s時的速度矢量圖，只是流體速度增大了，為4.5～16 m/s；當(dāng)速度為5 m/s時，速度矢量圖類似氣速速度為4 m/s時的速度矢量圖，整體速度增大，為8～23 m/s。結(jié)果分析，可以看出鼓泡塔氣速為4 m/s時，流體速度分布較好。

圖3 流體速度矢量圖Fig.3 The water velocity vector

圖4 所示為不同入口氣速下的鼓泡塔中速度云圖分析可以看出：當(dāng)速度為2 m/s時，僅反應(yīng)器氣體入口處及其附近有速度，為0.3～1.7 m/s；當(dāng)速度為3 m/s時，由于噴管入口處有氣體噴入，因此反應(yīng)器入口處速度較大，且速度由內(nèi)向外依次減小，由于反應(yīng)器出口處壓力為 0，反應(yīng)器內(nèi)壓力較大，氣體噴出時速度會增大，反應(yīng)器中部速度變小，反應(yīng)器內(nèi)整體速度大小為 0.2～5.5 m/s ;當(dāng)速度為4 m/s時，反應(yīng)器速度矢量圖類似氣速為3 m/s的速度矢量圖，反應(yīng)器中部速度達(dá)到6 m/s左右，整體速度為2～17 m/s；當(dāng)速度為5 m/s時，速度云圖類似氣速速度為4 m/s時的速度云圖，整體速度增大速度為 0.5～27 m/s,速度波動范圍過大。 結(jié)果分析可以看出鼓泡塔的入口氣速為 4 m/s時，流體速度分布效果較好。

圖4 速度云圖Fig.4 The air velocity contour

圖5 所示為不同入口氣速情況下的鼓泡塔中水的體積分?jǐn)?shù)云圖分析可以看出：當(dāng)速度為2 m/s時，反應(yīng)器出口處有一小部分部分水的含量很高，為0.9左右，但其它部分水與空氣的含量很均勻，水的含量為 0.4～0.5左右，反應(yīng)器內(nèi)水的含量整體分布均勻；當(dāng)速度為3 m/s時，反應(yīng)器出口處有一小部分部分水的含量很高，為0.9左右，由于噴氣管入口處進(jìn)氣，氣速增大了，所以氣體入口處水的含量較低，反應(yīng)器其它地方水與空氣含量較為均勻，水的含量為 0.4～0.5左右；當(dāng)速度為4 m/s時，由于噴氣管入口處進(jìn)氣，所以氣體入口處水的含量較低，水的體積分?jǐn)?shù)從反應(yīng)器內(nèi)部向外依次增大，整體水的含量稍大；當(dāng)速度為5 m/s時，由于噴氣管入口處進(jìn)氣速度偏大，所以氣體入口處水的含量偏低，水的含量在反應(yīng)器中從內(nèi)向外增加，且一邊水的含量偏大于另一邊，水分含量不太均勻。 結(jié)果分析可以看出鼓泡塔的入口氣速為2 m/s時，流體速度分布效果最佳，氣速為4 m/s時，流體速度分布效果次之。

圖5 水的體積分?jǐn)?shù)云圖Fig.5 The water volume fraction contour

4 不同氣速模擬結(jié)果分析

如圖6所示為不同氣速的軸心速度曲線，由圖可以看出：當(dāng)氣速為2 m/s時，整體氣速非常小，不利于氣液混合；當(dāng)氣速為3 m/s時，中心軸的整體速度比較小氣液混合不均勻；速度為4 m/s時速度最高可達(dá)12 m/s左右，氣液混合效果比較理想；速度為5 m/s時，速度減小，有點(diǎn)偏低氣液混合不理想。

圖6 不同氣速的軸心流體速度曲線Fig.6 The fluid velocity curve of axis

5 結(jié) 論

利用FLUENT計算軟件對旋流式鼓泡塔反應(yīng)器中氣速分別為2,3,4,5 m/s時的流場進(jìn)行了模擬，并且利用圖表及各種分析方法進(jìn)行了定性和定量的分析。得到了下面結(jié)果：

（1）鼓泡塔的入口氣速為4 m/s時的速度矢量和速度云圖都比其它入口速度的效果好些，入口氣速為2 m/s時的水的體積分?jǐn)?shù)云圖效果最佳，但氣速過小不利于氣液混合，入口氣速為4 m/s效果次之。從整體效果來看，取鼓泡塔入口氣速為4 m/s時最佳。

（2）從不同氣速的軸心流體速度曲線可以看出氣速為4 m/s時氣液混合效果比較理想。

［1］Sanyal J,Vasquez S,Roy S,et al. Numerical simulation of gas-liquid dynamics in cylindrical bubble column reactors［J］.Chem Eng Sci,1999,54(21):5071 － 5083.

［2］葉澤月, 張麗娟, 石康生,等. 分布器對鼓泡塔中氣泡直徑分布的影響[J].化工生產(chǎn)與技術(shù) , 2007:14(4) 43-46.

［3］王樹立,孫桂大,解 茂 昭.鼓泡塔內(nèi)氣液兩相流的研究[J].石油化工,1999,28 (10):676-681.

［4］林松,李良超,王嘉駿,等.鼓泡塔中氣泡尺寸分布和局部氣含率研究[J].化學(xué)工程,2008,36(2):21-24.

［5］樊亞超,靳海波,秦嶺,等.多級鼓泡塔內(nèi)液體速度分布的實(shí)驗研究,化學(xué)工程,2011,39(10):68-72.

［6］王樹立,趙會軍. 鼓泡塔內(nèi)氣液兩相湍流實(shí)驗研究[J].化工科技,2004,12(2):6-11.

［7］翟甜,郝惠娣,高利博,等.鼓泡塔氣液兩相流內(nèi)部流場的流體力學(xué)特性[J].化工進(jìn)展,2013,32（10）:2319-2323.

［8］李 光, 楊曉鋼,蔡清白,等. CFD 優(yōu)化大型淺層鼓泡塔管式氣體分布器結(jié)構(gòu)[J]. 中國科技論文在線,2008,3（12）:809-896.

［9］趙陸海波,廖波,王小泉,等. 不同分布器對鼓泡塔氣液兩相流影響的CFD 模擬[J].現(xiàn)代化工,2012,32(11):101-104.

［10］王福軍.計算流體動力學(xué)分析—CFD軟件原理與應(yīng)用 [M]. 北京:清華大學(xué)出版社,2004: 210-212.

［11］朱紅鑰. Fluent 12 流體分析及工程仿真[M]. 北京:清華大學(xué)出版社,2011:34-37.

Experimental Simulation of Inner Flow in the Swirling Bubbling Reactor

CAO Qiao-qiao，HAO Hui-di，WANG Jin
（School of Chemical Engineering, Northwest University, Shaanxi Xi’an 710069, China）

The experimental simulation of inner flow in the swirling bubbling reactor was carried out by CFD, effect of the inlet speed ( 2,3,4,5 m/s) on inner flow in the swirling bubbling reactor was investigated, and the velocity vector, the velocity contour and the water volume fraction contour were analyzed. The results show that the fluid mixed effect with 4 m/s is better.

New type of cyclone; Bubbling reactor; Fluid; CFD simulation

TQ 018

： A

： 1671-0460（2015）04-0815-03

陜西省工業(yè)攻關(guān)項目（2011k10—21）；陜西省教育廳專項基金項目（11Jk0619）。

2014-11-14

曹喬喬（1990-），女 ，陜西咸陽人，西北大學(xué)碩士研究生，研究方向：主要從事攪拌反應(yīng)設(shè)備研究。E-mail: 364967149@qq.com。