Front tracking applied to Rayleigh Taylor instability
SIAM Journal on Scientific and Statistical Computing
Sedimentation in ordered emulsions of drops at low Reynolds numbers
Zeitschrift für Angewandte Mathematik und Physik (ZAMP)
The bifurcation of tracked scalar waves
SIAM Journal on Scientific and Statistical Computing - Telecommunication Programs at U.S. Universities
Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations
Journal of Computational Physics
Journal of Computational Physics
A front-tracking method for viscous, incompressible, multi-fluid flows
Journal of Computational Physics
A continuum method for modeling surface tension
Journal of Computational Physics
Improved volume conservation in the computation of flows with immersed elastic boundaries
Journal of Computational Physics
A level set approach for computing solutions to incompressible two-phase flow
Journal of Computational Physics
Robust Computational Algorithms for Dynamic Interface Tracking in Three Dimensions
SIAM Journal on Scientific Computing
Interfacial dynamics for Stokes flow
Journal of Computational Physics
Boundary integral methods for multicomponent fluids and multiphase materials
Journal of Computational Physics
Journal of Computational Physics
Level set methods: an overview and some recent results
Journal of Computational Physics
A front-tracking method for the computations of multiphase flow
Journal of Computational Physics
Journal of Computational Physics
A moving mesh interface tracking method for 3D incompressible two-phase flows
Journal of Computational Physics
Numerical simulation of bubble rising in viscous liquid
Journal of Computational Physics
Journal of Computational Physics
Journal of Computational Physics
Connectivity-free front tracking method for multiphase flows with free surfaces
Journal of Computational Physics
Simulation of bubble dynamics in a microchannel using a front-tracking method
Computers & Mathematics with Applications
Hi-index | 31.46 |
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.