Showing 1–38 of 38 results for author: Fumagalli, A

Two Nitsche-based mixed finite element discretizations for the seepage problem in Richards' equation

Authors: Federico Gatti, Andrea Bressan, Alessio Fumagalli, Domenico Gallipoli, Leonardo Maria Lalicata, Simone Pittaluga, Lorenzo Tamellini

Abstract: This paper proposes two algorithms to impose seepage boundary conditions in the context of Richards' equation for groundwater flows in unsaturated media. Seepage conditions are non-linear boundary conditions, that can be formulated as a set of unilateral constraints on both the pressure head and the water flux at the ground surface, together with a complementarity condition: these conditions in pr… ▽ More This paper proposes two algorithms to impose seepage boundary conditions in the context of Richards' equation for groundwater flows in unsaturated media. Seepage conditions are non-linear boundary conditions, that can be formulated as a set of unilateral constraints on both the pressure head and the water flux at the ground surface, together with a complementarity condition: these conditions in practice require switching between Neumann and Dirichlet boundary conditions on unknown portions on the boundary. Upon realizing the similarities of these conditions with unilateral contact problems in mechanics, we take inspiration from that literature to propose two approaches: the first method relies on a strongly consistent penalization term, whereas the second one is obtained by an hybridization approach, in which the value of the pressure on the surface is treated as a separate set of unknowns. The flow problem is discretized in mixed form with div-conforming elements so that the water mass is preserved. Numerical experiments show the validity of the proposed strategy in handling the seepage boundary conditions on geometries with increasing complexity. △ Less

Submitted 10 July, 2024; originally announced July 2024.

A hybrid upwind scheme for two-phase flow in fractured porous media

Authors: Enrico Ballini, Luca Formaggia, Alessio Fumagalli, Eirik Keilegavlen, Anna Scotti

Abstract: Simulating the flow of two fluid phases in porous media is a challenging task, especially when fractures are included in the simulation. Fractures may have highly heterogeneous properties compared to the surrounding rock matrix, significantly affecting fluid flow, and at the same time hydraulic aperture that are much smaller than any other characteristic sizes in the domain. Generally, flow simula… ▽ More Simulating the flow of two fluid phases in porous media is a challenging task, especially when fractures are included in the simulation. Fractures may have highly heterogeneous properties compared to the surrounding rock matrix, significantly affecting fluid flow, and at the same time hydraulic aperture that are much smaller than any other characteristic sizes in the domain. Generally, flow simulators face difficulties with counter-current flow, generated by gravity and pressure gradients, which hinders the convergence of non-linear solvers (Newton). In this work, we model the fracture geometry with a mixed-dimensional discrete fracture network, thus lightening the computational burden associated to an equi-dimensional representation. We address the issue of counter-current flows with appropriate spatial discretization of the advective fluid fluxes, with the aim of improving the convergence speed of the non-linear solver. In particular, the extension of the hybrid upwinding to the mixed-dimensional framework, with the use of a phase potential upstreaming at the interfaces of subdomains. We test the method across several cases with different flow regimes and fracture network geometry. Results show robustness of the chosen discretization and a consistent improvements, in terms of Newton iterations, compared to use the phase potential upstreaming everywhere. △ Less

Submitted 6 July, 2024; originally announced July 2024.

Numerical validation of an adaptive model for the determination of nonlinear-flow regions in highly heterogeneous porous media

Authors: Alessio Fumagalli, Francesco S. Patacchini

Abstract: An adaptive model for the description of flows in highly heterogeneous porous media is developed in~\cite{FP21,FP23}. There, depending on the magnitude of the fluid's velocity, the constitutive law linking velocity and pressure gradient is selected between two possible options, one better adapted to slow motion and the other to fast motion. We propose here to validate further this adaptive approac… ▽ More An adaptive model for the description of flows in highly heterogeneous porous media is developed in~\cite{FP21,FP23}. There, depending on the magnitude of the fluid's velocity, the constitutive law linking velocity and pressure gradient is selected between two possible options, one better adapted to slow motion and the other to fast motion. We propose here to validate further this adaptive approach by means of more extensive numerical experiments, including a three-dimensional case, as well as to use such approach to determine a partition of the domain into slow- and fast-flow regions. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Application of Deep Learning Reduced-Order Modeling for Single-Phase Flow in Faulted Porous Media

Authors: Enrico Ballini, Luca Formaggia, Alessio Fumagalli, Anna Scotti, Paolo Zunino

Abstract: We apply reduced-order modeling (ROM) techniques to single-phase flow in faulted porous media, accounting for changing rock properties and fault geometry variations using a radial basis function mesh deformation method. This approach benefits from a mixed-dimensional framework that effectively manages the resulting non-conforming mesh. To streamline complex and repetitive calculations such as sens… ▽ More We apply reduced-order modeling (ROM) techniques to single-phase flow in faulted porous media, accounting for changing rock properties and fault geometry variations using a radial basis function mesh deformation method. This approach benefits from a mixed-dimensional framework that effectively manages the resulting non-conforming mesh. To streamline complex and repetitive calculations such as sensitivity analysis and solution of inverse problems, we utilize the Deep Learning Reduced Order Model (DL-ROM). This non-intrusive neural network-based technique is evaluated against the traditional Proper Orthogonal Decomposition (POD) method across various scenarios, demonstrating DL-ROM's capacity to expedite complex analyses with promising accuracy and efficiency. △ Less

Submitted 6 March, 2024; originally announced March 2024.

Deep learning based reduced order modeling of Darcy flow systems with local mass conservation

Authors: Wietse M. Boon, Nicola R. Franco, Alessio Fumagalli, Paolo Zunino

Abstract: We propose a new reduced order modeling strategy for tackling parametrized Partial Differential Equations (PDEs) with linear constraints, in particular Darcy flow systems in which the constraint is given by mass conservation. Our approach employs classical neural network architectures and supervised learning, but it is constructed in such a way that the resulting Reduced Order Model (ROM) is guara… ▽ More We propose a new reduced order modeling strategy for tackling parametrized Partial Differential Equations (PDEs) with linear constraints, in particular Darcy flow systems in which the constraint is given by mass conservation. Our approach employs classical neural network architectures and supervised learning, but it is constructed in such a way that the resulting Reduced Order Model (ROM) is guaranteed to satisfy the linear constraints exactly. The procedure is based on a splitting of the PDE solution into a particular solution satisfying the constraint and a homogenous solution. The homogeneous solution is approximated by mapping a suitable potential function, generated by a neural network model, onto the kernel of the constraint operator; for the particular solution, instead, we propose an efficient spanning tree algorithm. Starting from this paradigm, we present three approaches that follow this methodology, obtained by exploring different choices of the potential spaces: from empirical ones, derived via Proper Orthogonal Decomposition (POD), to more abstract ones based on differential complexes. All proposed approaches combine computational efficiency with rigorous mathematical interpretation, thus guaranteeing the explainability of the model outputs. To demonstrate the efficacy of the proposed strategies and to emphasize their advantages over vanilla black-box approaches, we present a series of numerical experiments on fluid flows in porous media, ranging from mixed-dimensional problems to nonlinear systems. This research lays the foundation for further exploration and development in the realm of model order reduction, potentially unlocking new capabilities and solutions in computational geosciences and beyond. △ Less

