閆紅杰，趙國建，劉柳，段珺皓(.中南大學(xué) 能源科學(xué)與工程學(xué)院，湖南 長沙，40083；.中南大學(xué) 粉末冶金研究院，湖南 長沙，40083)

靜止水中單氣泡形狀及上升規(guī)律的實(shí)驗(yàn)研究

閆紅杰1，趙國建1，劉柳1，段珺皓2
(1.中南大學(xué) 能源科學(xué)與工程學(xué)院，湖南 長沙，410083；
2.中南大學(xué) 粉末冶金研究院，湖南 長沙，410083)

為了研究靜止液體中單氣泡的形狀及運(yùn)動(dòng)規(guī)律，以靜止水中運(yùn)動(dòng)的單氣泡為研究對(duì)象，采用高速攝影法和數(shù)字圖像處理技術(shù)，通過實(shí)驗(yàn)研究單氣泡在水中生成及上升的過程。提出用氣泡圓形度C描述氣泡形狀，采用不同量綱一參數(shù)對(duì)氣泡圓形度進(jìn)行預(yù)測，并分析氣泡運(yùn)動(dòng)軌跡、終速度的變化規(guī)律。研究結(jié)果表明：采用We數(shù)或Re數(shù)對(duì)氣泡圓形度進(jìn)行預(yù)測效果最好；氣泡在運(yùn)動(dòng)過程中的軌跡呈直線型、之字型、螺旋型3種形態(tài)，且終速度隨氣泡變形的加劇而變化，氣泡形狀為球形時(shí)運(yùn)動(dòng)速度最小，為球帽形時(shí)速度最大。

高速攝影法；氣泡形狀；氣泡圓形度；氣泡運(yùn)動(dòng)軌跡；氣泡終速度

氣泡在水中的運(yùn)動(dòng)廣泛存在于工業(yè)過程中，如污水處理中的曝氣工藝、冶金過程中的浮選工藝等都需要考慮氣泡的生成及運(yùn)動(dòng)過程，而此類問題與氣泡的動(dòng)力學(xué)特性密切相關(guān)，因此，深入研究氣泡形狀及其運(yùn)動(dòng)過程的一般規(guī)律，對(duì)于解決許多實(shí)際工程中的問題具有重要的理論意義[1]。氣泡形狀是氣泡動(dòng)力學(xué)特性的重要參數(shù)之一，與液體物性參數(shù)、氣泡大小等密切相關(guān)。SURESH等[2]通過實(shí)驗(yàn)研究了氣泡形狀與氣體流量的關(guān)系，結(jié)果表明隨著氣體流量的增大，氣泡運(yùn)動(dòng)受表面張力影響變小，受慣性力影響變大，其形狀由球形向半球形、球帽形轉(zhuǎn)變；CLIFT等[3]研究了不同物性參數(shù)液體內(nèi)氣泡的形變，繪制出氣泡形狀與Re(雷諾數(shù))、Eo(奧托斯數(shù))和Mo(莫頓數(shù))等量綱一參數(shù)間關(guān)系的氣泡相圖。WELLEK等[4-9]通過將量綱一參數(shù)(Eo，We，Re，Ta(塔達(dá)基數(shù)))與氣泡縱橫比E(氣泡高度與寬度的比值)關(guān)聯(lián)，得到一系列經(jīng)驗(yàn)公式。SHEN等[10]在球形氣泡的測量方法中提出了“氣泡圓形度”(氣泡區(qū)域等面積圓的周長和該區(qū)域周長的比值)的概念，張東衡等[11-12]在氣泡形態(tài)參數(shù)檢測、特征參數(shù)提取的研究中采用氣泡圓形度C對(duì)氣泡形狀進(jìn)行了描述。以上研究表明，可采用氣泡縱橫比E與圓形度C對(duì)氣泡形狀定量描述，但與縱橫比E不同，當(dāng)前的研究并沒有將圓形度C與量綱一參數(shù)關(guān)聯(lián)，因此，對(duì)氣泡圓形度的預(yù)測是本文研究的重點(diǎn)。另外，氣泡形狀、液相物性參數(shù)等對(duì)氣泡的上升過程有很大的影響。LESAGE等[13]對(duì)氣泡脫離體積進(jìn)行了實(shí)驗(yàn)研究，得到氣泡脫離體積與氣體流量、噴嘴內(nèi)徑及氣泡形狀等參數(shù)密切相關(guān)；AMIRNIA等[14]采用高速攝影法對(duì)液體中氣泡的上升過程進(jìn)行了實(shí)驗(yàn)研究，結(jié)果表明直徑較小時(shí)氣泡呈直線上升，隨著直徑的增大，氣泡呈之字或螺旋上升；PETERS等[15]通過實(shí)驗(yàn)研究了自來水中運(yùn)動(dòng)的氣泡，得到氣泡運(yùn)動(dòng)速度與形狀及其表面振蕩有很大的關(guān)系。但氣泡形狀、液相物性參數(shù)等與氣泡運(yùn)動(dòng)過程間的規(guī)律尚無明確的結(jié)論，有待進(jìn)一步研究。大量研究表明[16-20]，高速攝影法具有非接觸、對(duì)流場無干擾且可直接獲得流體流動(dòng)的原始圖像等優(yōu)點(diǎn)，是當(dāng)前氣泡動(dòng)力學(xué)參數(shù)測量領(lǐng)域的重要手段。本文以水中運(yùn)動(dòng)的單氣泡為對(duì)象，采用高速攝影法及圖像處理算法提取氣泡特征參數(shù)，對(duì)氣泡形狀、運(yùn)動(dòng)軌跡及終速度進(jìn)行研究，采用氣泡圓形度C對(duì)氣泡形狀進(jìn)行定量描述，重點(diǎn)分析Eo，We，Re及Ta等參數(shù)與圓形度的關(guān)系，并給出圓形度的預(yù)測公式。

