Various linear system solvers with multi-physics analysis schemes are compared focusing on the near-field region considering thermal-hydraulic-chemical (THC) coupled multi-physics phenomena. APro, developed at KAERI for total system performance assessment (TSPA), performs a finite element analysis with COMSOL, for which the various combinations of linear system solvers and multi-physics analysis schemes should to be compared. The KBS-3 type disposal system proposed by Sweden is set as the target system and the near-field region, which accounts for most of the computational burden is considered. For comparison of numerical analysis methods, the computing time and memory requirement are the main concerns and thus the simulation time is set up to one year. With a single deposition hole problem, PARDISO and GMRESSSOR are selected as representative direct and iterative solvers respectively. The performance of representative linear system solvers is then examined through a problem with an increasing number of deposition holes and the GMRES-SSOR solver with a segregated scheme shows the best performance with respect to the computing time and memory requirement. The results of the comparative analysis are expected to provide a good guideline to choose better numerical analysis methods for TSPA.

Identifying plausible scenarios is necessary to evaluate the performance of the repository reliably over a very long period. All features, events, and processes (FEPs) expected in the repository should be comprehensively well-defined and structured into scenarios based on the relation analysis. A platform for the FEP DB management and relation analysis is needed to facilitate the efficient composition of the scenarios. For this purpose, the CYPRUS program was developed, but abandoned due to suspended FEPs and scenario research. Thus, it became necessary to build a new easy-tomaintain platform that inherits the legacy of CYPRUS and reflects the latest research. The data structure and user interface configuration were derived to develop a new platform. The new platform provides extensive data such as the assessment context, the FEP DB, the interaction between FEP contents, the relevance to other project FEPs, the influence on performance, the scenarios for the TSPA, the AMF, and the PA Data. The platform displays the long-term evolution FEPs developed by KAERI, the international and major project FEPs in table format. The correlation between FEP items is composed of a detailed interaction matrix and visualized as the chord diagram or arc diagram. The relevance and linkages between the project FEP items are mapped and presented in the form of network diagrams and network tables. The platform designed in this study will be used to manage the FEP DB, analyze and visualize the relationship between the FEP and scenarios, and finally construct the performance assessment scenarios. It is expected that the platform itself will be used as a part of the knowledge management system and facilitate efficient collaboration and knowledge exchange among experts.

Excavation Damaged Zone (EDZ) is created by the excavation of deposition holes and disposal tunnels at high-level radioactive waste repository that causes macro- and micro-fracturing in the surrounding rock. Since EDZ can significantly increase the hydraulic transmissivity in the rock and act as a major pathway of leaked radionuclides, consideration of EDZ in terms of safety assessment is very important. Moreover, long-term stress changes such as stress redistribution due to excavation of nearby deposition holes and disposal tunnels, thermal stress due to temperature rise, effective stress change due to pore pressure change, and swelling pressure of bentonite buffer can increase EDZ size and change in thermal-hydraulic-mechanical properties, and consequently, it can affect the transport of radionuclides. Therefore, in order to analyze the effect of long-term evolution of EDZ on radionuclide transport, it is essential to conduct numerical analysis considering the coupled Thermal-Hydraulic- Mechanical (THM) behavior in EDZ. In order to simulate the behavior of EDZ, coupled THM model was developed using the Adaptive Process-based total system performance assessment framework for a geological disposal system (APro) proposed by the Korea Atomic Energy Research Institute (KAERI). The concept of damage was introduced to demonstrate the jointed rock as a continuous medium. Among several damage models, Mazars damage model was applied in this study. Mazars damage model is the most well-known model for concrete which has similar behavior with rock as brittle material, and the input data of the model can be easily obtained through laboratory testing. If damage occurs due to the influence of thermal-hydraulic-mechanical coupled behavior at the bedrock, the properties change according to the degree of damage, and as a result, the migration of the radionuclide is affected. Based on this conceptual model, radionuclide transport model in the near field considering the long-term evolution of EDZ was developed. To investigate the effect of EDZ in terms of process-based performance assessment, the modeling results with and without EDZ were compared. Finally, by simulating the coupled THM behavior of EDZ with damage model, the effect of long-term evolution of EDZ on radionuclide transport was investigated.

