A comparative study of explicit high-resolution schemes for compositional simulations

In this paper, compositional flow of two- and three-phase
fluids in one dimensional porous media is studied numerically and a comparison is made between several upwind and central numerical schemes.
Implicit Pressure Explicit Composition (IMPEC) procedure is used for discretization of governing equations. The pressure equation is solved implicitly while the mass conservation equations are solved explicitly using different Upwind (UPW) and Central (CEN) numerical schemes. These include Classical Upwind (UPW-CLS), Flux-based Decomposition Upwind (UPW-FLX), Variable-based Decomposition Upwind (UPW-VAR), Roe's Upwind (UPW-ROE), Local Lax Friedrichs (CEN-LLF), Dominant
Wave (CEN-DW), Harten-Lax-van Leer (HLL), and newly proposed Modied Dominant Wave (CEN-MDW)
schemes. To achieve higher resolution, high-order data generated by either MUSCL or WENO reconstructions
are utilized.
It was found that the new CEN-MDW scheme can accurately solve multiphase compositional flow equations. This scheme utilizes most of the information in
flux function while it has a moderate computational cost as a consequence of using simple algebraic formula for the wave speed approximation. Moreover, numerically calculated wave structure is shown to be used as a tool for a priori estimation of problematic regions, i.e., degenerate, umbilic, and elliptic points, which require to apply correction procedures to produce physically acceptable (entropy) solutions.
This paper is concerned with one-dimensional study of compositional two- and three-phase flows in porous media. Temperature is assumed constant and the physical model accounts for miscibility and compressibility of fluids while gravity and capillary effects are neglected.
The proposed numerical scheme can be efficiently used for solving two- and three-phase compositional flows in porous media with a low computational cost which is especially useful when the number of chemical species increases.
A new central scheme is proposed that leads to improved accuracy and computational efficiency. Moreover, to the best of authors knowledge, this is the first time that the wave structure of compositional model is investigated numerically to determine the problematic situations during numerical solution and adopt appropriate correction techniques.