1　實(shí)驗(yàn)方案

通過分析氣泡運(yùn)動(dòng)的各種測量方式，確立基于高速攝影法的氣泡行為測量方案。實(shí)驗(yàn)裝置如圖1所示。測量系統(tǒng)由三大結(jié)構(gòu)即氣泡發(fā)生裝置、高速攝影測量系統(tǒng)、圖像處理系統(tǒng)組成。

實(shí)驗(yàn)裝置主體為1個(gè)長方體(長×寬×高為150 mm×150 mm×500 mm)有機(jī)玻璃容器，側(cè)面壁厚為8 mm，底面壁厚為10 mm。在容器底面中心開直徑為30 mm的圓孔，于圓孔中塞入橡皮塞，并在橡皮塞中心插入不同內(nèi)徑的平口不銹鋼注射針或玻璃噴嘴做氣體噴嘴，噴嘴高度為30 mm。根據(jù)不同的實(shí)驗(yàn)要求，在有機(jī)玻璃容器中加入水作為液相，每次實(shí)驗(yàn)的液相高度均為300 mm，液相溫度等于室溫后即可進(jìn)行實(shí)驗(yàn)。實(shí)驗(yàn)采用空氣作為氣相。實(shí)驗(yàn)過程如下：注射泵產(chǎn)生的空氣經(jīng)軟管與噴嘴相連，在噴嘴出口處生成氣泡并向上運(yùn)動(dòng)，高速攝影儀拍攝氣泡的運(yùn)動(dòng)圖像并傳入計(jì)算機(jī)，通過數(shù)字圖像處理得到氣泡的特征參數(shù)，進(jìn)而分析并討論氣泡的上升規(guī)律。

圖1　實(shí)驗(yàn)裝置示意圖Fig.1　Schematic diagram of experimental system

實(shí)驗(yàn)所用水的密度采用液體密度計(jì)測量，黏度采用DV-II型旋轉(zhuǎn)黏度計(jì)測量，表面張力采用鉑金環(huán)法測量。實(shí)驗(yàn)中所用水的物性參數(shù)見表1。

表1　液相物性參數(shù)Table 1　Phisical parameters of liquid phase