In the geological disposal system whose host rock is crystalline rock, fractures play a significant role in the safety assessment as they are the main pathway of the radionuclide migration. From the perspective of long-term safety assessment, the properties of fractures can be changed by tectonic movement such as earthquake, uplift, etc. In general, methods for simulating fractures include Discrete Fracture Network (DFN), which directly simulates the fracture surface, and Equivalent Continuous Porous Media (ECPM), which is equivalent to the ratio of the fractures in a certain rock volume. DFN is generally appropriate for deterministic fractures with large scale and high flow velocity, but ECPM may be more appropriate for small scale and sporadically distributed stochastic fractures because the flow velocity is slow and thus the rock matrix diffusion needs to be considered. In fact, several commercial software, such as FracMan, are already in use to convert DFN to ECPM. However, in order to consider the change in properties of fractures due to tectonic movement in the long-term safety assessment, a model that converts DFN to ECPM needs to be modularized and embedded into the safety assessment model. In this study, therefore, an in-house MATLAB code was developed to convert DFN to ECPM, which can be used as a submodule. The algorithm of converting from DFN to ECPM basically followed the Oda’s method. As the first step of the algorithm, in order to obtain the volume ratio of the fracture in a certain mesh element, the cross-sectional area of the fracture and the mesh element was calculated. Then, porosities of each mesh element were calculated as the volume fraction of fractures passing through the mesh element. Based on the Oda’s method, the permeability tensors of each mesh element were calculated by using an empirical fracture tensor which is weighted by the cross-sectional area and transmissivity of each fracture. Finally, the newly developed module was verified by a benchmark test, in which the ECPM results converted from a certain DFN data by using the numerical module developed in this study were compared with those by using FracMan. The newly developed module will be installed in the process-based total system performance assessment framework (APro) being developed by KAERI.

The analysis of uranium migration is crucial for the accurate safety assessment of high-level radioactive waste (HLW) repository. Previous studies showed that the migration of the uranium can be affected by various physical and chemical processes, such as groundwater flow, heat transfer, sorption/ desorption and, precipitation/dissolution. Therefore, a coupled Thermal-Hydrological-Chemical (THC) model is required to accurately simulate the uranium migration near the HLW repository. In this study, COMSOL-PHREEQC coupled model was used to simulate the uranium migration. In the model, groundwater flow, heat transfer, and non-reactive solute transport were calculated by COMSOL, and geo-chemical reaction was calculated by PHREEQC. Sorption was primarily considered as geo-chemical reaction in the model, using the concept of two-site protolysis nonelctrostatic surface complexation and cation exchange (2 SP NE SC/CE). A modified operator splitting method was used to couple the results of COMSOL and PHREEQC. Three benchmarks were done to assess the accuracy of the model: 1) 1D transport and cation exchange model, 2) cesium transport in the column experiment done by Steefel et al. (2002), and 3) the batch sorption experiment done by Fernandes et al. (2012), and Bradbury and Baeyens (2009). Three benchmark results showed reliable matching with results from the previous studies. After the validation, uranium 1D transport simulation on arbitrary porewater condition was conducted. From the results, the evolution of the uranium front with sequentially saturating sites was observed. Due to the limitation of operator splitting method, time step effect was observed, which caused the uranium to sorbed at further sites then it should. For further study, 3 main tasks were proposed. First, precipitation/ dissolution will be added to the reaction part. Second, multiphase flow will be considered instead of single phase Darcy flow. Last, the effect of redox potential will be considered.

