Compositional modeling by the combined discontinuous Galerkin and mixed methods

H. Hoteit*, A. Firoozabadi

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

119 Scopus citations


In this work, we present a numerical procedure that combines the mixed finite-element (MFE) and the discontinuous Galerkin (DG) methods. This numerical scheme is used to solve the highly nonlinear coupled equations that describe the flow processes in homogeneous and heterogeneous media with mass transfer between the phases. The MFE method is used to approximate the phase velocity based on the pressure (more precisely average pressure) at the interface between the nodes. This approach conserves the mass locally at the element level and guarantees the continuity of the total flux across the interfaces. The DG method is used to solve the mass-balance equations, which are generally convection-dominated. The DG method associated with suitable slope limiters can capture sharp gradients in the solution without creating spurious oscillations. We present several numerical examples in homogeneous and heterogeneous media that demonstrate the superiority of our method to the finite-difference (FD) approach. Our proposed MFE-DG method becomes orders of magnitude faster than the FD method for a desired accuracy in 2D.

Original languageEnglish (US)
Pages (from-to)19-34
Number of pages16
JournalSPE Journal
Issue number1
StatePublished - Mar 2006
Externally publishedYes

ASJC Scopus subject areas

  • Energy Engineering and Power Technology
  • Geotechnical Engineering and Engineering Geology


Dive into the research topics of 'Compositional modeling by the combined discontinuous Galerkin and mixed methods'. Together they form a unique fingerprint.

Cite this