We propose a novel second order in time numerical scheme for Cahn-Hilliard-Navier- Stokes phase field model with matched density. The scheme is based on second order convex-splitting for the Cahn-Hilliard equation and pressure-projection for the Navi
er-Stokes equation. We show that the scheme is mass-conservative, satisfies a modified energy law and is therefore unconditionally stable. Moreover, we prove that the scheme is uncondition- ally uniquely solvable at each time step by exploring the monotonicity associated with the scheme. Thanks to the weak coupling of the scheme, we design an efficient Picard iteration procedure to further decouple the computation of Cahn-Hilliard equation and Navier-Stokes equation. We implement the scheme by the mixed finite element method. Ample numerical experiments are performed to validate the accuracy and efficiency of the numerical scheme.
We study the well-posedness of a coupled Cahn-Hilliard-Stokes-Darcy system which is a diffuse-interface model for essentially immiscible two phase incompressible flows with matched density in a karstic geometry. Existence of finite energy weak soluti
on that is global in time is established in both 2D and 3D. Weak-strong uniqueness property of the weak solutions is provided as well.
