Efficient two-phase mass-conserving level set method for simulation of incompressible turbulent free surface flows with large density ratio

An efficient three-dimensional two-phase flow model for simulation of free surface flows is developed. The method is based in the solution of Navier-Stokes equation for quasi-incompressible flows. The equations in generalized coordinates are solved by extension of the fully explicit MacCormack scheme, second order in time and fourth order in space. The free surface is implicitly captured by the zero Level Set of a smooth function and by the Ghost Fluid Method to capture accurately shape discontinuities for properties in the vicinity of the interface. Finally, the Volume of Fluid method is used to ensure mass conservation. Turbulence is described by large eddy simulation, where only the large-scale eddies are solved, while the small scales are modeled by using the selective structure function subgrid-scale model. Boundary shapes are represented through the immersed boundary method on the Cartesian grid. The numerical model is validated by some free surface problems. Numerical predictions by all cases are in good agreement with experimental data. (C) 2016 Elsevier Ltd. All rights reserved.