It is expected that around 576,000 bundles of CANDU spent nuclear fuels (SNF) will be generated from the four CANDU reactors located at the Wolsong site, according to the 2nd National Plan for the management of High-Level radioactive Waste (HLW). The CANDU SNFs are currently stored at the dry storage facilities at the Wolsong site. The authors proposed KRS+ geological disposal system consisting of two different concepts, Swedish KBS-3V type and Canadian NWMO type, for the final management of CANDU SNF. Both the concepts were designed based on the geological data obtained from the KURT (KAERI Underground Research Tunnel). The NWMO type is an in-room horizontal placement method. In this study, we try to determine the reference concept among the two proposed concepts at 500 meters below the ground surface. Assuming 10,000 tU of CANDU SNF and the KURT site, we design two engineered barrier systems, that is disposal canisters and buffers. The copper disposal canister is designed with a copper thickness of 10 mm based on a cold spray coating technique for both the disposal concepts. The domestic Ca-bentonite is used for the compact bentonite buffer with dry density of 1.6 g/cm3. Two concepts are compared in terms of safety, economics of the engineered barriers, and environment-friendliness. Because the same amounts of CANDU SNF are disposed of at the same depth, the differences in the disposal area are neglected. For the comparison in terms of safety, the corrosion lifetimes of the disposal canisters of two disposal systems are quantitatively calculated, and the capacities for retarding radionuclide releases of the compacted bentonite buffers are assessed. A computer tool developed by the authors is used in order to assess the lifetime of a disposal canister. In this study, the case that corrosion of a copper canister by sulfide from groundwater through intact buffer is analyzed. The sulfide concentration in groundwater is assumed to be 3 ppm. The most important safety function of buffer is to retard the radionuclide release. Twelve long-lived radionuclides are selected to compare the capacities for retarding the radionuclide transport through the buffer using an analytical solution. The retention time by an engineered barrier consisting of a disposal canister and a buffer is compared with twenty times the half-life of each radionuclide for both the disposal systems. The selected reference concept will be compared with the alternative geological concepts through a further study.

With the increase of temporarily-stored radioactive waste in Korea, the disposal of radioactive waste in a deep geological repository, which is located in crystalline rock at a depth of hundreds of meters below the ground level, has received great attention nowadays. To ensure the permanent isolation of radionuclides from the human and surrounding ecosystems, the safety assessment for the high-level radioactive waste disposal facilities is essential. For the reliable safety assessment of fractured rock, it is especially important to input proper hydraulic properties of fractures such as aperture and hydraulic conductivity, which can directly affect the fluid flow and radionuclide transport. Meanwhile, it has become important to consider sudden fault behavior caused by an earthquake with the recent occurrence of high-intensity earthquakes in the Korean Peninsula. The sudden fault behavior can induce the changes of the hydraulic properties of fractures. Since the changes of the hydraulic properties directly affects to the radionuclide transport in the fractured rock, it is important to estimate the effect of earthquake-induced stress change on hydraulic properties of fractures in the perspective of long-term safety assessment. In this study, the effect of an earthquake on the hydraulic properties of fractures was explored by a numerical approach. The static Coulomb stress change after the earthquake was calculated using software ‘Coulomb 3’ developed by United States Geological Survey (USGS) with the assumption for several mechanical properties such as Young’s modulus, Poisson’s ratio and effective coefficient of friction. The final stress after earthquake occurrence was calculated as the sum of the initial stress and the stress change. Thereafter, the normalized transmissivity of fracture after the earthquake was calculated using the final stress from the stress-transmissivity relationship. Using the methodology for calculating fracture transmissivity change induced by the earthquake developed in this study, the effect of several factors, such as the earthquake magnitude and the distance between fracture and epicenter, was additionally explored. The newly developed methodology will be applied to the processbased total system performance assessment framework (APro) being developed by KAERI, and this study is expected to be helpful for the safety assessment considering long-term evolution phenomena including earthquakes.

Montmorillonite plays a key role in engineered barrier systems in the high-level radioactive waste repository because of its large sorption capacity and high swelling pressure. However, the sorption capacity of montmorillonite can be largely varied dependent on the surrounding environments. This study conducted the batch simulation for U(VI) sorption on Na-montmorillonite by utilizing the cation exchange and surface complexation coupled (2SP-NE-SC/CE) model and evaluated the effects of physicochemical properties (i.e., pH, temperature, competing cations, U(VI) concentration, and carbonate species) on U(VI) sorption. The simulation demonstrated that the U(VI) sorption was affected by physicochemical properties: the pH and temperature relate to aqueous U(VI) speciation, the competing cations relate to the cation exchange process and selectivity, the U(VI) concentration relates to saturation at sorption sites. For example, the Kd (L kg−1) of Na-montmorillonite represented the largest values of 2.7×105 L kg−1 at neutral pH condition and had significantly decreased at acidic pH<3, showing non-linear and diverse U(VI) sorption at the ranged pH from 2 to 11. Additionally, the U(VI) sorption on montmorillonite significantly decreased in presence of carbonate species. The U(VI) sorption for long-term in actual porewater chemistry and temperature of high-level radioactive waste repository represented that the sorption capacity of Na-montmorillonite was affected by various external properties such as concentration of competing cation, temperature, pH, and carbonate species. These results indicate that geochemical sorption capacity of bentonite should be evaluated by considering both geological and aquifer environments in the high-level radioactive waste repository.

