Conservative unconditionally stable decoupled numerical schemes for the Cahn-Hilliard-Navier-Stokes-Darcy-Boussinesq system

We propose two mass and heat energy conservative, unconditionally stable, decoupled numerical algorithms for solving the Cahn-Hilliard-Navier-Stokes-Darcy-Boussinesq system that models thermal convection of two-phase flows in superposed free flow and porous media. The schemes totally decouple the computation of the Cahn-Hilliard equation, the Darcy equations, the heat equation, the Navier-Stokes equations at each time step, and thus significantly reducing the computational cost. We rigorously show that the schemes are conservative and energy-law preserving. Numerical results are presented to demonstrate the accuracy and stability of the algorithms.