In Korea, two types of spent nuclear fuels (SNFs) are generated, pressurized light water reactor type (PWR) and pressurized heavy water reactor type (PHWR; CANDU), that differ greatly in size, decay heat, and radioactive characteristics. Technology development for the disposal of SNFs has mainly focused on PWR SNFs that are large in size and have extremely high decay heat and radioactivity. However, CANDU SNFs should be considered differently from PWR SNFs in deep geological disposal systems because their characteristics significantly differ from those of PWR SNFs in terms of their dimensions, number of SNF bundles, and handling systems in nuclear power plant sites. In this paper, after reviewing the status of the CANDU SNF disposal concept by Canada and Korea, concepts related to the direct geological disposal of CANDU SNFs were described, and two concepts were proposed based on the results of the development. The engineered barrier systems developed using these two concepts were comparatively analyzed in terms of disposal safety, disposal efficiency, and technical maturity. Based on the results of the comparative analyses, a vertical-type emplacement disposal concept was determined as a reference concept for the deep geological disposal of CANDU SNFs.

The safety of deep geological disposal systems has to be ensured to guarantee the isolation of radionuclides from human and related environments for over a million years. Over such a long timeframe, disposal systems can be influenced by climate change, leading to significant long-term impacts on the hydrogeological condition, including changes in temperature, precipitation and sea levels. These changes can affect groundwater flow, alter geochemical conditions, and directly/ indirectly impact the stability of the repository. Hence, it is essential to conduct a safety assessment that considers the long-term evolution induced by climate change. In this context, the Korea Atomic Energy Research Institute (KAERI) is developing the Adaptive Process-based total system performance assessment framework for a geological disposal system (APro). Currently, numerical modules for APro are under development to account for the longterm evolution that can influence groundwater flow and radionuclide transport in the far-field of the disposal system. This study focuses on the development of two numerical modules designed to model permafrost formation and buoyance force due to relative density changes. Permafrost is defined as a ground in which temperature remains below zero-isotherm (0°C) continuously for more than two consecutive years. In regions where permafrost forms, the relative permeability of porous media is significantly reduced. The changes in permeability due to permafrost formation are modelled by calculating the unfrozen fluid content within a porous medium. Meanwhile, buoyancy force can occur when there is a difference in density at the boundary of two distinct water groups, such as seawater (salt water) and freshwater. Sea level change associated with climate change can alter the boundary between seawater and freshwater, resulting in changes in groundwater flow. The buoyancy force due to relative density is modelled by adjusting concentration boundary conditions. Using the developed numerical modules, we evaluated the long-term evolution’s effects by analyzing radionuclide transport in the far-field of the disposal system. Incorporating permafrost and buoyancy force modelling into the APro framework will contribute valuable insights into the complex interactions between geological and climatic factors, enhancing our ability to ensure the secure isolation of radionuclides for extended periods.

Over the years, in the field of safety assessment of geological disposal system, system-level models have been widely employed, primarily due to considerations of computational efficiency and convenience. However, system-level models have their limitations when it comes to phenomenologically simulating the complex processes occurring within disposal systems, particularly when attempting to account for the coupled processes in the near-field. Therefore, this study investigates a machine learning-based methodology for incorporating phenomenological insights into system-level safety assessment models without compromising computational efficiency. The machine learning application targeted the calculation of waste degradation rates and the estimation of radionuclide flux around the deposition holes. To develop machine learning models for both degradation rates and radionuclide flux, key influencing factors or input parameters need to be identified. Subsequently, process models capable of computing degradation rates and radionuclide flux will be established. To facilitate the generation of machine learning data encompassing a wide range of input parameter combinations, Latin-hypercube sampling will be applied. Based on the predefined scenarios and input parameters, the machine learning models will generate time-series data for the degradation rates and radionuclide flux. The time-series data can subsequently be applied to the system-level safety assessment model as a time table format. The methodology presented in this study is expected to contribute to the enhancement of system-level safety assessment models when applied.

