周曉琳，田 震，王麗軍，成有為，李 希

浙江大學(xué)化學(xué)工程與生物工程學(xué)院，浙江 杭州 310027

鼓泡塔中氣泡群運(yùn)動(dòng)的實(shí)驗(yàn)與模擬

周曉琳，田 震，王麗軍，成有為，李 希

浙江大學(xué)化學(xué)工程與生物工程學(xué)院，浙江 杭州 310027

采用高速攝像法測(cè)量了均勻鼓泡流狀態(tài)下，水以及體積分?jǐn)?shù)分別為20%和40%的甘油-水體系中氣泡群的浮升運(yùn)動(dòng)，考察了氣含率、雷諾數(shù)和分布器孔徑對(duì)氣泡尺寸、形狀、浮升速率和曳力系數(shù)的影響。構(gòu)建了立方體單元胞模型，并根據(jù)雷諾數(shù)的不同選取層流和湍流模型，模擬得到氣泡的浮升速率和曳力系數(shù)與實(shí)驗(yàn)值吻合較好。結(jié)果表明：隨著氣含率及液體黏度的提高，氣泡群浮升阻力增大，浮升速率減少；隨著雷諾數(shù)和氣泡直徑的增加，曳力系數(shù)減小，氣泡浮升速率增大。單元胞模型能較好地反映氣泡群浮升過(guò)程中各因素的影響，是處理氣泡群運(yùn)動(dòng)的有效工具。

氣泡群 高速攝像 單元胞模型 曳力系數(shù)

鼓泡塔具有結(jié)構(gòu)簡(jiǎn)單、操作方便、無(wú)機(jī)械傳動(dòng)部件及相間接觸面積大等優(yōu)點(diǎn)，被廣泛地應(yīng)用于化工生產(chǎn)過(guò)程。鼓泡塔中氣泡群的浮升及其與液體間的相互作用決定了鼓泡塔的流動(dòng)特性，有關(guān)鼓泡塔內(nèi)氣泡群的研究較多。Martinez-Mercado等[1,2]使用雙電導(dǎo)探針?lè)ㄑ芯苛藲馀萑焊∩^(guò)程中氣泡直徑和相對(duì)滑移速率等參數(shù)。Rabha等[3]采用高速攝像技術(shù)研究了氣泡群在均勻分布和非均勻分布條件下的氣泡直徑和浮升速率隨氣含率的變化，高速攝像所獲得的圖像更直觀，可靠性較高，更適用于氣速較低的均勻鼓泡的情況。在氣泡群運(yùn)動(dòng)模擬方面，目前采用的模型主要為直接數(shù)值模擬（DNS）和單元胞模型。DNS能更好地描述氣泡群的真實(shí)運(yùn)動(dòng)，但模型的計(jì)算量較大。Bunner等[4]用DNS模擬了216個(gè)氣泡，計(jì)算耗時(shí)數(shù)十天。單元胞模型由Happel[5]提出，該模型的核心是將液體均分，通過(guò)計(jì)算典型單個(gè)氣泡，獲得氣泡群浮升運(yùn)動(dòng)規(guī)律。相對(duì)DNS模型，單元胞模型的計(jì)算量較小。Leclair等[6]用單元胞模型模擬了雷諾數(shù)（Re）小于1 000時(shí)的球形氣泡群，模擬與實(shí)驗(yàn)值吻合較好。Manjunath等[7]用單元胞模型研究了Re小于100時(shí)氣泡群曳力系數(shù)與氣含率和雷諾數(shù)的關(guān)系。毛在砂[8]模擬了Re為50時(shí)氣泡的運(yùn)動(dòng)，并考慮了氣泡在運(yùn)動(dòng)過(guò)程中的變形。本工作使用高速攝像方法拍攝3種不同黏度狀態(tài)下液體中氣泡群浮升過(guò)程，測(cè)量了氣泡大小、形態(tài)變化以及浮升速率等，并利用單元胞模型進(jìn)行氣泡群浮升速率和曳力計(jì)算，通過(guò)選取合適的單元胞結(jié)構(gòu)和流體方程，將單元胞模型應(yīng)用于更寬的雷諾數(shù)范圍。從實(shí)驗(yàn)和理論兩方面對(duì)鼓泡塔中氣泡群浮升運(yùn)動(dòng)規(guī)律進(jìn)行研究。

