A computational method for the hydrodynamics of fractured‐porous media

A numerical method based on the boundary‐fitted finite difference method (BFDM) is presented in this paper. The boundaries are external (the boundary of the physical domain) and internal (which corresponds to the fracture network). The difference between this approach and the usual one lies in the inclusion of discrete fractures in the volume that represents the porous medium. The numerical model has been used in the prediction of the flow pattern in several internationally recognized verification cases and applied to the solution of hypothetical problems of interest to us in the field of nuclear waste repository modelling. The results obtained show that the numerical approach considered gives accurate and reliable predictions of the hydrodynamics of fractured‐porous media, thus justifying its use for the above‐mentioned studies.