The safety assessment of a geological disposal system is performed over a period of hundreds of thousands of years, during which the activity of radionuclides in spent nuclear fuel decreases to natural radioactivity levels. During this period, the biosphere also experiences the long-term evolution of the surface environment including climate, terrain, and ecosystem changes. These changes cause changes in the water balance, which in turn change the pathways of radionuclides in the subsurface. Therefore, it is essential to consider these long-term changes in the surface environment for a reasonable biosphere safety assessment. For this purpose, this study developed the biosphere assessment module considering the long-term evolution of the surface environment, as a sub-module of APro (Adaptive process-based total system performance assessment framework). As a preceding study, the biosphere assessment module was previously developed using COMSOL for hydraulic and radionuclide transport processes, to simulate the pathway of radionuclides traveling from the shallow aquifer to the surface water body and soil. To consider the long-term evolution of the surface environment, the previous module needed to be improved to apply different water balances as boundary conditions of the module at each snapshot, which is a sub-time period divided based on the surface evolution data. To this end, this study utilized SWAT (Soil and Water Assessment Tool) which calculates the water balance using the surface environmental data including climate, terrain, land cover, and soil type. Conceptually, SWAT calculated annual water balance considering surface environmental changes, and certain components (i.e., groundwater recharge and hydraulic head of water bodies) of water balance were transferred to COMSOL as external data to simulate the pathway of radionuclide transport and spatio-temporal variability of radionuclides. At the current stage, the biosphere computational module has been developed to correspond to its conceptual model, and we plan to further test the applicability of the module using different surface environmental data.

Since it takes hundreds of thousands of years for the radiotoxicity of spent nuclear fuel to decrease to natural levels, interactions between each repository barrier, climate change, and geological evolutions are inevitable. These processes should be defined as the long-term evolution FEPs and considered in the performance assessment to ensure the long-term safety of the disposal system. The literature survey on geological characteristics and history of the Korean peninsula was conducted, and the list of A-KRS-FEPs which are directly or indirectly related to long-term evolutions was identified in this study. The ice age and geological change are the capital phenomena considered in the exceedingly long-term evolution before/after climate change. The historical data on ice sheets and permafrost were analyzed to investigate the effects of the ice ages on the Korean peninsula. The sealevel changes were investigated based on the research on the coastal terrace to identify the impact on uplift and shoreline change accompanying the ice age. Also, the survey on the geological history data was conducted from the perspective of tectonic activity, metamorphism, igneous activity, and seismic activities to consider the geodynamic evolution of the Korean peninsula. As results, it was suggested that 14 FEPs were directly related to climate change, 18 FEPs were directly related to geological evolution, and 47 FEPs were indirectly relevant to long-term geodynamics. The consent-based FEPs and scenarios for the long-term evolution will be developed shortly, including most of the critical long-term evolution phenomena defined in this study and which are highly probable in domestic disposal conditions. The evaluation and verification of the APro system for long-term safety will accomplish using these FEPs and scenarios.

An objective of a safety assessment for geological disposal is to evaluate the radiological impact by radionuclides release from radioactive wastes. Computational estimation of all radionuclides transport in the disposal system, however, is not neccessary because some radionuclides has negligible effect on radiological doses. For this reason, prioritization of radionuclides list is preceded before the safety assessment. The Korea Atomic Energy Research Institue (KAERI) has assessed the long-term safety of a disposal system for spent nculear fuels. Currently, thirty eight radionuclides and twenty three elements are considered in the safety assessment activity of the KAERI. Nevertheless, a screening process for radionulides selection has not been articulated yet. In this study, we reviewed radionuclides selection process in forign countries to re-establish screening criteria for the KAERI’s radionuclides list. Screeing models of the Swedish Nuclear Fuel and Waste Management Company (SKB), the Deparment of Eenrgy (US DOE), and the Japan Nuclear Cycle Development Istitute (JNC) were compared. We found that each country developed different screening model depending on scenarios of radionuclides release. Nonetheless, there were common properties that determines the importance of radionuclides. These properties for radionuclides include halflife, radiotoxicity (or specific activity), and mobility in underground medium. Based on the review results, we proposed radionuclides selection process to prioritize the importance of radionucldies in the KAERI safety assessment.

