肖翾，馬楊，沈陽(yáng)，程永攀，徐進(jìn)良

(華北電力大學(xué)能源動(dòng)力與機(jī)械工程學(xué)院 低品位能源多相流與傳熱北京市重點(diǎn)實(shí)驗(yàn)室，北京 102206)

氣泡在化工、生物、醫(yī)藥、動(dòng)力設(shè)備、核反應(yīng)堆、航天和熱能等多個(gè)領(lǐng)域有著廣泛的應(yīng)用，譬如：大型船舶的氣泡減阻、油氣輸運(yùn)、血管中氣泡的運(yùn)動(dòng)。船舶發(fā)動(dòng)機(jī)水下排氣會(huì)造成螺旋槳的氣蝕,利用氣泡可以加速反應(yīng)裝置中的物質(zhì)混合、熱量交換及化學(xué)反應(yīng)過(guò)程；在水處理方面，曝氣、氣浮等工藝的處理能耗和處理效率與氣泡運(yùn)動(dòng)特性密切相關(guān)。在沸騰傳熱中，氣泡脫離是沸騰傳熱的重要過(guò)程，氣泡從基板上脫離受多個(gè)因素的影響，包括氣泡的形狀、大小、加熱模式、基板表面的潤(rùn)濕性等。

Fritz和Ende[1]通過(guò)氣泡脫離時(shí)的受力分析，得到著名的Fritz公式，在常壓情況下，該公式的計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果較為吻合，但當(dāng)壓力過(guò)高或過(guò)低時(shí)有誤差。Cole[2]在低壓條件下進(jìn)行一系列的氣泡生長(zhǎng)和脫離的實(shí)驗(yàn)，并得到新的擬合關(guān)系式，該式中引入一個(gè)Jakob數(shù)，后來(lái)得到廣泛引用。Rohsenow[3]發(fā)現(xiàn)氣泡生長(zhǎng)需要的一定過(guò)熱度，其大小由氣液界面的表面張力決定。Lee和Nydahl[4]對(duì)氣泡的生長(zhǎng)和脫離進(jìn)行數(shù)值模擬，假設(shè)氣泡具有球缺形狀并且在整個(gè)氣泡生長(zhǎng)期間存在微層,他們發(fā)現(xiàn)，微層蒸發(fā)的熱量占整個(gè)核態(tài)沸騰換熱量的87%。Zeng等[5]提出，氣泡在上升與分離時(shí)主要由氣泡受到的浮力以及表面張力決定，據(jù)此獲得氣泡生長(zhǎng)和氣泡分離直徑的表達(dá)式。Son等[6]使用Level Set方法對(duì)水沸騰期間氣泡的生長(zhǎng)和離開進(jìn)行數(shù)值模擬。Fuchs等[7]利用數(shù)值模擬的方法，分析壁面導(dǎo)熱效應(yīng)和蒸發(fā)彎月面對(duì)氣泡生長(zhǎng)和脫離的影響。本文利用數(shù)值模擬的方法，對(duì)于氣泡脫離基板的機(jī)理進(jìn)行深入全面的研究。

1 模型和公式

本文利用Level Set方法建立二維軸對(duì)稱非穩(wěn)態(tài)數(shù)學(xué)模型，用有限差分法對(duì)控制方程進(jìn)行離散處理，在給定的邊界條件下，追蹤氣液界面的演化規(guī)律，獲得氣泡在浮力、黏性力、慣性力與表面張力共同作用下的運(yùn)動(dòng)規(guī)律，分析不同接觸角、Re數(shù)和Bo數(shù)對(duì)氣泡脫離基板的影響，并獲得氣泡從基板脫離時(shí)的臨界參數(shù)。

本文中，以單個(gè)氣泡為研究對(duì)象，氣泡初始形狀為球缺形。由于氣泡沿y軸方向是軸對(duì)稱的，我們僅對(duì)氣泡x軸正方向的部分進(jìn)行模擬，計(jì)算區(qū)域底部為無(wú)滑移邊界條件，其余邊界為自由邊界條件，如圖1所示。數(shù)值模擬結(jié)果表明，當(dāng)氣液的密度比和黏性比小于0.05時(shí)，該比值對(duì)于氣泡的運(yùn)動(dòng)規(guī)律幾乎沒有影響，因此選取氣液的密度比和黏性比為0.05。氣泡無(wú)初始速度，在脫離基板的過(guò)程中滿足質(zhì)量守恒。經(jīng)過(guò)校驗(yàn)后，網(wǎng)格數(shù)101×301，時(shí)間步長(zhǎng)為1.0×10-4，計(jì)算域?yàn)?×9。