Submitted 24 November, 2023; originally announced November 2023.

A mixed-dimensional model for direct current simulations in presence of a thin high-resistivity liner

Authors: Alessio Fumagalli, Lorenzo Panzeri, Luca Formaggia, Anna Scotti, Diego Arosio

Abstract: In this work we present a mixed-dimensional mathematical model to obtain the electric potential and current density in direct current simulations when a thin liner is included in the modelled domain. The liner is used in landfill management to prevent leakage of leachate from the waste body into the underground and is made of a highly-impermeable high-resistivity plastic material. The electrodes a… ▽ More In this work we present a mixed-dimensional mathematical model to obtain the electric potential and current density in direct current simulations when a thin liner is included in the modelled domain. The liner is used in landfill management to prevent leakage of leachate from the waste body into the underground and is made of a highly-impermeable high-resistivity plastic material. The electrodes and the liner have diameters and thickness respectively that are much smaller than their other dimensions, thus their numerical simulation might be too costly in an equi-dimensional setting. Our approach is to approximate them as objects of lower dimension and derive the corresponding equations. The obtained mixed-dimensional model is validated against laboratory experiments of increasing complexity showing the reliability of the proposed mathematical model. △ Less

Submitted 22 March, 2023; originally announced March 2023.

Mixed and multipoint finite element methods for rotation-based poroelasticity

Authors: Wietse M. Boon, Alessio Fumagalli, Anna Scotti

Abstract: This work proposes a mixed finite element method for the Biot poroelasticity equations that employs the lowest-order Raviart-Thomas finite element space for the solid displacement and piecewise constants for the fluid pressure. The method is based on the formulation of linearized elasticity as a weighted vector Laplace problem. By introducing the solid rotation and fluid flux as auxiliary variable… ▽ More This work proposes a mixed finite element method for the Biot poroelasticity equations that employs the lowest-order Raviart-Thomas finite element space for the solid displacement and piecewise constants for the fluid pressure. The method is based on the formulation of linearized elasticity as a weighted vector Laplace problem. By introducing the solid rotation and fluid flux as auxiliary variables, we form a four-field formulation of the Biot system, which is discretized using conforming mixed finite element spaces. The auxiliary variables are subsequently removed from the system in a local hybridization technique to obtain a multipoint rotation-flux mixed finite element method. Stability and convergence of the four-field and multipoint mixed finite element methods are shown in terms of weighted norms, which additionally leads to parameter-robust preconditioners. Numerical experiments confirm the theoretical results. △ Less

Submitted 23 December, 2022; originally announced December 2022.

MSC Class: 65N12; 65N22; 65N30

Mixed Virtual Element approximation of linear acoustic wave equation

Authors: Franco Dassi, Alessio Fumagalli, Ilario Mazzieri, Giuseppe Vacca

Abstract: We design a Mixed Virtual Element Method for the approximated solution to the first-order form of the acoustic wave equation. In absence of external load, the semi-discrete method exactly conserves the system energy. To integrate in time the semi-discrete problem we consider a classical theta-method scheme. We carry out the stability and convergence analysis in the energy norm for the semi-discret… ▽ More We design a Mixed Virtual Element Method for the approximated solution to the first-order form of the acoustic wave equation. In absence of external load, the semi-discrete method exactly conserves the system energy. To integrate in time the semi-discrete problem we consider a classical theta-method scheme. We carry out the stability and convergence analysis in the energy norm for the semi-discrete problem showing optimal rate of convergence with respect to the mesh size. We further study the property of energy conservation for the fully-discrete system. Finally, we present some verification tests as well as engineering application of the method. △ Less

Submitted 23 September, 2022; originally announced September 2022.

A multipoint vorticity mixed finite element method for incompressible Stokes flow

Authors: Wietse M. Boon, Alessio Fumagalli

Abstract: We propose a mixed finite element method for Stokes flow with one degree of freedom per element and facet of simplicial grids. The method is derived by considering the vorticity-velocity-pressure formulation and eliminating the vorticity locally through the use of a quadrature rule. The discrete solution is pointwise divergence-free and the method is pressure robust. The theoretically derived conv… ▽ More We propose a mixed finite element method for Stokes flow with one degree of freedom per element and facet of simplicial grids. The method is derived by considering the vorticity-velocity-pressure formulation and eliminating the vorticity locally through the use of a quadrature rule. The discrete solution is pointwise divergence-free and the method is pressure robust. The theoretically derived convergence rates are confirmed by numerical experiments. △ Less

Submitted 29 August, 2022; originally announced August 2022.

Uncertainty quantification for mineral precipitation and dissolution in fractured porous media

Authors: Michele Botti, Alessio Fumagalli, Anna Scotti

Abstract: In this work, we present an uncertainty quantification analysis to determine the influence and importance of some physical parameters in a reactive transport model in fractured porous media. An accurate description of flow and transport in the fractures is key to obtain reliable simulations, however, fractures geometry and physical characteristics pose several challenges from both the modeling and… ▽ More In this work, we present an uncertainty quantification analysis to determine the influence and importance of some physical parameters in a reactive transport model in fractured porous media. An accurate description of flow and transport in the fractures is key to obtain reliable simulations, however, fractures geometry and physical characteristics pose several challenges from both the modeling and implementation side. We adopt a mixed-dimensional approximation, where fractures and their intersections are represented as objects of lower dimension. To simplify the presentation, we consider only two chemical species: one solute, transported by water, and one precipitate attached to the solid skeleton. A global sensitivity analysis to uncertain input data is performed exploiting the Polynomial Chaos expansion along with spectral projection methods on sparse grids. △ Less

Submitted 13 July, 2022; originally announced July 2022.

