Bound‐preserving discontinuous Galerkin methods for compressible two‐phase flows in porous media

This paper presents a numerical study of immiscible, compressible two‐phase flows in porous media, that takes into account heterogeneity, gravity, anisotropy and injection/production wells. We formulate a fully implicit stable discontinuous Galerkin solver for this system that is accurate, that respects maximum principle for the approximation of saturation, and that is locally mass conservative. To completely eliminate the overshoot and undershoot phenomena, we construct a flux limiter that produces bound‐preserving elementwise average of the saturation. The addition of a slope limiter allows to recover a pointwise bound‐preserving discrete saturation. Numerical results show that both maximum principle and monotonicity of the solution are satisfied. The proposed flux limiter does not impact the local mass error and the number of nonlinear solver iterations.