1 實(shí)驗(yàn)部分

1.1 實(shí)驗(yàn)裝置

實(shí)驗(yàn)裝置如圖1所示，主要包括流體控制和檢測(cè)兩個(gè)部分。流體控制系統(tǒng)由空氣壓縮機(jī)、儲(chǔ)蓄罐、調(diào)壓器、測(cè)試塔和閥門(mén)等組成。鼓泡塔為長(zhǎng)方體透明有機(jī)玻璃水槽，尺寸為170 cm×15 cm×5 cm，厚度方向比較小，近似可作為二維鼓泡塔，塔底部采用分布器進(jìn)氣，進(jìn)氣口均勻分布于塔底，為考察進(jìn)氣孔徑的影響，選擇了0.2和0.6 mm兩種孔徑分布器。測(cè)量系統(tǒng)由壓差計(jì)和圖像采集裝置組成。壓差計(jì)用于測(cè)量塔內(nèi)平均氣含率。圖像裝置設(shè)備包括照明、高速攝像機(jī)和計(jì)算機(jī)等單元。采用德國(guó)Mikrotron公司的Cube7高速攝像機(jī)測(cè)量二維鼓泡塔中的氣泡群浮升運(yùn)動(dòng)，每秒1 000幅的拍攝速度，采用1 000 W的冷光源逆向照射，拍攝垂直方向氣泡的運(yùn)動(dòng)圖像，采集圖像分辨率為1 024×1 024，幀頻為500 fps。

圖1 實(shí)驗(yàn)裝置及典型氣泡示意Fig.1 Experimental set-up and typical bubble figure1,7-computer; 2-A/D conversion card; 3-pressure transmitter; 4-light; 5-bubble column; 6-camera; 8-gas distributor; 9-gas chamber; 10-flow meter; 11-gas buffer tank; 12-compressor

實(shí)驗(yàn)所用氣體為空氣，為考察雷諾數(shù)變化對(duì)氣泡浮升運(yùn)動(dòng)的影響，選擇了水及兩種配比的水-甘油混合物（室溫10.3 ℃）為液相流體，相關(guān)組成及物性列于表1。在此3種氣液體系下，當(dāng)氣體的表觀氣速小于0.05 m/s，鼓泡流型為均勻鼓泡流，氣泡群浮升過(guò)程中不存在氣泡聚并和破碎現(xiàn)象。

表1 液相組成及物理性質(zhì)Table 1 Liquid phase composition and physical properties

1.2 實(shí)驗(yàn)步驟

首先進(jìn)行尺寸標(biāo)記，向鼓泡塔中注入液體，并放入6 mm管，調(diào)節(jié)高速攝像機(jī)焦距以獲得最清晰的圖像，作為測(cè)量氣泡尺寸的基準(zhǔn)；然后打開(kāi)空氣壓縮機(jī)，調(diào)節(jié)流量計(jì)，使空氣通過(guò)分布器以鼓泡形式均勻進(jìn)入鼓泡塔中；待鼓泡塔運(yùn)行穩(wěn)定后，進(jìn)行圖像實(shí)時(shí)采集，不斷調(diào)節(jié)流量計(jì)改變表觀氣速，獲得不同表觀氣速實(shí)驗(yàn)數(shù)據(jù)。為考察分布器孔徑影響，更換不同孔徑分布器孔，再次重復(fù)以上實(shí)驗(yàn)步驟。

1.3 圖像處理

圖像分析采用Matlab編程實(shí)現(xiàn)，將采集的圖像轉(zhuǎn)化為灰度圖，如圖2所示。氣泡所對(duì)應(yīng)的灰度較高，因而可利用圖像中灰度的梯度變化識(shí)別每一個(gè)氣泡的邊緣信息，并通過(guò)全局閾值設(shè)置調(diào)節(jié)圖像的分辨率。為計(jì)算氣泡的尺寸，使用標(biāo)準(zhǔn)6 mm的實(shí)心管在測(cè)量位置上進(jìn)行標(biāo)定，通過(guò)拍攝到的氣泡大小與標(biāo)定管的尺寸比較，確定變形橢球氣泡長(zhǎng)軸dmax，短軸dmin和縱橫比x（x=dmax/dmin），得到氣泡的等效直徑。氣泡浮升速率采用顆粒跟蹤算法（PTV）計(jì)算，根據(jù)所獲得的動(dòng)態(tài)圖像，在較短的時(shí)間間隔內(nèi)，跟蹤每個(gè)氣泡的運(yùn)動(dòng)軌跡，計(jì)算出單個(gè)氣泡浮升速率的瞬時(shí)速率，進(jìn)而得到平均浮升速率（Vg），根據(jù)式（1）計(jì)算氣泡的曳力系數(shù)（CD）：

