馬 超，薄涵亮

(1.中國(guó)核動(dòng)力研究設(shè)計(jì)院 核反應(yīng)堆系統(tǒng)設(shè)計(jì)技術(shù)重點(diǎn)實(shí)驗(yàn)室，四川 成都 610041；2.清華大學(xué) 核能與新能源技術(shù)研究院，北京 100084)

單氣泡破裂產(chǎn)生膜液滴空間分布實(shí)驗(yàn)研究

馬 超1，薄涵亮2

(1.中國(guó)核動(dòng)力研究設(shè)計(jì)院 核反應(yīng)堆系統(tǒng)設(shè)計(jì)技術(shù)重點(diǎn)實(shí)驗(yàn)室，四川 成都 610041；2.清華大學(xué) 核能與新能源技術(shù)研究院，北京 100084)

本文研究了常溫、常壓條件下自由液面單氣泡破裂產(chǎn)生的膜液滴現(xiàn)象。利用染色液滴撞擊熟宣紙表面不洇散特性，加工宣紙紙筒進(jìn)行稀釋墨水中氣泡破裂產(chǎn)生膜液滴空間分布實(shí)驗(yàn)，并通過后續(xù)圖片掃描及編寫Matlab圖像處理程序獲得膜液滴在圓柱紙筒上的尺寸分布、空間分布信息。通過改變實(shí)驗(yàn)條件，討論分析氣泡尺寸及紙筒半徑對(duì)膜液滴空間分布結(jié)果的影響，獲得了常溫、常壓條件下自由液面單氣泡破裂產(chǎn)生膜液滴空間分布規(guī)律。

蒸汽發(fā)生器；汽水分離；單氣泡破裂；膜液滴；空間分布實(shí)驗(yàn)

在壓水堆蒸汽發(fā)生器二次側(cè)，由于氣泡破裂或機(jī)械打碎等原因產(chǎn)生大量細(xì)小的液滴，這些液滴被夾帶在蒸汽發(fā)生器二次側(cè)產(chǎn)生的飽和蒸汽中，如果這些小液滴進(jìn)入反應(yīng)堆二回路系統(tǒng)會(huì)造成極大的破壞，因此在蒸汽發(fā)生器上部設(shè)有汽水分離裝置，汽水分離裝置的作用至關(guān)重要。確定液滴來源信息是研究汽水分離機(jī)理工作的基礎(chǔ)。蒸汽發(fā)生器內(nèi)液滴按照其產(chǎn)生機(jī)理主要可劃分為兩種：氣泡破裂產(chǎn)生的液滴及機(jī)械打碎產(chǎn)生的液滴。機(jī)械打碎包括氣流撕破液膜產(chǎn)生卷吸液滴、原液滴聚并或斷裂產(chǎn)生新液滴、液滴撞擊壁面或液膜產(chǎn)生飛濺液滴。氣泡破裂產(chǎn)生液滴按照方式不同，又可分為兩種液滴：噴射液滴與膜液滴[1]。由于氣泡破裂產(chǎn)生液滴是蒸汽發(fā)生器內(nèi)液滴產(chǎn)生的主要方式之一[2]，研究氣泡破裂產(chǎn)生液滴的現(xiàn)象及機(jī)理，對(duì)汽水分離裝置研究設(shè)計(jì)具有重要意義。

由于氣泡破裂產(chǎn)生膜液滴現(xiàn)象幾何尺度小，時(shí)間間隔短，物理機(jī)理較為復(fù)雜，針對(duì)該現(xiàn)象的實(shí)驗(yàn)研究始終占據(jù)主導(dǎo)地位。對(duì)研究汽水分離機(jī)理而言，初始速度與位置是膜液滴在流場(chǎng)中運(yùn)動(dòng)的重要初始參數(shù)。前人實(shí)驗(yàn)主要通過建立高速攝像裝置研究該物理現(xiàn)象。由于遇到的實(shí)際工程問題不同及受限于高速攝像的二維限值，針對(duì)氣泡破裂產(chǎn)生膜液滴現(xiàn)象的實(shí)驗(yàn)研究主要關(guān)注膜液滴的產(chǎn)生數(shù)量、尺寸，對(duì)膜液滴的初始速度、初始位置等信息關(guān)注度較低[3-7]。目前雖不能通過直接的方法測(cè)量氣泡破裂產(chǎn)生膜液滴的初始速度及位置，但可通過間接的方法在靠近氣泡附近收集膜液滴產(chǎn)生并飛行較短距離后的空間分布信息。

