High-order bound-preserving finite difference methods for incompressible two-phase flow in porous media

Document Type

Article

Publication Date

12-1-2025

Department

Department of Mathematical Sciences

Abstract

In this paper, we develop high-order bound-preserving (BP) finite difference (FD) methods for solving the incompressible and immiscible two-phase flow problem with capillary pressure in porous media. We use the implicit pressure explicit saturation (IMPES) scheme to solve for the pressure, auxiliary variables, and saturations of each phase in the coupled system. The boundedness of the saturations of the two phases, Sw and Sn, between 0 and 1 is an important physical characteristic. Applying non-physical numerical approximations may lead to significant oscillations in the numerical results and cause instability in the simulation. We apply high-order FD method and BP technique to maintain the high-order accuracy and the boundary of saturations. In the BP technique, the main idea is to choose an appropriate time step and apply positivity-preserving (PP) technique to Sw and Sn, respectively, and ensure that Sw+Sn=1. In addition, the high-order accuracy is obtained by the parameterized flux limiter. Numerical examples are presented to demonstrate the high-order accuracy of the scheme and the effectiveness of the BP technique.

Publication Title

Journal of Computational and Applied Mathematics

Share

COinS