The effect of various physicochemical processes, such as seawater intrusion, on the performance of the engineered barrier should be closely analyzed to precisely assess the safety of high-level radioactive waste repository. In order to evaluate the impact of such processes on the performance of the engineered barrier, a thermal-hydrological-chemical model was developed by using COMSOL Multiphysics and PHREEQC. The coupling of two software was achieved through the application of a sequential non-iterative approach. Model verification was executed through a comparative analysis between the outcomes derived from the developed model and those obtained in prior investigations. Two data were in a good agreement, demonstrating the model is capable of simulating aqueous speciation, adsorption, precipitation, and dissolution. Using the developed model, the geochemical evolution of bentonite buffer under a general condition was simulated as a base case. The model domain consists of 0.5 m of bentonite and 49.5 m of granite. The uraninite (UO2) was assigned at the canister-bentonite interface as the potential source of uranium. Assuming the lifetime of canister as 1,000 years, the porewater mixing without uranium leakage was simulated for 1,000 years. After then, the uranium leakage through the dissolution of uraninite was initiated and simulated for additional 1,000 years. In the base case model, where the porewater mixing between the bentonite and granite was the only considered process, the gypsum tended to dissolve throughout the bentonite, while it precipitated in the vicinity of bentonite-granite boundary. However, the precipitation and dissolution of gypsum only showed a limited effect on the performance of the bentonite. Due to the low solubility of uraninite in the reduced environment, only infinitesimal amounts of uranium dissolved and transported through the bentonite. Additional cases considering various environmental processes, such as seawater or cement porewater intrusion, will be further investigated.

Effective containment and disposal of high-level radioactive waste is critical to ensure long-term environmental and human safety. Especially bentonite, which is widely used as a buffer material due to its favorable characteristics such as swelling ability and low permeability, plays an important role in preventing the migration of radioactive waste into the surrounding environment. However, the long-term performance of bentonite buffer remains an area of ongoing investigation, with particular attention focused on erosion mechanisms induced by swelling and groundwater flow. The erosion of the bentonite buffer can significantly impact the integrity of buffer and lead to the formation of colloids, which could potentially facilitate the transport of radionuclides through groundwater. Therefore, quantification of bentonite buffer erosion based on an understanding of the underlying mechanisms and factors that influence bentonite buffer erosion, is essential for the safety assessment of high-level radioactive waste repositories. In this study, we aimed to develop a bentonite buffer erosion model using the Adaptive Processbased total system performance assessment framework for a geological disposal system (APro) proposed by the Korea Atomic Energy Research Institute (KAERI). The impact of bentonite erosion on performance assessment can be broadly divided into bentonite property degradation by the penetration of the bentonite buffer into rock fractures and the formation of pseudocolloids. To simulate this phenomenon, Two-region model based on a dynamic bentonite diffusion model is adopted, which can quantify the extent of bentonite intrusion and loss by erosion. Using this Tworegion model, a numerical model was developed to simulate the degradation of bentonite properties based on the amount of bentonite intrusion, as well as to simulate the migration of pseudocolloids in the near-field by deriving the amount of pseudocolloid production based on the loss of bentonite and the sorption rate of radionuclides. To check the applicability of the developed numerical model, preliminary analysis was performed for the effect of bentonite erosion in terms of process-based performance assessment. It is anticipated that this comprehensive model developed in this study will contribute to the accurate and reliable assessment of the long-term performance and safety of high-level radioactive waste repositories.

The post-closure safety assessment of a repository is typically conducted over an extensive timescale from ten thousand to a million years. Considering that biosphere ecosystems may undergo significant changes over such lengthy periods, it is essential to incorporate the long-term evolution of the biosphere into the safety assessment. Climate change and landscape development are identified as critical drivers with the potential to impact the hydrogeological and hydrogeochemical characteristics of the biosphere. These changes can subsequently alter the migration patterns of radionuclides through the biosphere and influence human exposure doses. Therefore, this study formulates scenarios within the context of long-term biosphere evolution. We examine biosphere assessment processes employed in other countries and conduct a comparative study on scenario conditions. For example, biosphere assessment in Finland has identified sea-level changes and land-use alterations as significant factors in the long-term evolution of the biosphere. These factors are linked to Features, Events, and Processes (FEPs) associated with climate change and human activities. Sea-level changes are related to FEPs regarding climate change, land uplift, and shoreline displacement, while land-use changes are based on human activity-related FEPs (e.g., crop type, livestock and forest management, well construction, and demographics). Based on the literature review, this study has configured long-term evolution scenarios for the safety assessment of a deep geological repository for spent fuels.