在MgO涂層法基礎(chǔ)上[8-9]，本文利用熟宣紙加工的紙筒測(cè)量膜液滴的空間分布。

1 實(shí)驗(yàn)簡(jiǎn)介

氣泡破裂產(chǎn)生膜液滴空間分布的實(shí)驗(yàn)裝置示意圖如圖1所示。應(yīng)用氣泡發(fā)生裝置，包括空氣泵、針閥、質(zhì)量流量計(jì)、噴嘴等，在水槽底部形成氣泡。由前人理論分析可知，膜液滴產(chǎn)生現(xiàn)象主要同氣泡露出的液帽部分有關(guān)，因此水槽中液面高度設(shè)置為僅高于噴嘴5 mm，使得氣泡在產(chǎn)生、長(zhǎng)大、破裂過程中均不脫離噴嘴，以保證氣泡對(duì)中性。將實(shí)驗(yàn)臺(tái)架設(shè)置為對(duì)中支架，夾持制備好的宣紙紙筒。調(diào)節(jié)針閥，控制穩(wěn)定的氣體流量，產(chǎn)生的氣泡在液面處破裂產(chǎn)生膜液滴，紙筒測(cè)量記錄膜液滴空間分布的信息。實(shí)驗(yàn)中測(cè)量氣泡破裂產(chǎn)生膜液滴在紙筒上的空間分布。紙筒測(cè)量結(jié)束后，從夾持裝置中取下展平，利用掃描儀器將紙面信息掃描成圖片，使用Matlab圖像處理程序進(jìn)行進(jìn)一步分析處理，通過運(yùn)用邊緣檢測(cè)、連通域分析、特征識(shí)別等函數(shù)，得到掃描圖片中的膜液滴數(shù)量、直徑、空間位置等信息。

圖1 氣泡破裂產(chǎn)生膜液滴空間分布的實(shí)驗(yàn)裝置示意圖Fig.1 Schematic of space distribution experiment device of film drops produced by bubble bursting

染色液滴選擇稀釋的水性墨水，墨與水體積比為5∶100，該比例混合既可保證染色混合流體密度接近于純水，且紙面著墨視覺效果良好以便后續(xù)實(shí)驗(yàn)數(shù)據(jù)掃描處理。通過實(shí)驗(yàn)測(cè)量，得到該混合染色液體的物性參數(shù)為：密度ρ=960 kg/m3，黏性系數(shù)η=4.68×10-3Pa·s，表面張力系數(shù)σ=54.5 mN/m，染色液滴同宣紙之間靜態(tài)接觸角θ=131°。墨水中炭黑顆粒直徑在20～40 nm之間[10]，實(shí)驗(yàn)中氣泡曲率半徑大于20 mm，液膜厚度為10 μm以上，因此認(rèn)為墨水中微小炭黑顆粒對(duì)膜液滴產(chǎn)生方式無影響。

實(shí)驗(yàn)中紙筒夾持支架加工精度為9級(jí)，對(duì)中誤差范圍為±0.38 mm。該實(shí)驗(yàn)方法測(cè)量精度主要受后處理過程中圖像掃描設(shè)備分辨率限制。實(shí)驗(yàn)中掃描儀分辨率為400像素/英寸(1英寸=2.54 cm)，圖片中像素尺寸為64 μm，根據(jù)膜液滴撞擊壁面形成最大鋪展直徑關(guān)系式[11]，由圖片顯示尺寸換算成膜液滴原始尺寸。該實(shí)驗(yàn)方法在現(xiàn)有掃描設(shè)備精度條件下能夠測(cè)量的最小膜液滴直徑為34 μm。

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

根據(jù)Toba[12]的氣泡浮動(dòng)與破裂理論和實(shí)驗(yàn)結(jié)果，浮動(dòng)氣泡液帽為軸對(duì)稱結(jié)構(gòu)，因此推測(cè)破裂點(diǎn)周向方位角位置符合等概率均勻分布。此處應(yīng)用同一紙筒(紙筒半徑Rp=60 mm)收集相同大小氣泡(氣泡曲率半徑R=28.4 mm)多次破裂產(chǎn)生膜液滴的周向角度分布情況，如圖2所示。由圖2可見，隨著氣泡數(shù)量的增加，紙筒表面周向分布直方圖逐漸趨于平均，證明破裂點(diǎn)周向方位角位置等概率分布。