2　實(shí)驗(yàn)結(jié)果分析

2.1氣泡形狀

圖2中的氣泡圖像是噴嘴內(nèi)徑Do=1.283 mm，氣體流量q=5.8 mL/min時(shí)，水中1個(gè)氣泡從開始生成到脫離噴嘴再運(yùn)動(dòng)至穩(wěn)定上升的過程。由圖2可知：氣泡在成長階段(從開始生成至脫離噴嘴)呈球形；隨著時(shí)間推移，氣泡形狀向橢球形轉(zhuǎn)變，變形較小，體積變大；在上升階段(從脫離噴嘴至穩(wěn)定上升)，由于表面振蕩作用，氣泡變形加劇，從橢球形迅速向球帽狀、冠狀轉(zhuǎn)變，最后穩(wěn)定在冠狀。

圖3所示為水中不同形狀的氣泡圖像。從圖3可見：隨著氣泡直徑的增大及變形幅度的加劇，氣泡縱橫比整體呈越來越小的趨勢。但隨著氣泡直徑繼續(xù)增大，氣泡變形幅度較大，如圖3(d)~3(f)所示。氣泡表面出現(xiàn)不規(guī)則的變形，導(dǎo)致氣泡縱橫比與氣泡形狀的關(guān)聯(lián)性減弱，此時(shí)采用縱橫比對(duì)氣泡形狀進(jìn)行描述不再準(zhǔn)確。本文用氣泡圓形度C來定量描述氣泡形狀，并對(duì)圓形度C進(jìn)行預(yù)測。

2.2氣泡圓形度

氣泡圓形度C用來描述氣泡接近圓形的程度。針對(duì)經(jīng)圖像處理后得到的氣泡邊緣圖像(如圖4所示)，氣泡圓形度等于與氣泡區(qū)域等面積圓的周長和該區(qū)域周長的比值，計(jì)算式為

式中：S為區(qū)域面積，m2;P為區(qū)域周長， m。C范圍為0＜C≤1，C越大，氣泡形狀越接近圓形，當(dāng)C=1時(shí)，氣泡形狀為圓形。

2.2.1氣泡圓形度與量綱一參數(shù)的關(guān)系

氣泡在液體中運(yùn)動(dòng)時(shí)，其形狀主要受浮力、表面張力、黏性力和慣性力的作用，因此，可用Eo，We，Re及Ta等參數(shù)對(duì)氣泡的形狀進(jìn)行描述。Eo，We，Re 及Ta數(shù)的計(jì)算式為

圖2　水中氣泡生成及上升圖像Fig.2　Generation and rising image of bubble in water

圖3　水中氣泡形狀Fig.3　Bubble shapes in water

圖4　氣泡邊緣圖像Fig.4　Image of bubble edge

式中：g為重力加速度，m/s2； d為氣泡直徑， m；lρ為液體的密度，kg/m3；gρ為氣體的密度，kg/m3；σ為表面張力，N/m；Tv為氣泡終速度，m/s；lμ為液體黏度，kg/(s2·m)。

水中氣泡圓形度C與Eo，We，Re及Ta等參數(shù)間的關(guān)系如圖5所示。

圖5(a)所示為水中氣泡圓形度C與Eo數(shù)間的關(guān)系。從圖5(a)可以看出：當(dāng)Eo＜0.1時(shí)，氣泡圓形度保持或接近于1，氣泡形狀為球形或接近球形；當(dāng)Eo＞0.1時(shí)，隨著氣泡直徑的增大，氣泡變形加劇，隨著Eo數(shù)的增大，氣泡圓形度分布開始變得分散，即圓形度與Eo數(shù)不能很好地關(guān)聯(lián)。

