頂部浸沒豎直向下平口管口氣泡膨脹脫離特性

吳　晅，焦晶晶，梁盼龍，金光，武文斐

（1內蒙古科技大學能源與環境學院，內蒙古 包頭 014010；2內蒙古科技大學內蒙古自治區白云鄂博礦多金屬資源綜合利用重點實驗室，內蒙古 包頭014010）

?

頂部浸沒豎直向下平口管口氣泡膨脹脫離特性

吳晅1，2，焦晶晶1，梁盼龍1，金光1，武文斐2

（1內蒙古科技大學能源與環境學院，內蒙古 包頭 014010；2內蒙古科技大學內蒙古自治區白云鄂博礦多金屬資源綜合利用重點實驗室，內蒙古 包頭014010）

利用高速攝像機拍攝頂部浸沒豎直向下管口氣泡膨脹及脫離演變過程，對比分析管徑、平均氣速對氣泡尺寸、形成時間的影響規律。引入泡齡k表征管口氣泡位置特性，用橢球形氣泡長短軸比L表征氣泡形狀特征。研究發現，氣泡生成機理有間歇生成氣泡、含抽吸間歇生成氣泡和連續生成氣泡3種；氣泡脫離直徑隨平均氣速變化過程明顯存在氣泡脫離形態轉變點；氣泡膨脹脫離模式隨平均氣速變化分為單個氣泡形成、氣泡聚并形成兩種。平均氣速為3.1 m·s-1時管口處存在抽吸現象，此時最大泡齡kmax最小；脫離階段氣泡長短軸比L隨脫離時間變化過程明顯存在時間分界點。在分界點前，L值在2.0附近波動，氣泡主要以橢球形存在。分界點后，L值明顯增加，氣泡形狀向扁橢球體發展。

氣液兩相流；向下管口；氣泡；流體動力學

DOI：10.11949/j.issn.0438-1157.20151478

引　言

氣液兩相流動廣泛存在于能源、化工、石油、航空、冶金等工業領域中。氣泡作為流體動力學黏性流體中一個分散相，其運動是典型的兩相流動現象，運動特性極其復雜，在20世紀90年代就已開展了對氣泡生長行為及影響因素的研究［1-4］。氣泡多是從浸沒管口或孔生成，浸沒方式可分為管口（或孔）向上底部浸沒、管口（或孔）水平布置側面浸沒和管口（或孔）向下頂部浸沒3種［5］。

前兩種浸沒方式應用廣泛，研究成果資料較多且細致。劉柳［6］采用高速攝像機獲取不同溶液中垂直上升管氣泡運動過程，得到氣泡運動軌跡、脫離管口后形狀分布，但未涉及氣泡直徑變化及生成時間。徐炯等［7］采用高速相機獲取靜止水下單個氣泡上升過程，指出氣泡上升運動具有勻速特性，氣泡體積增加導致氣泡變形等，未對氣泡直徑變化進行研究。Montante等［8］對攪拌槽氣液兩相流動進行PIV實驗研究，采用兩臺CCD同時拍攝氣泡和液體流場的運動，分析得到氣泡尺寸分布和兩相速度場分布，關于氣泡形態分布未涉及。Zahedi等［9］對單個氣泡形成、脫離、上升等進行實驗和數值模擬，得到不同液體中管徑對氣泡直徑、形成時間的影響規律，數值模擬結果與實驗結果吻合較好。宋慶唐等［10］采用高速攝像法拍攝氣液鼓泡床內不同氣速下氣泡分布，指出氣泡尺寸分布具有正態分布特征，氣速可促進氣含率、氣泡聚并和破碎，但其氣泡尺寸是以氣泡個數表示的。Liu等［11］建立數學計算模型研究孔口方向和流體流型對氣泡生成的影響。

相比前兩種浸沒方式，管口向下頂部浸沒方式也存在于一些工業過程中。鋼鐵冶金氧氣底吹工藝中Savard-Lee雙重套管噴槍在底部噴槍頂端形成蘑菇頭狀氣泡確保工藝順利進行［12］。Huda等［13］采用CFD研究攪拌槽頂部浸沒氣體射流淹沒深度、氣體流量等對管口氣泡形成特性的影響，給出預測氣體噴射深度的經驗公式，但缺乏實驗結果的有力支撐。在中間包氣體噴吹中，詹樹華等［14］采用數值模擬方法研究氣液兩相鼓泡流特征，模擬結果與實驗結果一致。水煤漿氣化爐激冷室中將高溫合成氣經頂端浸沒下降管通入液池來實現高溫合成氣的粗洗滌［15］，吳晅等［16］采用VOF模型模擬再現下降管口氣泡生成、破碎等一系列運動過程，其模型中管口尺寸較大，文中氣泡為較大的氣塊。王沖等［17］對貧化電爐還原油槍口處氣液兩相運動進行實驗研究，得出不同流量下攪拌深度和攪拌半徑的分布規律，但其實驗中采用空氣-水的氣液兩相混合頂吹模型。Rivera-Salinas等［18］用CFD模擬研究冶金精煉過程頂部浸沒管口氣液固三相流動過程，得到管口氣泡形態和管口周圍流場分布特性，同時給出固體顆粒與反應效率關系，模擬結果與現有文獻吻合。楊濮亦等［19］采用尼康單反相機拍攝銅冶煉熔池中頂吹氣泡上升運動，給出不同浸沒深度下氣泡運動速率分布規律，但未涉及氣泡形狀特性的研究。