雖然大量氣泡破裂情況下周向角度符合軸對(duì)稱規(guī)律，但單個(gè)氣泡破裂過程中破裂點(diǎn)產(chǎn)生的隨機(jī)性較大，使得在紙筒上膜液滴周向分布峰值角度均會(huì)不同。為便于對(duì)比分析單氣泡破裂產(chǎn)生膜液滴空間分布統(tǒng)計(jì)規(guī)律，根據(jù)氣泡液帽軸對(duì)稱性，假設(shè)實(shí)驗(yàn)中氣泡破裂點(diǎn)周向方位角固定，對(duì)紙筒表面膜液滴的周向角度根據(jù)360°周期進(jìn)行對(duì)中處理，將分布曲線峰值(或直方圖峰值)調(diào)整至150°～180°區(qū)間。

雖然氣泡液帽部分是單點(diǎn)破裂，但由于破裂點(diǎn)軸向相對(duì)位置不固定，每次破裂產(chǎn)生膜液滴的空間分布都是變化的，因此只有通過統(tǒng)計(jì)同一尺寸氣泡多次破裂結(jié)果來考察膜液滴的空間分布規(guī)律。實(shí)驗(yàn)中變化工況，考察氣泡大小及紙筒半徑對(duì)膜液滴分布的影響。每個(gè)工況相應(yīng)地記錄20次氣泡重復(fù)破裂情況，同一曲率半徑(R=22.9 mm)下氣泡多次破裂產(chǎn)生膜液滴在不同半徑(分別為30、40、50、60 mm)紙筒上的空間分布譜如圖3所示。

圖2 多氣泡破裂產(chǎn)生膜液滴周向分布直方圖Fig.2 Angle distribution histogram of film drops produced by multi-bubble bursting

紙筒半徑：a——30 mm；b——40 mm；c——50 mm；d——60 mm

將該組實(shí)驗(yàn)數(shù)據(jù)繪制成膜液滴尺寸分布方框圖，如圖4所示。由圖4可見，隨著紙筒半徑的增大，膜液滴平均直徑逐漸變小，尺寸分布范圍逐漸縮小，方框圖形狀上寬下窄，平均值高于中間值，表明紙筒上撞擊的小尺寸膜液滴份額較多，大尺寸膜液滴份額較少。將實(shí)驗(yàn)數(shù)據(jù)畫成周向角度分布曲線，如圖5所示，分區(qū)統(tǒng)計(jì)間隔為30°。由圖5可見，各曲線變化趨勢(shì)均是中間角度區(qū)域高、邊緣角度區(qū)域低。隨著紙筒半徑的增大，周向分布曲線峰值增大，曲線趨勢(shì)由平緩逐漸變得陡峭，表明膜液滴由周向均勻分布逐漸向中間方位角區(qū)域集中。將實(shí)驗(yàn)數(shù)據(jù)畫成軸向高度分布曲線，如圖6所示，分區(qū)統(tǒng)計(jì)間隔為10 mm。由圖6可見，各曲線變化趨勢(shì)均是紙筒底部區(qū)域高、頂部區(qū)域低，隨著紙筒半徑的增大，軸向膜液滴的最大高度逐漸變大，軸向分布曲線峰值降低，曲線趨勢(shì)由陡峭逐漸變得平緩，表明膜液滴由集中在紙筒底部區(qū)域逐漸向紙筒頂部區(qū)域擴(kuò)散。

圖4 紙筒半徑對(duì)膜液滴尺寸分布的影響Fig.4 Effect of paper cylinder radius on size distribution of film drops

圖5 紙筒半徑對(duì)膜液滴周向角度分布的影響Fig.5 Effect of paper cylinder radius on angle distribution of film drops

圖6 紙筒半徑對(duì)膜液滴軸向高度分布的影響Fig.6 Effect of paper cylinder radius on height distribution of film drops