氣泡圓形度C與We數(shù)間關(guān)系如圖5(b)所示。由圖5(b)可知：當(dāng)We＜0.02時(shí)，氣泡圓形度保持為1，氣泡形狀為球形；當(dāng)We＞0.02時(shí)，圓形度隨We數(shù)的增大而變小，氣泡形狀轉(zhuǎn)變?yōu)闄E球形；當(dāng)We＞0.8時(shí)，氣泡變形加劇，圓形度隨We數(shù)的增大開始分散，但與圖5(a)中Eo數(shù)與圓形度的關(guān)系相比，分散程度小很多，圓形度與We數(shù)的關(guān)聯(lián)性很強(qiáng)，說明表面張力和慣性力對(duì)氣泡的變形有重要的影響。

圖5　水中氣泡圓形度C與Eo，We，Re和Ta數(shù)的關(guān)系Fig.5　Relationship between bubble circularity in water and Eo,We,Re and Ta number

圖5(c)所示為氣泡圓形度C與Re數(shù)間的關(guān)系。由圖5(c)可知：Re數(shù)對(duì)氣泡圓形度的影響與Eo數(shù)和We數(shù)的影響類似：在Re數(shù)大于200時(shí)，氣泡直徑增大，變形加劇，氣泡形狀從球形向橢球形、球帽形轉(zhuǎn)變，氣泡圓形度與Re數(shù)的關(guān)系開始分散，其分散程度說明Re數(shù)對(duì)氣泡變形的影響大于Eo數(shù)對(duì)其的影響。

氣泡圓形度C與Ta數(shù)間關(guān)系如圖5(d)所示。從圖5(d)可以看出：當(dāng)Ta數(shù)較小時(shí)，氣泡直徑較小，氣泡形狀保持為球形，氣泡圓形度基本等于1；隨著Ta數(shù)的逐漸增大，圓形度逐漸減小，在Ta達(dá)到0.7后，圓形度隨Ta數(shù)的增大開始出現(xiàn)分散，其分散程度較圖5(a)中的小，較圖5(b)中的大，說明在分析氣泡變形時(shí)應(yīng)綜合分析考慮黏性力、慣性力和表面張力的影響。

2.2.2氣泡圓形度的預(yù)測

從圖5可看出：隨著Eo，We，Re和Ta的增大，氣泡圓形度C總體上呈現(xiàn)減小的趨勢，即氣泡圓形度C與Eo，We，Re和Ta數(shù)呈反函數(shù)的關(guān)系，因此，可構(gòu)建圓形度與Eo，We，Re和Ta數(shù)之間的關(guān)聯(lián)式，通過對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行擬合，找出氣泡圓形度C的最佳預(yù)測模型，通過多次試驗(yàn)，構(gòu)建如下經(jīng)驗(yàn)關(guān)聯(lián)式對(duì)氣泡圓形度C進(jìn)行預(yù)測：

式中：a，b和c為待定系數(shù)；IL為量綱一參數(shù)。分別采用Eo，We，Re或Ta數(shù)作為式(3)中的參數(shù)IL，以實(shí)驗(yàn)數(shù)據(jù)為基礎(chǔ)，通過最小二乘法對(duì)式(3)進(jìn)行擬合，可得到待定系數(shù)(a，b和c)、預(yù)測值與實(shí)驗(yàn)值之間的絕對(duì)偏差平均值、標(biāo)準(zhǔn)差、殘差平方和、相關(guān)系數(shù)等參數(shù)如表2所示。表2中，S為平均偏差，Std為標(biāo)準(zhǔn)差，SSE為殘差平方和，r為相關(guān)系數(shù)。

從表2可以看出：采用Re數(shù)對(duì)氣泡圓形度進(jìn)行預(yù)測時(shí)，預(yù)測值和實(shí)驗(yàn)值的絕對(duì)偏差平均值、標(biāo)準(zhǔn)偏差、殘差平方和最小，相關(guān)系數(shù)最大。因此，在水中，Re數(shù)能較好地預(yù)測氣泡圓形度、描述氣泡的形狀，預(yù)測公式為

表2　液體水待定系數(shù)、平均偏差、標(biāo)準(zhǔn)差、殘差平方和與相關(guān)系數(shù)Table 2　Undeterminined coefficient,mean deviation,standard deviation,residual sum of squares and correlation coefficient of liquid water

