Multiscale finite-element method for linear elastic geomechanics

N Castelletto, Hadi Hajibeygi, HA Tchelepi

Research output: Contribution to journalArticleScientificpeer-review

27 Citations (Scopus)
30 Downloads (Pure)


The demand for accurate and efficient simulation of geomechanical effects is widely increasing in the geoscience community. High resolution characterizations of the mechanical properties of subsurface formations are essential for improving modeling predictions. Such detailed descriptions impose severe computational challenges and motivate the development of multiscale solution strategies. We propose a multiscale solution framework for the geomechanical equilibrium problem of heterogeneous porous media based on the finite-element method. After imposing a coarsescale grid on the given fine-scale problem, the coarse-scale basis functions are obtained by solving local equilibrium problems within coarse elements. These basis functions form the restriction and prolongation operators used to obtain the coarse-scale system for the displacement-vector. Then, a two-stage preconditioner that couples the multiscale system with a smoother is derived for the iterative solution of the fine-scale linear system. Various numerical experiments are presented to demonstrate accuracy and robustness of the method.
Original languageEnglish
Pages (from-to)337-356
Number of pages30
JournalJournal of Computational Physics
Publication statusPublished - 2016


  • Multiscale methods
  • multiscale finite-element method
  • geomechanics
  • reconditioning
  • porous media


Dive into the research topics of 'Multiscale finite-element method for linear elastic geomechanics'. Together they form a unique fingerprint.

Cite this