圖1 氣泡脫離基板示意圖Fig.1 Schematic of bubble detachment from substrate

由于氣泡是在浮升力的驅(qū)動(dòng)下脫離的，因此以浮升力來(lái)定義特征速度

(1)

引入無(wú)量綱參數(shù)：

Re=ρlUR/μl,

(2)

Bo=ρlgR2/σ.

(3)

分析得到氣泡脫離壁面的無(wú)量綱控制方程如下：

連續(xù)性方程：

(4)

動(dòng)量方程：

(5)

式中:ρ*=ρg/ρl，μ*=μg/μl分別為氣泡內(nèi)外的密度比和黏度比。

2 結(jié)果與討論

氣泡脫離基板的過(guò)程受眾多因素的影響，將控制方程進(jìn)行無(wú)量綱化以后，這些因素可以歸結(jié)為Re數(shù)、Bo數(shù)和接觸角的影響。

2.1 模型校核

為了驗(yàn)證模型的正確性，我們與Sussman 和Smereka[8]的結(jié)果進(jìn)行比較，如圖2所示。結(jié)果表明，我們的模擬結(jié)果和他們的結(jié)果非常符合。氣泡在浮力作用下，由靜止開始上升，圓形氣泡下部開始變形，進(jìn)而導(dǎo)致氣泡破裂。

圖2 Re=100和Bo=2時(shí)氣泡運(yùn)動(dòng)過(guò)程的形狀演化Fig.2 Bubble shape evolution at Re=100 and Bo=2

2.2 Re數(shù)對(duì)氣泡脫離基板的影響

首先通過(guò)改變氣泡脫離基板時(shí)的Re數(shù)(其中Bo=1，θ=60°)，研究Re數(shù)的變化對(duì)氣泡脫離基板的影響。由圖3可知，不同Re數(shù)下氣泡脫離基板的垂直速度明顯不同，當(dāng)Re=1時(shí)，由于黏性的影響，氣泡無(wú)法脫離，氣泡的垂直速度很小，可忽略不計(jì)；當(dāng)Re=10時(shí)，氣泡脫離基板的垂直速度隨時(shí)間增長(zhǎng)而增加，t=2.4時(shí)開始脫離基板，；當(dāng)Re=100時(shí)，氣泡脫離明顯提前，t=1.3時(shí)氣泡就能夠開始脫離，氣泡脫離基板的垂直速度也隨時(shí)間增長(zhǎng)而增加，且相比Re=10時(shí)，氣泡脫離基板的垂直速度較大。

圖3 不同Re數(shù)下氣泡脫離基板的垂直速度變化規(guī)律(Bo=1，θ=60°)Fig.3 Variations in vertical velocities of bubble detachment from substrate at different Reynolds numbers at Bo=1 and θ=60°

圖4是不同Re數(shù)下氣泡形狀隨時(shí)間的變化(Bo=1，θ=60°)。當(dāng)Re=1時(shí)，氣泡雖然不能脫離，但是由于浮力的作用，氣泡形狀略有變化，在豎直方向上被拉伸；當(dāng)Re=10時(shí)，氣泡能夠沿豎直方向脫離，氣泡先是在豎直方向上被拉伸，然后在t=2.4時(shí)刻氣泡開始脫離，氣泡在上升過(guò)程中基本為半球形，形狀變化不劇烈；當(dāng)Re=100時(shí)，氣泡在t=1.3時(shí)刻就開始脫離基板，且氣泡在上升過(guò)程中形狀變化較為劇烈，脫離之后氣泡變得更加扁平?？傮w來(lái)說(shuō)，當(dāng)Bo數(shù)和接觸角不變時(shí)，Re數(shù)越大，氣泡越容易脫離，氣泡脫離基板的垂直速度越大，形狀變化越劇烈。

圖4 不同Re數(shù)下氣泡形狀的變化(Bo=1，θ=60°)Fig.4 Evolution of bubble shapes at different Reynolds numbers at Bo=1 and θ=60°

2.3 Bo數(shù)對(duì)氣泡脫離基板的影響

圖5為不同Bo數(shù)下氣泡脫離基板的垂直速度的變化規(guī)律 (Re=10，θ=60°)?？梢钥闯觯跉馀菝撾x基板前，氣泡的垂直速度隨著時(shí)間推移而增加，但是Bo=1時(shí)的垂直速度較小；在Bo=5或Bo=10時(shí)，垂直速度隨著時(shí)間幾乎是線性的增加。而在氣泡脫離基板后，Bo=1時(shí)，氣泡速度持續(xù)增加，而Bo=5或Bo=10時(shí)，氣泡速度達(dá)到穩(wěn)定值。總體來(lái)說(shuō)，當(dāng)Re數(shù)和接觸角不變時(shí)，Bo數(shù)越大，氣泡脫離基板的垂直速度越大，但是Bo數(shù)到達(dá)一定程度時(shí)，其對(duì)氣泡的影響會(huì)很小。