Well-posedness and variational numerical scheme for an adaptive model in highly heterogeneous porous media

Authors: Alessio Fumagalli, Francesco Saverio Patacchini

Abstract: Mathematical modeling of fluid flow in a porous medium is usually described by a continuity equation and a chosen constitutive law. The latter, depending on the problem at hand, may be a nonlinear relation between the fluid's pressure gradient and velocity. The actual shape of this relation is normally chosen at the outset of the problem, even though, in practice, the fluid may experience velociti… ▽ More Mathematical modeling of fluid flow in a porous medium is usually described by a continuity equation and a chosen constitutive law. The latter, depending on the problem at hand, may be a nonlinear relation between the fluid's pressure gradient and velocity. The actual shape of this relation is normally chosen at the outset of the problem, even though, in practice, the fluid may experience velocities outside of its range of applicability. We propose here an adaptive model, so that the most appropriate law is locally selected depending on the computed velocity. From the analytical point of view, we show well-posedness of the problem when the law is monotone in velocity and show existence in one space dimension otherwise. From the computational point of view, we present a new approach based on regularizing via mollification the underlying dissipation, i.e., the power lost by the fluid to the porous medium through drag. The resulting regularization is shown to converge to the original problem using $Γ$-convergence on the dissipation in the monotone case. This approach gives rise to a variational numerical scheme which applies to very general problems and which we validate on three test cases. △ Less

Submitted 16 June, 2022; originally announced June 2022.

A Reduced Basis Method for Darcy flow systems that ensures local mass conservation by using exact discrete complexes

Authors: Wietse M. Boon, Alessio Fumagalli

Abstract: A solution technique is proposed for flows in porous media that guarantees local conservation of mass. We first compute a flux field to balance the mass source and then exploit exact co-chain complexes to generate a solenoidal correction. A reduced basis method based on proper orthogonal decomposition is employed to construct the correction and we show that mass balance is ensured regardless of th… ▽ More A solution technique is proposed for flows in porous media that guarantees local conservation of mass. We first compute a flux field to balance the mass source and then exploit exact co-chain complexes to generate a solenoidal correction. A reduced basis method based on proper orthogonal decomposition is employed to construct the correction and we show that mass balance is ensured regardless of the quality of the reduced basis approximation. The method is directly applicable to mixed finite and virtual element methods, among other structure-preserving discretization techniques, and we present the extension to Darcy flow in fractured porous media. △ Less

Submitted 31 May, 2022; originally announced May 2022.

MSC Class: 65N22; 76M10; 55U15; 65N30

Bend 3d Mixed Virtual Element Method for Elliptic Problems

Authors: Franco Dassi, Alessio Fumagalli, Anna Scotti, Giuseppe Vacca

Abstract: In this study, we propose a virtual element scheme to solve the Darcy problem in three physical dimensions. The main novelty, here proposed, is that curved elements are naturally handled without any degradation of the solution accuracy. In fact, in presence of curved boundaries, or internal interfaces, the geometrical error introduced by planar approximations may dominate the convergence rate limi… ▽ More In this study, we propose a virtual element scheme to solve the Darcy problem in three physical dimensions. The main novelty, here proposed, is that curved elements are naturally handled without any degradation of the solution accuracy. In fact, in presence of curved boundaries, or internal interfaces, the geometrical error introduced by planar approximations may dominate the convergence rate limiting the benefit of high-order approximations. We consider the Darcy problem in its mixed form to directly obtain, with our numerical scheme, accurate and mass conservative fluxes without any post-processing. An important step to derive this new scheme is the actual computation of polynomials over curved polyhedrons, here presented and discussed. Finally, we show the theoretical analysis of the scheme as well as several numerical examples to support our findings △ Less

Submitted 20 November, 2021; originally announced November 2021.

A Virtual Element Method for the wave equation on curved edges in two dimensions

Authors: Franco Dassi, Alessio Fumagalli, Ilario Mazzieri, Anna Scotti, Giuseppe Vacca

Abstract: In this work we present an extension of the Virtual Element Method with curved edges for the numerical approximation of the second order wave equation in a bidimensional setting. Curved elements are used to describe the domain boundary, as well as internal interfaces corresponding to the change of some mechanical parameters. As opposite to the classic and isoparametric Finite Element approaches, w… ▽ More In this work we present an extension of the Virtual Element Method with curved edges for the numerical approximation of the second order wave equation in a bidimensional setting. Curved elements are used to describe the domain boundary, as well as internal interfaces corresponding to the change of some mechanical parameters. As opposite to the classic and isoparametric Finite Element approaches, where the geometry of the domain is approximated respectively by piecewise straight lines and by higher order polynomial maps, in the proposed method the geometry is exactly represented, thus ensuring a highly accurate numerical solution. Indeed, if in the former approach the geometrical error might deteriorate the quality of the numerical solution, in the latter approach the curved interfaces/boundaries are approximated exactly guaranteeing the expected order of convergence for the numerical scheme. Theoretical results and numerical findings confirm the validity of the proposed approach. △ Less

Submitted 11 June, 2021; originally announced June 2021.

Model adaptation in a discrete fracture network: existence of solutions and numerical strategies

Authors: Alessio Fumagalli, Francesco Saverio Patacchini

Abstract: Fractures are normally present in the underground and are, for some physical processes, of paramount importance. Their accurate description is fundamental to obtain reliable numerical outcomes useful, e.g., for energy management. Depending on the physical and geometrical properties of the fractures, fluid flow can behave differently, going from a slow Darcian regime to more complicated Brinkman or… ▽ More Fractures are normally present in the underground and are, for some physical processes, of paramount importance. Their accurate description is fundamental to obtain reliable numerical outcomes useful, e.g., for energy management. Depending on the physical and geometrical properties of the fractures, fluid flow can behave differently, going from a slow Darcian regime to more complicated Brinkman or even Forchheimer regimes for high velocity. The main problem is to determine where in the fractures one regime is more adequate than others. In order to determine these low-speed and high-speed regions, this work proposes an adaptive strategy which is based on selecting the appropriate constitutive law linking velocity and pressure according to a threshold criterion on the magnitude of the fluid velocity itself. Both theoretical and numerical aspects are considered and investigated, showing the potentiality of the proposed approach. From the analytical viewpoint, we show existence of weak solutions to such model under reasonable hypotheses on the constitutive laws. To this end, we use a variational approach identifying solutions with minimizers of an underlying energy functional. From the numerical viewpoint, we propose a one-dimensional algorithm which tracks the interface between the low- and high-speed regions. By running numerical experiments using this algorithm, we illustrate some interesting behaviors of our adaptive model on a single fracture and small networks of intersecting fractures. △ Less

