Preconditioning Techniques for the Numerical Solution of Flow in Fractured Porous Media

This work deals with the efficient iterative solution of the system of equations stemming from mimetic finite difference discretization of a hybrid-dimensional mixed Darcy problem modeling flow in fractured porous media. We investigate the spectral properties of a mixed discrete formulation based on mimetic finite differences for flow in the bulk matrix and finite volumes for the fractures, and present an approximation of the factors in a set of approximate block factorization preconditioners that accelerates convergence of iterative solvers applied to the resulting discrete system. Numerical tests on significant three-dimensional cases have assessed the properties of the proposed preconditioners.

Comparison of the response to geometrical complexity of methods for unstationary simulations in discrete fracture networks with conforming, polygonal, and non-matching grids

The aim of this study is to compare numerical methods for the simulation of single-phase flow and transport in fractured media, described here by means of the discrete fracture network (DFN) model. A Darcy problem is solved to compute the advective field, then used in a subsequent time-dependent transport-diffusion-reaction problem. The numerical schemes are benchmarked in terms of flexibility in handling geometrical complexity, mass conservation, and stability issues for advection-dominated flow regimes. To this end, two benchmark cases, along with an additional one from a previous work, have been specifically designed and are here proposed and investigated, representing some of the most critical issues encountered in DFN simulations.