采用Re數(shù)和We數(shù)預(yù)測氣泡圓形度的預(yù)測值與實(shí)驗(yàn)值，如圖6(a)和6(b)所示。由圖6可知：形狀保持或接近圓形的氣泡，采用We數(shù)描述氣泡形狀相對(duì)采用Re數(shù)描述氣泡形狀效果較好，氣泡圓形度的預(yù)測值與實(shí)驗(yàn)值擬合度更高。這是因?yàn)榇藭r(shí)氣泡直徑較小，氣泡的運(yùn)動(dòng)過程受表面張力作用比較大，而僅僅采用Re數(shù)來預(yù)測圓形度，并沒有考慮到表面張力的影響，導(dǎo)致預(yù)測值與實(shí)驗(yàn)值相對(duì)誤差偏大。于是，在氣泡形狀的描述中加入表面張力的因素，綜合考慮黏性力、慣性力及表面張力等對(duì)氣泡變形的影響，在預(yù)測公式中將Re數(shù)和We數(shù)結(jié)合對(duì)預(yù)測模型進(jìn)行修正。因此，水中氣泡圓形度修正預(yù)測模型為

圖6　水中采用Re和We數(shù)計(jì)算氣泡圓形度C預(yù)測值與實(shí)驗(yàn)值對(duì)比Fig.6　Comparison of predicted values and experiment value of bubble circularity according to Re and We number in water

對(duì)式(5)進(jìn)行最小二乘擬合，得到相關(guān)參數(shù)如表3所示。

對(duì)比表3和表2中各項(xiàng)擬合參數(shù)可知：采用We數(shù)及Re數(shù)對(duì)氣泡圓形度C進(jìn)行預(yù)測時(shí)，預(yù)測值和實(shí)驗(yàn)值的平均偏差、標(biāo)準(zhǔn)差、殘差平方和以及相關(guān)系數(shù)都比單獨(dú)采用Eo，We，Re或Ta數(shù)的預(yù)測模型的小。

綜上所述，利用氣泡圓形度C來定量描述氣泡形狀，其預(yù)測公式為

此式適用范圍為：1＜Re＜3 041，0.000 2＜We＜8.6。

表3　液體水待定系數(shù)、平均偏差、標(biāo)準(zhǔn)偏差、殘差平方和與相關(guān)系數(shù)Table 3　Undeterminined coefficient,mean deviation,standard deviation,residual sum of squares and correlation coefficient

2.3氣泡運(yùn)動(dòng)軌跡與終速度

2.3.1氣泡運(yùn)動(dòng)軌跡

圖7　d=0.83 mm時(shí)氣泡運(yùn)動(dòng)軌跡Fig.7　Bubble rising path when d=0.83 mm

圖8　d=1.62 mm時(shí)氣泡運(yùn)動(dòng)軌跡Fig.8　Bubble rising path when d=1.62 mm

通過分析同一內(nèi)徑噴嘴(Do=7.7 mm)生成的氣泡在二維空間內(nèi)水平位置X、豎直速度vv和水平速度vh與氣泡上升高度Z的變化關(guān)系，推導(dǎo)出氣泡在三維空間內(nèi)運(yùn)動(dòng)軌跡主要有直線型、之字型及螺旋型3種。圖7~9所示分別為氣泡脫離噴嘴后運(yùn)動(dòng)軌跡為直線型、之字型及螺旋型的水平位置X、豎直速度vv和水平速度vh隨氣泡上升而變化的過程。由圖7可知：當(dāng)氣泡直徑較小時(shí)，氣泡形狀為球形或近似球形，運(yùn)動(dòng)軌跡為直線型，豎直速度vv和水平速度vh保持不變，水平位置X與豎直高度Z呈線性關(guān)系；隨著氣泡直徑增大，氣泡形狀向橢球形轉(zhuǎn)變，運(yùn)動(dòng)軌跡表現(xiàn)為之字型。由圖8可知：水平位置X隨上升高度Z增大呈周期性震蕩，豎直方向速度vv和水平方向速度vh也呈周期性震蕩，且vv的振蕩頻率是X與vh振蕩頻率的2倍；隨著氣泡直徑繼續(xù)增大，氣泡變形加劇，運(yùn)動(dòng)軌跡表現(xiàn)為螺旋型，如圖9所示，與之字型軌跡相比，豎直速度vh保持穩(wěn)定。

