Abstract
This work presents an extension of grid generation techniques for finite-volume discretizations of densitydriven flow in fractured porous media, in which fractures are considered as low-dimensional manifolds and are resolved by sides of grid elements. The proposed technique introduces additional degrees of freedom for the unknowns assigned to the fractures and thus allows to reconstruct jumps of the solution over a fracture. Through the concept of degenerated elements, the proposed technique can be used for arbitrary junctions of fractures but is sufficiently simple regarding the implementation and allows for the application of conventional numerical solvers. Numerical experiments presented at the end of the paper demonstrate the applicability of this technique in two and three dimensions for complicated fracture networks.