Numerical simulation for a two-phase porous medium flow problem with rate independent hysteresis

The paper is devoted to the numerical simulation of a multiphase flow in porous medium with a hysteretic relation between the capillary pressures and the saturations of the phases. The flow model we use is based on Darcy's law. The hysteretic relation between the capillary pressures and the saturations is described by a play-type hysteresis operator. We propose a numerical algorithm for treating the arising system of equations, discuss finite element schemes and present simulation results for the case of two phases.