The effect of numerical integration in mixed finite element approximation in the simulation of miscible displacement

We consider the effect of numerical integration in finite element procedures applied to a nonlinear system of two coupled partial differential equations describing the miscible displacement of one incompressible fluid by another in a porous meduim. We consider the use of the numerical quadrature scheme for approximating the pressure and velocity by a mixed method using Raviart - Thomas space of index and the concentration by a standard Galerkin method. We also give some sufficient conditions on the quadrature scheme to ensure that the order of convergence is unaltered in the presence of numerical integration. Optimal order estimates are derived when the imposed external flows are smoothly distributed.