A rotational pressure-correction discontinuous Galerkin scheme for the Cahn-Hilliard-Darcy-Stokes system

This paper is devoted to the numerical approximations of the Cahn-Hilliard-Darcy-Stokes system, which is a combination of the modified Cahn-Hilliard equation with the Darcy-Stokes equation. A novel discontinuous Galerkin pressure-correction scheme is proposed for solving the coupled system, which can achieve the desired level of linear, fully decoupled, and unconditionally energy stable. The developed scheme here is implemented by combining several effective techniques, including by adding an additional stabilization term artificially in Cahn-Hilliard equation for balancing the explicit treatment of the coupling term, the stabilizing strategy for the nonlinear energy potential, and a rotational pressure-correction scheme for the Darcy-Stokes equation. We rigorously prove the unique solvability, unconditional energy stability, and optimal error estimates of the proposed scheme. Finally, a number of numerical examples are provided to demonstrate numerically the efficiency of the present formulation.