Safety assessments for geological disposal systems extend over tens of thousands of years, taking into account the radiotoxicity decay period of spent nuclear fuel. During this extensive period, the biosphere experiences multiple glacial cycles, and fluctuations in seawater amounts, attributed to the formation and melting of glaciers, lead to global sea level changes known as eustacy. These sea level changes can directly influence the land-sea interface and groundwater flow dynamics, consequently affecting the pathways of radionuclide transport - an essential element of dose assessment. Therefore, this study aims to investigate how glacial cycles and sea level changes impact radionuclide transport within geological disposal systems, especially in the biosphere. To achieve this objective, we obtained climate evolution data including sea level changes for the Korean Peninsula over a 200,000-years, simulated by a General Circulation Model (GCM). These data were then employed to predict site and hydrology evolutions. The study site was conceptualized biosphere of Artificial Disposal System (ADioS), and we utilized the Soil and Water Assessment Tool (SWAT) to simulate hydrological evolution. These datasets, encompassing climate, site, and hydrology evolution, were collectively employed as inputs for the biosphere module of Adaptive Process-Based Total System Performance Assessment Framework (APro). Subsequently, the APro’s biosphere module calculated radionuclide transport in groundwater flow and its release into surface water bodies, considering the influences of glacial cycles and sea level changes. The results show that hydrologic changes due to sea level change are relatively minor, while the impact of sea level change on groundwater flow and discharge is significant. Additionally, we identified that among the water bodies within ADioS, including rivers, lakes, and oceans, the ocean exhibits the most substantial radionuclide outflow throughout the entire period. The spatiotemporal distributions of radionuclides computed within APro will be further processed into a grid format and used as input for the dose assessment module. Through this study, it was possible to determine the impact of long-term glacial cycles and sea level changes on radionuclide transport. Additionally, this module can serve as a valuable tool for providing the spatiotemporal variability of radionuclides required for enhanced dose assessments.

Understanding the long-term geochemical evolution of engineered barrier system is crucial for conducting safety assessment in high-level radioactive waste disposal repository. One critical scenario to consider is the intrusion of seawater into the engineered barrier system, which may occur due to global sea level rise. Seawater is characterized by its high ionic strength and abundant dissolved cations, including Na, K, and Mg. When seawater infiltrates an engineered barrier, such dissolved cations displace interlayer cations within the montmorillonite and affect to precipitation/ dissolution of accessory minerals in bentonite buffer. These geochemical reactions change the porewater chemistry of bentonite buffer and influence the reactive transport of radionuclides when it leaked from the canister. In this study, the adaptive process-based total system performance assessment framework (APro), developed by the Korea Atomic Energy Research Institute, was utilized to simulate the geochemical evolution of engineered barrier system resulting from seawater intrusion. Here, the APro simulated the geochemical evolution in bentonite porewater and mineral composition by considering various geochemical reactions such as mineral precipitation/dissolution, temperature, redox processes, cation exchange, and surface complexation mechanisms. The simulation results showed that the seawater intrusion led to the dissolution of gypsum and partial precipitation of calcite, dolomite, and siderite within the engineered barrier system. Additionally, the composition of interlayer cation in montmorillonite was changed, with an increase in Na, K, and Mg and a decrease in Ca, because the concentrations of Na, K, and Mg in seawater were 2-10 times higher than those in the initial bentonite porewater. Further studies will evaluate the geochemical sorption and transport of leaked uranium-238 and iodine-129 by applying TDB-based sorption model.