Submitted 2 March, 2021; originally announced March 2021.

The mixed virtual element method for grids with curved interfaces

Authors: Franco Dassi, Alessio Fumagalli, Davide Losapio, Stefano Scialò, Anna Scotti, Giuseppe Vacca

Abstract: In many applications the accurate representation of the computational domain is a key factor to obtain reliable and effective numerical solutions. Curved interfaces, which might be internal, related to physical data, or portions of the physical boundary, are often met in real applications. However, they are often approximated leading to a geometrical error that might become dominant and deteriorat… ▽ More In many applications the accurate representation of the computational domain is a key factor to obtain reliable and effective numerical solutions. Curved interfaces, which might be internal, related to physical data, or portions of the physical boundary, are often met in real applications. However, they are often approximated leading to a geometrical error that might become dominant and deteriorate the quality of the results. Underground problems often involve the motion of fluids where the fundamental governing equation is the Darcy law. High quality velocity fields are of paramount importance for the successful subsequent coupling with other physical phenomena such as transport. The virtual element method, as solution scheme, is known to be applicable in problems whose discretizations requires cells of general shape, and the mixed formulation is here preferred to obtain accurate velocity fields. To overcome the issues associated to the complex geometries and, at the same time, retaining the quality of the solutions, we present here the virtual element method to solve the Darcy problem, in mixed form, in presence of curved interfaces in two and three dimensions. The numerical scheme is presented in detail explaining the discrete setting with a focus on the treatment of curved interfaces. Examples, inspired from industrial applications, are presented showing the validity of the proposed approach. △ Less

Submitted 18 November, 2020; originally announced November 2020.

A multi-layer reactive transport model for fractured porous media

Authors: Luca Formaggia, Alessio Fumagalli, Anna Scotti

Abstract: An accurate modeling of reactive flows in fractured porous media is a key ingredient to obtain reliable numerical simulations of several industrial and environmental applications. For some values of the physical parameters we can observe the formation of a narrow region or layer around the fractures where chemical reactions are focused. Here the transported solute may precipitate and form a salt,… ▽ More An accurate modeling of reactive flows in fractured porous media is a key ingredient to obtain reliable numerical simulations of several industrial and environmental applications. For some values of the physical parameters we can observe the formation of a narrow region or layer around the fractures where chemical reactions are focused. Here the transported solute may precipitate and form a salt, or vice-versa. This phenomenon has been observed and reported in real outcrops. By changing its physical properties this layer might substantially alter the global flow response of the system and thus the actual transport of solute: the problem is thus non-linear and fully coupled. The aim of this work is to propose a new mathematical model for reactive flow in fractured porous media, by approximating both the fracture and these surrounding layers via a reduced model. In particular, our main goal is to describe the layer thickness evolution with a new mathematical model, and compare it to a fully resolved equidimensional model for validation. As concerns numerical approximation we extend an operator splitting scheme in time to solve sequentially, at each time step, each physical process thus avoiding the need for a non-linear monolithic solver, which might be challenging due to the non-smoothness of the reaction rate. We consider bi- and tridimensional numerical test cases to asses the accuracy and benefit of the proposed model in realistic scenarios. △ Less

Submitted 9 November, 2020; originally announced November 2020.

The Mixed Virtual Element Method on curved edges in two dimensions

Authors: Franco Dassi, Alessio Fumagalli, Davide Losapio, Stefano Scialò, Anna Scotti, Giuseppe Vacca

Abstract: In this work, we propose an extension of the mixed Virtual Element Method (VEM) for bi-dimensional computational grids with curvilinear edge elements. The approximation by means of rectilinear edges of a domain with curvilinear geometrical feature, such as a portion of domain boundary or an internal interface, may introduce a geometrical error that degrades the expected order of convergence of the… ▽ More In this work, we propose an extension of the mixed Virtual Element Method (VEM) for bi-dimensional computational grids with curvilinear edge elements. The approximation by means of rectilinear edges of a domain with curvilinear geometrical feature, such as a portion of domain boundary or an internal interface, may introduce a geometrical error that degrades the expected order of convergence of the scheme. In the present work a suitable VEM approximation space is proposed to consistently handle curvilinear geometrical objects, thus recovering optimal convergence rates. The resulting numerical scheme is presented along with its theoretical analysis and several numerical test cases to validate the proposed approach. △ Less

Submitted 27 July, 2020; originally announced July 2020.

A mathematical model for thermal single-phase flow and reactive transport in fractured porous media

Authors: Alessio Fumagalli, Anna Scotti

Abstract: In this paper we present a mathematical model and a numerical workflow for the simulation of a thermal single-phase flow with reactive transport in porous media, in the presence of fractures. The latter are thin regions which might behave as high or low permeability channels depending on their physical parameters, and are thus of paramount importance in underground flow problems. Chemical reaction… ▽ More In this paper we present a mathematical model and a numerical workflow for the simulation of a thermal single-phase flow with reactive transport in porous media, in the presence of fractures. The latter are thin regions which might behave as high or low permeability channels depending on their physical parameters, and are thus of paramount importance in underground flow problems. Chemical reactions may alter the local properties of the porous media as well as the fracture walls, changing the flow path and possibly occluding some portions of the fractures or zones in the porous media. To solve numerically the coupled problem we propose a temporal splitting scheme so that the equations describing each physical process are solved sequentially. Numerical tests shows the accuracy of the proposed model and the ability to capture complex phenomena, where one or multiple fractures are present. △ Less

Submitted 19 May, 2020; originally announced May 2020.

Numerical Methods for Flow in Fractured Porous Media

Authors: Luca Formaggia, Anna Scotti, Alessio Fumagalli

Abstract: In this work we present the mathematical models for single-phase flow in fractured porous media. An overview of the most common approaches is considered, which includes continuous fracture models and discrete fracture models. For the latter, we discuss strategies that are developed in literature for its numerical solution mainly related to the geometrical relation between the fractures and porous… ▽ More In this work we present the mathematical models for single-phase flow in fractured porous media. An overview of the most common approaches is considered, which includes continuous fracture models and discrete fracture models. For the latter, we discuss strategies that are developed in literature for its numerical solution mainly related to the geometrical relation between the fractures and porous media grids. △ Less

Submitted 31 March, 2020; originally announced March 2020.

