Numerical modeling of radionuclide migration in water-saturated planar fracture: study of performance of bentonite in the far-field region

The analysis of radionuclide migration in fractured porous media is an important part of the safety assessment of a deep geologic disposal for high level radioactive wastes. In this work, numerical solution for simple geometry was developed to study radionuclide migration, including decay chain from a hypothetical repository, whereas the initial region of fracture is filled with bentonite which expanded from EBS (Engineered Barrier System). The following cases were considered: convective transport with constant velocity along the fracture, longitudinal hydrodynamic dispersion in the fracture along the fracture axis, molecular diffusion from fracture into rock matrix, sorption within rock matrix, sorption onto the surface of the fracture, radioactive decay, decay chain, and diffusion in bentonite extrusion region. For conservative analysis, the porous matrix adjacent to the fracture was subdivided into two different subdomains, each with different set of parameters and considering that the radionuclides were available for migration in the solubility limit, at fracture inlet, from the initial time. The partial differential equations that govern the physical system were discretized by finite differences, by using the Implicit Euler Method with forward scheme in the convective term. In this study, numerical simulation was performed for 100, 1000 and 10000 years, with and without bentonite extrusion, in order to compare the migration retardation obtained by bentonite located at the beginning of the fracture in saturated environment. The numerical simulation results showed the importance of extruded area in the far field region of the fractured host rock. (author)