Multiscale finite volume method for finite-volume-based simulation of poroelasticity

Irina Sokolova, Muhammad Gusti Bastisya, Hadi Hajibeygi

Research output: Contribution to journalArticleScientificpeer-review

23 Citations (Scopus)
1 Downloads (Pure)


We propose a multiscale finite volume method (MSFV) for simulation of coupled flow-deformation in heterogeneous porous media under elastic deformation (i.e., poroelastic model). The fine-scale fully resolved system of equations is obtained based on a conservative finite-volume method in which the displacement and pore pressure unknowns are located in a staggered configuration. The coupling is treated through a fully-coupled fully-implicit formulation. On this fully-coupled finite-volume system, coarse-scale grids for flow and deformation are imposed. Local basis functions for scalar pore pressure and vectorial displacement unknowns are then solved over their respective local domains at the beginning of the simulation, and reused for the rest of the time-dependent simulations. These local basis functions are then clustered to form the prolongation operator. As for the finite-volume nature of the proposed multiscale system, finite-volume restriction operators for poroelastic systems are utilised. Once the coarse-scale system is solved, its solution is prolonged back to the original fine-scale resolution, providing approximate fine-scale solution. The finite-volume multiscale formulation provides conservative stress and mass flux both at fine and coarse scale. Several numerical test cases are provided first to validate the fine-scale finite-volume discrete fully-implicit simulation, and then to investigate the accuracy of the proposed multiscale formulation. Moreover we compare our fully implicit MSFV method with hybrid multiscale Finite Element-Finite Volume (h-MSFE-FV). Our multiscale method allows for quantification of the elastic geomechanical behaviour with using only a fraction of the fine-scale grid cells, even for highly heterogeneous time-dependent models. As such, it casts a promising approach for field-scale quantification of the mechanical deformation and stress field due to injection and production in a subsurface formation.

Original languageEnglish
Pages (from-to)309-324
Number of pages16
JournalJournal of Computational Physics
Publication statusPublished - 2019


  • Algebraic multiscale solver
  • Finite volume method
  • Geomechanics
  • Multiscale finite-volume method
  • Poromechanics
  • Porous media


Dive into the research topics of 'Multiscale finite volume method for finite-volume-based simulation of poroelasticity'. Together they form a unique fingerprint.

Cite this