Second order linear decoupled energy dissipation rate preserving schemes for the Cahn-Hilliard-Extended-Darcy model

Abstract

The thermodynamically consistent Cahn-Hilliard-Extended-Darcy (CHED) model has been used to describe the transient motion of a binary incompressible fluid flow in porous media. In this paper, we develop a series of linear, second-order, energy-dissipation-rate preserving numerical algorithms for the CHED model based on the energy quadratization strategy. We first extend the incompressible CHED model into a weakly compressible, thermodynamically consistent one using the generalized Onsager principle. Guided by the weakly compressible model, we then devise a series of linear, second-order, decoupled, semi-discrete, temporal algorithms in the spirit of the projection and the energy quadratization (EQ) method. The fully discrete algorithms are obtained by the use of the second-order finite difference method on staggered grids in space. We show that the obtained numerical algorithms respect the energy-dissipation-rate at the discrete level without any restrictions on the time step, making them unconditionally energy stable. Mesh refinement tests, coarsening dynamics of binary fluids, and the buoyancy-driven binary fluid motion in porous media are investigated numerically. The numerical results compare well with the results in the literature.