The virtual element method for discrete fracture network simulations

Abstract In this work, an optimization based approach presented in Berrone et al. (2013, 2014) [10–12] for Discrete Fracture Network simulations is coupled with the Virtual Element Method (VEM) for the space discretization of the underlying Darcy law. The great flexibility of the VEM in handling rather general polygonal elements allows, in a natural way, for an effective description of irregular solutions starting from an arbitrary triangulation, which is built independently of the mesh on other fractures. Only partial conformity is in fact obtained with this approach. Numerical results performed on several DFN configurations confirm the viability and efficiency of the resulting method.

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

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

[3]  W. Marsden I and J , 2012 .

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

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

[6]  Wenqing Wang,et al.  Geometric modelling and object-oriented software concepts applied to a heterogeneous fractured network from the Grimsel rock laboratory , 2007 .

[7]  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..

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

[9]  M. Saunders,et al.  Solution of Sparse Indefinite Systems of Linear Equations , 1975 .

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

[11]  Franco Brezzi,et al.  Virtual Element Methods for plate bending problems , 2013 .

[12]  Jean-Raynald de Dreuzy,et al.  Hydraulic properties of two‐dimensional random fracture networks following a power law length distribution: 2. Permeability of networks based on lognormal distribution of apertures , 2001 .

[13]  T. Belytschko,et al.  The extended/generalized finite element method: An overview of the method and its applications , 2010 .

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

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

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

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

[18]  G. Marsily,et al.  Modeling fracture flow with a stochastic discrete fracture network: calibration and validation: 1. The flow model , 1990 .