李少白,徐雙,范俊賡,黃思源

(沈陽航空航天大學 能源與環(huán)境學院,沈陽 110136)

能源與環(huán)境工程

非牛頓流體中在線雙氣泡相互作用的數(shù)值模擬

李少白,徐雙,范俊賡,黃思源

(沈陽航空航天大學 能源與環(huán)境學院,沈陽 110136)

利用Fluent軟件采用SIMPLE算法,數(shù)值模擬了非牛頓流體中在線雙氣泡之間的相互作用。先行氣泡和跟隨氣泡之間的相互作用與兩氣泡之間的中心間距有關,無因次距離S*越小,先行氣泡和跟隨氣泡聚并融合為一個氣泡的時間越短,當無因次距離S*增大到一定程度時,兩氣泡之間不會發(fā)生聚并,最終會以各自不同的速度上升。隨著氣泡之間的無因次距離S*增大,先行氣泡與跟隨氣泡速度之比不斷減小,直至先行氣泡和跟隨氣泡速度之比為1,跟隨氣泡不受先行氣泡的尾流影響。液相的流變性質不同,先行氣泡和跟隨氣泡速度之比為1時的無因次距離S*也不同,非牛頓性越強,先行氣泡尾流部分的剪切變稀作用越強,先行氣泡對跟隨氣泡的影響的距離越大,先行氣泡與跟隨氣泡速度之比減小的越快。

非牛頓流體;氣泡;相互作用;數(shù)值模擬

非牛頓流體中的氣泡運動廣泛存在于能源、環(huán)境、化工等工業(yè)過程中，如沼氣液、污泥或高分子溶液等都是典型的非牛頓流體[1],因此不管是理論基礎研究還是工程實際應用,研究非牛頓流體中的氣泡運動行為都具有較高的學術價值[2-3]。在這些過程中,氣泡之間的相互作用會影響其周圍的液相湍動狀況、氣液兩相的界面面積,進而影響到相間的傳熱、傳質速度[4-6]。為此,近年來許多學者針對非牛頓流體中的相互作用展開了研究。Kee D等人[7]利用高速攝像技術對非牛頓流體中的氣泡形狀及聚并行為進行了觀察,考察了流變性質、表面張力對氣泡性狀及聚并的影響。Sousa等人[8]利用PIV技術對Taylor氣泡在非牛頓流體中的相互作用時的流場信息進行了觀察,認為非牛頓流體中的氣泡相互作用距離大于牛頓流體。Fan等人[9]則利用PIV觀察了非牛頓流體中平行上升雙氣泡的流場特性。事實上,非牛頓流體中相互作用的氣泡之間有平行、在線及呈一定角度等多種方式,其中雙氣泡的在線相互作用是最簡單,也是最基礎的一種方式。因此,本文通過數(shù)值模擬方法對非牛頓流體中垂直在線的相互作用進行了研究,探討了氣泡間的距離和流體流變性質對氣泡相互作用的影響。

1 計算模型

1.1 控制方程

氣泡在非牛頓流體中的運動所涉及到的控制方程主要有連續(xù)性方程、動量方程,本過程為常溫常壓過程,液相可看做不可壓縮流體,其連續(xù)性方程為

(1)

其中,u表示速度矢量。

不可壓縮流體的動量方程如式(2)所示。

(2)

式中,ρ表示密度,F為質量力,μ表示流體的粘度,p表示流體內應力。

混合流體的密度和黏度計算方法如式(3)、式(4)所示。

ρ(F)=ρl(F)+ρg(1-F)

(3)

μ(F)=μl(F)+μg(1-F)

(4)

對于非牛頓流體而言,其粘度隨著剪切速率的變化而變化,可用冪律模型表示為

(5)

(6)

采用VOF方法追蹤氣液界面,采用該方法時,計算域內的流體的含量通過體積分數(shù)函數(shù)表示。其中F=0表示計算域內為純氣相,F=1表示純液相,0

(7)

追蹤界面通過求解體積分率函數(shù)的對流方程來實現(xiàn),如式(8)所示。

(8)

其中,對于氣液兩相的體積分率而言,應滿足歸一化,如式(9)所示。

Fl+Fg=1

(9)