APro, a modularized framework of the process-based total system performance assessment, has been developed by KAERI to simulate the radionuclide transport in geological disposal system considering multi-physics phenomena. However, the target problem including more than 10,000 boreholes and over 100,000 years of simulation time is computationally challenging to deal with numerical solvers provided by COMSOL Multiphysics constituting APro. To alleviate the computational burden, machine learning (ML) techniques have been studied to develop a surrogate model replacing the heavy computation part. In recent studies, attempts have been made to integrate the knowledge of physics and numerical methods into the ML model for partial differential equations (PDEs). Unlike conventional ML approaches solely relying on data-driven method, the integration can help to make the ML model more specialized for solving PDEs. The hybrid neural network (NN) solver method is one of the strategies to develop more efficient PDE solver by interleaving NN with numerical solvers like finite element method (FEM). The hybrid NN model on the premise of numerical solver is easier to train and more stable than the purely data-driven model. For example, one previous study has used the hybrid NN model as a corrector for an incomplete numerical solver for the advection-diffusion problem. In every time step of simulation, NN corrects the error of incomplete solution obtained by a relaxed numerical solver with coarse meshing. The simulation in the next time step starts from the corrected solution, so NN interacts with the numerical solver iteratively. If the corrector is successfully trained, the incomplete but fast solver with corrector can provide reliable results comparable to the original massive solver. This study adopts the hybrid concept to develop a surrogate model for the near-field region, which is the heavy computation part in the simulation of geological disposal system. Various incomplete models such as coarse meshing or emptying the borehole domain are studied to construct a hybrid NN solver. This study also covers how to embed the hybrid NN in COMSOL Multiphysics to train and use it during the simulation.

Domain decomposition method (DDM) has been widely employed for the numerical analysis of large-scale problems due to its applicability to parallel computing. DDM divides the modeling domain into a set of subdomains and obtains the entire solution iteratively until the values of each subdomain which are shared with other subdomains, such as boundary values, are converged. Therefore, in general, DDM is a memory-efficient iterative algorithm with inherent parallelism on the geometric level. APro, the process-based total system performance assessment model, aims for simulating the radionuclide transport considering coupled multi-physics phenomena occurring in large-scale geological disposal system, which are inevitably accompanied by huge memory burden. Therefore, DDM is applicable for the large-scale problem of APro and its performance in parallel computing needs to be examined. The DDM solvers provided by COMSOL which constitute APro can be classified into two methods. One is the overlapping Schwarz method that each subdomain overlaps its neighboring domains and the other is the Schur complement method that subdomains are non-overlapping and separated by boundary domains. For the Schwarz method, the additive, hybrid, multiplicative and symmetric methods can be selected according to the solution update scheme. And for the Schur method, the additive and multiplicative ordering options can be chosen for solving Schur complement system. In this study, the calculation efficiency of the DDM solvers in COMSOL and the applicability to the cluster environment were examined. In aspect of efficiency, the memory requirements with different number of subdomains and calculation schemes were compared in a single node. Then, the memory requirements with increasing number of disposal tunnels and deposition holes were investigated in multiple nodes. As a result, on the cluster environment, with the help of distributed memory architecture which enables efficient memory usage, the applicability of DDM solvers to the large-scale problem of APro was confirmed.

APro, developed by KAERI as a process-based total system performance assessment model, can simulate the radionuclide transport affected by thermal, hydraulic, mechanical and geochemical changes that may occurs in the engineering and natural barriers of a geological disposal system. APro targets a large-scale and heterogeneous 3D system that includes more than 10,000 boreholes located about 500 m underground and hundreds of fractures of different sizes distributed within an area of several km2. Simulating transport and reaction phenomena for such a system through the global implicit approach (GIA) may require considerable computational resources or be intractable in some cases. Therefore, APro adopts the sequential non-iterative approach (SNIA), one of the operator splitting (OS) methods, to separate the mass transport and reaction phenomena into independent problems. By using SNIA, the parallel computation performance in APro with multiple cores is expected to be improved. In this study, the effect of SNIA on the parallel computation performance was analyzed through a simple 1D reactive transport problem. Without SNIA, finite difference equations, discretized from the partial differential equations (PDEs) describing the reactive transport problem, have to be solved at once because all dependent variables are nonlinearly and spatially interconnected through reaction and mass transport terms. When the reaction and mass transport terms are separated through SNIA, the mass transport problem can be converted into independent linear equations for each chemical and the efficient linear system solver can be applied to each linear equation. In particular, since the reaction problem is changed to independent nonlinear equations for each node, the parallel computation performance can be greatly improved. To verify this, the 1D reactive transport problem was implemented in MATLAB, and SNIA and GIA were applied to solve the problem. As a result, there was no significant difference in results between SNIA and GIA for proper spatial and temporal discretization, which verified the accuracy of SNIA. In order to see the parallel computation performance, the calculation times for SNIA and GIA with increasing number of cores were measured and compared. As the number of cores increased, the SNIA calculation speed became faster than that of GIA, which verified that SNIA could improve parallel computation performance in APro. In the future, the effect of SNIA on the parallel computation performance will be verified for the numerical analysis of large-scale geological disposal systems.