圖9　d=1.96 mm時(shí)氣泡運(yùn)動(dòng)軌跡Fig.9　Bubble rising path when d=1.96 mm

2.3.2氣泡終速度

氣泡終速度是指氣泡上升過程中最終穩(wěn)定的速度，圖10所示為氣泡終速度隨氣泡直徑增大的變化趨勢圖。由圖10可知：當(dāng)d＜0.83 mm時(shí)，終速度隨氣泡直徑的增大呈線性增大，氣泡運(yùn)動(dòng)主要受黏性阻力和浮力控制；隨著氣泡直徑繼續(xù)增大(0.83＜d＜6.00 mm)，終速度出現(xiàn)多值化，氣泡運(yùn)動(dòng)主要受表面張力和慣性力控制；當(dāng)d＞6.00 mm時(shí)，終速度隨氣泡直徑增大而增大并趨于穩(wěn)定，其運(yùn)動(dòng)主要受慣性力控制。

2.3.3氣泡形狀對(duì)運(yùn)動(dòng)軌跡和終速度的影響

氣泡形狀、液體物性參數(shù)等對(duì)氣泡在水中的運(yùn)動(dòng)行為均有不同程度的影響，其中，氣泡形狀的影響可用圓形度C來定量分析。

水中氣泡運(yùn)動(dòng)軌跡、終速度與圓形度間關(guān)系如圖11所示。由圖11可知：隨著圓形度減小即隨著氣泡變形加劇，氣泡運(yùn)動(dòng)軌跡從直線型向之字型、螺旋型轉(zhuǎn)變；當(dāng)0.99＜C＜1時(shí)，形狀從球形向橢球形變化，此時(shí)氣泡運(yùn)動(dòng)主要受黏性阻力和浮力控制，其表面振蕩不明顯，形狀相對(duì)穩(wěn)定，運(yùn)動(dòng)軌跡保持為直線型；隨著氣泡變形的加劇，當(dāng)0.92＜C＜0.99時(shí)，氣泡形狀向球帽形轉(zhuǎn)變，其運(yùn)動(dòng)主要受表面張力和慣性力控制，此時(shí)隨著氣泡直徑的增大，氣泡表面振蕩加強(qiáng)，氣液交界面不穩(wěn)定，運(yùn)動(dòng)軌跡從直線型向之字型轉(zhuǎn)變；隨著圓形度繼續(xù)變小，氣泡變形較大，氣泡表面振蕩強(qiáng)烈，形狀為橢球形或球帽形，其運(yùn)動(dòng)主要受慣性力控制，運(yùn)動(dòng)軌跡為之字型或螺旋型。

圖10　氣泡終速度與直徑的關(guān)系Fig.10　Relationship between bubble final velocity and diameter

圖11　氣泡運(yùn)動(dòng)軌跡、終速度與圓形度C的關(guān)系Fig.11　Relationship among bubble trajectory,final velocity and circularity

隨著氣泡形狀從球形、橢球形向球帽形變化，氣泡運(yùn)動(dòng)終速度出現(xiàn)相應(yīng)變化：當(dāng)圓形度接近于1，即氣泡形狀為球形、橢球形時(shí)，氣泡直徑較小，其運(yùn)動(dòng)受黏性阻力影響較大，氣液交界面的振蕩對(duì)氣泡運(yùn)動(dòng)影響較弱，終速度隨圓形度的減小有呈線性減小的趨勢；隨著氣泡變形加劇即隨著圓形度繼續(xù)減小，氣泡形狀向球帽形轉(zhuǎn)變，氣泡終速度與圓形度呈多值關(guān)系。這是因?yàn)殡S著氣泡直徑的增大，氣泡表面振蕩加強(qiáng)，表面張力與慣性力同時(shí)影響氣泡的運(yùn)動(dòng)，相同直徑的氣泡其表面振蕩幅度不一，圓形度可大可小，從而導(dǎo)致終速度不同。當(dāng) C＜0.92時(shí)，氣泡形狀不穩(wěn)定，直徑較大，表面振蕩劇烈，其運(yùn)動(dòng)受慣性力影響較大，終速度隨氣泡表面的振蕩出現(xiàn)細(xì)微變化，整體上趨于穩(wěn)定。

