A hybrid numerical scheme for multi-phase flow in heterogeneous porous media

This work describes a numerical scheme for simulation of two-phase flow in heterogeneous porous media. The system of governing equations is transformed into an equivalent problem, as proposed by Peaceman (1977). For the numerical solution, the operator split technique is used, separating the advection-diffusion equation in a hyperbolic and a parabolic part. The Finite Element and Finite Volume Methods are used for the spatial discretization of the problem. The reliability of the scheme is demonstrated for a typical problem of Petroleum Engineering, showing better performance than a commercial software.