Numerical study of wall effects on buoyant gas-bubble rise in a liquid-filled finite cylinder

被引:64
作者
Mukundakrishnan, Karthik
Quan, Shaoping
Eckmann, David M.
Ayyaswamy, Portonovo S. [1 ]
机构
[1] Univ Penn, Dept Mech Engn & Appl Mech, Philadelphia, PA 19104 USA
[2] Univ Penn, Dept Anesthesiol & Crit Care, Philadelphia, PA 19104 USA
来源
PHYSICAL REVIEW E | 2007年 / 76卷 / 03期
关键词
D O I
10.1103/PhysRevE.76.036308
中图分类号
O35 [流体力学]; O53 [等离子体物理学];
学科分类号
070204 ; 080103 ; 080704 ;
摘要
The wall effects on the axisymmetric rise and deformation of an initially spherical gas bubble released from rest in a liquid-filled, finite circular cylinder are numerically investigated. The bulk and gas phases are considered incompressible and immiscible. The bubble motion and deformation are characterized by the Morton number (Mo), Eotvos number (Eo), Reynolds number (Re), Weber number (We), density ratio, viscosity ratio, the ratios of the cylinder height and the cylinder radius to the diameter of the initially spherical bubble (H-*=H/d(0), R-*=R/d(0)). Bubble rise in liquids described by Eo and Mo combinations ranging from (1,0.01) to (277.5,0.092), as appropriate to various terminal state Reynolds numbers (Re-T) and shapes have been studied. The range of terminal state Reynolds numbers includes 0.02 < Re-T < 70. Bubble shapes at terminal states vary from spherical to intermediate spherical-cap-skirted. The numerical procedure employs a front tracking finite difference method coupled with a level contour reconstruction of the front. This procedure ensures a smooth distribution of the front points and conserves the bubble volume. For the wide range of Eo and Mo examined, bubble motion in cylinders of height H-*=8 and R-*>= 3, is noted to correspond to the rise in an infinite medium, both in terms of Reynolds number and shape at terminal state. In a thin cylindrical vessel (small R-*), the motion of the bubble is retarded due to increased total drag and the bubble achieves terminal conditions within a short distance from release. The wake effects on bubble rise are reduced, and elongated bubbles may occur at appropriate conditions. For a fixed volume of the bubble, increasing the cylinder radius may result in the formation of well-defined rear recirculatory wakes that are associated with lateral bulging and skirt formation. The paper includes figures of bubble shape regimes for various values of R-*, Eo, Mo, and Re-T. Our predictions agree with existing results reported in the literature.
引用
收藏
页数:15
相关论文
共 49 条
[1]   A numerical method for the incompressible Navier-Stokes equations based on an approximate projection [J].
Almgren, AS ;
Bell, JB ;
Szymczak, WG .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 1996, 17 (02) :358-369
[2]   Numerical simulation of behavior of gas bubbles using a 3-D front-tracking method [J].
Annaland, MV ;
Dijkhuizen, W ;
Deen, NG ;
Kuipers, JAM .
AICHE JOURNAL, 2006, 52 (01) :99-110
[3]  
Batchelor GK, 1968, INTRO FLUID MECH
[4]   A 2ND-ORDER PROJECTION METHOD FOR THE INCOMPRESSIBLE NAVIER STOKES EQUATIONS [J].
BELL, JB ;
COLELLA, P ;
GLAZ, HM .
JOURNAL OF COMPUTATIONAL PHYSICS, 1989, 85 (02) :257-283
[5]   A 2ND-ORDER PROJECTION METHOD FOR VARIABLE-DENSITY FLOWS [J].
BELL, JB ;
MARCUS, DL .
JOURNAL OF COMPUTATIONAL PHYSICS, 1992, 101 (02) :334-348
[6]   BUBBLES IN VISCOUS-LIQUIDS - SHAPES, WAKES AND VELOCITIES [J].
BHAGA, D ;
WEBER, ME .
JOURNAL OF FLUID MECHANICS, 1981, 105 (APR) :61-85
[7]   Transition from spherical cap to toroidal bubbles [J].
Bonometti, T ;
Magnaudet, J .
PHYSICS OF FLUIDS, 2006, 18 (05)
[8]   A CONTINUUM METHOD FOR MODELING SURFACE-TENSION [J].
BRACKBILL, JU ;
KOTHE, DB ;
ZEMACH, C .
JOURNAL OF COMPUTATIONAL PHYSICS, 1992, 100 (02) :335-354
[9]   Interfacial dynamics of stationary gas bubbles in flows in inclined tubes [J].
Cavanagh, DP ;
Eckmann, DM .
JOURNAL OF FLUID MECHANICS, 1999, 398 :225-244
[10]   The development of a bubble rising in a viscous liquid [J].
Chen, L ;
Garimella, SV ;
Reizes, JA ;
Leonardi, E .
JOURNAL OF FLUID MECHANICS, 1999, 387 :61-96