3　結(jié)論

1)在靜止水中，氣泡形狀從球形向橢球形、球帽形轉(zhuǎn)變，得到圓形度的預(yù)測公式。

2)氣泡在上升過程中，運(yùn)動(dòng)軌跡呈直線型、之字型或螺旋型3種形態(tài)，且隨著氣泡變形的加劇，其運(yùn)動(dòng)軌跡從直線型向之字型、螺旋型轉(zhuǎn)變。

3)氣泡終速度隨著氣泡形狀的變化而變化，氣泡形狀為球形時(shí)終速度最小，為球帽形時(shí)終速度最大。

[1]郭烈錦.兩相與多相流動(dòng)力學(xué)[M].西安:西安交通大學(xué)出版社,2002:1-13. GUO Liejin.Fluid dynamics of two phase and multiphase system[M].Xi’an:Xi’an Jiaotong University Press,2002:1-13.

[2]SURESHM,MANI A.Experimentalstudiesonbubble characteristicsforR134a–DMFbubbleabsorber[J]. Experimental Thermal and Fluid Science,2012,39(39):79-89.

[3]CLIFT R,GRACE J R,WEBER M E.Bubbles,drops,and particles[M].New York,America:Courier Dover Publications, 2005:16-28.

[4]WELLEK R M,AGRAWAL A K,SKELLAND A H P.Shape of liquid drops moving in liquid media[J].AIChE Journal,1966, 12(5):854-862.

[5]OKAWA T,TANAKA T,KATAOKA I,et al.Temperature effect on single bubble rise characteristics in stagnant distilled water[J]. International Journal of Heat and Mass Transfer,2003,46(5): 903-913.

[6]CELATA GP,D’ANNIBALEF,DI MARCOP,etal. Measurements of rising velocity of a small bubble in a stagnant fluidinone-andtwo-componentsystems[J].Experimental Thermal and Fluid Science,2007,31(6):609-623.

[7]DONG H,WANG X,LIU L,et al.The rise and deformation of a single bubble in ionic liquids[J].Chemical Engineering Science, 2010,65(10):3240-3248.

[8]LIU L,YAN H,ZHAO G.Experimental studies on the shape and motion of air bubbles in viscous liquids[J].Experimental Thermal and Fluid Science,2015,62:109-121.

[9]TADAKI T,MAEDA S.On the shape and velocity of single air bubbles rising in various liquids[J].Kagaku Kogaku Ronbunshu, 1961,25(4):254-264.

[10]SHEN X,NAKAMURA H.Spherical-bubble-based four-sensor probesignalprocessingalgorithmfortwo-phaseflow measurement[J].International Journal of Multiphase Flow,2014, 60(2):11-29.

[11]張東衡,唐志航,葉鴻明,等.一種氣液兩相流氣相參數(shù)圖像檢測方法[J].計(jì)算機(jī)測量與控制,2006,14(5):597-599. ZHANG Dongheng,TANG Zhihang,YE Hongming,et al. Measurement method of gas-liquid two-phase flow parameters based on digital image processing[J].Computer Measurement& Control,2006,14(5):597-599.

[12]薛婷,孟欣東,張濤.氣液兩相流中氣泡形態(tài)及運(yùn)動(dòng)特征參數(shù)提取[J].光電子·激光,2010,21(8):1218-1221. XUE Ting,MENG Xindong,ZHANG Tao.Extraction of bubble shape and motion feature parameters in the gas-liquid two-phase flow[J].JournalofOptoelectronicsLaser,2010,21(8): 1218-1221.

