Skip to Main content Skip to Navigation
Journal articles

A generalized mixed hybrid mortar method for solving flow in stochastic discrete fracture networks

Abstract : The simulation of flow in fractured media requires handling both a large number of fractures and a complex interconnecting network of these fractures. Networks considered in this paper are three-dimensional domains made up of two-dimensional fractures intersecting each other and randomly generated. Due to the stochastic generation of fractures, intersections can be highly intricate. The numerical method must generate a mesh and define a discrete problem for any discrete fracture network (DFN). A first approach [Erhel, de Dreuzy, and Poirriez, SIAM J. Sci. Comput., 31 (2009), pp. 2688-2705] is to generate a conforming mesh and to apply a mixed hybrid finite element method. However, the resulting linear system becomes very large when the network contains many fractures. Hence a second approach [Pichot, Erhel, and de Dreuzy, Appl. Anal., 89 (2010), pp. 1629-1643] is to generate a nonconforming mesh, using an independent mesh generation for each fracture. Then a Mortar technique applied to the mixed hybrid finite element method deals with the nonmatching grids. When intersections do not cross or overlap, pairwise Mortar relations for each intersection are efficient [Pichot, Erhel, and de Dreuzy, 2010]. But for most random networks, discretized intersections involve more than two fractures. In this paper, we design a new method generalizing the previous one and that is applicable for stochastic networks. The main idea is to combine pairwise Mortar relations with additional relations for the overlapping part. This method still ensures the continuity of fluxes and heads and still yields a symmetric positive definite linear system. Numerical experiments show the efficiency of the method applied to complex stochastic fracture networks. We also study numerical convergence when reducing the mesh step. This method makes it easy to perform mesh optimization and appears to be a very promising tool to simulate flow in multiscale fracture networks.
Document type :
Journal articles
Complete list of metadatas

Cited literature [23 references]  Display  Hide  Download
Contributor : Isabelle Dubigeon <>
Submitted on : Thursday, October 10, 2013 - 11:40:48 AM
Last modification on : Friday, March 6, 2020 - 1:28:21 AM
Long-term archiving on: : Saturday, January 11, 2014 - 2:30:25 AM


Publisher files allowed on an open archive



Géraldine Pichot, Jocelyne Erhel, Jean-Raynald de Dreuzy. A generalized mixed hybrid mortar method for solving flow in stochastic discrete fracture networks. SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2012, 34 (1), pp.B86-B105. ⟨10.1137/100804383⟩. ⟨insu-00681662⟩



Record views


Files downloads