High-Order Bound-Preserving Local Discontinuous Galerkin Methods for Incompressible and Immiscible Two-Phase Flows in Porous Media

Document Type


Publication Date



Department of Mathematical Sciences


In this paper, we develop high-order bound-preserving (BP) local discontinuous Galerkin methods for incompressible and immiscible two-phase flows in porous media, and employ implicit pressure explicit saturation (IMPES) methods for time discretization, which is locally mass conservative for both phases. Physically, the saturations of the two phases, Sw and Sn, should belong to the range of [0, 1]. Nonphysical numerical approximations may result in instability of the simulation. Therefore, it is necessary to construct a BP technique to obtain physically relevant numerical approximations. However, the saturation does not satisfy the maximum principle, so most of the existing BP techniques cannot be applied directly. The main idea is to apply the positivity-preserving techniques to both Sw and Sn, respectively, and enforce Sw+Sn=1 simultaneously. Numerical examples are given to demonstrate the high-order accuracy of the scheme and effectiveness of the BP technique.

Publisher's Statement

© The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature 2024. Publisher’s version of record: https://doi.org/10.1007/s10915-024-02532-2">https://doi.org/10.1007/s10915-024-02532-2

Publication Title

Journal of Scientific Computing