This script illustrates the basic use of the DFM module. We create a simple network of two fractures in a Cartesian 2D grid and add the fracuters to the grid with a hybrid representation. Transmissibilities are then computed with a two-point flux expression, and solve a tracer transport problem.
A somewhat advanced example, where a fracture network is gridded with triangles using a fairly simple approach. The flow equations are then discretized with both two- and multi-point flux expressions, and a two-phase transport problem is solved.
The module provides control volume discretizations (two- and multi-point flux approximations) of fractured media by representing fractures as lower-dimensional objects that still are represented as cells in the computational grid. The discretization is described by Karimi-Fard et al. [1] and Sandve et al. [2]. The approach assumes that the fractures are represented as faces in the grid, thus the module also provide some utility functions for triangulating fractured domains. Finally, the transport solvers (both explicit and implicit) from the incomp module are modified to account for fracture transport. MRST files that are modified have a suffix _DFM to distinguish them from the original files.
The module is mainly developed by Tor Harald Sandve (at IRIS from Jan 2013), minor modifications are introduced by Eirik Keilegavlen (University of Bergen). Some utility files are partly written at Uni Research. Technical details can be found in the implementation note.
T.H. Sandve, and I. Berre, J. M. Nordbotten. An efficient multi-point flux approximation based approach for discrete fracture matrix simulation. Journal of Computational Physics, 231(9), 3784-3800, 2012. DOI:10.1016/j.jcp.2012.01.023.
Published December 8, 2016