結(jié)合高速攝像實(shí)驗(yàn)[3,13]分析以上尺寸、空間分布的結(jié)果。由于膜液滴飛行速度及方向由產(chǎn)生時(shí)刻射流所處球面方位角決定：中間角度膜液滴產(chǎn)生初始時(shí)刻沿液帽切向向上飛行，隨著紙筒半徑增大，撞擊高度逐漸增大；邊緣角度膜液滴產(chǎn)生初始時(shí)刻沿著液帽切向向下飛行，隨著紙筒半徑增加，撞擊高度逐漸降低，最終撞擊水面消失。在真實(shí)情況下，由于液膜厚度受到重力影響，氣泡球面液膜中間薄邊緣厚，邊緣方位角產(chǎn)生的膜液滴的直徑大于中間方位角產(chǎn)生的膜液滴的直徑，飛行過程中較大直徑液滴在重力的作用下在到達(dá)紙筒之前撞擊水面。綜上，隨著紙筒半徑的增大，逐漸形成膜液滴平均直徑變小、周向角度分布變窄、軸向高度分布變寬的空間圖譜。

同一半徑(Rp=50 mm)紙筒上不同大小氣泡(曲率半徑分別為21.5、22.9、24.6、28.4 mm)下多次破裂產(chǎn)生膜液滴的空間分布譜如圖7所示。

將該組實(shí)驗(yàn)數(shù)據(jù)繪制成膜液滴尺寸分布方框圖，如圖8所示。由圖8可見，方框圖形狀上寬下窄，平均值高于中間值，表明紙筒上撞擊的小尺寸膜液滴份額較多，大尺寸膜液滴份額較少。隨著氣泡曲率半徑的增加，膜液滴平均直徑逐漸增大，膜液滴尺寸分布范圍逐漸變大。將實(shí)驗(yàn)數(shù)據(jù)繪制成周向角度分布曲線，如圖9所示，分區(qū)統(tǒng)計(jì)間隔為30°。由圖9可見，各曲線變化趨勢(shì)均是中間角度區(qū)域高、邊緣角度區(qū)域低，隨著氣泡曲率半徑的增加，周向分布曲線峰值降低，曲線趨勢(shì)由陡峭逐漸變得平緩，表明膜液滴由集中在中間方位角區(qū)域逐漸向邊緣方位角擴(kuò)散。將實(shí)驗(yàn)數(shù)據(jù)繪制成軸向高度分布曲線，如圖10所示，分區(qū)統(tǒng)計(jì)間隔為10 mm。由圖10可見，各曲線變化趨勢(shì)均是紙筒底部區(qū)域高、頂部區(qū)域低，隨著氣泡曲率半徑的增加，軸向上膜液滴最大高度逐漸變大，軸向分布曲線峰值降低，曲線趨勢(shì)由陡峭逐漸變得平緩，表明膜液滴由集中在紙筒底部區(qū)域逐漸向頂部區(qū)域擴(kuò)散。

曲率半徑：a——21.5 mm；b——22.9 mm；c——24.6 mm；d——28.4 mm

圖8 氣泡曲率半徑對(duì)膜液滴尺寸分布的影響Fig.8 Effect of different bubble radiion size distribution of film drops

結(jié)合此前高速攝像實(shí)驗(yàn)結(jié)果，分析以上尺寸、空間分布結(jié)果。膜液滴尺寸主要同氣泡大小有關(guān)，隨氣泡曲率半徑增大，膜液滴平均直徑逐漸增大，膜液滴尺寸分布范圍逐漸變寬。與此同時(shí)，液帽所對(duì)應(yīng)的球心仰角逐漸增大，不穩(wěn)定射流頂端速度仰角逐漸增大，使得膜液滴初始最大速度仰角(膜液滴速度同水面交角)增大，因此在同一半徑紙筒表面撞擊高度逐漸增大，周向角度、軸向高度分布逐漸變寬。綜上，隨著氣泡曲率半徑增大，逐漸形成膜液滴平均直徑變大，周向角度分布變寬，軸向高度分布變寬的空間譜。

圖9 氣泡曲率半徑對(duì)膜液滴周向角度分布的影響Fig.9 Effect of bubble radius on angle distribution of film drops

圖10 氣泡曲率半徑對(duì)膜液滴軸向高度分布的影響Fig.10 Effect of bubble radius on height distribution of film drops

3 結(jié)束語

本文利用實(shí)驗(yàn)方法研究了紙筒半徑及氣泡大小對(duì)膜液滴空間分布的影響，得到以下結(jié)論。隨紙筒半徑的增大，膜液滴尺寸分布范圍逐漸減小，膜液滴平均直徑減小，中間角度區(qū)域膜液滴份額提高，邊緣角度區(qū)域膜液滴份額降低，紙筒底部區(qū)域膜液滴份額降低，頂部區(qū)域膜液滴份額提高。隨氣泡曲率半徑的增大，膜液滴尺寸分布范圍逐漸增大，膜液滴平均直徑增大，中間角度區(qū)域膜液滴份額降低，邊緣角度區(qū)域膜液滴份額提高，紙筒底部區(qū)域膜液滴份額降低，頂部區(qū)域膜液滴份額提高。

