采用界面追踪方法(FTM)研究了狭窄通道内受浮力上升的单气泡经过突然扩张通道时的运动形态，考察了狭窄通道和宽阔通道的通道宽度比及Eotvos数(Eo)和Morton数(Mo)对气泡形态和上升速度的影响. 结果表明，在狭窄通道、扩张部位和宽阔通道内气泡形态呈不同特性，气泡在离开扩张处时高宽比最大. 随通道宽度比增大，气泡最大高宽比减小. Eo数较小时，减小Mo数，气泡高宽比减小，气泡上升速度增大；增大Eo数，气泡上升时形变更严重. 改变初始时刻气泡相对于狭窄通道的位置对气泡形态有影响，初始位置靠近单侧壁面时，气泡横向变形较严重，在两侧壁面之间摆动上升.
A single bubble driven by buoyance rising through a sudden expansive channel was simulated by Front Tracking Method (FTM). Firstly, the influence of the width ratio of the narrow channel and the wide channel on shape and rising velocity of bubble was investigated. It was found that the morphology of bubble interface presented different characteristics in the three parts, i.e., the narrow channel, the expansion part and the wide channel, respectively. When the bubble left the expansion, the surrounding fluid flowed from both sides of the bubble into the bottom area, causing both sides of the bubble to be squeezed, the bubble reached the highest aspect ratio. With the increase of the channel width ratio, the flow field of the bubble was less inhibited by both sides of the wall, and the maximum aspect ratio of the bubble decreased. When Wn/We reached over 4/6, the influence of expansion segment can be ignored. Then the influence of Eotvos number Eo and Mo number were analyzed. At low Eo, when the Mo number reduced, the bubble was more likely to deform, the aspect ratio of the bubble decreased and the bubble?s rising speed increased. When the Eo number increased, the deformation became more serious when the bubble rised. Changing the position of the bubble relative to the narrow channel at the initial moment also affected the shape of the bubble. When the initial position of the bubble was close to the single side wall, due to the low flow velocity of the surrounding fluid near the wall, the bubbles would be stretched horizontally while moving, and the bubble deformed badly in the horizontal direction. Asymmetrical pressure distributions occurred on the left and right sides of the bubble level. The pressure on the side of the bubble near the wall surface was lower, making the movement path of the bubble in the narrow channel appear to swing left and right. After leaving the expansion part, the movement of bubbles gradually stabilized due to the absence of wall constraint on the bubbles. When the bubble was far away from the expansion part, the influence of the wall surface decreased. After 2.3 times from the expansion part, it can be considered that the bubble enters a stable rising state.
[1]陶然，权晓波，徐建中.微尺度流动研究中的几个问题[J].工程热力学报, 2001, 22(5):575-577
[2]陈光文，袁权.微化工技术[J].化工学报, 2003, 54(4):427-439
[3] Mehmet Y G, Ahmet I K, Hasan D.Analytical prediction of pressure loss through a sudden-expansion in two-phase pneumatic conveying linesAdvanced Powder Technology, 2009, 20: 48-54.
[4]YIN Z G, XIE S H, CHENG D S, Long X Y.Oxygen transfer characteristics of water and bubble mixture pipe flow through two sudden contractions and expansions[J].Journal of Hydrodynamics, 2014, 26(5):745-750
[5]Grace J R.Shapes and velocities of bubbles rising in infinite liquids[J].Trans. Inst. Chem. Eng, 1973, 51(2):116-120
[6] D.Bhaga, M.E, Weber. Bubbles in viscous liquids: shapes, wakes and velocities [J]. Journal of Fluid Mechanics, 1981, 105: 61-86.
[7] Chen I Y, Tseng C Y, Lin Y T, Wang C C.Two-phase flow pressure change subject to sudden contraction in small rectangular channels [J]. International Journal of Multiphase Flow, 2009, 35: 297-306.
[8] A.Ahmadpour, S.M.A, Noori Rahim Abadi, R. Kouhikamali. Numerical simulation of two-phase gas-liquid flow through gradual expansions/contractions [J]. International Journal of Multiphase Flow, 2016, 79: 31-49.
[9] Erfan Kadivar, Mojtaba Farrokhbin.A numerical procedure for scaling droplet deformation in a microfluidic expansion channel [J]. Physica A, 2017, 479: 449-459.
[10] Stephen Ambrose, Ian S.Lowndes, David M. Hargreaves, Barry Azzopardi. Numerical modelling of the rise of Taylor bubbles through a change in pipe diameter [J]. Computers and Fluids, 2017, 148: 10-25.
[11]R H Chen, W.X. Tian,GH. Su,S.Z. Qiu,Yuki Ishiwatari,Yoshiaki Oka. Numerical investigation on coalescence of bubble pairs rising in a stagnant liquid[J].Chemical Engineering Science, 2011, 66(21):5055-5063
[12] Ryu Seungyeob, Ko Sungho.A comparative study of lattice Boltzmann and volume of fluid method for two-dimensional multiphase flows [J]. Nuclear Engineering and Technology, 2012, 44(6), 623– 638.
[13]陈斌.静止水中单个上升气泡的直接数值模拟[J].工程热物理学报, 2005, 26(2):980-982
[14] 陈斌.高粘度流体中上升气泡的直接数值模拟 [J]. 工程热物理学报, 2006. 27(2): 255-258.
[15]王焕然，李彦鹏，杨栋，席光.黏性液体中单个气泡上升的形状特性[J].工程热物理学报, 2009, 30(9):1492-1494
[16] Néstor Balcázar, Oriol Lehmkuhl, Lluís Jofre, Assensi Oliva.Level-set simulations of buoyancy-driven motion of single and multiple bubbles [J]. International Journal of Heat and Fluid Flow, 2015, 56: 91-107.
