An energy stable linear numerical method for thermodynamically consistent modeling of two-phase incompressible flow in porous media

In this paper, we consider numerical approximation of a thermodynamically consistent model of two-phase flow in porous media, which obeys an intrinsic energy dissipation law. The model under consideration is newly-developed, so there is no energy stable numerical scheme proposed for it at present. This model consists of two nonlinear degenerate parabolic equations and a saturation constraint, but lacking an independent equation for the pressure, so for the purpose of designing efficient numerical scheme, we reformulate the model forms as well as the free energy function, and further prove the corresponding energy dissipation inequality. Based on the alternative reformulations, using the invariant energy quadratization approach and subtle semi-implicit treatments for the pressure and saturation, we for the first time propose a linear and energy stable numerical method for this model. The fully discrete scheme is devised combining the upwind approach for the phase mobilities and the cell-centered finite difference method. The unique solvability of numerical solutions and unconditional energy stability are rigorously proved for both the semi-discrete time scheme and the fully discrete scheme. Moreover, the scheme can guarantee the local mass conservation for both phases. We also show that the upwind mobility approach plays an essential role in preserving energy stability of the fully discrete scheme. Numerical results are presented to demonstrate the performance of the proposed scheme.