Performances of the mixed virtual element method on complex grids for underground flow

Authors: Alessio Fumagalli, Anna Scotti, Luca Formaggia

Abstract: The numerical simulation of physical processes in the underground frequently entails challenges related to the geometry and/or data. The former are mainly due to the shape of sedimentary layers and the presence of fractures and faults, while the latter are connected to the properties of the rock matrix which might vary abruptly in space. The development of approximation schemes has recently focuse… ▽ More The numerical simulation of physical processes in the underground frequently entails challenges related to the geometry and/or data. The former are mainly due to the shape of sedimentary layers and the presence of fractures and faults, while the latter are connected to the properties of the rock matrix which might vary abruptly in space. The development of approximation schemes has recently focused on the overcoming of such difficulties with the objective of obtaining numerical schemes with good approximation properties. In this work we carry out a numerical study on the performances of the Mixed Virtual Element Method (MVEM) for the solution of a single-phase flow model in fractured porous media. This method is able to handle grid cells of polytopal type and treat hybrid dimensional problems. It has been proven to be robust with respect to the variation of the permeability field and of the shape of the elements. Our numerical experiments focus on two test cases that cover several of the aforementioned critical aspects. △ Less

Submitted 27 February, 2020; originally announced February 2020.

Verification benchmarks for single-phase flow in three-dimensional fractured porous media

Authors: Inga Berre, Wietse M. Boon, Bernd Flemisch, Alessio Fumagalli, Dennis Gläser, Eirik Keilegavlen, Anna Scotti, Ivar Stefansson, Alexandru Tatomir, Konstantin Brenner, Samuel Burbulla, Philippe Devloo, Omar Duran, Marco Favino, Julian Hennicker, I-Hsien Lee, Konstantin Lipnikov, Roland Masson, Klaus Mosthaf, Maria Giuseppina Chiara Nestola, Chuen-Fa Ni, Kirill Nikitin, Philipp Schädle, Daniil Svyatskiy, Ruslan Yanbarisov , et al. (1 additional authors not shown)

Abstract: Flow in fractured porous media occurs in the earth's subsurface, in biological tissues, and in man-made materials. Fractures have a dominating influence on flow processes, and the last decade has seen an extensive development of models and numerical methods that explicitly account for their presence. To support these developments, we present a portfolio of four benchmark cases for single-phase flo… ▽ More Flow in fractured porous media occurs in the earth's subsurface, in biological tissues, and in man-made materials. Fractures have a dominating influence on flow processes, and the last decade has seen an extensive development of models and numerical methods that explicitly account for their presence. To support these developments, we present a portfolio of four benchmark cases for single-phase flow in three-dimensional fractured porous media. The cases are specifically designed to test the methods' capabilities in handling various complexities common to the geometrical structures of fracture networks. Based on an open call for participation, results obtained with 17 numerical methods were collected. This paper presents the underlying mathematical model, an overview of the features of the participating numerical methods, and their performance in solving the benchmark cases. △ Less

Submitted 17 February, 2020; originally announced February 2020.

Reactive flow in fractured porous media

Authors: Alessio Fumagalli, Anna Scotti

Abstract: In this work we present a model reduction procedure to derive a hybrid-dimensional framework for the mathematical modeling of reactive transport in fractured porous media. Fractures are essential pathways in the underground which allow fast circulation of the fluids present in the rock matrix, often characterized by low permeability. However, due to infilling processes fractures may change their h… ▽ More In this work we present a model reduction procedure to derive a hybrid-dimensional framework for the mathematical modeling of reactive transport in fractured porous media. Fractures are essential pathways in the underground which allow fast circulation of the fluids present in the rock matrix, often characterized by low permeability. However, due to infilling processes fractures may change their hydraulic properties and become barriers for the flow creating impervious blocks. The geometrical as well as the physical properties of the fractures require a special treatment to allow the subsequent numerical discretization to be affordable and accurate. The aim of this work is to introduce a simple yet complete mathematical model to account for such diagenetic effects where chemical reactions will occlude or empty portions of the porous media and, in particular, fractures. △ Less

Submitted 15 January, 2020; originally announced January 2020.

PorePy: An Open-Source Software for Simulation of Multiphysics Processes in Fractured Porous Media

Authors: Eirik Keilegavlen, Runar Berge, Alessio Fumagalli, Michele Starnoni, Ivar Stefansson, Jhabriel Varela, Inga Berre

Abstract: Development of models and dedicated numerical methods for dynamics in fractured rocks is an active research field, with research moving towards increasingly advanced process couplings and complex fracture networks. The inclusion of coupled processes in simulation models is challenged by the high aspect ratio of the fractures, the complex geometry of fracture networks and the crucial impact of proc… ▽ More Development of models and dedicated numerical methods for dynamics in fractured rocks is an active research field, with research moving towards increasingly advanced process couplings and complex fracture networks. The inclusion of coupled processes in simulation models is challenged by the high aspect ratio of the fractures, the complex geometry of fracture networks and the crucial impact of processes that completely change characteristics on the fracture-rock interface. This paper provides a general discussion of design principles for introducing fractures in simulators, and defines a framework for integrated modeling, discretization and computer implementation. The framework is implemented in the simulation software PorePy, which can serve as a flexible prototyping tool or multiphysics problems in fractured rocks. Based on a representation of the fractures and their intersections as lower-dimensional objects, we discuss data structures for mixed-dimensional meshes, formulation of multiphysics problems and discretizations that utilize existing software. We further present the implementation of these concepts in the PorePy open-source software tool, which is aimed at coupled simulation of flow and transport in three-dimensional fractured reservoirs as well as deformation of fractures and the reservoir in general. We present validation by benchmarks for flow, poroelasticity and fracture deformation in fractured porous media. The flexibility of the framework is then illustrated by simulations of fully coupled flow and transport and of injection driven deformation of fractures. All results reported herein can be reproduced by openly available simulation scripts. △ Less

Submitted 26 August, 2019; originally announced August 2019.

Analysis of conforming, non-matching, and polygonal methods for Darcy and advection-diffusion-reaction simulations in discrete fracture networks

Authors: Andrea Borio, Alessio Fumagalli, Stefano Scialò