最后,本文采用SIMPLE算法進行雙氣泡相互作用過程的求解。

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

幾何模型如圖1所示,選取的模型的邊界尺寸為50 mm×150 mm,由于氣泡直徑一般選擇2～10 mm,所以在該矩形區(qū)域內可以忽略壁面影響。上邊界為出口,邊界條件設為壓力出口,其他

圖1 二維計算模型示意圖

幾個壁面為光滑無滑動壁面。根據(jù)網(wǎng)格的獨立性,選取邊長為0.2 mm、0.3 mm、0.4 mm三種網(wǎng)格系統(tǒng)。三種網(wǎng)格尺寸下,單氣泡上升速度如圖2所示,從圖中可以看出,網(wǎng)格尺寸對單氣泡上升速度沒有什么影響,考慮到網(wǎng)格數(shù)越多,計算時間越長,對計算機的要求也越高;又考慮到計算結果的精確性,所以選擇網(wǎng)格邊長為0.3 mm的網(wǎng)格系統(tǒng)。

圖2 網(wǎng)格對氣泡上升速度的影響

2 結果與分析

非牛頓流體中的氣泡運動主要受到流體的流變性質、氣泡尺寸以及氣泡的排列方式等因素的影響。本文主要模擬了三種不同流變性質的流體,其物理性質如表1所示。

表1 模擬流體的物理及流變性質

未驗證非牛頓流體中氣泡運動模擬結果的可靠性,最具說服力的方法就是將模擬結果與實驗結果相對比。但是,對于具體實驗而言,氣泡在線相互作用的數(shù)據(jù)獲得十分困難,為此,本研究將對模擬過程中穩(wěn)定上升后的單氣泡運動的實驗結果與數(shù)值模擬結果進行對比,并說明數(shù)值模擬的可靠性。圖3為非牛頓流體單氣泡模擬結果與實驗結果的對比,從圖中可以看出,總體上模擬結果和實驗結果吻合良好。國內外許多研究表明[10-11],與牛頓流體相比,在適當?shù)臈l件下,非牛頓流體中氣泡會出現(xiàn)尖端的尾部,由圖3可以看出,本研究的模擬也出現(xiàn)了該現(xiàn)象。

氣泡的不同排列方式對氣泡之間相互作用的影響各不相同,當兩個氣泡在同一條直線上時,先行氣泡的運動會對跟隨氣泡上升速度和運動狀態(tài)有一定的影響。在研究氣泡之間的間距與氣泡上升速度之間的關系時,為了方便研究,氣泡之間距離無因次變量S*,其定義為

S*=L/D

(10)

其中,L為氣泡中心間距,D為氣泡直徑。

圖3 實驗所得的單氣泡的形狀與模擬結果的對比

在線氣泡上升過程如圖4所示,從中可以看出,兩個體積相同的氣泡在線上升運動過程,先行氣泡和跟隨氣泡一起上升,在上升的過程中兩氣泡逐漸變形,先行氣泡的形狀在水平方向上拉伸,跟隨氣泡則在豎直方向上拉伸。圖4給出氣泡體積相同,但氣泡之間初始間距不相同的三種情況下氣泡相互作用模擬結果,從圖中可以看出無因次距離S*越小,先行氣泡和跟隨氣泡聚并融合為一個氣泡的時間越短,兩氣泡聚并為一個氣泡后,新的氣泡體積是原來氣泡體積的兩倍,速度相對于在同等條件下單氣泡要略大。當無因次距離S*增大到一定程度時,兩氣泡之間不會發(fā)生聚并,最終以各自最終速度上升。先行氣泡在上升過程,其尾部會產(chǎn)生一個尾流區(qū)域,在該區(qū)域非牛頓流體剪切變稀,在先行氣泡尾流的作用下,跟隨氣泡上升過程中加速度比先行氣泡要大,經(jīng)過一段時間后,當跟隨氣泡最終追趕上先行氣泡,與先行氣泡發(fā)生相互作用,最終兩氣泡聚并為一個氣泡。模擬所得的非牛頓流體中的相互作用與文獻中報道的非牛頓流體中氣泡相互作用過程十分相符[12-13],表明了模擬結果的準確性和可靠性。

