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.
Yunzhou BAO Chenhan WU Ying ZHANG
. Numerical simulation of interface evolution of single bubble rising through sudden expansion channel[J]. The Chinese Journal of Process Engineering, 2018
, 18(6)
: 1178
-1186
.
DOI: 10.12034/j.issn.1009-606X.217435
[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.