Abstract: The aim of this study is to compare numerical methods for the simulation of single-phase flow and transport in fractured media, described here by means of the Discrete Fracture Network (DFN) model. A Darcy problem is solved to compute the advective field, then used in a subsequent time dependent transport-diffusion-reaction problem. The numerical schemes are benchmarked in terms of flexibility in… ▽ More The aim of this study is to compare numerical methods for the simulation of single-phase flow and transport in fractured media, described here by means of the Discrete Fracture Network (DFN) model. A Darcy problem is solved to compute the advective field, then used in a subsequent time dependent transport-diffusion-reaction problem. The numerical schemes are benchmarked in terms of flexibility in handling geometrical complexity, mass conservation and stability issues for advection dominated flow regimes. To this end, three benchmark cases have been specifically designed and are here proposed, representing some of the most critical issues encountered in DFN simulations. △ Less

Submitted 29 July, 2019; originally announced July 2019.

Robust linear domain decomposition schemes for reduced non-linear fracture flow models

Authors: Elyes Ahmed, Alessio Fumagalli, Ana Budiša, Eirik Keilegavlen, Jan Martin Nordbotten, Florin Adrian Radu

Abstract: In this work, we consider compressible single-phase flow problems in a porous media containing a fracture. In the latter, a non-linear pressure-velocity relation is prescribed. Using a non-overlapping domain decomposition procedure, we reformulate the global problem into a non-linear interface problem. We then introduce two new algorithms that are able to efficiently handle the non-linearity and t… ▽ More In this work, we consider compressible single-phase flow problems in a porous media containing a fracture. In the latter, a non-linear pressure-velocity relation is prescribed. Using a non-overlapping domain decomposition procedure, we reformulate the global problem into a non-linear interface problem. We then introduce two new algorithms that are able to efficiently handle the non-linearity and the coupling between the fracture and the matrix, both based on linearization by the so-called L-scheme. The first algorithm, named MoLDD, uses the L-scheme to resolve the non-linearity, requiring at each iteration to solve the dimensional coupling via a domain decomposition approach. The second algorithm, called ItLDD, uses a sequential approach in which the dimensional coupling is part of the linearization iterations. For both algorithms, the computations are reduced only to the fracture by pre-computing, in an offline phase, a multiscale flux basis (the linear Robin-to-Neumann co-dimensional map), that represent the flux exchange between the fracture and the matrix. We present extensive theoretical findings and in particular, the stability and the convergence of both schemes are obtained, where user given parameters are optimized to minimise the number of iterations. Examples on two important fracture models are computed with the library PorePy and agree with the developed theory. △ Less

Submitted 1 October, 2020; v1 submitted 13 June, 2019; originally announced June 2019.

A multi-layer reduced model for flow in porous media with a fault and surrounding damage zones

Authors: Alessio Fumagalli, Anna Scotti

Abstract: In this work we present a new conceptual model to describe fluid flow in a porous media system in presence of a large fault. Geological faults are often modeled simply as interfaces in the rock matrix, but they are complex structure where the high strain core is surrounded by the so called damage zones, characterized by the presence of smaller fractures which enhance the permeability of the medium… ▽ More In this work we present a new conceptual model to describe fluid flow in a porous media system in presence of a large fault. Geological faults are often modeled simply as interfaces in the rock matrix, but they are complex structure where the high strain core is surrounded by the so called damage zones, characterized by the presence of smaller fractures which enhance the permeability of the medium. To obtain reliable simulation outcomes these damage zone, as well as the fault, have to be accurately described. The new model proposed in this work considers both these two regions as lower dimensional and embedded in the rock matrix. The model is presented, analyzed, and tested in several configurations to prove its robustness and ability to capture many important features, such as hight contrast and heterogeneity of permeability. △ Less

Submitted 4 March, 2019; originally announced March 2019.

A multiscale flux basis for mortar mixed discretizations of reduced Darcy-Forchheimer fracture models

Authors: Elyes Ahmed, Alessio Fumagalli, Ana Budiša

Abstract: In this paper, a multiscale flux basis algorithm is developed to efficiently solve a flow problem in fractured porous media. Here, we take into account a mixed-dimensional setting of the discrete fracture matrix model, where the fracture network is represented as lower-dimensional object. We assume the linear Darcy model in the rock matrix and the non-linear Forchheimer model in the fractures. In… ▽ More In this paper, a multiscale flux basis algorithm is developed to efficiently solve a flow problem in fractured porous media. Here, we take into account a mixed-dimensional setting of the discrete fracture matrix model, where the fracture network is represented as lower-dimensional object. We assume the linear Darcy model in the rock matrix and the non-linear Forchheimer model in the fractures. In our formulation, we are able to reformulate the matrix-fracture problem to only the fracture network problem and, therefore, significantly reduce the computational cost. The resulting problem is then a non-linear interface problem that can be solved using a fixed-point or Newton-Krylov methods, which in each iteration require several solves of Robin problems in the surrounding rock matrices. To achieve this, the flux exchange (a linear Robin-to-Neumann co-dimensional mapping) between the porous medium and the fracture network is done offline by pre-computing a multiscale flux basis that consists of the flux response from each degree of freedom on the fracture network. This delivers a conserve for the basis that handles the solutions in the rock matrices for each degree of freedom in the fractures pressure space. Then, any Robin sub-domain problems are replaced by linear combinations of the multiscale flux basis during the interface iteration. The proposed approach is, thus, agnostic to the physical model in the fracture network. Numerical experiments demonstrate the computational gains of pre-computing the flux exchange between the porous medium and the fracture network against standard non-linear domain decomposition approaches. △ Less

Submitted 30 October, 2018; originally announced October 2018.

Journal ref: Computer Methods in Applied Mechanics and Engineering 354 (2019)16-36

Call for participation: Verification benchmarks for single-phase flow in three-dimensional fractured porous media

Authors: Inga Berre, Wietse Boon, Bernd Flemisch, Alessio Fumagalli, Dennis Gläser, Eirik Keilegavlen, Anna Scotti, Ivar Stefansson, Alexandru Tatomir

Abstract: This call for participation proposes four benchmark tests to verify and compare numerical schemes to solve single-phase flow in fractured porous media. With this, the two-dimensional suite of benchmark tests presented by Flemisch et al. 2018 is extended to include three-dimensional problems. Moreover, transport simulations are included as a means to compare discretization methods for flow. With th… ▽ More This call for participation proposes four benchmark tests to verify and compare numerical schemes to solve single-phase flow in fractured porous media. With this, the two-dimensional suite of benchmark tests presented by Flemisch et al. 2018 is extended to include three-dimensional problems. Moreover, transport simulations are included as a means to compare discretization methods for flow. With this publication, we invite researchers to contribute to the study by providing results to the test cases based on their applied discretization methods. △ Less

