Mixed Discontinuous Finite Element Methods for Multiphase Flow in Porous Media

In this paper we present an application of mixed discontinuous finite element methods to the simulation of two phase immiscible flow in porous media. The partial differential system describing this flow is written in terms of an elliptic equation for a global pressure and a parabolic equation for a saturation. Both of these two equations are solved using the mixed discontinuous finite element methods. Numerical results are presented for P0, P1, and P2 discontinuous elements. The P0 results are very close to thoseb generated by the standard finite difference method, while the P1 and P2 results seem more accurate.