Numerical study of wall effects on buoyant gas-bubble rise in a liquid-filled finite cylinder.
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), Eötvös 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/d0, R*=R/d0). 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 (ReT) and shapes have been studied. The range of terminal state Reynolds numbers includes 0.02<ReT<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*> or =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 ReT. Our predictions agree with existing results reported in the literature.
Mukundakrishnan
Karthik
K
Department of Anesthesiology and Critical Care, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA.
Quan
Shaoping
S
Eckmann
David M
DM
Ayyaswamy
Portonovo S
PS
D003198:Computer Simulation; D005740:Gases / Q000737:chemistry*; D006112:Gravitation; D044085:Microfluidics / Q000379:methods*; D008962:Models, Theoretical; D012996:Solutions / Q000737:chemistry*; D013499:Surface Properties
D005740:Gases; D012996:Solutions