Enriched Galerkin discretization for modeling poroelasticity and permeability alteration in heterogeneous porous media

JOURNAL OF COMPUTATIONAL PHYSICS

卷:427

Article

Kadeethum, T.^1,2  Nick, H. M.¹  Lee, S.³  Ballarin, F.⁴ 

[1] Tech Univ Denmark, Lyngby, Denmark

[2] Cornell Univ, Ithaca, NY 14850 USA

[3] Florida State Univ, Tallahassee, FL 32306 USA

[4] SISSA, mathLab, Math Area, Siena, Italy

关键词: Deformable porous media;    Poroelastic effects;    Biot's system;    Enriched Galerkin;    Finite element method;    Heterogeneity;   

DOI  :  10.1016/j.jcp.2020.110030

来源: Elsevier

【 摘 要 】

In this paper, we utilize the enriched Galerkin (EG) finite element method for the flow equation in Biot's system, which provides a robust locally conservative flux in heterogeneous porous media. The computational algorithm to solve the coupled system with the permeability alteration is presented with the linearization and Picard's iterative scheme. The block structure is utilized for the linear system in numerical discretization, and the computer code is shared in the open-source platform. In the numerical experiments, we compare the proposed EG method with the classical continuous Galerkin (CG) and discontinuous Galerkin (DG) finite element methods in different scenarios, including the North sea reservoirs setup. While DG and EG methods provide similar approximations for the pressure solutions, the CG method produces spurious oscillations in fluid pressure and volumetric strain solutions near the material interfaces, especially for the soft materials. The difference of flux approximation between EG and DG methods is insignificant; still, the EG method demands approximately two and three times fewer degrees of freedom than the DG method for twoand three-dimensional geometries. (C) 2020 Elsevier Inc. All rights reserved.

【 授权许可】

Free