Analysis of conforming, non-matching, and polygonal methods for Darcy and advection-diffusion-reaction simulations in discrete fracture networks
Source code and examples for the paper
"Analysis of conforming, non-matching, and polygonal methods for Darcy and advection-diffusion-reaction simulations in discrete fracture networks" by Andrea Borio, Alessio Fumagalli, Stefano Scialò. See arXiv pre-print.
Runscripts for all test cases of the work available here.
Note that you may have to revert to an older version of PorePy to run the examples.
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, three benchmark cases have been specifically designed and are here proposed, representing some of the most critical issues encountered in DFN simulations.
If you use this work in your research, we ask you to cite the following publication arXiv:1907.12514 [math.NA].
If you want to run the code you need to install PorePy and revert to commit f82a91c1b6ef83b6954ecd26723f8b7443880891
Newer versions of PorePy may not be compatible with this repository.
See license.