CIRM, Luminy, Marseille
Date(s) : 05/12/2022 - 09/12/2022 iCal
0 h 00 min
CIRM – Jean-Morlet Chair
Martin GANDER & Florence HUBERT
Numerical Methods for PDEs: Discretization, Iterative Solution and Parallelization
Méthodes numériques pour les EDP: discrétisation, solution itérative et parallélisation
Discrete Fracture Networks
Date: 5-9 décembre 2022
Discrete fracture networks are models to simulate flow behavior in the underground in the presence of fractures. In the fractures, the flow behavior is very different from the flow in the so-called surrounding matrix. The simulation of one fracture thus requires both a different flow model in the fracture from the flow model in the matrix, and coupling conditions between the thin fracture and the matrix. In addition, the mesh, in general, needs to be extremely fine in the fracture compared to the matrix. While this is possible for the simulation of one (or a few fractures), it becomes prohibitive when many fractures are present. In the limit of a very large number, one could then use homogenization techniques, but we are interested here in the intermediate, hard regime, where the fractures still need to be modeled individually, and there are many of them. A common method to establish such models consists in integrating the fracture equations over the fracture width and using some ad hoc approximations for the coupling conditions [2, 5, 1, 4]. We are interested here in a systematic construction of reduced order elliptic PDE models for such fractured domains. In , we derived coupling conditions between matrix subdomains, when the fracture is represented as a hypersurface embedded in the surrounded rock matrix. These coupling conditions were obtained by analytically eliminating the fracture and encoding the information of the fracture in the coupling conditions via a continuous variant of a Schur complement. We then derived local approximations of these coupling conditions using asymptotic expansions when 1 the fracture width goes to zero. A complete analysis of this approach can be found in . In the present project, we are interested in a different asymptotic regime, namely a fixed small fracture width, but when the mesh size of the numerical approximation is going to zero. The advantage of this second approach is that it will be suitable also for fractures with a more realistic fracture width, and at the same time optimized for the numerical method used. We will follow in our approach the analogy we know from domain decomposition methods, in particular optimized Schwarz methods, where for transport dominated flows, optimized transmission conditions were derived both for the physical parameter represented by the Peclet number becoming large, or the numerical mesh parameter becoming small.
 P. Angot, F. Boyer, and F. Hubert. Asymptotic and numerical modelling of flows in fractured porous media. ESAIM: Mathematical Modelling and Numerical Analysis-Modélisation Mathématique et Analyse Numérique, 43 (2):239–275, 2009.  E. Flauraud, F. Nataf, I. Faille, and R. Masson. Domain decomposition for an asymptotic geological fault modeling. Comptes Rendus Mécanique, 331(12):849–855, 2003.  M. J. Gander, J. Hennicker, and R. Masson. Asymptotic analysis for the coupling between subdomains in discrete fracture matrix models. In Domain Decomposition Methods in Science and Engineering XXV, LNCSE. Springer Verlag, 2020. [
4] M. Lesinigo, C. D’Angelo, and A. Quarteroni. A multiscale darcy–brinkman model for fluid flow in fractured porous media. Numerische Mathematik, 117(4):717–752, 2011. V. Martin, J. Jaffré, and J. E. Roberts. Modeling fractures and barriers as interfaces for flow in porous media. SIAM Journal on Scientific Computing, 26(5):1667–1691, 2005.