Radionuclide Transport Calculations from High-level Long-lived Radioactive Waste Disposal in Deep Clayey Geologic Formation toward Adjacent Aquifers

In the context of high-level nuclear waste repository safety calculations, the modeling of radionuclide migration is of first importance. Three dimensional radionuclide transport calculations in geological repository need to describe objects of the meter scale embedded in geologic layer formations of kilometer extension. A complete and refined spatial description would end up with at least meshes of hundreds of millions to billions elements. The resolution of this kind of problem is today not reachable with classical computers due to resources limitations. Although parallelized computation appears as potential tool to handle multiscale calculations, to our knowledge no attempt have been yet performed. One emerging solution for repository safety calculations on very large cells meshes consists in using a domain decomposition approach linked to massive parallelized computer calculation. In this approach, the repository domain is divided in small elementary domains and transport calculation are performed independently on different processor for each elementary domain. Before to develop this possible solution, we performed some preliminary test in order to access the order of magnitude of cells needed to perform converged calculation on one elementary disposal domain and to check if Finite Volume (FV) based on Multi Point Flux Approximation (MPFA) spatial scheme or more classical Mixed Hybrid Finite Element (MHFE) spatial scheme were adapted for those calculations in highly heterogeneous porous media. Our preliminary results point out that MHFE and VF schemes applied on non-parralellepipedic hexahedral cells for flow and transport calculations in highly heterogeneous media gave satisfactory results. Nevertheless further investigations and additional calculations are needed in order to exhibit the mesh discretization level needed to perform converged calculations.