A volume of fluid method algorithm for simulation of surface tension dominant two-phase flows

Abstract

In this article, a two-phase flow algorithm that uses an explicit Adams–Bashforth scheme coupled with volume of fluid (VOF) method in a uniform staggered Cartesian grid for surface tension dominant simulation is presented. The interface reconstruction procedure is implemented using the Weighed Linear Interface Calculation (WLIC) algorithm. The level set (LS) function instead of VOF function is adopted for calculation of not only the interface normal vectors in the WLIC algorithm but also the surface tension terms in the Navier–Stokes equations. This VOF method is compared with LS method and other experiment results in the literature for vortex deforming, milk crown, droplet impacting into deep liquid pool and bubble bursting at a free surface problems to show whether mass can be well conserved and breaking and merging phenomena can be accurately simulated.