劉國強，孫立成,2，閻昌琪，田道貴

(1.哈爾濱工程大學 核安全與仿真技術國防重點學科實驗室，黑龍江 哈爾濱 150001;2.四川大學 水利水電學院，四川 成都 610065)

管內(nèi)向下兩相流動現(xiàn)象在核電、化工生產(chǎn)、新型換熱設備的設計及石油開采等工業(yè)領域都會遇到[1-3]。例如：在新一代輕水反應堆核電站的非能動余熱排出系統(tǒng)中，常用的一種設計是利用自然循環(huán)將堆芯或蒸汽發(fā)生器產(chǎn)生的蒸汽引入到豎直布置于高位水箱中的冷凝換熱器內(nèi)，蒸汽在豎直布置的冷凝管內(nèi)不斷冷凝，從而將堆芯熱量導出，此時冷凝管內(nèi)為向下兩相流；在石油開采中，為了保持油藏的壓力，使油井能靠自身的壓力產(chǎn)油，通常向油藏中注入氣水混合物，這是提高石油采收率的一個重要手段；另外，在一些類型的鍋爐中，也會出現(xiàn)向下氣液兩相流動。

兩相流動中界面結(jié)構(gòu)對相間的界面輸運過程至關重要，而截面含氣率和界面面積濃度是表征界面特征最重要的兩個參數(shù)，在兩相流熱工水力特性分析中至關重要。因此，獲得空泡份額、界面面積濃度及氣泡頻率等局部參數(shù)的分布特性，是進一步研究向下兩相流動中界面輸運過程的必要條件。

目前，國內(nèi)外學者對豎直向上兩相流動的研究已開展多年，取得了很多研究成果[4-6]，但對于向下兩相流動的研究還不多，特別是對界面輸運過程的研究還很不充分。本文借助自研的雙探頭光纖探針及光學測量系統(tǒng)，對管內(nèi)向下流動泡狀流界面參數(shù)分布特性進行研究，旨為向下兩相流動的界面輸運過程研究奠定基礎。

1 實驗裝置

1.1 實驗回路

實驗回路主要由水回路和氣回路組成，如圖1所示。實驗段選用內(nèi)徑為50 mm的透明有機玻璃圓管，水通過離心泵的驅(qū)動進入混合腔，而空氣通過空壓機加壓貯存于高壓氣罐內(nèi)，并由高壓儲氣罐供給。兩者相繼進入混合腔，水和空氣在混合腔內(nèi)進行充分混合，然后進入豎直的可視實驗段。通過實驗段后，空氣、水經(jīng)回水管路匯入水箱，水在水箱內(nèi)循環(huán)使用，而空氣則經(jīng)自然分離釋放到大氣中。通過在回水管路設置調(diào)節(jié)閥使實驗段內(nèi)壓力為正壓，并借助排水閥排除由混合腔進入氣回路內(nèi)的液體。

1.2 氣液混合腔

實驗采用文丘里管結(jié)構(gòu)的氣液混合腔，氣液混合腔由套管及文丘里管組成(圖2)。借助文丘里管突縮突擴的獨特結(jié)構(gòu)，使液相流經(jīng)混合腔時在喉管內(nèi)形成負壓，從而將空氣通過文丘里管與套管之間的小孔吸入喉管，并在漸擴區(qū)域被強烈的紊流作用打碎成細小氣泡，形成均勻的兩相流動。

圖1 實驗回路簡圖

圖2 混合腔剖面圖

水流量和氣流量全部由質(zhì)量流量計測量，測量精度分別為0.2級和1級。在實驗段L/D(L為實驗段長度，D為圓管內(nèi)徑)為24及34處安裝有精度為0.1級的壓力傳感器。探針距實驗段入口1.2 m，沿半徑方向從中心線向壁面設置11個測點，與通道中心的相對距離ri/R(i=0，…，11；R為通道半徑)分別為0、0.16、0.32、0.44、0.56、0.68、0.76、0.84、0.88、0.92和0.96。液相表觀速度范圍為0.43～0.71 m/s，氣相表觀速度范圍為0.004～0.077 m/s。

