Numerical methods for a model for compressible miscible displacement in porous media

Abstract

A nonlinear parabolic system is derived to describe compressible miscible displacement in a porous medium. The system is consistent with the usual model for incompressible miscible displacement. Two finite element procedures are introduced to approximate the concentration of one of the fluids and the pressure of the mixture. The concentration is treated by a Galerkin method in both procedures, while the pressure is treated by either a Galerkin method or by a parabolic mixed finite element method. Optimal order estimates in <inline-formula content-type="math/mathml"> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" alttext="upper L squared"> <mml:semantics> <mml:mrow class="MJX-TeXAtom-ORD"> <mml:msup> <mml:mi>L</mml:mi> <mml:mn>2</mml:mn> </mml:msup> </mml:mrow> <mml:annotation encoding="application/x-tex">{L^2}</mml:annotation> </mml:semantics> </mml:math> </inline-formula> and essentially optimal order estimates in <inline-formula content-type="math/mathml"> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" alttext="upper L Superscript normal infinity"> <mml:semantics> <mml:mrow class="MJX-TeXAtom-ORD"> <mml:msup> <mml:mi>L</mml:mi> <mml:mi mathvariant="normal"> ∞ </mml:mi> </mml:msup> </mml:mrow> <mml:annotation encoding="application/x-tex">{L^\infty }</mml:annotation> </mml:semantics> </mml:math> </inline-formula> are derived for the errors in the approximate solutions for both methods.