The timescale for the post-closure safety assessment of a deep geological repository ranges from ten thousand to a million year. In such a long period of time, the biosphere inevitably undergoes changes. Therefore, the long-term evolution of a biosphere is recognized as an important issue in the post-closure safety assessment of a deep geological repository for spent fuels. In this study, we reviewed the approaches to address the long-term evolution of a biosphere. The major drivers of longterm evolution of a biosphere are the climate change and the resulting landscape development. They can affect the hydrogeological and hydrogeochemical characteristics of a biosphere, and then the radionuclide migration through the biosphere followed by the exposure doses for the critical groups. In addition, human activities and the social developments can affect the climate change resulting in the long-term evolution of a biosphere. To make a biosphere assessment, the long-term evolution scenarios for the biosphere should be formulated considering these climate change, landscape development, and human activities. In addition, features, events, and processes (FEPs) that affect the long-term evolution of a biosphere should be used. According to the Safety Case reports of Finland, the major long-term evolution scenario drivers of a biosphere are local sea-level change due to climate change and land use related to crop type, irrigation procedures, livestock, forest management, construction of a well, and demographics. The climate change causing the local sea-level change can be simulated using various earth system models such as CLIMBER-2, MPI/UW, and UVic and an icesheet model such as SICOPOLIS. The review results of this study and FEPs related to the climate change, the landscape development, and human activities will be used to formulate long-term evolution scenarios for the safety assessment of a deep geological repository for spent fuels.

The design of the high-level radioactive waste (HLW) repository is made for isolating the HLW from the groundwater system by using artificial and natural barriers. Granite is usually considered to be a great natural barrier for the HLW repository in various countries including Sweden, Canada, and Korea due to its low hydraulic permeability. However, many fractures that can act as conduits for groundwater and radionuclides exist in granite. Furthermore, the decay heat generated by the HLW can induce groundwater acceleration through the fracture. Since the direction, magnitude, and lasting time of the heat-induced groundwater flow can be differed depending on the fracture geometry, the effect of fracture geometry on the groundwater flow around the repository should be carefully analyzed. In this study, groundwater models were conducted with various fracture geometries to quantify the effect of various properties of fractures (or fracture networks) on the heat-induced groundwater flow. In all models, the pressure around the repository only lasted for a short period after it peaked at 0.1 years. In contrast, the temperature lasted for 10,000 years after the disposal inducing the convective groundwater flow. Single fracture models with different orientations were conducted to evaluate the variations in groundwater velocities around the repository depending on the fracture slope. According to the results, the groundwater velocity on the fracture was the fastest when the regional groundwater flow direction and the fracture direction coincided. In double fracture models, various inclined fractures were added to the horizontal fracture. Due to the intersecting, the groundwater flow velocity showed a discontinuous change at the intersecting point. Lastly, the discrete fracture network models were conducted with different fracture densities, length distributions, and orientations. According to the modeling results, the groundwater flow was significantly accelerated when the fracture network density increased, or the average fracture length increased. However, the effect of the fracture orientation was not significant compared to the other two network properties.