式中：Vrel=Vg-Vl，在本實(shí)驗(yàn)所考慮的均勻鼓泡流范圍內(nèi)，Vl≈0，Vrel≈Vg。

2 單元胞數(shù)學(xué)模型

單元胞模型是處理鼓泡塔內(nèi)氣泡群浮升運(yùn)動(dòng)的重要方法，其假設(shè)氣泡群中的氣泡均勻地分布在液體中，并且每個(gè)氣泡表面攜帶的液體量相等，單元胞由位于中心的氣泡和它周?chē)囊后w所構(gòu)成，單元胞中氣液間的體積比例與鼓泡塔的平均氣含率一致。

2.1 單元胞模型

Happel等[5-7]采用球形單元胞，球形的對(duì)稱(chēng)性好且理論推導(dǎo)方便。Sirignano等[9]認(rèn)為更合理的單元胞形狀是立方體形，不足之處是立方體單元胞屬于三維問(wèn)題，求解計(jì)算量較大，但隨著計(jì)算機(jī)能力的提高，求解3維計(jì)算流體力學(xué)（CFD）問(wèn)題目前已不困難，因而本工作采用立方體單元胞模型模擬氣泡群運(yùn)動(dòng)。

圖2 單元胞模型示意Fig.2 Diagram of cell model

圖3 單氣泡繞流示意Fig.3 Diagram of flow around the single bubble

立方體單元胞如圖2所示，假設(shè)氣泡以相同的間隔均勻地分布在液體里，每個(gè)氣泡分屬于各自的立方體內(nèi)，研究一個(gè)立方體中的流場(chǎng)即可掌握整個(gè)氣泡群的運(yùn)動(dòng)規(guī)律。以氣泡作為運(yùn)動(dòng)參考系，圖3給出了立方體單元胞中的繞流場(chǎng)的速度分布。

2.2 流體力學(xué)方程

有關(guān)單元胞流體力學(xué)模型均采用層流方程，適用雷諾數(shù)較低的情況[7,8,10]，本實(shí)驗(yàn)在水體系下，氣泡的浮升雷諾數(shù)接近1 000，應(yīng)該使用湍流模型進(jìn)行模擬，而對(duì)于水-甘油混合體系，液體黏度較高，雷諾數(shù)為50～400，可以采用Navier-Stokes方程計(jì)算，相關(guān)流體力學(xué)方程列于表2。

表2 流體力學(xué)方程Table 2 Fluid mechanics equations

單元胞中心氣泡為橢球形，橢球長(zhǎng)軸和短軸大小通過(guò)實(shí)驗(yàn)拍照確定。橢球所占立方體單元胞體積分?jǐn)?shù)等于全塔氣含率，故立方體單元胞邊長(zhǎng)（a）為：

邊界條件為：（1）單元胞之間的分界壁為對(duì)稱(chēng)邊界條件，二維鼓泡塔中單元胞與容器壁面邊界條件采用無(wú)滑移邊界條件，湍流模型使用標(biāo)準(zhǔn)壁函數(shù)；（2）氣泡表面采用部分滑移邊界條件，參照文獻(xiàn)[10]，設(shè)置污染角為90°；（3）入口和出口定義為充分發(fā)展邊界，以保證進(jìn)出口流體分布規(guī)律相同，給定入口及出口邊界的質(zhì)量流率。

2.3 求解方法

