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 two- and three-dimensional geometries.

Enriched Galerkin discretization for modeling poroelasticity and permeability alteration in heterogeneous porous media / Kadeethum, T.; Nick, H. M.; Lee, S.; Ballarin, F.. - In: JOURNAL OF COMPUTATIONAL PHYSICS. - ISSN 0021-9991. - 427:(2021). [10.1016/j.jcp.2020.110030]

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

Ballarin, F.
2021-01-01

Abstract

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 two- and three-dimensional geometries.
2021
427
110030
https://arxiv.org/abs/2010.06653
Kadeethum, T.; Nick, H. M.; Lee, S.; Ballarin, F.
File in questo prodotto:
File Dimensione Formato  
2010.06653.pdf

non disponibili

Descrizione: arXiv:2010.06653
Tipologia: Documento in Post-print
Licenza: Non specificato
Dimensione 6.98 MB
Formato Adobe PDF
6.98 MB Adobe PDF   Visualizza/Apri   Richiedi una copia

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/20.500.11767/117775
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 12
  • ???jsp.display-item.citation.isi??? 9
social impact