實驗中，利用探針獲取的光信號，經(jīng)放大轉(zhuǎn)化系統(tǒng)處理后，變成電壓數(shù)字信號，再由NI數(shù)據(jù)采集系統(tǒng)進行采集，并存于計算機中。采集時間和采集頻率分別為120 s和10 kHz，以保證得到足夠的數(shù)據(jù)量，確保數(shù)據(jù)統(tǒng)計結(jié)果的準確。

2 實驗數(shù)據(jù)處理

實驗采用自研雙探頭光纖探針對界面面積濃度、空泡份額、氣泡頻率等局部界面參數(shù)進行測量。文獻[7]詳細介紹了雙探頭光纖探針的結(jié)構(gòu)及其對氣液兩相局部界面參數(shù)的測量原理。文獻[8]分別采用高速攝影儀和光纖探針兩種方法獲得空泡份額，并對兩種方法進行了比較，兩者的相對偏差在10%以內(nèi)，這進一步說明應用光纖探針對空泡份額進行測量是可行的。本實驗采用的光纖探針及數(shù)據(jù)處理方法與文獻[8]相同。

局部界面面積濃度是指單位體積內(nèi)兩相界面的面積大小，采用Kataoka等提出的模型進行計算：

tan 0.5α0ln(sin 0.5α0))

(1)

(2)

對于標準偏差σz，首先通過多普勒技術得到viz(界面速度在z軸方向上的分量)的光譜，再由σz與viz的關系進行求解，由下式計算：

(3)

局部氣泡頻率指單位時間內(nèi)通過局部測點的氣泡數(shù)目，可通過統(tǒng)計單位時間內(nèi)前端探頭信號出現(xiàn)高電勢的次數(shù)得到。

3 實驗結(jié)果分析

3.1 局部界面參數(shù)徑向分布

豎直向下兩相流動中空泡份額、界面面積濃度、氣泡頻率3個局部界面參數(shù)的徑向分布如圖3所示。可發(fā)現(xiàn)，各局部界面參數(shù)徑向分布特點相似。根據(jù)目前國內(nèi)外的研究結(jié)果，向上兩相流動中局部界面參數(shù)徑向分布主要呈“壁峰”或“核峰”兩種類型分布，而向下流動局部界面參數(shù)徑向分布呈“喇叭”型、“偏離核峰”或“核峰”型分布[9-13]。本實驗條件下，由于液相速度較低，各局部界面參數(shù)出現(xiàn)“壁峰”及“寬峰”型分布。“寬峰”型分布與前文提到的“喇叭”型、“偏離核峰”型類似。由圖3可看出：液相表觀速度jf=0.43 m/s且氣流量較高時，局部界面參數(shù)在圓管中心區(qū)域較均勻，在近壁區(qū)呈現(xiàn)峰值，即呈“壁峰”型分布；氣流量較低且00.8時，局部界面參數(shù)逐漸減小到一定值并出現(xiàn)一拐點，拐點后呈平滑分布，即為“寬峰”型分布。當0.57 m/s≤jf≤0.71 m/s，氣流量較低時，局部界面參數(shù)徑向分布均呈“寬峰”型分布；氣流量較高時，隨著水流量的增加，局部界面參數(shù)逐漸由“壁峰”轉(zhuǎn)變?yōu)椤皩挿濉毙头植迹欢S著水流量的不斷加大，空泡份額、界面面積濃度不斷減少，同時沿徑向分布變得更為均勻，這是由于湍流擴張力隨著水流速度的增加而增大，使氣泡的徑向分布更加均勻。隨著水流量的不斷增加，氣泡頻率逐漸增大，這主要是氣泡尺寸隨著液流量的增加而變小所引起的。

圖3 局部界面參數(shù)徑向分布

3.2 豎直向上、向下兩相流動局部參數(shù)對比

