Modeling of Multicomponent Diffusions and Natural Convection in Unfractured and Fractured Media by Discontinuous Galerkin and Mixed Methods

Hussein Hoteit, Abbas Firoozabadi

Research output: Contribution to journalArticlepeer-review

4 Scopus citations

Abstract

Computation of the distribution of species in hydrocarbon reservoirs from diffusions (thermal, molecular, and pressure) and natural convection is an important step in reservoir initialization. Current methods, which are mainly based on the conventional finite difference approach, may not be numerically efficient in fractured and other media with complex heterogeneities. In this work, the discontinuous Galerkin (DG) method combined with the mixed finite element (MFE) method is used for the calculation of compositional variation in fractured hydrocarbon reservoirs. The use of unstructured gridding allows efficient computations for fractured media when the crossflow equilibrium concept is invoked. The DG method has less numerical dispersion than the upwind finite difference (FD) methods. The MFE method ensures continuity of fluxes at the interface of the grid elements. We also use the local discontinuous Galerkin (LDG) method instead of the MFE calculate the diffusion fluxes. Results from several numerical examples are presented to demonstrate the efficiency, robustness, and accuracy of the model. Various features of convection and diffusion in homogeneous, layered, and fractured media are also discussed.
Original languageEnglish (US)
Pages (from-to)535-556
Number of pages22
JournalInternational Journal for Numerical Methods in Engineering
Volume114
Issue number5
DOIs
StatePublished - Feb 21 2018

Bibliographical note

KAUST Repository Item: Exported on 2020-10-01

Fingerprint Dive into the research topics of 'Modeling of Multicomponent Diffusions and Natural Convection in Unfractured and Fractured Media by Discontinuous Galerkin and Mixed Methods'. Together they form a unique fingerprint.

Cite this