Simultaneous resolution of reactive radioactive decay, non-isothermal flow, and migration with application to the performance assessment for HLW repositories

Abstract

Radioactive decay is an important subject to take into account when studying the thermo-hydro-dynamic behavior of the buffer clay material used in the containment of radioactive waste. The modern concepts for the multibarrier design of a repository of high level waste in deep geologic formations consider that once canisters have failed, the buffer clay material must ensure the retention and/or delay of radionuclides within the time framework given in the assessment studies. Within the clay buffer, different chemical species are retarded/fixed according to several physicochemical processes (ion exchange, surface complexation, precipitation, matrix diffusion, ...) but typical approaches do not consider the eventuality that radioactive species change their chemical nature (i.e. phase) thereby affecting their reactive behavior.

The radioactive decay of an element takes place independently of the phase (aqueous, solid or gaseous) to which it belongs. This means that, in terms of radionuclide fixation, some geochemical processes will be effective scavengers (for instance mineral precipitation of crystal growth) while others will not (for instance ion exchange and/or sorption).

In this contribution we present a reactive radioactive decay model of any number of chemical components including those that belong to decay series. The model, which is named FLOW-DECAY, also takes into account flow and isotopic migration and it has been applied considering a hypothetical model scenario provided by the project ENRESA 2000 and direct comparison with the results generated by the probabilistic code GoldSim. Results indicate that FLOW-DECAY may simulate the decay processes in a similar way that GoldSim, being the differences related to factors associated to code architecture.