液相表觀速度為0.57 m/s時，豎直向上、向下兩種條件下局部界面面積濃度和空泡份額沿徑向分布特性示于圖4?？砂l(fā)現(xiàn)，豎直向上流動時，兩者呈現(xiàn)“壁峰”型分布，而豎直向下流動時則呈“寬峰”型分布，且在主流區(qū)域參數(shù)值均比向上流動時的大，僅在壁面附近參數(shù)值比向上時的?。幌蛳铝鲃訒r空泡份額截面平均值比向上流動時大119.6%～145.0%，而界面面積濃度截面平均值比向上流動時大18.8%～82.5%。這主要是由于在氣液兩相向下流動中，浮力與氣泡運動方向相反，降低了氣泡的運動速度，導致通道內(nèi)氣泡數(shù)量增加，從而使界面面積濃度和空泡份額均隨之增加。

3.3 局部界面參數(shù)徑向分布形成原因

通過以上分析可發(fā)現(xiàn)：本實驗條件下，在豎直向下兩相流動中，局部界面參數(shù)徑向分布主要呈“寬峰”型或“壁峰”型分布。根據(jù)國內(nèi)外對局部界面參數(shù)徑向分布的研究結(jié)果得知，氣泡徑向受力是決定局部界面參數(shù)徑向分布的主要原因。其中，Tomiyama提出的橫向升力模型[14]為：

圖4 局部界面參數(shù)徑向分布

(4)

式中：CT、ρL、db、VG、VL分別為升力系數(shù)、液體密度、氣泡直徑、氣相軸向絕對速度及液相軸向絕對速度；rotVL為液相速度旋度。

根據(jù)式(4)可知 ，升力方向與升力系數(shù)、相對速度有關，而升力系數(shù)主要與氣泡尺寸有關。當兩者均為正值時，升力方向指向管道壁面，從而驅(qū)使氣泡向壁面運動；當兩參數(shù)為一正一負時，升力方向指向圓管中心，驅(qū)使氣泡向管道中心區(qū)域運動。本實驗條件下，氣泡尺寸均小于5 mm，故升力系數(shù)為正值。氣相速度由于受到浮力的抑制作用使其在圓管中心區(qū)域明顯低于對應的液相速度；在壁面附近，不同于向上兩相流動中液相速度徑向分布在近壁區(qū)逐漸降低，向下流動液相速度在壁面附近速度較高[9]，而相應位置的氣相速度則較低，故向下流動中氣液相對速度沿直徑方向均為負值。結(jié)合以上分析可發(fā)現(xiàn)，由于相對速度項為負值，使得升力為負值，從而驅(qū)使氣泡向圓管中心區(qū)域運動，而又由于湍流擴散力具有使管內(nèi)氣泡均勻分布的作用，故兩力的共同作用致使局部界面參數(shù)徑向呈“寬峰”型分布。而根據(jù)Rouhani[15]的研究結(jié)果可知，由于管道內(nèi)湍流的作用產(chǎn)生離心力，驅(qū)使氣泡向壁面運動，而由于壁面力抑制氣泡向壁面運動從而使局部參數(shù)呈“壁峰”分布。

4 結(jié)論

通過自研的雙探頭光纖探針，獲得了豎直向下氣液兩相流的局部空泡份額、氣泡頻率及界面面積濃度實驗數(shù)據(jù)及分布規(guī)律，主要研究結(jié)論如下：

1) 空泡份額、氣泡頻率及界面面積濃度徑向分布相類似，不同工況下均呈“寬峰”或“壁峰”型分布。

2) 氣泡頻率隨液相表觀速度的增加逐漸增大，而空泡份額、界面面積濃度則逐漸減小。

3) 向下氣液兩相流動中空泡份額、界面面積濃度截面平均值較向上兩相流動中的高。在液相表觀速度為0.57 m/s時，向下流動時空泡份額截面平均值均比向上流動時大119.6%～145.0%，界面面積濃度截面平均值比向上流動時大18.8%～82.5%。

參考文獻：

[1] 車得福. 多相流及其應用[M]. 西安：西安交通大學出版社，2007：2-10.

