Analysis of incompressible miscible displacement in porous media by characteristics collocation method

Miscible displacement of one incompressible fluid by another in a porous medium is modelled by a coupled system of two partial differential equations. The pressure equation is elliptic, whereas the concentration equation is parabolic but normally convection-dominated. In this article, the collocation scheme is used to approximate the pressure equation and another characteristics collocation scheme to treat concentration equation. Existence and uniqueness of solutions of the algorithm are proved. Optimal order error estimate is demonstrated. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006