模型采用Star-CCM+10.0軟件求解，采用三維穩(wěn)態(tài)模擬，使用分離求解器，收斂標(biāo)準(zhǔn)為相對(duì)殘差小于10-3。利用幾何結(jié)構(gòu)的對(duì)稱(chēng)性，僅選取一半的物理區(qū)域進(jìn)行計(jì)算，以減少計(jì)算負(fù)荷。使用Star CCM+特色的網(wǎng)格劃分工具劃分多面體網(wǎng)格，所需內(nèi)存較少，求解速率較快。網(wǎng)格密度設(shè)置為：當(dāng)αg小于2%時(shí)，網(wǎng)格基本尺寸為3×10-4m；當(dāng)αg大于2%時(shí)，網(wǎng)格基本尺寸為1.5×10-4m，該網(wǎng)格密度通過(guò)網(wǎng)格無(wú)關(guān)性檢驗(yàn)。

3 結(jié)果與討論

3.1 平均氣含率

圖4 不同表觀氣速下全塔氣含率Fig.4 Gas holdup of whole tower as a function of superficial gas velocity

全塔平均氣含率隨表觀氣速的變化如圖4所示。由圖可知，隨著表觀氣速的增大，平均氣含率單調(diào)遞增，且3種液體中氣含率隨表觀氣速變化規(guī)律相近。在本實(shí)驗(yàn)范圍內(nèi)，氣含率與表觀氣速滿足線性關(guān)系。

3.2 氣泡尺寸和形狀

在相同表觀氣速下3種不同體系中拍攝得到的典型氣泡群浮升圖像如圖5所示。圖5（a）～（c）采用的分布器孔口直徑為0.2 mm，圖5（d）～（f）使用的分布器孔口直徑為0.6 mm。由圖可知，氣泡為橢球形，水體系中氣泡等效直徑大于水-甘油混合體系，并且隨著氣含率的增大氣泡變形程度增大。在0.6 mm分布器條件下，水體系中氣泡變形高于水-甘油混合體系。

圖5 典型氣泡群運(yùn)動(dòng)Fig.5 Motion diagram of typical bubble swarma, d-water; b, e-20% glycerol-water; c, f-40% glycerol-water

圖6為在不同氣含率下3個(gè)體系中的氣泡等效直徑。由圖可知，隨著氣含率的增大，氣泡等效直徑增大，純水中氣泡直徑明顯大于水-甘油混合體系。這是由于降低液體表面張力可以減小氣泡尺寸，而水的表面張力最大。對(duì)比0.2和0.6 mm孔徑分布器，大孔徑分布器所得到的氣泡直徑也較大，說(shuō)明在均勻鼓泡狀態(tài)下，氣泡尺寸受分布器影響較大。

圖6 不同氣含率下氣泡等效直徑Fig.6 Equivalent diameter of bubble as a function of gas holdup

3.3 氣泡浮升速率

圖7為實(shí)驗(yàn)和模擬得到的氣泡浮升速率。由圖可以看出，氣泡浮升速率均隨著氣含率增大而減小，這是由于氣泡群隨著氣含率的增大，氣泡間距減少，氣泡浮升過(guò)程中氣泡間隙處的氣液相對(duì)滑移速率增大，使得氣泡所受阻力增大，導(dǎo)致氣泡浮升速率減小。由圖還可以看到，隨著甘油濃度的增加，氣泡浮升速率減少，這是由于甘油含量的增大使得液體黏度增大，摩擦阻力增大。對(duì)比不同分布器孔徑結(jié)果，分布器孔徑增加，氣泡浮升速率增大，這是由于大孔分布器產(chǎn)生的氣泡等效直徑大，所受浮力更大。

圖7 不同氣含率下氣泡浮升速率Fig.7 Bubble velocity as a function of gas holdup

對(duì)比模擬與實(shí)驗(yàn)結(jié)果可知，在0.2 mm分布器條件下，模擬和實(shí)驗(yàn)吻合結(jié)果更好，0.6 mm分布器條件下偏差增大，這可能是因?yàn)榇罂讖椒植计鳟a(chǎn)生的氣泡在浮升過(guò)程中動(dòng)態(tài)變形較大，單元胞模型很難描述這種現(xiàn)象。此外，隨著體系黏度增加，模擬與實(shí)驗(yàn)值間誤差也增大，實(shí)驗(yàn)發(fā)現(xiàn)黏度提高，氣泡團(tuán)聚有所增加，導(dǎo)致實(shí)驗(yàn)值與模擬值的偏差。

