A new physics-preserving IMPES scheme for incompressible and immiscible two-phase flow in heterogeneous porous media

TL;DR

This paper introduces a novel physics-preserving IMPES scheme for simulating incompressible, immiscible two-phase flow in heterogeneous porous media, ensuring local mass conservation and bounded saturations with improved stability.

Contribution

The paper develops a new IMPES scheme that is inherently physics-preserving, mass conservative, and unbiased for two-phase flow in heterogeneous media, with enhanced stability and robustness.

Findings

The scheme maintains local mass conservation for both phases.

Saturations remain within physical bounds under certain time step conditions.

The method demonstrates high efficiency and robustness in numerical examples.

Abstract

In this work we consider a new efficient IMplicit Pressure Explicit Saturation (IMPES) scheme for the simulation of incompressible and immiscible two-phase flow in heterogeneous porous media with capillary pressure. Compared with the conventional IMPES schemes, the new IMPES scheme is inherently physics-preserving, namely, the new algorithm is locally mass conservative for both phases and it also enjoys another appealing feature that the total velocity is continuous in the normal direction. Moreover, the new scheme is unbiased with regard to the two phases and the saturations of both phases are bounds-preserving if the time step size is smaller than a certain value. The key ideas in the new scheme include that the Darcy flows for both phases are rewritten in the formulation based on the total velocity and an auxiliary velocity referring to as the capillary potential gradient, and the total discretized conservation equation is obtained by the summation of the discretized conservation equation for each phase. The upwind strategy is applied to update the saturations explicitly, and the upwind mixed finite element methods are used to solve the pressure-velocity systems which can be decoupled. We also present some interesting examples to demonstrate the efficiency and robustness of the new algorithm.