A numerical method for two-phase flow in fractured porous media with non-matching grids

We propose a novel computational method for the efficient simulation of two-phase flow in fractured porous media. Instead of refining the grid to capture the flow along the faults or fractures, we represent the latter as immersed interfaces, using a reduced model for the flow and suitable coupling conditions. We allow for non matching grids between the porous matrix and the fractures to increase the flexibility of the method in realistic cases. We employ the extended finite element method for the Darcy problem and a finite volume method that is able to handle cut cells and matrix-fracture interactions for the saturation equation. Moreover, we address through numerical experiments the problem of the choice of a suitable numerical flux in the case of a discontinuous flux function at the interface between the fracture and the porous matrix. A wrong approximate solution of the Riemann problem can yield unphysical solutions even in simple cases.