在線氣泡上升過程,上升速度隨時間變化如圖5所示,圖中給出三種不同濃度聚丙烯酰胺-PAA溶液單氣泡與雙氣泡速度之間的關系。從圖中可以看出,氣泡體積相同的情況下,在線雙氣泡的速度比單獨條件下的單氣泡速度大。雙氣泡在線運動過程中,跟隨氣泡的速度要大于先行氣泡,這是由于先行氣泡在運動過程中,其氣泡尾部有一個“尾流”區(qū)域[14]。在這個先行氣泡的“尾流”區(qū)域內,由于非牛頓流體的剪切變稀的性質,跟隨氣泡在進入“尾流”區(qū)域進行加速運動,所以跟隨氣泡的速度要大于先行氣泡的速度。此外,先行氣泡尾流部分的流場速度要大于周圍流體的流速。溶液濃度越低,“尾流”區(qū)域內流體的粘度越小,因此流體1中先行氣泡與跟隨氣泡速度要大于流體2和流體3的氣泡運動速度。

圖4 在線雙氣泡上升過程模擬結果圖

圖5 氣泡上升速度隨時間變化關系

先行氣泡和跟隨氣泡是否發(fā)生聚并對氣液兩相流之間的相互作用很重要。初始距離是影響氣泡是否發(fā)生聚并的一個很重要的因素,為了研究在線氣泡運動過程中兩氣泡的臨界間距對聚并的影響,在模擬過程中選取了距離的無因次變量S*=2、3、4、5、6、7、8、9、10。圖6給出了流體流變性質及兩氣泡速度比與無因次距離S*之間的關系,從圖中可以看出,隨著氣泡之間的間距S*增大,先行氣泡和跟隨氣泡的速度之比接近1,這是由于兩氣泡之間的間距越小,跟隨氣泡越容易受到先行氣泡的影響,當兩氣泡之間的間距增大到一定程度,跟隨氣泡不受先行氣泡尾流的影響。從圖6中還可以看出,液相的流變性質不同,先行氣泡和跟隨氣泡速度之比為1時氣泡的間距S*也不同,流體1中雙氣泡的臨界S*約為4.5左右,流體2中的臨界S*約為5左右,流體3中的臨界S*約為5.5左右。這是由于非牛頓性越強,先行氣泡尾流部分的剪切變稀作用越強,先行氣泡對跟隨氣泡的影響的距離越大。先行氣泡與跟隨氣泡速度之比減小的越快,這與Li[15]等人的實驗結果類似。

圖6 氣泡間距與兩氣泡速度之比的關系

3 結論

采用數(shù)值模擬手段對非牛頓流體中在線雙氣泡之間的相互作用進行了研究,主要得到如下主要結論。

(1)先行氣泡和跟隨氣泡之間的相互作用與兩氣泡之間的中心間距有關,無因次距離S*越小,先行氣泡和跟隨氣泡聚并融合為一個氣泡的時間越短,當無因次距離S*增大到一定程度時,兩氣泡之間不會發(fā)生聚并,最終會以各自不同的速度上升。

(2)隨著氣泡之間的無因次距離S*增大,先行氣泡與跟隨氣泡速度之比不斷減小,直至先行氣泡和跟隨氣泡速度之比為1時,跟隨氣泡不受先行氣泡的尾流影響。

(3)液相的流變性質不同,先行氣泡和跟隨氣泡速度之比為1時的無因次距離S*也不同,非牛頓性越強,先行氣泡尾流部分的剪切變稀作用越強,先行氣泡對跟隨氣泡的影響的距離越大。先行氣泡與跟隨氣泡速度之比減小得越快。

[1] CHHABRA R P.Bubbles,Drops and Particles in Non-Newtonian Fluids[M].2nded CRC Press:Boca Raton,FL,2007.

[2] FAN W,SUN Y,CHEN H.Bubble Volume and Aspect Ratio Generated in Non‐Newtonian Fluids[J].Chemical Engineering & Technology,2015,37(9):1566-1574.

[3] 李少白,徐雙,范俊賡,等.粘彈性流體中氣泡上升過程的軌跡[J].沈陽航空航天大學學報,2016,33(5):88-92.

