A Novel Energy Stable Numerical Scheme for Navier-Stokes-Cahn-Hilliard Two-Phase Flow Model with Variable Densities and Viscosities

Abstract

A novel numerical scheme including time and spatial discretization is offered for coupled Cahn-Hilliard and Navier-Stokes governing equation system in this paper. Variable densities and viscosities are considered in the numerical scheme. By introducing an intermediate velocity in both Cahn-Hilliard equation and momentum equation, the scheme can keep discrete energy law. A decouple approach based on pressure stabilization is implemented to solve the Navier-Stokes part, while the stabilization or convex splitting method is adopted for the Cahn-Hilliard part. This novel scheme is totally decoupled, linear, unconditionally energy stable for incompressible two-phase flow diffuse interface model. Numerical results demonstrate the validation, accuracy, robustness and discrete energy law of the proposed scheme in this paper.