Comparison and application of non-conforming mesh models for flow in fractured porous media using dual Lagrange multipliers
- Zulian, Patrick ORCID Euler Institute (EUL), Università della Svizzera italiana, Switzerland
- Schädle, Philipp ORCID Institute of Geophysics, Department of Earth Sciences, ETH, Zürich, Switzerland
- Karagyaur, Liudmila ORCID Euler Institute (EUL), Università della Svizzera italiana, Switzerland
- Nestola, Maria Giuseppina Chiara ORCID Euler Institute (EUL), Università della Svizzera italiana, Switzerland
- 2022
Published in:
- Journal of Computational Physics. - 2022, vol. 449
Discrete fracture model
Non-matching meshes
Flow in fractured porous media
Finite element method
Non-conforming meshes
Lagrange multiplier method
English
Geological settings, such as reservoirs, include fractures with different material properties and geometric features. Hence, numerical simulations in applied geophysics demands for computational frameworks which efficiently allow us to integrate various fracture geometries in a porous medium matrix. This study focuses on a modeling approach for single-phase flow in fractured porous media and its application to different types of non-conforming mesh models. We propose a combination of the Lagrange multiplier method with variational transfer techniques for simulating flow through fractured porous media by employing complex non-conforming geometries as well as hybrid-and equi-dimensional models and discretizations. The variational transfer is based on the L2-projection and enables an accurate and highly efficient parallel projection of fields between non-conforming meshes (e.g., between fracture and porous matrix domain).
We present the different techniques as a unified mathematical framework with a practical perspective. By means of numerical examples we discuss both, performance and applicability of the particular strategies. Comparisons of finite element simulation results to widely adopted 2D benchmark cases show good agreement and the dual Lagrange multiplier spaces show good performance. In an extension to 3D fracture network, we first provide complementary results to a recently developed benchmark case and afterwards we explore a complex scenario which leverages the different types of fracture meshes. Complex and highly conductive fracture networks are found more suitable in combination with embedded hybrid-dimensional fractures. However, thick and blocking fractures are better approximated by equi-dimensional embedded fractures and the equi-dimensional mortar method, respectively.
We present the different techniques as a unified mathematical framework with a practical perspective. By means of numerical examples we discuss both, performance and applicability of the particular strategies. Comparisons of finite element simulation results to widely adopted 2D benchmark cases show good agreement and the dual Lagrange multiplier spaces show good performance. In an extension to 3D fracture network, we first provide complementary results to a recently developed benchmark case and afterwards we explore a complex scenario which leverages the different types of fracture meshes. Complex and highly conductive fracture networks are found more suitable in combination with embedded hybrid-dimensional fractures. However, thick and blocking fractures are better approximated by equi-dimensional embedded fractures and the equi-dimensional mortar method, respectively.
- Collections
- Language
-
- English
- Classification
- Computer science and technology
- License
- CC BY
- Open access status
- green
- Identifiers
-
- DOI 10.1016/j.jcp.2021.110773
- ARK ark:/12658/srd1322377
- Persistent URL
- https://n2t.net/ark:/12658/srd1322377
Statistics
Document views: 18
File downloads:
- Zulian_2021_Else_jcp.pdf: 30