Thermo-hydrodynamic-mechanical multiphase flow model for CO2 sequestration in fracturing porous media

Mehdi Musivand Arzanfudi*, Rafid Al-Khoury

*Corresponding author for this work

Research output: Contribution to journalArticleScientificpeer-review

9 Citations (Scopus)


In this paper, we introduce a fully coupled thermo-hydrodynamic-mechanical computational model for multiphase flow in a deformable porous solid, exhibiting crack propagation due to fluid dynamics, with focus on CO2 geosequestration. The geometry is described by a matrix domain, a fracture domain, and a matrix-fracture domain. The fluid flow in the matrix domain is governed by Darcy's law and that in the crack is governed by the Navier-Stokes equations. At the matrix-fracture domain, the fluid flow is governed by a leakage term derived from Darcy's law. Upon crack propagation, the conservation of mass and energy of the crack fluid is constrained by the isentropic process. We utilize the representative elementary volume-averaging theory to formulate the mathematical model of the porous matrix, and the drift flux model to formulate the fluid dynamics in the fracture. The numerical solution is conducted using a mixed finite element discretization scheme. The standard Galerkin finite element method is utilized to discretize the diffusive dominant field equations, and the extended finite element method is utilized to discretize the crack propagation, and the fluid leakage at the boundaries between layers of different physical properties. A numerical example is given to demonstrate the computational capability of the model. It shows that the model, despite the relatively large number of degrees of freedom of different physical nature per node, is computationally efficient, and geometry and effectively mesh independent.

Original languageEnglish
Pages (from-to)635–674
JournalInternational Journal for Numerical Methods in Fluids
Issue number11
Publication statusPublished - 2017


  • Deformable porous media
  • Fluid dynamics
  • Fracturing porous media
  • Mixed discretization finite element scheme
  • Multiphase flow
  • Phase change


Dive into the research topics of 'Thermo-hydrodynamic-mechanical multiphase flow model for CO2 sequestration in fracturing porous media'. Together they form a unique fingerprint.

Cite this