An arbitrary Lagrangian‐Eulerian framework with exact mass conservation for the numerical simulation of 2D rising bubble problem

Abstract

SummaryAn arbitrary Lagrangian‐Eulerian framework, which combines the advantages of both Lagrangian and Eulerian methods, is presented to solve incompressible multiphase flow problems. The incompressible Navier‐Stokes equations are discretized using the side‐centered unstructured finite volume method, where the velocity vector components are defined at the midpoint of each cell face, while the pressure term is defined at element centroids. The pressure field is treated to be discontinuous across the interface with the discontinuous treatment of density and viscosity. The surface tension term at the interface is treated as a force tangent to the interface and computed with several different approaches including the use of Legendre polynomials. In addition, the several different discretizations of interface kinematic boundary conditions are investigated. For the application of the interface kinematic boundary condition, a special attention is given to satisfy both local and global discrete geometric conservation law to conserve the total mass of both species at machine precision. The mesh vertices are deformed by solving the linear elasticity equations due to the normal displacement of interface. The resulting algebraic equations are solved in a fully coupled manner, and a one‐level restricted additive Schwarz preconditioner with a block‐incomplete factorization within each partitioned subdomain is used for the resulting fully coupled system. The method is validated by simulating the classical benchmark problem of a single rising bubble in a viscous fluid due to buoyancy. The results of numerical simulations are found out to be in an excellent agreement with the earlier results in the literature. The mass of the bubble is conserved, and discontinuous pressure field is obtained to avoid errors due to the incompressibility condition in the vicinity of the interface, where the density and viscosity jumps occur.