When simulating high Reynolds number two-phase flow, boundary layers develop at the interface, which are much thinner compared to the capillary length-scales that are of interest. Resolving such an interface layer is expensive and therefore it is often not resolved in a simulation. Numerically such an underresolved interface layer results in a velocity discontinuity tangential to the interface. We propose to include such tangential velocity discontinuities in our numerical model. This results in a sharp two-fluid model for two-phase flow, where only the interface-normal component of the velocity field is smooth. This condition is implicitly enforced via a new jump condition on the pressure gradient, which we discretize using a multidimensional variant of the ghost fluid method . Results are shown of breaking waves  as well as (breaking) waves impacting a solid wall  where we compare to state-of-the-art methods [3, 4]. We show that our proposed method is able to accurately simulate high Reynolds number two-phase flow without the need for resolving, or artificially thickening, of the interface layer.
Are you one of the authors of this document?