本研究將以靜止液體中頂部浸沒豎直向下平口管口氣液兩相流單個氣泡為對象。前期工作［20］已采用數值模擬的方法對氣泡受力、形成尺寸及影響因素進行分析研究，建立兩階段氣泡運動方程。而本研究將搭建氣液兩相實驗平臺，通過調節不同的進氣量和平口管管徑來改變實驗工況。采用高速攝像機拍攝得到同一浸沒深度不同直徑下管口氣泡膨脹及脫離演變過程，獲取各進氣量下氣泡生長關鍵幀圖片，并提取氣泡數據，得到平口管管徑和平口管內平均氣速對氣泡尺寸和形成時間的影響。引入泡齡表征管口氣泡位置運動特性，并利用長短軸比表征氣泡形狀分布特征。

1　實　驗

1.1實驗材料

所采用的實驗介質為壓縮空氣和水，壓縮空氣由空壓機制取，壓縮空氣儲存在儲氣罐經氣體管路供給實驗系統，水采用去離子的靜止水，實驗為常溫、常壓。

實驗所用到的平口細管為內外精密拋光的不銹鋼平口細管，其尺寸規格見表1，其中Do表示平口管外徑，Di表示平口管內徑。

1.2實驗設備與分析儀器

1.2.1實驗系統組成實驗測試系統如圖1所示，該實驗系統主要包括3個模塊，其中A為氣體供給及計量模塊，B為氣泡制備及生成模塊，C為氣泡采集及分析處理模塊。氣體流量由玻璃轉子流量計精確控制，氣體壓力由U形管量測。水箱采用透明有機玻璃制成，為提高測量結果的準確度和合理性，測試區域采用黑色背景。

表1　R不銹鋼平口細管尺寸Table 1　Size of stainless steel flat pipe

圖1　頂部浸沒管口氣泡形成特性實驗測試系統Fig.1　Experimental test system of bubble formation characteristics from downward nozzle1—gasholder； 2—valve； 3—flowmeter； 4—U-shape tube； 5—water tank；6—pool； 7—stainless steel pipe； 8—tank export； 9—light source； 10—high-speed camera； 11—computer processing system

1.2.2氣源及氣泡制取設備實驗中由靜音無油空壓機制取壓縮空氣，制得的空氣儲存在儲氣罐中。該空壓機采用無油設計，同時添加干燥器，使實驗結果不受空氣品質的影響，減小實驗誤差。

實驗中氣泡觀測過程均在矩形透明有機玻璃水箱（48 cm×48 cm×120 cm）中完成，水箱中存有去離子水，不銹鋼平口管管口以頂部浸沒方式浸沒在靜態水面以下0.1 m處。

1.2.3實驗采集設備實驗采集系統主要由高速攝像機、拍攝光源、計算機及連線組成。采用Revealer高速攝像機對氣泡生成過程進行采集。高速攝像機采用CameraLink接口，該接口傳輸速度最高可達400萬像素。高速攝像采集的分辨率為1080×720，曝光時間為600 μs，采集周期為1 μs，幀率為827 幀/秒，能夠完整記錄氣泡在膨脹脫離過程瞬間形態。鏡頭采用佳能微距鏡頭，光源采用LED太陽燈。

實驗采集過程由Revealer CL 3.0高速圖像采集系統精確控制，載有氣泡動力特性的視頻直接被存儲在電腦主機硬盤中。利用圖像處理技術分析處理視頻，獲取頂部浸沒豎直向下平口管口氣泡膨脹及脫離過程特征參數數據信息。

1.3實驗方法與實驗數據獲取

1.3.1氣泡形成膨脹及脫離實驗對水下頂部浸沒管口氣泡膨脹及脫離運動特性進行實驗研究，氣體經平口管頂部管口以鼓泡形式冒出，經出口排出水箱。實驗拍攝過程中，高速攝像機、光源和觀測目標處于同一水平面上，高速攝像機與電腦主機相連，將不同管徑不銹鋼平口管頂部浸沒在水面以下0.1 m，調節不同的平口管進氣量，在實驗工況穩定的條件下，通過圖像軟件實時獲取同一浸沒深度不同管徑下，平口管頂部浸沒向下管口氣泡生成過程視頻。