[13]LESAGE F J,MAROIS F.Experimental and numerical analysis ofquasi-staticbubblesizeandshapecharacteristicsat detachment[J].International Journal of Heat and Mass Transfer, 2013,64(3):53-69.

[14]AMIRNIA S,DE BRUYN J R,BERGOUGNOU M A,et al. Continuous rise velocity of air bubbles in non-Newtonian biopolymer solutions[J].Chemical Engineering Science,2013, 94(5):60-68.

[15]PETERS F,ELS C.An experimental study on slow and fast bubbles in tap water[J].Chemical Engineering Science,2012, 82(1):194-199.

[16]李娜.氣泡運(yùn)動(dòng)的數(shù)字圖像測量方法研究[D].北京:首都經(jīng)濟(jì)貿(mào)易大學(xué)安全與環(huán)境工程學(xué)院,2011:7-32. LI Na.Research on digital image measurement method of bubble motion[D].Beijing:Capital University of Economics and Business.School of Safety and Environmental Engineering, 2011:7-32.

[17]HU Y,MENG X,SUN W.A novel quantificational assessment method of frothers effect on bubble characteristics[J].Journal of CentralSouthUniversityofTechnology,2011,18(6): 1910-1916.

[18]AVACHAT S,ZHOU M.High-speed digital imaging and computationalmodelingofdynamicfailureincomposite structuressubjectedtounderwaterimpulsiveloads[J]. International Journal of Impact Engineering,2015,77:147-165.

[19]ZHAO D,LV M,WANG P,et al.Can the plant area index of a submergedvegetationcanopybeestimatedusingdigital hemisphericalphotography?[J].AgriculturalandForest Meteorology,2014,192(4):69-77.

[20]唐朝暉,劉金平,桂衛(wèi)華,等.基于數(shù)字圖像處理的浮選泡沫速度特征提取及分析[J].中南大學(xué)學(xué)報(bào)(自然科學(xué)版),2009, 40(6):1616-1622. TANG Chaohui,LIU Jinping,GUI Weihua,et al.Froth bubbles speed characteristic extraction and analysis based on digital imageprocessing[J].JournalofCentralSouthUniversity (Science and Technology),2009,40(6):1616-1621.

(編輯陳燦華)

Experimental study on shape and rising behavior of single bubble in stagnant water

YAN Hongjie1,ZHAO Guojian1,LIU Liu1,DUAN Junhao2
(1.School of Energy Science and Engineering,Central South University,Changsha 410083,China;
2.Powder Metallurgy Research Institute,Central South University,Changsha 410083,China)

In order to study the shape and rising behavior of single bubble in stagnant liquid,the formation and rising process of single bubble in water were experimentally studied based on high-speed photography and digital image processing technology.The bubble circularity C which was used to describe the bubble shape was correlated by using different dimensionless parameters,and the bubble rising path and terminal velocity were studied.The results show that bubble circularity can be correlated well by using the dimensionless parameter We or Re number,the bubble rising paths are linear,zig-zag and helical,and the terminal velocity changes with the aggravation of bubble deformation.The velocity of the spherical bubble is smaller,and the velocity of the cap bubble is larger.

high-speed photography;bubble shape;bubble circularity;bubble rising path;bubble terminal velocity

閆紅杰，博士，教授，從事多相流流動(dòng)、熱工設(shè)備的仿真與優(yōu)化研究；E-mail:s-rfy@csu.edu.cn

O359

A

1672-7207(2016)07-2513-08

10.11817/j.issn.1672-7207.2016.07.045

2015-07-18；

2015-09-22

國家自然科學(xué)基金創(chuàng)新群體項(xiàng)目(61321003)；中南大學(xué)中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金資助項(xiàng)目(2015zzts205)；湖南省研究生科研創(chuàng)新項(xiàng)目(CX2014B063)(Project(61321003)supported by the Foundation for Innovative Research Groups of the National Natural Science Foundation of China;Project(2015zzts205)supported by the Fundamental Research Funds for the Central Universities of Central South University;Project(CX2014B063)supported by the Hunan Provincial Innovation Foundation for Postgraduate)