A Cartesian Method with Second-Order Pressure Resolution for Incompressible Flows with Large Density Ratios

An Eulerian method to numerically solve incompressible bifluid problems with high density ratio is presented. This method can be considered as an improvement of the Ghost Fluid method, with the specificity of a sharp second-order numerical scheme for the spatial resolution of the discontinuous elliptic problem for the pressure. The Navier–Stokes equations are integrated in time with a fractional step method based on the Chorin scheme and discretized in space on a Cartesian mesh. The bifluid interface is implicitly represented using a level-set function. The advantage of this method is its simplicity to implement in a standard monofluid Navier–Stokes solver while being more accurate and conservative than other simple classical bifluid methods. The numerical tests highlight the improvements obtained with this sharp method compared to the reference standard first-order methods