Multiscale finite-element method for linear elastic geomechanics

JOURNAL OF COMPUTATIONAL PHYSICS

卷:331

Article

Castelletto, Nicola¹  Hajibeygi, Hadi²  Tchelepi, Hamdi A.¹ 

[1] Stanford Univ, Dept Energy Resources Engn, Stanford, CA 94305 USA

[2] Delft Univ Technol, Fac Civil Engn & Geosci, Dept Geosci & Engn, Delft, Netherlands

关键词: Multiscale methods;    Multiscale finite-element method;    Geomechanics;    Preconditioning;    Porous media;   

DOI  :  10.1016/j.jcp.2016.11.044

来源: Elsevier

【 摘 要 】

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 coarse-scale 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. (C) 2016 Elsevier Inc. All rights reserved.

【 授权许可】

Free