Numerical study of rising bubbles with path instability using conservative level-set and adaptive mesh refinement

• Published in: Computers & fluids Vol. 187; pp. 83 - 97
• Main Authors: Antepara, Oscar, Balcázar, Néstor, Rigola, Joaquim, Oliva, Assensi
• Format: Journal Article Publication
• Language: English
• Published: Amsterdam Elsevier Ltd 15-06-2019; Elsevier BV; Elsevier
• Subjects: Adaptive mesh refinement; Bombolles; Bubbles; Computational fluid dynamics; Computer simulation; Conservative level-set; Correlation analysis; Dinàmica; Direct numerical simulation; Discretization; Drag coefficients; Dynamics; Enginyeria mecànica; Finite element method; Fluid flow; Grid refinement (mathematics); Mecànica de fluids; Path instability; Reynolds number; Simulació per ordinador; Stability; Velocity coupling; Wobbling bubble; Àrees temàtiques de la UPC; Conservative level-set; Path instability; Adaptive mesh refinement; Direct numerical simulation; Wobbling bubble
• ISSN: 0045-7930; 1879-0747
• DOI: 10.1016/j.compfluid.2019.04.013

•Numerical study of rising bubbles with path instability.•Combination of CLS and AMR methods.•Bubble path, bubble shape and vortical structures at high Reynolds number.•Effect of the dimensionless number in the wobbling regime.
This paper focuses on three-dimensional direct numerical simulations of rising bubbles in the wobbling regime, and the study of its dynamical behavior for Eötvös number 1  ≤  Eo  ≤  10 and Morton number 1e−11  ≤ M ≤  1e−9. The computational methodology is based on a mass Conservative Level-Set method, whereas the spatial discretization of the computational domain employs an Adaptive Mesh Refinement strategy for the reduction of computational resources. The Navier–Stokes equations are discretized using the finite-volume approach on a collocated unstructured mesh; the pressure-velocity coupling is solved using a classical fractional-step projection method. This methodology is applied to a series of verification and validation tests, which are compared with experiments and numerical results from the literature. Finally, buoyancy bubbles rising in the wobbling regime are researched at moderate to high Reynolds numbers (100 < Re < 3000). Terminal Reynolds number, drag coefficient and frequency of path oscillations are compared with empirical correlations and numerical studies from the literature. Results show the discharge of alternate oppositely-oriented hairpin vortex structures. Moreover, depending on the characteristics numbers of the system, different path features, bubble shape, and vortical structures in the wake are reported.