[1] 張謹(jǐn)奕. 三維流場(chǎng)中單液滴運(yùn)動(dòng)模型和應(yīng)用研究[D]. 北京：清華大學(xué)核能與新能源技術(shù)研究院，2012.

[2] 李宋恒. 電站鍋爐汽水分離裝置的原理和設(shè)計(jì)[M]. 北京：水利水電出版社，1986.

[3] LHUISSIER H, VILLERMAUX E. Bursting bubble aerosols[J]. Journal of Fluid Mechanics, 2012, 696: 5-44.

[4] BLANCHARD D C, SYZDEK L D. Film drop production as a function of bubble-size[J]. Jour-nal of Geophysical Research-Oceans, 1988, 93(C4): 3 649-3 654.

[5] RESCH F, AFETI G. Film drop distributions from bubbles bursting in seawater[J]. Journal of Geophysical Research-Oceans, 1991, 96(C6): 10 681-10 688.

[6] RESCH F, AFETI G. Submicron film drop production by bubbles in seawater[J]. Journal of Geophysical Research-Oceans, 1992, 97(C3): 3 679-3 683.

[7] RESCH F, DARROZES J, AFETI G. Marine liquid aerosol production from bursting of air bubbles[J]. Journal of Geophysical Research-Oceans, 1986, 91(C1): 1 019-1 029.

[8] SPIEL D E. A hypothesis concerning the peak in film drop production as a function of bubble size[J]. Journal of Geophysical Research-Oceans, 1997, 102(C1): 1 153-1 161.

[9] SPIEL D E. On the births of film drops from bubbles bursting on seawater surfaces[J]. Journal of Geophysical Research-Oceans, 1998, 103(C11): 24 907-24 918.

[10]張君紅，羅春榮，趙曉鵬. 用于電子墨水的納米炭黑顆粒表面電荷改性[J]. 功能材料,2011,42(9):1 559-1 562,1 568.

ZHANG Junhong，LUO Chunrong，ZHAO Xiao-peng.Modification of nano carbon black particles to obtain surface charge for electronic ink[J]. Journal of Functional Materials, 2011, 42(9): 1 559-1 562, 1 568(in Chinese).

[11]MAO T, KUHN D C S, TRAN H. Spread and rebound of liquid droplets upon impact on flat surfaces[J]. Aiche Journal, 1997, 43(9): 2 169-2 179.

[12]TOBA Y. Drop production by bursting of air bubbles on the sea surface, ii: Theoretical study on the shape of floating bubbles[J]. OceanogrSocJapan, 1959, 15: 121-130.

[13]馬超. 自由液面單氣泡破裂產(chǎn)生膜液滴現(xiàn)象實(shí)驗(yàn)與理論研究[D]. 北京：清華大學(xué)核能與新能源技術(shù)研究院，2014.

Experimental Research on Space-distribution of Film Drops Produced by Single Bubble Bursting

MA Chao1, BO Han-liang2

(1.ScienceandTechnologyonReactorSystemDesignTechnologyLaboratory,NuclearPowerInstituteofChina,Chengdu610041,China;2.InstituteofNuclearandNewEnergyTechnology,TsinghuaUniversity,Beijing100084,China)

The phenomenon of the film drops produced by the single bubble bursting at a free water surface was researched. The rice paper cylinder, which was made under special technique process and did not soak the ink at the surface, was used to conduct the space distribution experiment of the film drops produced by bubble bursting in the ink water. The size and space distributions of film drops on the different radius cylinders were obtained by scanning the picture and using the image processing code programmed in Matlab. With experiment conditions changing, the impact of the bubble size and the radius of the rice paper cylinder on space distributions of film drops was discussed. The rule of the space distribution of film drops produced by single bubble bursting at a free water surface was obtained under the normal temperature and pressure condition.

steam generator; steam water separation; single bubble bursting; film drops; space distribution experiment

2014-06-27;

2015-02-01

馬 超(1987—)，男，遼寧鐵嶺人，工程師，博士，核科學(xué)與技術(shù)專業(yè)

TL333

A

1000-6931(2015)10-1766-06

10.7538/yzk.2015.49.10.1766