Numerical investigation on the influence of surface tension and viscous force on the bubble dynamics with a CLSVOF method

• Published in: Journal of mechanical science and technology Vol. 30; no. 6; pp. 2547 - 2556
• Main Authors: Wang, Zhiying, Li, Yikai, Huang, Biao, Gao, Deming
• Format: Journal Article
• Language: English
• Published: Seoul Korean Society of Mechanical Engineers 01.06.2016; Springer Nature B.V; 대한기계학회
• Subjects: Algebra; Axisymmetric; Axisymmetric flow; Bond number; Bubbles; Computational fluid dynamics; Computer simulation; Conservation; Control; Deformation; Dynamical Systems; Engineering; Fluid flow; Incompressible flow; Industrial and Production Engineering; Investigations; Mathematical analysis; Mathematical models; Mechanical Engineering; Navier-Stokes equations; Reynolds number; Surface tension; Two phase flow; Vibration; 기계공학; Bubble shape and motion; CLSVOF; Surface tension; Viscous force
• ISSN: 1738-494X; 1976-3824
• DOI: 10.1007/s12206-016-0516-8

We numerically investigated the rising of bubbles in a quiescent liquid layer. The numerical simulation is performed by solving the incompressible, multiphase Navier-Stokes equations via computational code in axisymmetric coordinates using a Coupled level-set and volume-of-fluid (CLSVOF) method. The numerical results show that the CLSVOF method with a novel algebraic relation between F and ∅ for axisymmetric two-phase flows not only can predict the bubble surface accurately, but also overcome the deficiency in preserving volume conservation. The effects of the Reynolds number Re and the Bond number Bo on the bubble deformation and its motion are investigated. The results show that with the increasing of Re (10 < Re < 150), the bubble shape transfers from oblate ellipsoidal cap to toroidal when Bo = 116. With the increasing of Bo (10 < Bo < 700), the bubble shape transfers from oblate ellipsoidal to toroidal when Re = 30. Although the toroidal bubble shapes are reached in these two cases, the transition modes are different. For the case Bo = 116, the bubble front is pierced by an upward jet from the rear of the bubble. While for the case Re = 30, the rear of the bubble is pierced by a downward jet from the front part.