Submitted 18 September, 2018; originally announced September 2018.

Conforming, non-conforming and non-matching discretization couplings in discrete fracture network simulations

Authors: Alessio Fumagalli, Eirik Keilegavlen, Stefano Scialò

Abstract: Simulations of fluid flow in naturally fractured rocks have implications for several subsurface applications, including energy storage and extraction, and waste storage. We are interested in flow in discrete fracture networks, which explicitly represent flow in fracture surfaces, but ignore the impact of the surrounding host rock. Fracture networks, generated from observations or stochastic simula… ▽ More Simulations of fluid flow in naturally fractured rocks have implications for several subsurface applications, including energy storage and extraction, and waste storage. We are interested in flow in discrete fracture networks, which explicitly represent flow in fracture surfaces, but ignore the impact of the surrounding host rock. Fracture networks, generated from observations or stochastic simulations, will contain intersections of arbitrary length, and intersection lines can further cross, forming a highly complex geometry. As the flow exchange between fractures, thus in the network, takes place in these intersections, an adequate representation of the geometry is critical for simulation accuracy. In practice, the intersection dynamics must be handled by a combination of the simulation grid, which may or may not resolve the intersection lines, and the numerical methods applied on the grid. In this work, we review different classes of numerical approaches proposed in recent years, covering both methods that conform to the grid, and non-matching cases. Specific methods considered herein include finite element, mixed and virtual finite elements and control volume methods. We expose our methods to an extensive set of test cases, ranging from artificial geometries designed to test difficult configurations, to a network extruded from a real fracture outcrop. The main outcome is guidances for choice of simulation models and numerical discretization with a trade off on the computational cost and solution accuracy. △ Less

Submitted 6 March, 2018; v1 submitted 5 March, 2018; originally announced March 2018.

Unified approach to discretization of flow in fractured porous media

Authors: Jan M. Nordbotten, Wietse M. Boon, Alessio Fumagalli, Eirik Keilegavlen

Abstract: In this paper, we introduce a mortar-based approach to discretizing flow in fractured porous media, which we term the mixed-dimensional flux coupling scheme. Our formulation is agnostic to the discretizations used to discretize the fluid flow equations in the porous medium and in the fractures, and as such it represents a unified approach to integrated fractured geometries into any existing discre… ▽ More In this paper, we introduce a mortar-based approach to discretizing flow in fractured porous media, which we term the mixed-dimensional flux coupling scheme. Our formulation is agnostic to the discretizations used to discretize the fluid flow equations in the porous medium and in the fractures, and as such it represents a unified approach to integrated fractured geometries into any existing discretization framework. In particular, several existing discretization approaches for fractured porous media can be seen as special instances of the approach proposed herein. We provide an abstract stability theory for our approach, which provides explicit guidance into the grids used to discretize the fractures and the porous medium, as dependent on discretization methods chosen for the respective domains. The theoretical results are sustained by numerical examples, wherein we utilize our framework to simulate flow in 2D and 3D fractured media using control volume methods (both two-point and multi-point flux), Lagrangian finite element methods, mixed finite element methods, and virtual element methods. As expected, regardless of the ambient methods chosen, our approach leads to stable and convergent discretizations for the fractured problems considered, within the limits of the discretization schemes. △ Less

Submitted 9 August, 2018; v1 submitted 16 February, 2018; originally announced February 2018.

Implementation of mixed-dimensional models for flow in fractured porous media

Authors: Eirik Keilegavlen, Alessio Fumagalli, Runar Berge, Ivar Stefansson

Abstract: Models that involve coupled dynamics in a mixed-dimensional geometry are of increasing interest in several applications. Here, we describe the development of a simulation model for flow in fractured porous media, where the fractures and their intersections form a hierarchy of interacting subdomains. We discuss the implementation of a simulation framework, with an emphasis on reuse of existing disc… ▽ More Models that involve coupled dynamics in a mixed-dimensional geometry are of increasing interest in several applications. Here, we describe the development of a simulation model for flow in fractured porous media, where the fractures and their intersections form a hierarchy of interacting subdomains. We discuss the implementation of a simulation framework, with an emphasis on reuse of existing discretization tools for mono-dimensional problems. The key ingredients are the representation of the mixed-dimensional geometry as a graph, which allows for convenient discretization and data storage, and a non-intrusive coupling of dimensions via boundary conditions and source terms. This approach is applicable for a wide class of mixed-dimensional problems. We show simulation results for a flow problem in a three-dimensional fracture geometry, applying both finite volume and virtual finite element discretizations. △ Less

Submitted 20 December, 2017; originally announced December 2017.

PorePy: An Open-Source Simulation Tool for Flow and Transport in Deformable Fractured Rocks

Authors: Eirik Keilegavlen, Alessio Fumagalli, Runar Berge, Ivar Stefansson, Inga Berre

Abstract: Fractures are ubiquitous in the subsurface and strongly affect flow and deformation. The physical shape of the fractures, they are long and thin objects, puts strong limitations on how the effect of this dynamics can be incorporated into standard reservoir simulation tools. This paper reports the development of an open-source software framework, termed PorePy, which is aimed at simulation of flow… ▽ More Fractures are ubiquitous in the subsurface and strongly affect flow and deformation. The physical shape of the fractures, they are long and thin objects, puts strong limitations on how the effect of this dynamics can be incorporated into standard reservoir simulation tools. This paper reports the development of an open-source software framework, termed PorePy, which is aimed at simulation of flow and transport in three-dimensional fractured reservoirs, as well as deformation of the reservoir due to shearing along fracture and fault planes. Starting from a description of fractures as polygons embedded in a 3D domain, PorePy provides semi-automatic gridding to construct a discrete-fracture-matrix model, which forms the basis for subsequent simulations. PorePy allows for flow and transport in all lower-dimensional objects, including planes (2D) representing fractures, and lines (1D) and points (0D), representing fracture intersections. Interaction between processes in neighboring domains of different dimension is implemented as a sequence of couplings of objects one dimension apart. This readily allows for handling of complex fracture geometries compared to capabilities of existing software. In addition to flow and transport, PorePy provides models for rock mechanics, poro-elasticity and coupling with fracture deformation models. The software is fully open, and can serve as a framework for transparency and reproducibility of simulations. We describe the design principles of PorePy from a user perspective, with focus on possibilities within gridding, covered physical processes and available discretizations. The power of the framework is illustrated with two sets of simulations; involving respectively coupled flow and transport in a fractured porous medium, and low-pressure stimulation of a geothermal reservoir. △ Less

