Simulation of the steady-state flow in discrete fracture networks with non-conforming meshes and extended finite elements

In this paper a numerical method for the simulation of the steady-state fluid flow in discrete fracture networks is described. It is based on the use of non-conforming meshes, enrichment functions and an optimization procedure. The meshing process is performed on each fracture independently of the other fractures, i.e. without geometrical conformity at the intersections (traces). The slope discontinuities due to the flux exchange at the traces are then captured with the enrichment functions of the extended finite elements, and finally a functional is minimized by resorting to an optimization procedure. The method can be easily implemented for parallel computers being based on many small independent problems. In order to show the effectiveness of the method and the quality of the results, simulations of fluid flow in simple networks are illustrated.