Numerical simulation of 3D bubbles rising in viscous liquids using a front tracking method

• Published in: Journal of computational physics Vol. 227; no. 6; pp. 3358 - 3382
• Main Authors: Hua, Jinsong, Stene, Jan F., Lin, Ping
• Format: Journal Article
• Language: English
• Published: Amsterdam Elsevier Inc 01.03.2008; Elsevier
• Subjects: Adaptive mesh refinement; Bubble rising; Computational fluid dynamics; Computational techniques; Exact sciences and technology; Front tracking method; Incompressible flow; Mathematical methods in physics; Moving reference frame; Multiphase flow; Physics; SIMPLE algorithm; Incompressible flow; Computational fluid dynamics; Multiphase flow; SIMPLE algorithm; Front tracking method; Adaptive mesh refinement; Moving reference frame; Bubble rising; Viscosity; Viscous liquid; Sensitivity analysis; High density; Reynolds number; Digital simulation; Calculation methods; Large Reynolds number; Newtonian flow; Three dimensional model; Algorithms; Referential; Dynamics; Two-phase flow; Calculation; Navier-Stokes equations
• ISSN: 0021-9991; 1090-2716; 1090-2716
• DOI: 10.1016/j.jcp.2007.12.002

The rise of bubbles in viscous liquids is not only a very common process in many industrial applications, but also an important fundamental problem in fluid physics. An improved numerical algorithm based on the front tracking method, originally proposed by Tryggvason and his co-workers, has been validated against experiments over a wide range of intermediate Reynolds and Bond numbers using an axisymmetric model [J. Hua, J. Lou, Numerical simulation of bubble rising in viscous liquid, J. Comput. Phys. 22 (2007) 769–795]. In the current paper, this numerical algorithm is further extended to simulate 3D bubbles rising in viscous liquids with high Reynolds and Bond numbers and with large density and viscosity ratios representative of the common air–water two-phase flow system. To facilitate the 3D front tracking simulation, mesh adaptation is implemented for both the front mesh on the bubble surface and the background mesh. On the latter mesh, the governing Navier–Stokes equations for incompressible, Newtonian flow are solved in a moving reference frame attached to the rising bubble. Specifically, the equations are solved using a finite volume scheme based on the Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) algorithm, and it appears to be robust even for high Reynolds numbers and high density and viscosity ratios. The 3D bubble surface is tracked explicitly using an adaptive, unstructured triangular mesh. The numerical model is integrated with the software package PARAMESH, a block-based adaptive mesh refinement (AMR) tool developed for parallel computing. PARAMESH allows background mesh adaptation as well as the solution of the governing equations in parallel on a supercomputer. Further, Peskin distribution function is applied to interpolate the variable values between the front and the background meshes. Detailed sensitivity analysis about the numerical modeling algorithm has been performed. The current model has also been applied to simulate a number of cases of 3D gas bubbles rising in viscous liquids, e.g. air bubbles rising in water. Simulation results are compared with experimental observations both in aspect of terminal bubble shapes and terminal bubble velocities. In addition, we applied this model to simulate the interaction between two bubbles rising in a liquid, which illustrated the model’s capability in predicting the interaction dynamics of rising bubbles.