1.3.2氣泡特征參數的獲取與分析利用圖像處理軟件分析氣泡膨脹及脫離運動過程視頻，采用Matlab軟件提取視頻文件中的幀圖片，并進行簡單預處理，使之便于特征提取。再通過Photoshop軟件中尺寸測量工具獲取氣泡尺寸信息。Photoshop軟件獲取單個氣泡尺寸原理如圖2所示。

不銹鋼平口細管實際直徑D已知，可以利用Photoshop軟件標尺工具獲取幀圖片中平口細管直徑d和橢球形氣泡的長軸長a和短軸長b，由比例關系可以獲取橢球形氣泡實際長軸長A和短軸長B，其原理公式如下

圖2　單個氣泡尺寸信息獲取原理Fig.2　Principle diagram of single bubble size informationa— spheroidicity bubble's long axis in photo； b— spheroidicity bubble's short axis in photo； d— diameter of flat pipe in photo

由式（1）可計算得到

頂部浸沒平口管口氣泡生成形狀為橢球形，為便于表征氣泡尺寸大小，將測量所得橢球形氣泡尺寸信息與圓進行等效，得到等效球形氣泡脫離直徑dg，經計算簡化得到

式中，A、B分別為橢球形氣泡實際長軸長和短軸長；D為平口管實際直徑。

2　結果與分析

2.1頂部浸沒豎直向下平口管口氣泡膨脹及脫離演變過程

圖3為25號豎直向下平口管在平口管內平均氣速u=1.7 m·s-1下高速攝像機獲取的頂部浸沒管口氣泡形成演變過程。由圖3可知，頂部浸沒豎直向下平口管口氣泡形成演變過程可分成兩個階段：膨脹階段和脫離階段。設定氣泡開始形成時刻為0時刻，在膨脹階段，如圖3（a）～（d）所示，管口氣泡呈現不斷膨脹的狀態，在16.9 ms以前，氣泡主要進行體積膨脹的徑向運動，該過程中遠離管口側氣泡外邊緣為圓弧狀。在16.9 ms時，由圖3（d）可知氣泡近似為半橢圓形和圓臺形的組合，該點稱為氣泡膨脹臨界點，該點之后氣泡上邊緣開始高于管口。究其原因是頂部浸沒管口氣泡在膨脹過程中主要受到浮力、表面張力、曳力以及脫離瞬間的慣性力等作用，如圖4（a）所示。氣泡的膨脹主要由氣泡內外側壓力差決定，當氣泡內側壓力大于外側壓力時，氣泡將進行體積膨脹的徑向運動，呈現圖3（a）～（d）所示的狀態。

在16.9 ms以后，管口氣泡體積仍會繼續膨脹。但當氣泡膨脹到一定體積時，由于浮力作用遠大于氣液界面表面張力等的作用，氣泡開始進行向上的脫離運動，此時脫離階段開始。脫離階段氣泡受力狀況如圖4（b）所示。豎直向上的浮力是氣泡脫離管口的主要因素。在脫離階段，氣泡主要進行體積膨脹的徑向運動和向上的脫離運動，如圖3（e）～（h）所示。在64.1 ms時，脫離階段結束。氣泡脫離管口時刻狀態如圖3（h）所示。在64.1 ms以后氣泡脫離管口，不再接收管口氣體的供給。由圖3還可以得出，25號平口管u=1.7 m·s-1時，管口氣泡脫離階段經歷的時間（47.2 ms）遠比膨脹階段經歷的時間（16.9 ms）長。

圖3給出了氣泡形成、膨脹及脫離的整個演變過程。可以看出，氣泡是以橢球形狀存在的。但為了便于描述氣泡尺寸特征，本文提出了等效球形氣泡脫離直徑dg。本文依據實驗獲得的氣泡演變過程，給出了氣泡膨脹及脫離階段各時刻等效球形氣泡物理模型如圖5所示，其分別與圖3中各時刻變化圖相對應，圖5中d、h分別為膨脹臨界點時的氣泡狀態和脫離管口時的氣泡狀態。

圖3　25號管頂部浸沒豎直向下平口管口氣泡形成演變過程（u=1.7 m·s-1）Fig.3　Bubble formation evolution process of downward nozzle for No.25（u=1.7 m·s-1）

圖4　氣泡形成過程受力分析Fig.4　Schematic diagram of stress analysis for bubble formation process

圖5　頂部浸沒豎直向下平口管口氣泡形成演變物理模型Fig.5　Physical model diagram of bubble formation evolution process from downward nozzle of vertical flat pipe