3.4 氣泡曳力系數(shù)

表3為實(shí)驗(yàn)和模擬得到的曳力系數(shù)相關(guān)數(shù)據(jù)匯總，氣泡的曳力系數(shù)實(shí)驗(yàn)值和模擬值的相對(duì)誤差在±10%左右，二者符合較好，說(shuō)明單元胞模型能較好地描述均勻鼓泡流中氣泡群所受到的阻力。由表可知，隨著氣含率的增大氣泡的變形程度增大，浮升速率減小，曳力系數(shù)增加。

使用立方體單元胞模型計(jì)算得到雷諾數(shù)對(duì)氣泡曳力系數(shù)的影響結(jié)果見(jiàn)圖8。由圖可知，隨著雷諾數(shù)的增加，氣泡曳力系數(shù)減小，在低氣含率層流條件下，模型計(jì)算結(jié)果與Tomiyama等[11]給出的單氣泡曳力公式符合較好。在高雷諾數(shù)情況下，湍流模型與層流模型計(jì)算結(jié)果差異較大，在水體系中湍流模型與實(shí)驗(yàn)結(jié)果吻合較好，說(shuō)明高雷諾數(shù)下，采用湍流模型更合理。

表3 曳力系數(shù)相關(guān)數(shù)據(jù)匯總Table 3 Summary of drag coefficient data

圖8 不同雷諾數(shù)下曳力系數(shù)值Fig.8 Drag coefficient as a function of Reynolds number

Manjunath等[7,8]采用單元胞模型研究了Re為50時(shí)，氣泡曳力系數(shù)與氣含率的關(guān)系，相同條件下采用本實(shí)驗(yàn)?zāi)Ｐ瓦M(jìn)行了模擬計(jì)算，結(jié)果見(jiàn)圖9（a）。由圖可知，隨著氣含率的增大，氣泡曳力系數(shù)增大，文獻(xiàn)[7]的結(jié)果明顯偏小。這主要是由于其沒(méi)有考慮氣泡變形，本實(shí)驗(yàn)?zāi)Ｐ陀?jì)算結(jié)果與文獻(xiàn)[8]的結(jié)果略有偏差，可能是由于本實(shí)驗(yàn)采用立方體單元胞，而文獻(xiàn)[8]采用球形單元胞。圖9（b）進(jìn)一步給出了在較高雷諾數(shù)情況下，預(yù)測(cè)氣泡曳力系數(shù)與氣含率的關(guān)系，隨著氣含率的增大，氣泡曳力系數(shù)隨之增大。

圖9 不同氣含率下曳力系數(shù)值Fig.9 Drag coefficient as a function of gas holdup

4 結(jié) 論

采用高速攝像法測(cè)量了水、20%甘油-水和40%甘油-水3種體系中均勻鼓泡流狀態(tài)下的氣泡群浮升運(yùn)動(dòng)，構(gòu)建了立方體結(jié)構(gòu)的單元胞數(shù)學(xué)模型，實(shí)驗(yàn)及模擬結(jié)果表明：

a）隨著氣含率的增大，氣泡等效直徑增加，氣泡變形程度增大，氣泡間隙處的流體運(yùn)動(dòng)受到阻礙，氣泡的曳力系數(shù)增大，浮升速率減小。

b）隨著液體黏度增大，氣泡等效直徑減小，并且0變形程度較小，黏度的增大使氣泡所受摩擦阻力增大，氣泡浮升速率減小，曳力系數(shù)增大。

c）氣泡尺寸與分布器孔徑有關(guān)，分布器孔徑對(duì)氣泡浮升速率有較大影響。

d）隨著雷諾數(shù)的增加，氣泡的曳力系數(shù)減小，水體系中采用湍流模型更合理。

e）計(jì)算得到的氣泡浮升速率及曳力系數(shù)與實(shí)驗(yàn)值符合較好，本工作所提出的單元胞模型能較好地預(yù)估均勻鼓泡流中氣泡群浮升的運(yùn)動(dòng)行為。

符號(hào)說(shuō)明

參考文獻(xiàn)：

