A fully mass-conservative iterative IMPEC method for multicomponent compressible flow in porous media

Abstract

In this paper we consider efficient and fully mass-conservative numerical methods for the multicomponent compressible single-phase Darcy flow in porous media. Compared with the classical IMplicit Pressure Explicit Concentration (IMPEC) scheme by which one of the components may be not massconservative, the new scheme enjoys an appealing feature that the conservation of mass is retained for each of the components. The pressure-velocity system is obtained by the summation of the discrete conservation equation for each component multiplying an unknown parameter which is nonlinearly dependent of the molar concentrations. This approach is quite different from the conventional method which is used in the classical IMPEC scheme. We utilize a fully mass-conservative iterative IMPEC method to solve the nonlinear system for molar concentration, pressure and velocity fields. The upwind mixed finite element methods are used to solve the pressure-velocity system. Although the Peng-Robinson equation of state (EOS) is utilized to describe the pressure as a function of the molar concentrations, our method is suitable for any type of EOS. Under some reasonable conditions, the iterative scheme can be proved to be convergent, and the molar concentration of each component is positivity-preserving. Several interesting examples of multicomponent compressible flow in porous media are presented to demonstrate the robustness of the new algorithm.2010 Mathematics Subject Classification. 65N30, 49S05. Key words and phrases. Multicomponent compressible flow, Peng-Robinson equation of state, upwind mixed finite element methods, fully mass-conservative iterative IMPEC method.

Nonlinearly preconditioned constraint-preserving algorithms for subsurface three-phase flow with capillarity

“…In comparison, the semiimplicit scheme relaxes the constraint on time stepping and meanwhile achieves enhanced stability, which makes it more extensively used than the fully explicit scheme in practice. Considering that the pressure information propagates more rapidly than the saturation information, the implicit pressure and explicit saturation (IMPES) scheme [2,3,4,5,6,7,8] was proposed by solving the linear pressure equation implicitly and nonlinear saturation equations explicitly. Analogous to the IMPES method, the implicit pressure and explicit concentration (IMPEC) scheme adopts a similar strategy to solve compositional flow problems [9].…”