2.2豎直向下平口管管徑對氣泡脫離直徑的影響

圖6為氣泡脫離管口前不同平口管管徑下氣泡脫離直徑隨平口管內平均氣速的變化規律。從圖中可以看出，在同一平均氣速下，氣泡尺寸隨管徑的增加而逐漸增大。對于同一管徑，氣泡脫離直徑隨平均氣速的增加先呈近似線性增加，之后平緩增加。氣泡尺寸在隨著平均氣速增加而變化的過程中明顯存在一個氣泡脫離形態轉變點，并且氣泡脫離形態轉變點所對應的氣泡脫離形態轉變氣速隨平口管管徑增加而增大。如圖6所示，各曲線氣泡脫離形態轉變點所對應的氣泡脫離形態轉變氣速分別為12.8、11.2、9.2、7.3 m·s-1。

如圖6所示，同一管徑下氣泡脫離管口時，在氣泡脫離形態轉變點之前，氣泡以單個橢球形態脫離管口。在氣泡脫離形態轉變點之后，氣泡則是以上下相連的兩個氣泡脫離管口。氣泡脫離形態轉變點之后，氣泡脫離管口時氣泡底端會形成二次氣泡，把管口脫離時刻大氣泡稱為母氣泡，把氣泡底端形成的小氣泡稱為子氣泡。圖6中左圖為18號管氣泡脫離形態轉變點之前平均氣速為4.4 m·s-1氣泡脫離管口時的形態圖，氣泡以單個橢球形態脫離管口。圖6中右圖為18號管氣泡脫離形態轉變點之后平均氣速為10.9 m·s-1時氣泡脫離管口時的形態圖，氣泡脫離管口時底端有子氣泡生成。

圖7為16號管和12號管內不同平均氣速下的氣泡膨脹脫離模式。從圖中可以看出，在同一管徑下的氣泡膨脹脫離模式可分成兩種：單個氣泡形成模式和氣泡聚并形成模式。在平均氣速較小時，管口以單個氣泡進行膨脹及脫離。而在平均氣速相對較大時，氣泡膨脹脫離過程存在聚并現象。分析圖4中氣泡膨脹脫離過程受力狀況可知，氣泡的膨脹脫離由氣泡內外側的壓力差決定。在單個氣泡形成模式下，氣泡內外側壓力差較小，氣泡膨脹脫離獲得的加速度也較小。在該模式下，由于管口前后兩個氣泡之間距離較大，使得氣泡之間相互影響較小，故呈現如圖7中左圖所示單個氣泡生成現象。而在氣泡聚并形成模式下，平均氣速較大，氣泡內外側壓力差較大，氣泡膨脹脫離所獲加速度較大，氣泡脫離管口速度也較大。后一個氣泡在脫離管口時會受到前一個氣泡尾流作用的影響，其速度明顯加快，隨后兩氣泡發生聚并現象，如圖7中右圖所示。在氣泡聚并形成模式下，氣泡的聚并是由于兩氣泡之間相互作用的結果。圖7中不同管徑下兩種模式的平均氣速分界點分別為6.9、17.2 m·s-1。

2.3豎直向下平口管膨脹氣泡位置隨時間變化分布特性

引入泡齡k來表征平口管頂部膨脹氣泡位置隨時間的變化規律，如圖8所示。泡齡k為管口底端到膨脹氣泡底端邊緣的垂直距離。規定膨脹氣泡底端低于管口底端所在平面時，k為正值；膨脹氣泡底端和管口底端在同一平面內，k為0；膨脹氣泡底端高于管口底端所在平面，則k為負值。分析圖5、圖8可知，在氣泡形成演變過程中，當k非負時，氣泡處于膨脹或脫離狀態。當k為負值時，氣泡已完成脫離并進入周圍水體中。本實驗主要研究管口氣泡k值為非負的情況。

圖6　不同管徑下氣泡脫離直徑隨平均氣速變化規律Fig.6　Change of bubble size with average gas velocity under different pipe diameter

圖7　不同管徑下氣泡膨脹脫離模式分布Fig.7　Distribution diagram of bubble's expansion and detachment model under different pipe diameter

圖8　泡齡k提取Fig.8　Schematic diagram of extracting age of bubble k

圖9給出了20號平口管頂部浸沒膨脹氣泡泡齡k隨時間t的變化規律，圖9中k0線為泡齡k=0時所有點的連線。k0線上方泡齡k值為正；k0線下方泡齡k值為負。設定氣體沖出平口管管口開始形成氣泡的時刻為0時刻，圖中觀測時間為300 ms。在平均氣速為2.7 m·s-1時，曲線斷成兩部分，每段曲線呈拋物線分布，整個曲線存在兩個波峰。從98 ms開始，前一段曲線延伸到k0線以下，k為負值，氣泡已經脫離管口。分析原因可知，曲線延伸到k0線以下并與后一段曲線斷開，說明前一氣泡已脫離管口而后一氣泡還未形成。而出現兩次波峰表明在300 ms這個時間段，管口處生成了兩個氣泡。此外，從平均氣速為3.5 m·s-1的泡齡變化曲線可知，在同一測量時間，管口形成了5個完整氣泡，而第6個氣泡正處于膨脹階段，整條曲線沒有斷點。