[1] Martínez-Mercado J, Palacios-Morales C A, Zenit R. Measurement of pseudoturbulence intensity in monodispersed bubbly liquids for 10＜Re＜ 500[J]. Physics of Fluids, 2007, 19(10): 103302.

[2] Simonnet M, Gentric C, Olmos E, et al. Experimental determination of the drag coefficient in a swarm of bubbles[J]. Chemical Engineering Science, 2007, 62(3): 858-866.

[3] Rabha S S, Buwa V V. Experimental investigations of rise behavior of monodispersed/polydispersed bubbly flows in quiescent liquids[J]. Industrial & Engineering Chemistry Research, 2010, 49(21): 10615-10626.

[4] Bunner B, Tryggvason G. Dynamics of homogeneous bubbly flows part 1: rise velocity and microstructure of the bubbles[J]. Journal of Fluid Mechanics, 2002, 466(1): 17-52.

[5] Happel J. Viscous flow in multiparticle systems: slow motion of fluids relative to beds of spherical particles[J]. AIChE Journal, 1958, 4(2): 197-201.

[6] Leclair B P, Hamielec A E. Viscous flow through particle assemblages at intermediate Reynolds numbers: a cell model for transport in bubble swarms[J]. The Canadian Journal of Chemical Engineering, 1971, 49(6): 713-720.

[7] Manjunath M, Tripathi A, Chhabra R P, et al. Numerical simulation of the drag on a swarm of bubbles[J]. International Journal of Engineering Science, 1994, 32(6): 927-933.

[8] 毛在砂. 用單元胞模型數(shù)值模擬氣泡群中氣泡的運(yùn)動(dòng)[J]. 化工學(xué)報(bào), 2007, 58(5): 1155-1162. Mao Zaisha. Cell model approach to motion of bubbles in swarm[J]. Journal of Chemical Industry and Engineering, 2007, 58(5): 1155-1162.

[9] Sirignano W A, Tal R. Cylindrical cell model for the hydrodynamics of particle assemblages at intermediate Reynolds numbers[J]. AIChE Journal, 1982, 28(2): 233-237.

[10] Nalajala V S, Kishore N, Chhabra R P. Effect of contamination on rise velocity of bubble swarms at moderate Reynolds numbers[J]. Chemical Engineering Research and Design, 2014, 92(6): 1016-1026.

[11] Tomiyama A, Kataoka I, Zun I, et al. Drag coefficients of single bubbles under normal and micro gravity conditions[J]. JSME International Journal Series B: Fluids and Thermal Engineering, 1998, 41(2): 472-479.

Experiment and Numerical Simulation of the Motion of Bubble Swarm in a Bubble Column

Zhou Xiaolin, Tian Zhen, Wang Lijun, Cheng Youwei, Li Xi
Department of Chemical and Biological Engineering, Zhejiang University, Hangzhou 310027, China

High-speed digital video processing was used to measure the motion of bubble swarm in water, 20% glycerol-water and 40% glycerol-water mixture under the uniform bubble flow state, respectively. The effects of gas holdup, Reynolds number and distributor aperture on the bubble size and shape, bubble velocity and drag coefficient were investigated. The cubic cell model was established and the laminar flow model and turbulent flow model were selected according to the different Reynolds number to simulate the motion of the bubble swarm. The results showed that the numerical rise velocity and the drag coefficient of the bubbles were in good agreement with the experimental data. The rise resistance of the bubble swarm increased and the rise velocity decreased with the increase of gas holdup and liquid viscosity. With the increase of the Reynolds number and the bubble equivalent diameter, the bubble swarm drag coefficient decreased and the rise velocity increased. It is suggested that the cell model can well reflect the influences of various factors in the process of the bubble swarm motion and it is an effective tool to deal with the movement of the bubble swarm.

bubble swarm; high-speed video; cell model; drag coefficient

TQ018；TQ021.1

：A

1001—7631 ( 2016 ) 06—0481—09

2016-03-21；

：2016-09-27。

周曉琳（1991—），女，碩士研究生；王麗軍（1976—），男，博士，副研究員，通訊聯(lián)系人。E-mail: wang_lijun@zju.edu.cn。

國(guó)家自然科學(xué)基金重大研究計(jì)劃（91334105）。