Anti-diffusion method for interface steepening in two-phase incompressible flow

In this paper, we present a method for obtaining sharp interfaces in two-phase incompressible flows by an anti-diffusion correction, that is applicable in a straight-forward fashion for the improvement of two-phase flow solution schemes typically employed in practical applications. The underlying discretization is based on the volume-of-fluid (VOF) interface-capturing method on unstructured meshes. The key idea is to steepen the interface, independently of the underlying volume-fraction transport equation, by solving a diffusion equation with reverse time, i.e. an anti-diffusion equation, after each advection time step of the volume fraction. As the solution of the anti-diffusion equation requires regularization, a limiter based on the directional derivative is developed for calculating the gradient of the volume fraction. This limiter ensures the boundedness of the volume fraction. In order to control the amount of anti-diffusion introduced by the correction algorithm we propose a suitable stopping criterion for interface steepening. The formulation of the limiter and the algorithm for solving the anti-diffusion equation are applicable to 3-dimensional unstructured meshes. Validation computations are performed for passive advection of an interface, for 2-dimensional and 3-dimensional rising-bubbles, and for a rising drop in a periodically constricted channel. The results demonstrate that sharp interfaces can be recovered reliably. They show that the accuracy is similar to or even better than that of level-set methods using comparable discretizations for the flow and the level-set evolution. Also, we observe a good agreement with experimental results for the rising drop where proper interface evolution requires accurate mass conservation.