圖9　不同平均氣速下泡齡k隨時間變化規律Fig.9　Change of age of bubble k with time under different average gas velocity

在脫離時間t大于58 ms后，平均氣速為3.5 m·s-1的泡齡變化曲線分別向k0線正方向和負方向延伸。泡齡變化曲線向k0線負方向延伸，k為負值，這表明氣泡已經脫離管口。同時泡齡變化曲線又向k0線上方延伸，k為正值，這說明在前一氣泡脫離管口時又生成了另一個氣泡，整個過程不存在間歇現象，如圖10（c）前4圖所示。平均氣速為2.7和3.5 m·s-1時管口處氣泡膨脹脫離整個演變過程如圖10（a）、（c）所示。圖10（a）展示了平均氣速為2.7 m·s-1時管口兩個氣泡生成過程。平均氣速為2.7 m·s-1時管口處氣泡生成存在間歇現象，如圖10（a）第4圖所示；圖10（c）給出了平均氣速為3.5 m·s-1時管口處氣泡連續生成的過程，氣泡生成不存在間歇現象。

圖10　20號管口頂端氣泡膨脹和脫離運動過程Fig.10　Bubble expansion and detachment movement of No.20 nozzle

此外如圖9所示，在平均氣速為3.1 m·s-1時，泡齡隨時間變化曲線與前兩條曲線的變化規律截然不同。曲線中k值先是出現一個較大的波峰，然后降至0以下，同時在k0線上方又有一個小波峰出現。究其原因是管口氣泡生成存在抽吸現象，如圖10（b）所示。前一氣泡在管口膨脹、脫離管口后又形成1個較小氣泡，但由于此時管內氣體壓力較小，氣泡外側受水壓等影響，該小氣泡不能繼續膨脹，又被抽吸進入管內，該過程如圖10（b）第3、4圖所示。但隨著氣體不斷供給，氣管內部壓力增大，使得氣泡又會在管口生成，如圖10（b）第5、6圖所示狀態。在該平均氣速下，管口氣泡的生成重復上述過程。其整個氣泡變化過程如圖10（b）所示。

另外分析圖9中3種平均氣速下泡齡-時間變化曲線可知，u=2.7 m·s-1時，最大泡齡kmax1值分布在8.1左右；u=3.1 m·s-1時，kmax2值在7.5左右；u=3.5 m·s-1時，kmax3值在8.5左右。在3.1 m·s-1平均氣速下（存在抽吸現象）的最大泡齡最小，這說明u=3.1 m·s-1時氣泡在管口體積膨脹位置較淺，抽吸現象不利于氣泡在管口體積膨脹的運動。從圖9和圖10綜合分析可知，單個氣泡膨脹及脫離時間分別為98、7和58 ms，這說明管內氣體平均氣速越大管口單個氣泡膨脹及脫離時間越短。

2.4平口管管徑對單個氣泡膨脹及脫離時間的影響

圖11　不同管徑單個氣泡膨脹及脫離時間變化規律Fig.11　Change of single bubble expansion and detachment time under different nozzle diameter

圖11給出了不同平口管管徑下單個氣泡膨脹及脫離時間隨平口管內平均氣速變化的規律。其中單個氣泡膨脹及脫離時間tg為氣泡在管口開始形成、膨脹再到氣泡最后脫離管口所經歷的總時間。從圖中可以看出，對于同一管徑，隨平均氣速的增加，單個氣泡膨脹及脫離時間先是減小，當平均氣速超過某一值后，tg開始增加。因此tg隨平均氣速u變化的過程存在最小值，并且tg的最小值隨管徑的增加而增大。圖11中各曲線單個氣泡膨脹及脫離時間最小值分別為35.1、32.6、28.9 ms。

2.5脫離階段氣泡圓度分布特性

圖12給出了平均氣速分別為1.1和2.3 m·s-1下25號管脫離階段橢球形氣泡長短軸比、等效球形氣泡脫離直徑dg隨時間的變化規律。本文利用橢球形氣泡長短軸比L來表征氣泡脫離階段各時刻形狀特征

式中，A、B分別為橢圓形氣泡的長軸和短軸。

圖12　脫離階段氣泡圓度變化分布（25號管）Fig.12　Bubble circular degree change distribution in detachment stage （No. 25 pipe）