[2] 陳聽寬，楊魯偉. 壓水堆冷凝回流特性的研究[J]. 工程熱物理學報，2000，21(4)：491-496.

CHEN Tingkuan, YANG Luwei. An investigation of PWR reflux condensation characteristics[J]. Journal of Engineering Thermophysics, 2000, 21(4): 491-496(in Chinese).

[3] HIBIKI T, GODA H, KIM S, et al. Axial development of interfacial structure of vertical downward bubbly flow[J]. International Journal of Heat and Mass Transfer, 2004, 48(5): 749-764.

[4] HIBIKI T, ISHII M. Distribution parameter and drift velocity of drift-flux model in bubbly flow[J]. International Journal of Heat and Mass Transfer, 2002, 45(4): 707-721.

[5] 孫波，孫立成，田道貴，等. 光纖探針方法測量界面面積濃度實驗研究[J]. 原子能科學技術，2013，47(3)：432-436.

SUN Bo, SUN Licheng, TIAN Daogui, et al. Experimental study of interfacial area concentration by optical fiber probe method[J]. Atomic Energy Science and Technology, 2013, 47(3): 432-436(in Chinese).

[6] SHEN X, MATSUI R, MISHIMA K, et al. Distribution parameter and drift velocity for two-phase flow in large diameter pipe[J]. Nuclear Engineering and Design, 2010, 240(7): 3 991-4 000.

[7] 田道貴，孫立成，高菲，等. 光學探針在氣液兩相流局部參數(shù)測量中應用研究[J]. 實驗流體力學，2012，26(6)：91-95.

TIAN Daogui, SUN Licheng, GAO Fei, et al. A study on application of optical probes for the measurement of local parameters in two-phase flow[J]. Journal of Experiments in Fluid Mechanics (China), 2012, 26(6): 91-95(in Chinese).

[8] 孫波，孫立成，幸奠川，等. 豎直大圓管內(nèi)界面面積濃度分布特性[J]. 化工學報，2012，63(6)：1 810-1 815.

SUN Bo, SUN Licheng, XING Dianchuan, et al. Distribution profile of interfacial area concentration in vertical and large circular tubes[J]. Journal of Chemical Industry and Engineering (China), 2012, 63(6): 1 810-1 815(in Chinese).

[9] ISHII M, PARANJAPE S S, KIM S, et al. Interfacial structures and interfacial area transport in downward two-phase bubbly flow[J]. International Journal of Multiphase Flow, 2004, 30(2): 779-801.

[10] BHAGWAT S M, GHAJAR A J. Similarities and differences in the flow patterns and void fraction in vertical upward and downward two phase flow[J]. Experimental Thermal and Fluid Science, 2012, 39(12): 213-227.

[11] HIBIKI T, GODA H, KIM S, et al. Bubble structure of vertical downward bubbly flow[J]. International Journal of Heat and Mass Transfer, 2004, 47(9): 1 847-1 863.

[12] 劉崗，李會雄，薛玉卿，等. 垂直下降管內(nèi)空氣-水兩相流流型特性研究[C]∥中國工程熱物理學會多相流學術年會論文集. 西安：中國工程熱物理學會，2012.

[13] 賈斗南，喻真烷，王正杰，等. 低流速氣液兩相流垂直向上和向下時空泡份額的研究[J]. 核科學與工程，1985，5(2)：181-184.

JIA Dounan, YU Zhenwan, WANG Zhengjie, et al. Experimental study of void fraction in vertically upward and downward liquid-gas two-phase flow at low velocity[J]. Chinese Journal of Nuclear Science and Engineering, 1985, 5(2): 181-184(in Chinese).

[14] TOMIYAMA A, TAMAI H, ZUN I, et al. Transverse migration of single bubbles in simple shear flows[J]. Chemical Engineering Science, 2002, 57(10): 1 849-1 858.

[15] ROUHANI Z. Effect of wall friction and vortex generation on the radial distribution of different phases[J]. International Journal of Multiphase Flow, 1976, 3(1): 35-50.