曹喬喬 郝惠娣 王 瑾

(西北大學化工學院)

鼓泡塔在工業(yè)領(lǐng)域應(yīng)用廣泛，如加氫、費托合成、污水處理及煤液化等[1]。迄今為止，大量學者對鼓泡塔的流體特性進行了模擬研究[2～8]。近年來，隨著計算流體力學技術(shù)和軟件的迅速發(fā)展，鼓泡塔內(nèi)的流場等詳細信息可以通過 CFD 模擬的方法獲取[9,10]。筆者利用Ansys CFD軟件對鼓泡塔內(nèi)的氣、液兩相進行混合性能的數(shù)值模擬，考察了通氣速度為0.02、0.05、0.08、0.11m/s時鼓泡塔內(nèi)流場的特性，并對其計算結(jié)果進行了分析。

1 幾何模型及網(wǎng)格劃分

筆者研究采用帶有分布器的鼓泡塔(圖1)，塔外徑為300mm，內(nèi)徑為295mm，高為1 000mm，分布器位于塔的中部，有37個小圓孔，孔徑為20mm，分布器模型如圖2所示。

圖1 鼓泡塔結(jié)構(gòu)模型

圖2 分布器模型

網(wǎng)格劃分[11,12]采用有限體積法，為結(jié)構(gòu)化四面體網(wǎng)格，最小尺寸為0.2mm，最大面尺寸為10mm，最大體尺寸為20mm。網(wǎng)格總節(jié)點數(shù)為141 716個，網(wǎng)格單元數(shù)為776 446個，壁面Inflation網(wǎng)格細化，邊界層定義為5層。

2 初始條件和邊界條件

初始時，流體區(qū)域全是水，故設(shè)置水的體積分數(shù)為1，空氣的體積分數(shù)為0，液體湍動形式選擇Medium Intensity and Eddy Viscosity Ratio。液體的壁面設(shè)置成無滑移的固體壁面邊界，氣體設(shè)置為自由滑移固體壁面邊界，質(zhì)量和動量選擇Fluid Dependent。邊界的條件參數(shù)，氣體入口為速度入口，氣相體積分數(shù)為1。

3 模擬結(jié)果及討論

3.1鼓泡塔內(nèi)氣速的計算結(jié)果

由圖3可知：入口附近鼓泡塔的底端兩個角處的液體幾乎是不運動的，運動區(qū)域呈現(xiàn)U型，為了避免容器內(nèi)部出現(xiàn)死角，建議在下端設(shè)計選型時選用U型或V型。容器中心部分像V字一樣帶動液體，隨著氣體速度Ug的增加，中心部位的V型也就越明顯，距離中心處越遠速度越小，對壁面附近流體的運動而言，速度大不易帶動容器壁面附近流體的流動；對于容器中心而言，速度大時的帶動作用比速度小時的強；對于整個塔器而言，速度大時的帶動作用比較顯著；分布器對不同氣速的流體分布作用范圍相差不是很大，但是分布器對高速流體的阻力明顯比低速流體的大，特別是在中心部位。分布器對于高速流體的分布效果比較好，速度在容器內(nèi)部分布比較均勻，對于消除或減弱中心趨勢效果比較好。

圖3 不同通氣速度下氣體的速度分布云圖

3.2鼓泡塔內(nèi)徑向氣含率的計算結(jié)果

圖4為Ug=0.05m/s時不同高度下氣含率的分布圖。由圖4a、b可以看出，中心處氣含率最高，從中心到壁面氣含率逐漸降低；隨著高度h的增加，整體氣含率逐漸降低，但是氣體的分布范圍越來越廣；圖4c為h=0.53m(即分布器)處氣含率的分布情況，此時氣體分布量較大，中心處氣體分布較多，距離壁面處氣體分布較少；由圖4d可知h=0.60m(分布器上方)時，氣含率迅速降低，但整體分布均勻，分布器的分布作用顯著；對比圖4d、e可知，隨著高度h的增大，氣體又有向中心處匯聚的趨勢；由圖4f可以看出，氣體在壁面處的分布比較理想，由于出口比較小，當氣體與上端壁面相撞后，氣體改變方向向四周擴散。

圖4 不同高度氣含率的分布

3.3速度矢量圖的結(jié)果分析

Ug=0.05m/s時，氣、液表觀速度計算結(jié)果如圖5、6所示。由圖5、6可知，在分布器出口處，氣體和液體有返混現(xiàn)象，相比之下，氣體的返混現(xiàn)象比較弱，如圖5所示，在壁面附近處才出現(xiàn)氣體少量逆流現(xiàn)象，并且速度也特別?。蝗鐖D6所示，液體在分布器和出口處返混特別嚴重，從中心沿著直徑方向液體的表觀速度逐漸降低，呈開口向下的拋物線狀，直至減小到0，然后又反向逐漸增加。

圖5 分布器上部氣體表觀速度矢量圖

圖6 整體截面液體表觀速度矢量圖

4 結(jié)論

4.1當氣體通過分布器以后，氣體呈倒U型上升，經(jīng)過分布器后氣體的徑向擴散速度明顯增加，比分布器下端氣體的徑向擴散速度要大。設(shè)計選型時為了避免死角應(yīng)選用U型、V型，或者橢圓型和圓型。

4.2氣體在上升過程中氣含率逐漸向四周擴散，當氣體經(jīng)過分布器后整體分布相對均勻，故分布器的分布作用特別明顯。

4.3在分布器出口處，氣體和液體有返混現(xiàn)象。氣體的返混現(xiàn)象比較弱，在壁面附近處才出現(xiàn)氣體少量逆流現(xiàn)象，并且速度也特別?。灰后w在分布器和出口處返混特別嚴重，從中心沿著直徑方向液體的表觀速度逐漸降低。

[1] Sanyal J，Vasquez S，Roy S，et al.Numerical Simulation of Gas- Liquid Dynamics in Cylindrical Bubble Column Reactors[J].Chemical Engineering Science，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].化學工程，2008，36(2)：21～24.

[5] 樊亞超，靳海波，秦嶺，等.多級鼓泡塔內(nèi)液體速度分布的實驗研究[J].化學工程，2011，39(10)：68～72.

[6] 王樹立，趙會軍.鼓泡塔內(nèi)氣液兩相湍流實驗研究[J].化工科技，2004，12(2)：6～11.

[7] 翟甜，郝惠娣，高利博，等.鼓泡塔氣液兩相流內(nèi)部流場的流體力學特性[J].化工進展，2013，32(10)：2319～2323.

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

[9] Luo X，Liu P L，Luo H.Modeling of Liquid Circulation in Gas- Liquid Bubble Columns[J].Journal of Chemical Industry and Engineering，2006，57(7)：1565～1569.

[10] Wang T F，Wang J F，Jin Y. A CFD- PBM Coupled Model for Gas- Liquid Flows[J].AIChE Journal，2006，52 (1)： 125～140.

[11] 謝龍漢，趙新宇，張炯明.ANSYS CFX流體分析及仿真[M].北京：電子工業(yè)出版社，2012.

[12] 紀兵兵，陳金瓶.ANSYS ICEM CFD網(wǎng)格劃分技術(shù)實例詳解[M].北京：中國水利水電出版社，2012.