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

Subsurface flows are influenced by the presence of faults and large fractures which act as preferential paths or barriers for the flow. In literature models were proposed to handle fractures in a porous medium as objects of codimension 1. In this work we consider the case of a network of intersecting fractures, with the aim of deriving physically consistent and effective interface conditions to impose at the intersection between fractures. This new model accounts for the angle between fractures at the intersections and allows for jumps of pressure across intersections. This fact permits to describe the flow when fractures are characterized by different properties more accurately with respect to other models that impose pressure continuity. The main mathematical properties of the model, derived in the two-dimensional setting, are analyzed. As concerns the numerical discretization we allow the grids of the fractures to be independent, thus in general non-matching at the intersection, by means of the extended finite element method (XFEM). This increases the flexibility of the method in the case of complex geometries characterized by a high number of fractures.

Abstract.We consider an incompressible flow problem in a N -dimensional fractured porous domain (Darcy's problem). The fracture is represented by a (N −1)-dimensional interface, exchanging fluid with the surrounding media. In this paper we consider the lowest-order (RT0, P0) Raviart-Thomas mixed finite element method for the approximation of the coupled Darcy's flows in the porous media and within the fracture, with independent meshes for the respective domains. This is achieved thanks to an enrichment with discontinuous basis functions on triangles crossed by the fracture and a weak imposition of interface conditions. First, we study the stability and convergence properties of the resulting numerical scheme in the uncoupled case, when the known solution of the fracture problem provides an immersed boundary condition. We detail the implementation issues and discuss the algebraic properties of the associated linear system. Next, we focus on the coupled problem and propose an iterative porous domain/fracture domain iterative method to solve for fluid flow in both the porous media and the fracture and compare the results with those of a traditional monolithic approach. Numerical results are provided confirming convergence rates and algebraic properties predicted by the theory. In particular, we discuss preconditioning and equilibration techniques to make the condition number of the discrete problem independent of the position of the immersed interface. Finally, two and three dimensional simulations of Darcy's flow in different configurations (highly and poorly permeable fracture) are analyzed and discussed.Mathematics Subject Classification. 76S05, 35Q86, 65L60.

No abstract

Abstract.We present a possible framework for the numerical simulation of flow in fractured porous media that couples mimetic finite differences for the porous matrix with a finite volume scheme for the flow in the fractures. The resulting method is theoretically analyzed in the case of a single fracture. Moreover, several numerical experiments show the capability of the method to deal also with complicated networks of fractures. Thanks to the implementation of rather general coupling conditions, it encompasses both "conductive fractures", i.e., fractures with high permeability and "sealed fractures", i.e., fractures with low permeability which act as a flow barrier.Mathematics Subject Classification. 65N30, 35Q86, 76S05.