圖12中dg（u=1.1 m·s-1）、L（u=1.1 m·s-1）和dg（u=2.3 m·s-1）、L（u=2.3 m·s-1）分別為平均氣速為1.1和2.3 m·s-1下等效球形氣泡脫離直徑和長短軸比。從圖中等效球形氣泡脫離直徑-時間變化曲線可以發現，平均氣速越小，氣泡脫離階段經歷的時間越長，兩種氣速下的脫離階段時間分別為69和56 ms。平均氣速較大時的管口氣泡脫離直徑也較大。分析橢球形氣泡長短軸比-脫離時間變化曲線可知，氣泡長短軸比L隨時間變化過程中存在時間分界點。時間分界點之前，L值在2.0附近波動，說明氣泡主要以橢圓形存在。在時間分界點之后到氣泡脫離階段結束時，L值明顯增加，氣泡向扁橢球體發展。圖12中兩種平均氣速下橢球形氣泡長短軸比L變化曲線的時間分界點分別為56和43 ms。

3　結　論

（1）頂端浸沒豎直向下平口管口氣泡形成演變過程可分成兩個階段：膨脹階段和脫離階段。本實驗條件下管口處氣泡脫離階段經歷的時間遠比膨脹階段經歷的時間長，如直徑為2.5 mm平口管在平均氣速為1.7 m·s-1時管口氣泡脫離階段經歷的時間為47.2 ms，膨脹階段經歷的時間為16.9 ms。

（2）同一平均氣速下氣泡脫離直徑隨管徑增加而逐漸增大。同一管徑下氣泡脫離直徑隨著平均氣速增加而增加的過程中明顯存在氣泡脫離形態轉變點，并且氣泡脫離形態轉變氣速隨管徑增加而減小。在氣泡脫離形態轉變點之前，氣泡脫離直徑呈近似線性增加，在氣泡脫離形態轉變點之后，氣泡脫離直徑平緩增加。同一管徑下氣泡膨脹脫離模式可分為單個氣泡形成和氣泡聚并形成兩種。

（3）氣泡生成機理有間歇生成氣泡、含抽吸間歇生成氣泡和連續生成氣泡3種。同一管徑下單個氣泡膨脹及脫離時間隨平均氣速增加而減小，對應的氣泡生成頻率增加。本實驗條件下直徑為2 mm的平口管在平均氣速為3.1 m·s-1時，管口處存在抽吸現象，最大泡齡kmax為7.5，小于其他平均氣速下的kmax值。

（4）脫離階段氣泡經歷的時間隨平均氣速增加而減小，氣泡脫離直徑隨平均氣速增加呈近似線性增加，且平均氣速越大增加值越大。本實驗條件下直徑為2.5 mm平口管兩種平均氣速下氣泡脫離階段經歷的時間分別為56和69 ms。脫離階段氣泡長短軸比-脫離時間變化曲線存在時間分界點，時間分界點前，橢球形氣泡長短軸比L值在2.0附近波動，氣泡主要以橢球形存在；時間分界點后，L值明顯增加，氣泡向扁橢球體發展。

符號說明

A ——橢球形氣泡實際長軸長，mm

a ——橢球形氣泡幀圖片中長軸長，mm

B ——橢球形氣泡實際短軸長，mm

b ——橢球形氣泡幀圖片中短軸長，mm

D ——平口細管實際直徑，mm

Di——平口細管實際內徑，mm

Do——平口細管實際外徑，mm

d ——平口細管幀圖片中直徑，mm

dg——等效球形氣泡脫離直徑，mm k ——泡齡

kmax——最大泡齡

k0——泡齡為0時所有點的連線

L——橢球形氣泡長短軸比

t——氣泡膨脹脫離時間，ms

tg——單個氣泡膨脹及脫離時間，ms

u——平口管內平均氣速，m·s-1

ug——氣泡脫離形態轉變氣速，m·s-1

下角標

g——氣相

i——內部

max——最大值

o——外部

0——所有值為0

References

［1］ 李佟茗， 賈紹義. 表面黏度對小氣泡聚并的影響［J］. 化工學報，1995， 46（5）： 532-538. LI D M， JIA S Y. Effects of surface viscosities on coalescence between two small bubbles ［J］. Journal of Chemical Industry and Engineering （China）， 1995， 46（5）： 532-538.

［2］ 王煥君， 費黎明， 陳定旺， 等. 三相外環流反應器的流動特性［J］.化學反應工程與工藝， 1992， 8（1）： 62-70. WANG H J， FEI L M， CHEN D W， et al. Flow Characteristics of three-phase external-airlift loop reactor ［J］. Chemical Reaction Engineering and Technology， 1992， 8（1）： 62-70.

