An efficient multi-point flux approximation method for Discrete Fracture–Matrix simulations

We consider a control volume discretization with a multi-point flux approximation to model Discrete Fracture–Matrix systems for anisotropic and fractured porous media in two and three spatial dimensions. Inspired by a recently introduced approach based on a two-point flux approximation, we explicitl...

Full description

Saved in:
Bibliographic Details
Published in:Journal of computational physics 2012-05, Vol.231 (9), p.3784-3800
Main Authors: Sandve, T.H., Berre, I., Nordbotten, J.M.
Format: Article
Language:English
Subjects:
Anisotropy
Approximation
Computational techniques
DFM
Discretization
Exact sciences and technology
Flux
Fracture mechanics
Fractured porous media
Intersections
Mathematical analysis
Mathematical methods in physics
Mathematical models
MPFA
Physics
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:We consider a control volume discretization with a multi-point flux approximation to model Discrete Fracture–Matrix systems for anisotropic and fractured porous media in two and three spatial dimensions. Inspired by a recently introduced approach based on a two-point flux approximation, we explicitly account for the fractures by representing them as hybrid cells between the matrix cells. As well as simplifying the grid generation, our hybrid approach excludes small cells in the intersection of the fractures and hence avoids severe time-step restrictions associated with small cells. Excluding the small cells also reduces the condition number of the discretization matrix. For examples involving realistic anisotropy ratios in the permeability, numerical results show significant improvement compared to existing methods based on two-point flux approximations. We also investigate the hybrid method by studying the convergence rates for different apertures and fracture/matrix permeability ratios. Finally, the effect of removing the cells in the intersections of the fractures are studied. Together, these examples demonstrate the efficiency, flexibility and robustness of our new approach.
ISSN:0021-9991
1090-2716
DOI:10.1016/j.jcp.2012.01.023