Unified approach to discretization of flow in fractured porous media

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