圖5 不同Bo數(shù)下氣泡脫離基板的垂直速度變化規(guī)律(Re=10，θ=60°)Fig.5 Variations in vertical velocities of bubble detachment from substrate at different Bond numbers at Re=10 and θ=60°

圖6為不同Bo數(shù)下時(shí)氣泡脫離基板過(guò)程中的形狀變化(Re=10，θ=60°)?？梢钥闯觯诓煌腂o數(shù)下，氣泡均能夠從基板脫離，且脫離的時(shí)間都在t=2.4左右。在氣泡脫離基板之前，氣泡在浮力的作用下被逐漸拉伸，在Bo=1時(shí)，氣泡能夠完全脫離基板，而在Bo=5和Bo=10時(shí)，由于氣泡表面張力比Bo=1時(shí)要低，在氣泡周圍液體慣性力的作用下，氣泡的底部會(huì)斷裂，從而使小部分氣體附著在基板上，而其余部分脫離基板，這也是圖4中Bo=5和10的時(shí)候，氣泡的垂直速度與Bo=1時(shí)不同的原因。

圖6 不同Bo數(shù)下氣泡形狀的變化(Re=10，θ=60°)Fig.6 Evolution of bubble shapes at different Bond numbers at Re=10 and θ=60°

2.4 接觸角大小對(duì)氣泡脫離基板的影響

圖7為不同接觸角下氣泡脫離基板的垂直速度變化曲線(Re=10,Bo=1)?？梢钥闯?，接觸角的大小對(duì)氣泡脫離基板有很大的影響。當(dāng)接觸角θ=120°時(shí)，氣泡無(wú)法成功脫離基板；當(dāng)θ=90°時(shí)，氣泡能夠脫離基板，但是需要較長(zhǎng)的時(shí)間，在t=9.7才能開始脫離，剛開始?xì)馀菝撾x基板的垂直速度的增長(zhǎng)較為緩慢；而當(dāng)θ=60°時(shí)，氣泡脫離基板的垂直速度隨時(shí)間的增長(zhǎng)迅速增大，在t=2.4時(shí)就開始脫離基板。

圖7 不同接觸角下氣泡脫離基板的垂直速度變化規(guī)律(Re=10，Bo=1)Fig.7 Variations in vertical velocities of bubble detachment from substrate at different contact angles at Re=10 and Bo=1

圖8 不同接觸角下氣泡形狀的變化(Re=10，Bo=1)Fig.8 Evolution of bubble shapes at different contact angles at Re=10 and Bo=1

圖8為不同接觸角下氣泡形狀的變化規(guī)律(Re=10，Bo=1)。可以看出，接觸角θ=120°時(shí)，氣泡無(wú)法成功脫離，氣泡形狀也基本沒有變化，豎直方向并沒有明顯地被拉伸；接觸角θ=90°時(shí)，氣泡在脫離基板的過(guò)程中由于浮力的作用逐漸被拉伸，而后逐漸脫離基板；接觸角θ=60°時(shí)，由于氣泡與基板接觸面積的減小，氣泡更容易脫離。

3 總結(jié)與展望

本文利用Level Set方法，針對(duì)氣泡脫離基板的運(yùn)動(dòng)建立瞬態(tài)二維軸對(duì)稱模型，分別研究Re數(shù)、Bo數(shù)和接觸角的大小對(duì)氣泡脫離基板規(guī)律的影響，所得結(jié)論如下：

1) 隨著Re數(shù)的增大，氣泡所受的慣性力越大，氣泡越容易脫離基板，即氣泡脫離基板的垂直速度增長(zhǎng)越快，開始脫離的時(shí)間越早。

2) 隨著Bo數(shù)的增大，在氣泡脫離之前，氣泡的垂直速度增加得越快，但是對(duì)于氣泡脫離的時(shí)間影響不大。Bo數(shù)增大時(shí)，氣泡從完全脫離變成部分脫離。

3) 氣泡與基板的接觸角越大，與基板的接觸面積越大，與基板底部和附近液體的作用力就越大，因而越難以脫離。這也進(jìn)而說(shuō)明基板的潤(rùn)濕性對(duì)于氣泡脫離有很大的影響。