［3］ 李佟茗. 小氣泡或小液滴之間的聚并［J］. 化工學報， 1994， 45（1）：38-44. LI D M. Coalescence between small bubbles or drops ［J］. Journal of Chemical Industry and Engineering （China）， 1994， 45（1）： 38-44.

［4］ 李新海， 陳新民， 莫鼎成， 等. 液-液界面上氣泡生長行為［J］. 中南礦冶學院學報， 1992， 23（3）： 359-363. LI X H， CHEN X M， MO D C， et al. The growth of gas bubbles on liquid-liquid interface ［J］. Journal of Central South University（Science and Technology）， 1992， 23（3）： 359-363.

［5］ LIOW J L. Quasi-equilibrium bubble formation during top-submerged gas injection ［J］. Chemical Engineering Science， 2000，55（20）： 4515-4524. DOI： 10.1016/S0009-2509（00）00095-6.

［6］ 劉柳. 垂直上升管中氣泡動力學特性實驗研究［D］. 長沙： 中南大學， 2013. LIU L. Experimental study of bubble dynamis characteristics in vertical upward pipe ［D］. Changsha： Central South University， 2013.

［7］ 徐炯， 王彤， 楊波， 等. 靜止水下氣泡運動特性的測試與分析［J］.水動力學研究與進展A輯， 2008， 23（6）： 709-714. XU J， WANG T， YANG B， et al. The measurement and analysis of motion behavior of bubbles in calm water ［J］. Chinese Journal of Hydrodynamics， 2008， 23（6）： 709-714.

［8］ MONTANTE G， HORN D， PAGLIANTI A. Gas-liquid flow and bubble size distribution in stirred tanks ［J］. Chemical Engineering Science， 2008， 63（8）： 2107-2118. DOI：10.1016/j.ces.2008.01.005.

［9］ ZAHEDI P， SALEH R， MORENO-ATANASIO R， et al. Influence of fluid properties on bubble formation， detachment， rising and collapse；Investigation using volume of fluid method ［J］. Korean Journal of Chemical Engineering， 2014， 31（8）： 1349-1361. DOI： 10.1007/ s11814-014-0063-x.

［10］ 宋慶唐， 金家琪， 吳桂英， 等. 二維鼓泡床內氣泡尺寸分布的實驗與CFD模擬［J］. 化工學報， 2008， 59（2）： 335-340. SONG Q T， JIN J Q， WU G Y， et al. Experiment and CFD simulation of bubble size distribution in 2D gas-liquid bubble column ［J］. Journal of Chemical Industry and Engineering （China）， 2008， 59（2）：335-340.

［11］ LIU C J， LIANG B， TANG S W， et al. Effects of orifice orientation and gas-liquid flow pattern on initial bubble size ［J］. Chinese Journal of Chemical Engineering， 2013， 21（11）： 1206-1215. DOI：10.1016/S1004-9541（13）60630-8.

［12］ 李秀軍. Savard-Lee雙重套管噴嘴在鋼鐵冶金和有色冶金行業應用發展綜述［J］. 中國有色冶金， 2015， （2）： 1-5， 29. DOI： 10.3969/ j.issn.1672-6103.2015.02.001. LI X J. The Savard-Lee shrouded injector： a review of its adoption and adaptation from ferrous to non-ferrous pyrometallurgy ［J］. China Nonferrous Metallurgy， 2015， （2）： 1-5， 29. DOI： 10.3969/j.issn. 1672-6103.2015.02.001.

［13］ HUDA N， NASER J， BROOKS G， et al. CFD Modeling of swirl and nonswirl gas injections into liquid baths using top submerged lances［J］. Metallurgical and Materials Transactions， 2010， 41（1）： 35-50. DOI： 10.1007/s11663-009-9316-1.

［14］ 詹樹華， 歐儉平， 蕭澤強. 底吹氣連鑄中間包內氣液兩相流的數值模擬［J］. 過程工程學報， 2005， 5（3）： 233-240. DOI： 10.3321/j.issn：1009-606X.2005.03.001. ZHAN S H， OU J P， XIAO Z Q. Numerical simulation on gas-liquid two phase flow in continuous casting tundish with bottom gas blowing ［J］. The Chinese Journal of Process Engineering， 2005， 5（3）：233-240. DOI： 10.3321/j.issn：1009-606X.2005.03.001.

［15］ 王輔臣， 于廣鎖， 龔欣， 等. 大型煤氣化技術的研究與發展［J］. 化工進展， 2009， 28（2）： 173-180. DOI： 10.16085/j.issn.1000-6613. 2009.02.017. WANG F C， YU G S， GONG X， et al. Research and development of large-scale coal gasification technology ［J］. Chemical Industry and Engineering Progress， 2009， 28（2）： 173-180. DOI： 10.16085/j.issn. 1000-6613.2009.02.017.

