Polymer flooding is an important stage of enhanced oil recovery in petroleum reservoir engineering. A model of this process is based on the study of multicomponent viscous flow in porous media with adsorption. This model can be expressed as a Brinkman-based model of flow in porous media coupled to a non-strictly hyperbolic system of conservation laws for multiple components (water and polymers) that form the aqueous phase. The discretisation proposed for this coupled flow-transport problem combines an H(div)-conforming discontinuous Galerkin (DG) method for the Brinkman flow problem with a classical DG method for the transport equations. The DG formulation of the transport problem is based on discontinuous numerical fluxes. An invariant region property is proved under the (mild) assumption that the underlying mesh is a B-triangulation [B. Cockburn, S. Hou, and C.-W. Shu, Math. Comp., 54 (1990), pp. 545-581]. This property states that only physically relevant (bounded and non-negative) saturation and concentration values are generated by the scheme. Numerical tests illustrate the accuracy and stability of the proposed method.

No abstract

No abstract

No abstract