[4] ABOULHASANZADEH B,TRYGGVASON G.Effect of bubble interactions on mass transfer in bubbly flow[J].International Journal of Heat & Mass Transfer,2014,79:390-396.

[5] GRACE S P S J R.Effect of Bubble Interaction on Interphase Mass Transfer in Gas Fluidized Beds[J].Chemical Engineering Science,1981,36(2):327-335.

[6] ABOULHASANXADEH B,TRYGGVASON G.Effect of bubble-bubble interaction on mass transfer in bubbly flow using a multi-scale approach[C].65th Annual Meeting of the APS Division of Fluid Dynamics.American Physical Society,San Diego,2012.

[7] KEE D D,CHHABRA R P.A photographic study of shapes of bubbles and coalescence in non-Newtonian polymer solutions[J].Rheologica Acta,1988,27(6):656-660.

[8] SOUSA R G,PINTO A M F R,CAMPOS J B L M.Interaction between Taylor bubbles rising in stagnant non-Newtonian fluids[J].International Journal of Multiphase Flow,2007,33(9):970-986.

[9] FAN W,YOUGUANG M A,XIAOLEI L I,et al.Study on the Flow Field around Two Parallel Moving Bubbles and Interaction Between Bubbles Rising in CMC Solutions by PIV[J].Chinese Journal of Chemical Engineering,2009,17(6):904-913.

[10]LIU Y J,LIAO T Y,JOSEPH D D.A two-dimensional cusp at the trailing edge of an air bubble rising in a viscoelastic liquid[J].Journal of Fluid Mechanics,1995,304(304):321-342.

[11]BELMONTE A.Self-oscillations of a cusped bubble rising through a micellar solution[J].Rheologica Acta,2000,39(6):554-559.

[12]LIN T J,LIN G M.The Mechanisms of bubble coalescence in a non-Newtonian fluid[J].Canadian Journal of Chemical Engineering,2003,81(3-4):476-482.

[13]LIN T J,LIN G M.Mechanisms of in-line coalescence of two-unequal bubbles in a non-Newtonian fluid[J].Chemical Engineering Journal,2009,155(3):750-756.

[14]SOUSA R G,RIETHMULLER M L,PINTO A M F R,et al.Flow around individual Taylor bubbles rising in stagnant polyacrylamide (PAA) solutions[J].Journal of Non-Newtonian Fluid Mechanics,2006,135(1):16-31.

[15]LI H Z,MOULINE Y,FUNFSCHILLING D,et al.Evidence for in-line bubble interactions in non-Newtonian fluids[J].Chemical Engineering Science,1998,53(12):2219-2230.

(責任編輯：劉劃 英文審校：趙歡)

Numericalsimulationofinteractionbetweenin-linetwobubblesinnon-Newtonianfluids

LI Shao-bai，XU Shuang，F(xiàn)AN Jun-geng，HUANG Si-yuan

(College of Energy and Environment,Shenyang Aerospace University,Shenyang 110136,China)

Interaction between in-line two bubbles in non-Newtonian fluids was studied numerically by Fluent software and SIMPLE algorithm.The results showed that the interaction between leading bubble and following bubble was related to center distance between the two bubbles.The smaller the dimensionless distanceS*,the shorter the time that the two bubbles coalesce each other.However,the coalescence will not occur when theS*increased and exceeded a certain value.Finally the two bubbles rose separately with different speed.Moreover,the velocity ratio between the two bubbles decreased with the increase ofS*,and then following bubble was independent on the leading bubble as the ratio approached to 1.The rheological property of liquid phase impacted theS*value when the ratio was equal to 1.The bigger the rheological property of liquid phase,the stronger the shear action of stern for the leading bubble and the effects of the leading bubbles on the following bubbles,which in turn resulted in the smaller velocity ratio.

non-Newtonian fluids;bubble;interaction;numerical simulation

2017-06-15

國家自然科學基金(項目編號：21406141)

李少白(1982-)，男，遼寧朝陽人，博士，主要研究方向:復雜的氣液兩相體系的傳遞特性，E-mail: lishaobai1982@163.com。

2095-1248(2017)04-0063-06

TQ021.1

： A

10.3969/j.issn.2095-1248.2017.04.008