The Korea Atomic Energy Research Institute (KAERI) is currently developing a process-based performance assessment model known as APro. Distinguished from the previous system-level safety assessment model developed by KAERI, APro exhibits the capacity to encompass a threedimensional biosphere domain, evolving over the long term. In this study, we elucidate the methodology employed in developing the dose assessment module of APro and present the module’s functionalities. The procedural steps underlying radiation dose calculations within the APro framework can be succinctly outlined as follows: 1) Definition of a landscape model, utilizing information derived from a specified snapshot period provided by the APro biosphere transport module; 2) Generation of unit biotope objects spanning the landscape; 3) Evaluation of radionuclide transfer within the soil medium; 4) Calculation of activity concentration for flora and fauna groups; 5) Assessment of the distribution of effective dose among representative human groups; 6) Progressing through successive time steps. The APro dose calculation module exhibits notable capabilities that encompass: 1) Accounting for radionuclide decay and ingrowth; 2) Facilitating transfer through unsaturated porous media; 3) Considering sorption effects; 4) Addressing the inheritance of radioactivity between various landscape models; 5) Offering customizable ecosystem parameters; 6) Providing flexibility for user-defined exposure pathways. Leveraging these functionalities of the dose assessment module, APro is proficient in evaluating the distribution of radiological doses and associated risks for representative population groups, all while accounting for the dynamic, long-term evolution of the biosphere, including alterations in land cover.

Conducting a TSPA (Total System Performance Assessment) of the entire spent nuclear fuel disposal system, which includes thousands of disposal holes and their geological surroundings over many thousands of years, is a challenging task. Typically, the TSPA relies on significant efforts involving numerous parts and finite elements, making it computationally demanding. To streamline this process and enhance efficiency, our study introduces a surrogate model built upon the widely recognized U-network machine learning framework. This surrogate model serves as a bridge, correcting the results from a detailed numerical model with a large number of small-sized elements into a simplified one with fewer and large-sized elements. This approach will significantly cut down on computation time while preserving accuracy comparable to those achieved through the detailed numerical model.

APro, a process-based total system performance assessment (TSPA) tool for a geological disposal system, has a framework for simulating the radionuclide transport affected by thermal, hydraulic, mechanical or geochemical changes occurred in the disposal system. APro aims to be applied for the TSPA to long-term (> 100,000) evolution scenarios in real-world repository having more than 10,000 boreholes. In this large-scale TSPA, it is important not only to develop a high-performance numerical approach, but also to apply an efficient post-processing approach to massive spatiotemporal data. The post-processing refers to validating numerical analysis results, analyzing and evaluating target systems through data processing or visualization. Since APro uses COMSOL interface, the postprocessing function in COMSOL can be used. However, when the data size increases due to largescale numerical analysis, the time for the COMSOL post-processing increases, resulting in a problem that the analysis and evaluation are not performed effectively. In this case, it is possible to extract necessary data using the COMSOL exporting function and importing it into an external postprocessing program for the analysis and evaluation. In this study, the efficiency of external post-processing with extracted data from COMSOL was reviewed. And, we derived a proper data extraction approach (format and structure) that can increase efficiency of external post-processing.

APro, developed in KAERI for the process-based total system performance assessment (TSPA) of deep geological disposal systems, performs finite element method (FEM)-based multiphysics analysis. In the FEM-based analysis, the mesh element quality influences the numerical solution accuracy, memory requirement, and computation time. Therefore, an appropriate mesh structure should be constructed before the mesh stability analysis to achieve an accurate and efficient process-based TSPA. A generic reference case of DECOVALEX-2023 Task F, which has been proposed for simulating stationary groundwater flow and time-dependent conservative transport of two tracers, was used in this study for mesh stability analysis. The relative differences in tracer concentration varying mesh structures were determined by comparing with the results for the finest mesh structure. For calculation efficiency, the memory requirements and computation time were compared. Based on the mesh stability analysis, an approach based on adaptive mesh refinement was developed to resolve the error in the early stage of the simulation time-period. It was observed that the relative difference in the tracer concentration significantly decreased with high calculation efficiency.