Unified approach to discretization of flow in fractured porous media

TL;DR

This paper presents a unified mortar-based discretization scheme for flow in fractured porous media, compatible with various existing methods, ensuring stability and convergence across different discretization approaches.

Contribution

It introduces a general, discretization-agnostic framework for modeling flow in fractured media, unifying multiple existing methods under a single approach.

Findings

Framework is compatible with multiple discretization methods.

Ensures stability and convergence in numerical simulations.

Validated through 2D and 3D numerical examples.

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