The Virtual Element Method for Underground Flow Simulations in Fractured Media

We focus on the problem of performing underground flow simulations in fractured media. The medium is modelled by means of the so-called Discrete Fracture Network (DFN) model. Within this framework, we discuss about the use of the Virtual Element Method (VEM) in performing simulations, and about its role in facilitating the meshing process. DFN models are characterized by a large number of planar fractures stochastically generated and placed in a 3D reference system starting from given probabilistic distribution of space position, dimension, aspect ratio, orientation and hydrogeological properties. For these reasons, realistic DFNs are usually characterized by the presence of very complex geometrical configurations on which it is very difficult to obtain a good quality Finite Element mesh. In particular, in order to circumvent these difficulties, we consider two different VEM approaches: an approach based on a totally conforming polygonal mesh obtained starting from a triangular mesh independently generated on the different fractures of the system and an approach based on a partially conforming mesh in which the connection between the fracture meshes is imposed through a mortaring approach borrowed from the classical domain decomposition methods. Concerning these two approaches, we discuss their relative advantages and present numerical results aimed at comparing them.

[1]  Jean-Raynald de Dreuzy,et al.  A Generalized Mixed Hybrid Mortar Method for Solving Flow in Stochastic Discrete Fracture Networks , 2012, SIAM J. Sci. Comput..

[2]  Stefano Berrone,et al.  A Posteriori Error Estimate for a PDE-Constrained Optimization Formulation for the Flow in DFNs , 2016, SIAM J. Numer. Anal..

[3]  L. Beirao da Veiga,et al.  Mixed Virtual Element Methods for general second order elliptic problems on polygonal meshes , 2014 .

[4]  Stefano Berrone,et al.  A globally conforming method for solving flow in discrete fracture networks using the Virtual Element Method , 2016 .

[5]  Barbara Wohlmuth.  Discretization Techniques Based on Domain Decomposition , 2001 .

[6]  Franco Brezzi,et al.  The Hitchhiker's Guide to the Virtual Element Method , 2014 .

[7]  J. Erhel,et al.  A mixed hybrid Mortar method for solving flow in discrete fracture networks , 2010 .

[8]  Stefano Berrone,et al.  A Parallel Solver for Large Scale DFN Flow Simulations , 2015, SIAM J. Sci. Comput..

[9]  Jeffrey D. Hyman,et al.  Conforming Delaunay Triangulation of Stochastically Generated Three Dimensional Discrete Fracture Networks: A Feature Rejection Algorithm for Meshing Strategy , 2014, SIAM J. Sci. Comput..

[10]  Claudio Canuto,et al.  Uncertainty quantification in Discrete Fracture Network models: Stochastic fracture transmissivity , 2015, Comput. Math. Appl..

[11]  Ahmed Alsaedi,et al.  Equivalent projectors for virtual element methods , 2013, Comput. Math. Appl..

[12]  J. Jaffré,et al.  Modeling flow in porous media with fractures; Discrete fracture models with matrix-fracture exchange , 2012 .

[13]  Stefano Berrone,et al.  Towards effective flow simulations in realistic discrete fracture networks , 2016, J. Comput. Phys..

[14]  Stefano Berrone,et al.  An optimization approach for large scale simulations of discrete fracture network flows , 2014, J. Comput. Phys..

[15]  C. Bernardi,et al.  A New Nonconforming Approach to Domain Decomposition : The Mortar Element Method , 1994 .

[16]  Martin Vohralík,et al.  Mixed and nonconforming finite element methods on a system of polygons , 2007 .

[17]  Lourenço Beirão da Veiga,et al.  Virtual Elements for Linear Elasticity Problems , 2013, SIAM J. Numer. Anal..

[18]  F. Brezzi,et al.  Basic principles of Virtual Element Methods , 2013 .

[19]  Richard S. Falk,et al.  Basic principles of mixed Virtual Element Methods , 2014 .

[20]  J. Erhel,et al.  A Mortar BDD Method for Solving Flow in Stochastic Discrete Fracture Networks , 2014 .

[21]  Jean-Raynald de Dreuzy,et al.  Synthetic benchmark for modeling flow in 3D fractured media , 2013, Comput. Geosci..

[22]  C. Fidelibus,et al.  Derivation of equivalent pipe network analogues for three‐dimensional discrete fracture networks by the boundary element method , 1999 .

[23]  Alessandro Russo,et al.  Mixed Virtual Element Methods for general second order elliptic problems on polygonal meshes , 2014, 1506.07328.

[24]  Giovanni Battista Barla,et al.  The Hydro-Mechanically Coupled Response of Rock Fractures , 2007 .

[25]  Stefano Berrone,et al.  A hybrid mortar virtual element method for discrete fracture network simulations , 2016, J. Comput. Phys..

[26]  Benoit Noetinger,et al.  A quasi steady state method for solving transient Darcy flow in complex 3D fractured networks , 2012, J. Comput. Phys..

[27]  F. Brezzi.  On the existence, uniqueness and approximation of saddle-point problems arising from lagrangian multipliers , 1974 .

[28]  Pierre M. Adler,et al.  Fractures and Fracture Networks , 1999 .

[29]  Alessio Fumagalli,et al.  A reduced model for Darcy’s problem in networks of fractures , 2014 .

[30]  Stefano Berrone,et al.  A PDE-Constrained Optimization Formulation for Discrete Fracture Network Flows , 2013, SIAM J. Sci. Comput..

[31]  L. Beirao da Veiga,et al.  A Virtual Element Method for elastic and inelastic problems on polytope meshes , 2015, 1503.02042.

[32]  Axel Klawonn,et al.  FETI Domain Decomposition Methods for Second Order Elliptic Partial Differential Equations , 2006 .

[33]  Stephen J. Wright,et al.  Numerical Optimization , 2018, Fundamental Statistical Inference.

[34]  Francesca Rapetti,et al.  Basics and some applications of the mortar element method , 2005 .

[35]  Barbara I. Wohlmuth,et al.  A Mortar Finite Element Method Using Dual Spaces for the Lagrange Multiplier , 2000, SIAM J. Numer. Anal..

[36]  B. Nœtinger.  A quasi steady state method for solving transient Darcy flow in complex 3D fractured networks accounting for matrix to fracture flow , 2015 .

[37]  Stefano Berrone,et al.  On Simulations of Discrete Fracture Network Flows with an Optimization-Based Extended Finite Element Method , 2013, SIAM J. Sci. Comput..

[38]  Sandra Pieraccini,et al.  On a PDE-Constrained Optimization Approach for Flow Simulations in Fractured Media , 2016 .

[39]  P. Raviart,et al.  Primal hybrid finite element methods for 2nd order elliptic equations , 1977 .

[40]  Stefano Berrone,et al.  The virtual element method for discrete fracture network simulations , 2014 .

[41]  Anna Scotti,et al.  Optimal Techniques to Simulate Flow in Fractured Reservoir , 2014 .

[42]  Faker Ben Belgacem,et al.  The Mortar finite element method with Lagrange multipliers , 1999, Numerische Mathematik.

[43]  Alessio Fumagalli,et al.  A numerical method for two-phase flow in fractured porous media with non-matching grids , 2013 .