To further study bubble coalescence and develop a highly efficient tracking interface program
an unsplit algorithm based on the PLIC interface reconstruction algorithm was applied to precisely calculate the fluid transport from one cell to its neighboring 26 cells
and the fluid transport of each cell was calculated using a unified expression which greatly reduces the programming difficulty and improves the calculation accuracy. The Navier-Stokes equations for gas-liquid two-phase flow were solved by combining an advection equation for liquid volume fraction to simulate the rise of the single bubble. The terminal velocities obtained by the numerical simulation for the single bubble with the diameter of 0.8-14 mm were in good agreement with the experimental data from literature. The co-axial coalescence of two gas bubbles in the initially quiescent liquid was then simulated by the unsplit algorithm. The numerical results show that when the surface tension is very high
the two bubbles rise and deform independently. The bubbles tend to break up when the surface tension is very small
and when the bubble breakup is serious
the following bubble would pass through the leading bubble. Additionally
the effect of liquid viscosity and surface tension on coalescence time is non-monotonous.
关键词
Keywords
references
WU M M, GHARIB M. Experimental studies on the shape and path of small air bubbles rising in clean water [J]. Physics of Fluids, 2002, 14(7): 49-52.
HE Dan, LI Yanpeng, LIU Yanyan. Effect of initial bubble shape on dynamics of a buoyancy-driven bubble [J]. Journal of Xi'an Jiaotong University, 2011, 45(1): 43-47.
INAMURO T, OGATA T, OGINO F. Numerical simulation of bubble flows by the lattice Boltzmann method [J]. FGCS, 2004, 20(6): 959-964.
YU Z, FAN L S. Direct simulation of the buoyant rise of bubbles in infinite liquid using level set method [J]. Canadian Journal of Chemical Engineering, 2008, 86(6): 267-275.
VAN SINT ANNALAND M, DEEN N G, KUIPERS J A M. Numerical simulation of gas bubbles behaviors using a three-dimensional volume of fluid method [J]. Chemical Engineering Science, 2005, 60(11): 2999-3011.
PU L, LI H X, LV X. Numerical simulation of bubble dynamics in microgravity [J]. Microgravity Sci Technol, 2008, 20(3/4): 247-251.
LI Yanpeng, ZHANG Qianlong, BAI Bofeng. Direct simulation of the rise of a bubble pair in a vertical channel [J]. Journal of Engineering for Thermal Energy and Power, 2007, 22(4): 375-379.
KOEBE M, BOTHE D, PRUESS J, et al. 3D direct numerical simulation of air bubbles in water at high Reynolds number [C]∥Proceeding of the 2002 ASME Fluids Engineering Division Summer Meeting. New York, USA: ASME, 2002: 14-18.
LIU Chun, MA Tianbao, NING Jianguo. Unsplit transport method in Eulerian method [J]. Transactions of Beijing Institute of Technology, 2008, 28(10): 847-850.
CLIFT R, GRACE J R, WEBER M E. Bubbles, drops and particles [M]. New York, USA: Academic Press, 1978.
JUNE K M, IL S P. Numerical study for laminar wavy motions of liquid film flow on vertical wall [J]. International Journal of Heat and Mass Transfer, 2011, 54(15/16): 3256-3266.
LORSTAD D, FUCHS L. High-order surface tension VOF-model for 3D bubble flows with high density ratio [J]. Journal of Computational Physics, 2004, 200(1): 153-176.
MENARD T, TANGUY S, BERLEMONT A. Coupling level set/VOF/ghost fluid methods: validation and application to 3D simulation of the primary break-up of a liquid jet [J]. International Journal of Multiphase Flow, 2007, 33(5): 510-524.
PASSANDIDEH-FARD M, FARHANGI M M. A numerical study on bubble rise and interaction in a viscous liquid [C]∥Proceeding of the 5th International Conference on Transport Phenomena in Multiphase Systems. Bialystok, Poland: Bialystok Technical University, 2008: 1-7.
BRERETON G, KOROTNEY D. Coaxial and oblique coalescence of two rising bubbles [C]∥The ASME Applied Mechanics Conference. New York, USA: ASME, 1991: 1-16.
Study of the Mechanism of Electrostatic Accumulation and Current Density Distribution Characteristics During Liquid Hydrogen Pipeline Transportation
Numerical Simulation of Space Charge Distribution Characteristics of Polyimide Materials under X-Ray Irradiation
Numerical Simulation Study on Geothermal Extraction Characteristics of Supercritical CO2 in Heterogeneous Fractures
Action Mechanism of Heat Transfer between the Working Medium and the Impeller on Aerodynamic Performance and Flow Characteristics of Radial Inflow Turbines with Supercritical Carbon Dioxide
Analysis of the Attached Vortex Model for the Flow Characteristics of Wing Near-Field Wakes
Related Author
SUN Wenhao
LIU Hongbao
WANG Lei
QING Ziyou
LI Zhuolun
MA Yuan
LI Yanzhong
ZHONG Hui
Related Institution
School of Energy and Power Engineering,Xi’an Jiaotong University
Aerospace System Engineering Shanghai
School of Electrical Engineering,Xi’an Jiaotong University
Northwest Institute of Nuclear Technology
Institute of Electronic Engineering,China Academy of Engineering Physics