［16］ 吳晅， 高慧杰， 王麗芳， 等. 氣化爐內氣體穿越液池過程中氣泡特性的數值模擬［J］. 動力工程學報， 2013， 33（12）： 962-968. WU X， GAO H J， WANG L F， et al. Bubble characteristics of syngas during passing through cistern in gasifier ［J］. Journal of Chinese Society of Power Engineering， 2013， 33（12）： 962-968.

［17］ 王沖， 熊靚， 朱道飛， 等. 貧化電爐氣液混合頂吹的氣泡群形態及攪拌效果［J］. 冶金能源， 2014， 33（2）： 30-33. DOI： 10.3969/j.issn. 1001-1617.2014.02.007. WANG C， XIONG L， ZHU D F， et al. Form of bubble group and stirring effect of top blowing with depleted furnace gas-liquid mixture［J］. Energy for Metallurgical Industry， 2014， 33（2）： 30-33. DOI：10.3969/j.issn.1001-1617.2014.02.007.

［18］ RIVERA-SALINAS J E， GUTIéRREZ-PéREZ V H， VARGASRAMíREZ M， et al. Application of computational fluid dynamic in aluminum refining through pneumatic injection of powders ［J］. Materials Research， 2014， 17（6）： 1550-1562. DOI： 10.1590/1516-1439. 275414.

［19］ 楊濮亦， 王沖， 王仕博， 等. 頂吹氣泡在兩相間運動的形變過程對熔池攪拌效果的影響［J］. 化工進展， 2014， 33（3）： 617-622. DOI：10.3969/j.issn.1000-6613.2014.03.016. YANG P Y， WANG C， WANG S B， et al. Experimental study on the influence for stirring effect of the bubbles deformation through two phases in top blowing bath ［J］. Chemical Industry and Engineering Progress， 2014， 33（3）： 617-622. DOI： 10.3969/j.issn. 1000-6613. 2014.03.016.

［20］ 吳晅， 李鐵， 袁竹林， 等. 靜止液體中頂部浸沒管口處氣泡形成的數值預測［J］. 東南大學學報（自然科學版）， 2008， 38（1）： 86-91. DOI：10.3969/j.issn.1001-0505.2008.01.017. WU X， LI T， YUAN Z L， et al. Numerical prediction of bubble formation from top-submerged nozzle in stagnant liquid ［J］. Journal of Southeast University （Natural Science Edition）， 2008， 38（1）： 86-91. DOI： 10.3969/j.issn.1001-0505.2008.01.017.

Expansion and detachment characteristics of bubble at downward nozzle of vertical flat pipe

WU Xuan1，2， JIAO Jingjing1， LIANG Panlong1， JIN Guang1， WU Wenfei2
（1School of Energy and Environment， Inner Mongolia University of Science and Technology， Baotou 014010， Inner Mongolia， China；2Key Laboratory of Integrated Exploitation of Bayan Obo Multi-Metal Resources， Inner Mongolia University of Science and Technology， Baotou 014010， Inner Mongolia， China）

Gas bubbles in liquid influence efficiency of mass and heat transfer. Much efforts on the observation of expansion and detachment of bubble at upward nozzle have been conducted， while that of downward nozzle is rarely known. In this article， the expansion and detachment of bubble at the downward nozzle of vertical flat pipe was observed by a high speed camera. Effect of the nozzle diameter and average gas velocity on bubble size and time of formation was investigated. The bubble age， k and the diameter ratio of long axis to short axis， L describe bubble position and bubble shape， respectively. It showed that three kinds of mechanism of intermittent formation，intermittent formation with suction and continued formation exist in bubble formation. With the variation of detachment diameter of bubble as a function of average gas velocity， a critical point representing the evolution of bubble detachment exists. The expansion and detachment of bubble along with average gas velocity， follow two types of single bubble formation and coalescence formation. The suction phenomenon occurred at 3.1 m·s-1of average gas velocity， of which kmaxis minimized. In the detachment stage， a critical point of detachment time exists in the evolution of L as a function of time， namely， L allocates around 2.0 with mostly ellipsoidal bubblebefore the critical point， however， L increases significantly with the formation of flat ellipsoidal bubble after the critical point.

date： 2015-09-21.

Prof. JIN Guang， wlzhebei@126.com

supported by the National Natural Science Foundation of China （51166010） and the Natural Science Foundation of Inner Mongolia（2014MS0530）.

gas-liquid flow； downward nozzle； bubble； hydrodynamics

TQ 027.3

A

0438—1157（2016）05—1868—10

2015-09-21收到初稿，2015-12-22收到修改稿。

聯系人：金光。第一作者：吳晅（1976—），男，博士，副教授。

國家自然科學基金項目（51166010）；內蒙古自治區自然科學基金項目（2014MS0530）。