Submitted 21 December, 2017; v1 submitted 1 December, 2017; originally announced December 2017.

Dual Virtual Element Methods for Discrete Fracture Matrix Models

Authors: Alessio Fumagalli, Eirik Keilegavlen

Abstract: The accurate description of fluid flow and transport in fractured porous media is of paramount importance to capture the macroscopic behaviour of an oil reservoir, a geothermal system, or a CO2 sequestration site, to name few applications. The construction of accurate simulation model for flow in fractures is challenging due to the high ratios between a fracture's length and width, which makes mod… ▽ More The accurate description of fluid flow and transport in fractured porous media is of paramount importance to capture the macroscopic behaviour of an oil reservoir, a geothermal system, or a CO2 sequestration site, to name few applications. The construction of accurate simulation model for flow in fractures is challenging due to the high ratios between a fracture's length and width, which makes modeling by lower-dimensional manifolds a natural option. In this paper we present a mixed-dimensional Darcy problem able to describe pressure and Darcy velocity in all the dimensions, i.e. in the rock matrix, in the fractures, and in their intersections. Moreover, we present a mixed-dimensional transport problem which, given the Darcy velocity, describes coupled advection and diffusion of a passive scalar into the fractured porous media. The approach can handle both conducting and blocking fractures. Our computational grids are created by coarsening of simplex tessellations that conform to the fractures surfaces. An accurate choice of the discrete approximation of the previous model, by virtual finite element and finite volume, allows us to simulate complex problem with a good balance in term of accuracy and computational cost. We illustrate the performance of our method by comparing to benchmark studies for two-dimensional fractured porous media, as well as a complex three-dimensional fracture geometry. △ Less

Submitted 6 November, 2017; originally announced November 2017.

A double-layer reduced model for fault flow on slipping domains with hybrid finite volume scheme

Authors: Alessio Fumagalli, Isabelle Faille

Abstract: In this work we are interested in dealing with single-phase flows in fractured porous media for underground processes. We focus our attention on domains where the presence of faults, with thickness several orders of magnitude smaller than other characteristic sizes, can allow one part of the domain to slide past to the other. We propose a mathematical scheme where a reduced model for the fault flo… ▽ More In this work we are interested in dealing with single-phase flows in fractured porous media for underground processes. We focus our attention on domains where the presence of faults, with thickness several orders of magnitude smaller than other characteristic sizes, can allow one part of the domain to slide past to the other. We propose a mathematical scheme where a reduced model for the fault flows is employed yielding a problem of co-dimension one. The hybrid finite volume method is used to obtain the discretized problem, which employs two different meshes for each fault, one associated with the porous-medium domain on each side of the fault. These two meshes can move with the corresponding domain, resulting in non-matching grids between the two parts of the fault. In an earlier paper a mathematical scheme was proposed where the numerical discretization considers the hybrid finite volume method. In this paper we focus on the well-posedness of the continuous problem, the convergence of the discretized problem, and with several numerical tests we support the theoretical findings. △ Less

Submitted 18 September, 2017; originally announced September 2017.

Dual virtual element method in presence of an inclusion

Authors: Alessio Fumagalli

Abstract: We consider a Darcy problem for saturated porous media written in dual formulation in presence of a fully immersed inclusion. The lowest order virtual element method is employ to derive the discrete approximation. In the present work we study the effect of cells with cuts on the numerical solution, able to geometrically handle in a more natural way the inclusion tips. The numerical results show th… ▽ More We consider a Darcy problem for saturated porous media written in dual formulation in presence of a fully immersed inclusion. The lowest order virtual element method is employ to derive the discrete approximation. In the present work we study the effect of cells with cuts on the numerical solution, able to geometrically handle in a more natural way the inclusion tips. The numerical results show the validity of the proposed approach. △ Less

Submitted 11 September, 2017; originally announced September 2017.

Benchmarks for single-phase flow in fractured porous media

Authors: Bernd Flemisch, Inga Berre, Wietse Boon, Alessio Fumagalli, Nicolas Schwenck, Anna Scotti, Ivar Stefansson, Alexandru Tatomir

Abstract: This paper presents several test cases intended to be benchmarks for numerical schemes for single-phase fluid flow in fractured porous media. A number of solution strategies are compared, including a vertex and a cell-centered finite volume method, a non-conforming embedded discrete fracture model, a primal and a dual extended finite element formulation, and a mortar discrete fracture model. The p… ▽ More This paper presents several test cases intended to be benchmarks for numerical schemes for single-phase fluid flow in fractured porous media. A number of solution strategies are compared, including a vertex and a cell-centered finite volume method, a non-conforming embedded discrete fracture model, a primal and a dual extended finite element formulation, and a mortar discrete fracture model. The proposed benchmarks test the schemes by increasing the difficulties in terms of network geometry, e.g. intersecting fractures, and physical parameters, e.g. low and high fracture-matrix permeability ratio as well as heterogeneous fracture permeabilities. For each problem, the results presented by the participants are the number of unknowns, the approximation errors in the porous matrix and in the fractures with respect to a reference solution, and the sparsity and condition number of the discretized linear system. All data and meshes used in this study are publicly available for further comparisons. △ Less

Submitted 5 January, 2017; originally announced January 2017.

Dual virtual element method for discrete fractures networks

Authors: Alessio Fumagalli, Eirik Keilegavlen

Abstract: Discrete fracture networks is a key ingredient in the simulation of physical processes which involve fluid flow in the underground, when the surrounding rock matrix is considered impervious. In this paper we present two different models to compute the pressure field and Darcy velocity in the system. The first allows a normal flow out of a fracture at the intersections, while the second grants also… ▽ More Discrete fracture networks is a key ingredient in the simulation of physical processes which involve fluid flow in the underground, when the surrounding rock matrix is considered impervious. In this paper we present two different models to compute the pressure field and Darcy velocity in the system. The first allows a normal flow out of a fracture at the intersections, while the second grants also a tangential flow along the intersections. For the numerical discretization, we use the mixed virtual finite element method as it is known to handle grid elements of, almost, any arbitrary shape. The flexibility of the discretization allows us to loosen the requirements on grid construction, and thus significantly simplify the flow discretization compared to traditional discrete fracture network models. A coarsening algorithm, from the algebraic multigrid literature, is also considered to further speed up the computation. The performance of the method is validated by numerical experiments. △ Less

Submitted 16 August, 2017; v1 submitted 10 October, 2016; originally announced October 2016.