Published on

High Performance Computing as a prognostic tool for Deep Geological Repositories for spent nuclear fuel



  • Name
    Hedieh Ebrahimi

The safety of Deep Geological Repositories (DGR) for spent nuclear fuel is a highly important issue, and should be assessed using the Best Available Technology (BAT). With this motivation in mind, Amphos 21 has recently developed a numerical tool, denoted as iDP (interface between DarcyTools and PFLOTRAN; Molinero et al., 2016), which provides aJuqueen unique High Performance Computing(HPC)-based framework for the simulation of multicomponent reactive transport problems in fractured crystalline media. iDP combines two powerful state-of-the-art software: DarcyTools (Svensson and Ferry, 2014), which is used to compute groundwater flow in sparsely fractured rocks using a stochastic continuum approach, and the massively parallel reactive transport code PFLOTRAN (Hammond and Lichtner, 2010). iDP has been tested during the JUQUEEN Extreme Scaling Workshop 2016 (Brömmel et al., 2016) showing excellent scalability up to 131,072 processes (8 racks).

In cooperation with the Forschungszentrum Jülich and the Swedish Nuclear Fuel and Waste Management Company (SKB), Amphos 21 has used iDP to carry out two challenging reactive transport simulations. The first set of calculations has focused on evaluating the potential ingress of oxygenated glacial-meltwater into initially anoxic fractured media. This is an issue of concern for the safety of deep geological Pocos de Caldasrepositories for spent nuclear fuel as oxidizing conditions would increase the solubility and mobility of many radionuclides and O2-rich waters would corrode canister materials containing the nuclear waste.
The reactive transport models have pointed out the strong buffering capacity of a typical crystalline medium (e.g. Forsmark in Sweden), with the infiltration of oxygen being limited to few highly transmissive deformation zones. These findings are in agreement with data from the natural analogue site of Poços de Caldas (Brazil) (Nordstrom et al., 1992). For more details about these reactive transport simulations, the reader is referred to a paper recently published in Journal of Contaminant Hydrology (Trinchero et al, 2017a).

Another important issue for the safety of a DGR is related to the potential of the host rock to retain those radionuclides that are released due to a possible failure of the engineered barriers.Total CS This potential, in turn, depends on two coupled processes: diffusion into the rock matrix and sorption onto the mineral surfaces of the available reactive minerals. A typical modelling assumption is to use an equivalent (homogeneous) description of sorption sites. This assumption is here assessed using a very detailed micro-scale model of the rock matrix, which is built upon a numerical grid consisting of more than 50 million elements. The results of the calculations, in which the distribution of sorption sites is described at the grain scale, show that mineralogical heterogeneity has a significant impact on the retention of radionuclides, namely cesium. Not considering that sorption sites are sparse and unevenly distributed results in an important overestimation of the rock retention capacity. Further details can be found in a paper recently published in Transport in Porous Media (Trinchero et al., 2017b).


This work has been funded by the Swedish Nuclear Fuel and Waste Management Company (SKB). The authors gratefully acknowledge the computing time granted by the JARA-HPC Vergabegremium and provided on the JARA-HPC Partition part of the supercomputer JUQUEEN at Forschungszentrum Jülich. The authors also thank Ignasi Puigdomenech, Björn Gylling, Guido Deissmann, Urban Svensson, Glenn Hammond, Peter Lichtner and the PFLOTRAN development group for their help and valuable suggestions during the project.


Brömmel, D., Wylie, B. J., Frings, W., 2016. JUQUEEN extreme scaling workshop 2016. Tech. Rep. FZJ-JSC-IB-2016-01, Jülich Supercomputing Centre.

Hammond, G. E., Lichtner, P. C., 2010. Field-scale model for the natural attenuation of uranium at the Hanford 300 area using high-performance computing. Water Resources Research 46 (9), W09527, 131.

Molinero, J., Trinchero, P., Ebrahimi, H., de Vries, L., Luna, M., Svensson, U., Lichtner, P. 2016. The BRIDGE project: development, testing and application of a High Performance Computing framework for reactive transport modelling in crystalline roks (iDP). R-15-17, Svensk Kärnbränslehantering AB, Stockholm, Sweden.

Nordstrom, D., McNutt, R., Puigdomènech, I., Smellie, J. A., Wolf, M. 1992. Ground water chemistry and geochemical modeling of water-rock interactions at the Osamu Utsumi mine and the Morro do Ferro analogue study sites, Poços de Caldas, Minas Gerais, Brazil. Journal of Geochemical Exploration 45 (1-3), 249-287

Svensson, U., Ferry, M., 2014. Darcytools: A computer code for hydrogeological analysis of nuclear waste repositories in fractured rock. Journal of Applied Mathematics and Physics 2 (06), 365.

Trinchero, P., Molinero, J., Deissmann, G., Svensson, U., Gylling, B., Ebrahimi, H., Hammond, G., Bosbach, D., Puigdomenech, I., 2017b. Implications of grain-scale mineralogical heterogeneity for radionuclide transport in fractured media. Transport in Porous Media 116: 73. doi:10.1007/s11242-016-0765-0

Trinchero, P., Puigdomenech, I., Molinero, J., Ebrahimi, H., Gylling, B., Svensson, U., Bosbach, D., Deissmann, G., 2017a. Continuum-based DFN-consistent numerical framework for the simulation of oxygen infiltration into fractured crystalline rocks. Journal of Contaminant Hydrology DOI: 10.1016/j.jconhyd.2017.04.001