SPLITTING SCHEMES FOR A NAVIER-STOKES-CAHN-HILLIARD MODEL FOR TWO FLUIDS WITH DIFFERENT DENSITIES

In this work, we focus on designing efficient numerical schemes to approximate a thermodynamically consistent Navier-Stokes/Cahn-Hilliard problem given in [3] modeling the mixture of two incompressible fluids with different densities. The model is based on a diffuse-interface phase-field approach that is able to describe topological transitions like droplet coalescense or droplet break-up in a natural way. We present a splitting scheme, decoupling computations of the Navier-Stokes part from the Cahn-Hilliard one, which is unconditionally energy-stable up to the choice of the potential approximation. Some numerical experiments are carried out to validate the correctness and the accuracy of the scheme, and to study the sensitivity of the scheme with respect to different physical parameters.