In recent years, the importance of the thermo-hydraulic-mechanical-chemical coupled processes is increasing in the performance assessment (PA) of the high-level radioactive waste repository. In the case of mechanical behavior, it is very important because it can affect fluid flow and radionuclide transport by changing the porosity and permeability of the medium. In particular, Excavation Damaged Zone (EDZ) should be considered essential in PA because the migration of radionuclide is affected by the enhanced hydraulic transmissivity and altered geomechanical behavior of EDZ. Furthermore, due to various thermo-hydraulic behaviors such as decay heat generated from radioactive waste, pore water pressure increase, and swelling pressure of bentonite buffer material, mechanical evolution is occurred which may change the size and physical properties of EDZ. Therefore, to solve this problem, analysis of coupled thermal-hydraulic-mechanical (THM) processes with the effect of long-term evolution of EDZ due to the mechanical behavior should be accompanied. In this study, numerical model for the long-term evolution due to mechanical behavior considering EDZ using the Adaptive Process-based total system performance analysis framework for a geological disposal system (APro) proposed by the Korea Atomic Energy Research Institute (KAERI). In the case of EDZ, the concept of Mazars’ damage evolution model was applied to simulate the behavior using the continuum model, and the change in hydraulic properties according to the degree of damage was considered. To investigate the importance of mechanical behavior in PA, the results were compared by performing numerical analysis according to the presence or absence of mechanical analysis. Finally, numerical analysis considering the mechanical evolution of EDZ was conducted using the model developed in this study to investigate the effect of EDZ.

With the increase of temporarily-stored spent radioactive fuels, there is an increasing necessity for the safe disposal of high-level radioactive waste (HLW). Among various methods for the disposal of HLW, a deep geological disposal system is adapted as a HLW disposal strategy in many countries. Before the construction of a repository in deep geological condition, a performance assessment, which means the use of numerical models to simulate the long-term behavior of a multi-barrier system in HLW repository, has been widely performed to ensure the isolation of radionuclides from human and related environments for more than a million years. Meanwhile, Korea Atomic Energy Research Institute (KAERI) is developing a process-based total system performance assessment framework for a geological disposal system (APro). To improve the reliability of APro, KAERI is participating in DECOVALEX-2023 Task F, which is the international joint program for the comparison of the models and methods used in deep geological performance assessment. As a final goal of Task F, the reference case for a generic repository in fractured crystalline rock is described. The three-dimensional generic repository is located in a domain of 5 km in length, 2 km in width, and 1 km in depth, and contains an engineering barrier system with 2,500 deposition holes in fractured crystalline rock. In this study, a numerical simulation of the reference case is performed with COMSOL Multiphysics as a part of Task F. The fractured crystalline rock is described with the discrete fracture matrix (DFM) model, which expresses major deterministic fractures explicitly in the domain and minor stochastic fractures implicitly with upscaled quantities. As an output of the numerical simulation, fluid flow at steady-state and radionuclide transport are evaluated for ~106 years. The result shows that fractures dominate the transport of radionuclides due to much higher hydraulic properties than rock matrix. The numerical modeling approaches used in this study are expected to provide a basis for performance assessment of nuclear waste disposal repository located in fractured crystalline rock.

Multiple sorptive sites on natural illitic clays (e.g., frayed edge [FES], type II [TS], and planar sites [PS]) play an important role to diverse 137Cs immobilization in soil and aquifer environments. This study investigated the Cs sorption capabilities of 10 natural illitic clays at ranged Cs concentrations (1 ×10−7 to 1×10−3 mol·L−1) under various competing potassium concentration (distilled water to 1×10−1 mol·L−1). Additionally, multisite cation exchange model was performed to evaluate the best-fit sorption model and optimize the sorption capacities and affinities of multiple sorptive sites for Cs. Here, the experimental Cs sorption isotherms varied among 10 illtic clays, indicating different sorption capacities of Cs on illitic clays. The best-fit sorption model exhibited that variable Cs sorption of 10 illitic clays was significantly related to the sorption capacities at the FES (1.76 × 10−5 to 1.12×10−4 eq·kg−1), TS (1.59×10−3 to 9.76×10−3 eq·kg−1), and PS (2.14×10−2 to 1.51×10−1 eq·kg−1), respectively. The FES predominantly contributed to Cs sorption at low aqueous concentrations, whereas the TS and PS sorbed Cs at high concentrations. These sorption capabilities of multiple sorptive sites were correlated to illite contents and crystallinity of illitic clays, implicating that such parameters could be key factors to predict the Cs sorption for natural illitic clays in soil and aquifer environments. Finally, 1-D transport simulations represented that the severe Cs retardation occurred at low Cs concentration, implying that the FES predominantly affected to Cs transport in actual radioactive contamination sites (i.e., where low Cs concentration prevails), compared to the TS and/or PS.