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.

Copper, mainly used as a material for outer canister, generates various corrosion products under aerobic and anaerobic conditions in the operational and/or post-closure phases of the deep geological repository. These products could affect performance of engineering barrier system (EBS) through interaction with surrounding bentonite that makes up the buffer and backfill materials. Accordingly, in this study, we suggested research items to be conducted to minimize degradation of EBS due to copper corrosion products, based on the phenomenological review results for copper corrosion mechanisms and interaction between resultant product and bentonite in the deep geological disposal environment. During the post-closure phase, condition in the disposal facility changes form aerobic to anaerobic over time, and thereby, causes and products of copper corrosion vary. Under aerobic condition, copper corrosion is mainly induced by oxygen (O2) in the repository, chloride (Cl-) and carbonate (CO3 2-) ions from groundwater flowing into the facility, resulting in corrosion products such as cuprite (Cu2O), tenorite (CuO), atacamite (CuCl2·3Cu(OH)2) and malachite (Cu2CO3(OH)2). And, copper corrosion under anaerobic condition is primarily due to hydrogen sulfide (H2S) and sulfate (SO4 2-) in groundwater flowing into the facility, leading to formation of chalcocite (Cu2S) and covellite (CuS) as corrosion products. Depending on environment of the disposal facility, copper corrosion products are dissolved and ionized to Cu2+ in groundwater, and subsequently adsorbed on the nearby smectite. Then, it causes a cation exchange reaction with exchangeable cations in the interlayer of smectite. As a result of reviewing the previous experiments, it was confirmed that Cu2+-exchanged bentonite has a slightly reduced basal spacing and swelling capacity. From the results as above, there is a possibility that performance of EBS may be degraded due to copper corrosion products. To minimize its effect of degradation in the domestic facility, items to be further studied are as follows: (a) Method for reducing copper corrosion such as selection of appropriate material and structure for the canister, and (b) How to control dissolution of copper canister product into groundwater through predicting type and ionization process. The results of this study could be directly used to developing design concept of EBS for the domestic disposal facility and to establishing roadmap of future R&D programs.

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.

The buffer is installed around the disposal canister, subjected to heating due to decay heat while simultaneously experiencing expansion influenced by groundwater inflow from the surrounding rock. The engineering barrier system for deep geological disposal require the evaluation of longterm evolution based on the verification of individual component performance and the interactions among components within the disposal environment. Thus, it is crucial to identify the thermalhydro- mechanical-chemical (THMC) processes of the buffer and assess its long- and short-term stability based on these interactions. Therefore, we conducted experimental evaluations of saturationswelling, dry heating, gas transport, and mineralogical alterations that the buffer may undergo in the heated-hydration environment. We simulated a 310 mm-thick buffer material in a cylindrical form, simulating the domestic disposal system concept of KRS+ (the improved KAERI reference disposal system for spent nuclear fuel), and subjected it to the disposal environment using heating cartridges and a hydration system. To monitor the thermal-hydro-mechanical behavior within the buffer material, load cells were installed in the hydration section, and both of thermal couples and relative humidity sensors were placed at regular intervals from the heat source. After 140 days of heating and hydration, we dismantled the experimental cell and conducted post-mortem analyses of the samples. In this post-mortem analysis, we performed functions of distance from the water contents, heat source, wet density, dry density, saturation, and X-Ray diffraction analysis (XRD). The results showed that after 140 days in the heated-hydration environment, the samples exhibited a significant decrease water contents and saturation near the heat source, along with very low wet and dry densities. XRD Quantitative Analysis did not indicate mineralogical changes. The findings from this study are expected to be useful for input parameters and THMC interaction assessments for the long-term stability evaluation of buffer in deep geological disposal.

Rock discontinuities in underground rock behave as weak planes and affect the safety of underground structures, such as high-level radioactive waste disposal and underground research facilities. In particular, rock discontinuities can be a main flow path of groundwater and induce large deformation caused by stress disturbance or earthquakes. Therefore, it is essential to investigate the characteristics of rock discontinuities considering in-situ conditions when constructing highlevel radioactive waste disposal, which needs to assure the long-term safety of the structure. We prepared Hwang-Deung granite rock block specimens, including a saw-cut rock surface, to perform multi-stage direct shear tests as a preliminary study. In the multi-stage direct shear tests, we can exclude possible errors induced by different specimens for obtaining a full failure envelope by using an identical specimen. We applied the initial normal stress of 3 MPa on the specimen and increased the normal stress to 5 and 10 MPa step by step after peak shear stress observation. We obtained the mechanical properties of saw-cut rock surfaces from the experiments, including friction coefficient and cohesion. Additionally, we investigated the effect of filling material between rock discontinuities, assuming the erosion and piping phenomenon in the buffer material of the engineering barrier system. When the filling material existed in the rock surfaces, the shear characteristics deteriorated, and the effect of bentonite was dominant on the shear behavior.

The engineered barrier system (EBS), composed of spent nuclear fuel, canister, buffer and backfill material, and near-field rock, plays a crucial role in the deep geological repository for high-level radioactive waste. Understanding the interactions between components in a thermo-hydro-mechanical -chemical (THMC) environment is necessary for ensuring the long-term performance of a disposal facility. Alongside the research project at KAERI, a comprehensive experimental facility has been established to elucidate the comprehensive performance of EBS components. The EBS performance demonstration laboratory, which installed in a 1,000 m2, consists of nine experimental modules pertaining to rock mechanics, gas migration, THMC characteristics, buffer-rock interaction, buffer & backfill development, canister corrosion, canister welding, canister performance, and structure monitoring & diagnostics. This facility is still conducting research on the engineering properties and complex interactions of EBS components under coupled THMC condition. It is expected to serve as an important laboratory for the development of the key technologies for assessing the long-term stability of engineered barriers

In order to ensure the long-term safety of a deep geological repository, the performance assessment of the Engineered Barrier System (EBS) considering a thermal process should be performed. The maximum temperature at the side wall of a disposal canister for the technical design requirement should not exceed 100°C. In this study, the thermal modelling was conducted to analyze the effects of the thermal process from a disposal canister to the surrounding near-field host rock using the PFLOTRAN code. The mesh was generated using the LaGriT code and the material properties were assigned by applying the FracMan code. Initial conditions were set as the average geothermal gradient (25.7°C/km) and an average surface temperature (14.7°C) in Korea. The highest temperature was observed at the middle of the canister side wall. The temperature of the buffer was lower than that of the canister, and the temperature increase of the deposition tunnel and the host rock was insignificant due to the lower effect of the heat source. The result of the thermal evolution of the EBS represented the highest thermal effects in the vicinity of the canister. In addition, the thermal effects were largely decreased after 10 years of the entire simulation period. It demonstrated that the model took 3 years to heat up the buffer around the canister. The temperature at the canister side wall increased until 3 years and then decreased after that time. This is because that the radioactive decay heat from the heat source was emitted enough to raise the overall temperature of the EBS by 3 years. However, the decay heat rate of the canister decreased exponentially with the disposal time and then its decay heat was not emitted enough after 3 years. In conclusion, the peak temperature results of the EBS were lower than 70°C to meet the technical design requirement.

Due to the necessity of isolating spent nuclear fuel (SNF) from the human life zone for a minimum of 106 years, deep geological disposal (DGD) has emerged as a prominent solution for SNF management in numerous countries. Consequently, the resilience of disposal canisters to corrosion over such an extended storage period becomes paramount. While copper exhibits a relatively low corrosion rate, typically measured in millimeters per million years, in geological environment, special attention must be directed towards verifying the corrosion resistance of copper canister welds. This validation becomes inevitable during the sealing of the disposal canister once SNFs are loaded, primarily because the weld zone presents a discontinuous microstructure, which can accelerate both uniform and localized corrosion processes. In this research, we conducted an in-depth analysis of the microstructural characteristics of copper welds manufactured by TIG-based wire are additive manufacturing, which is ideal for welding relatively large structures such as a disposal canister. To simulate the welds of copper canister, a 12 mm thick oxygen-free plate was prepared and Y and V grooves were applied to perform overlay welding. Both copper welding zones were very uniform, with negligible defects (i.e., void and cracks), and contained relatively large grains with columnar structure regardless of groove types. For improving microstructures at welds with better corrosion resistance, the effect of preheat temperature also investigated up to 600°C.

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.

The Barcelona Basic Model (BBM) is an elasto-plastic model used to describe the coupled thermo-hydro-mechanical behaviors of unsaturated soil. BBM is frequently adopted to model the unique swelling behavior of bentonite, which is generally considered as the buffer material between the host rock and the canister containing high-level radioactive waste in deep geological repositories, under the changing thermal, hydraulic, mechanical and chemical conditions during the lifetime of repository. Therefore, a variety of the continuum-based numerical methods tried to add the BBM for modelling the multi barrier systems of geological repository and succeeded to describe the elasto-plastic deformation of bentonite. However, to demonstrate the entire barrier systems the host rock should be modelled simultaneously with the buffer materials, and the continuum-based methods may be limited in their ability to reflect the fracture networks in the host rock which could be the major flow channels of groundwater. This research applies BBM in 3DEC, a three-dimensional block-based discrete element method, and validates the model by comparing the change of specific volume and mean effective stress during three numerical test cases. Discontinuum-based numerical methods with BBM can be extended to describe the coupled thermo-hydro-mechanical processes of multi-barrier systems in geological repositories, with a focus on the interaction between the host rock and bentonite.

The natural barrier system surrounding the geological repository for high-level radioactive waste plays a crucial role in preventing or delaying the leakage of radionuclides. Therefore, the natural barrier should ensure low permeability to prevent groundwater flow into the engineered barrier system throughout the repository’s lifetime. Crystalline rock, considered as the host rock for the geological repository in Korea, exhibits low intact rock permeability, but the crystalline rock often contains the multiple discontinuities due to its high brittleness that can allow the unexpected fluid flow. Therefore, the long-term hydraulic behavior of the discontinuity should be characterized while considering additional thermal, mechanical, and chemical effects. In comparison to thermal, hydraulic, and mechanical processes, the chemical processes on the discontinuities progress relatively slowly, resulting in limited researches to include these chemical processes. This research introduces mechanisms the involving coupled thermal-hydraulic-mechanicalchemical processes focusing on the rough fracture surfaces and asperities. The chemically-induced changes in mechanical and hydraulic properties are described based on pressure solution and precipitation concepts. A comprehensive review of laboratory tests, field tests, and numerical simulations is conducted related to the chemically-induced coupled processes in fractured rock. Laboratory tests, in particular, concentrate on microscopic changes in fracture asperities induced by pressure solution to analyze chemically-induced aperture changes. The TOUGHREACT, an integral finite difference method program for thermal-hydraulic-chemical simulations, is generally employed to model the chemical response of pressure solution and precipitation on fracture surfaces. The TOUGHREACT includes a module to describe effective porosity and permeability changes based on the modified cubic law, so the real-time change of the fracture permeability can be reflected during the flow simulation. Considering the coupled thermal-hydraulic-mechanicalchemical processes of discontinuity, it becomes evident that the chemical processes under repository conditions (long-term, high temperature, and high pressure) can disturb the hydraulic performance of the natural barrier, so further research is required to characterize the chemically-induced coupled processes for assessing the long-term performance of the natural barrier system.

The compacted bentonite buffer is a key component of the engineered barrier system in deep geological repositories for high-level radioactive waste disposal. Groundwater infiltration into the deep geological repository leads to the saturation of the bentonite buffer. Bentonite saturation results in bentonite swelling, gelation and intrusion into the nearby rock discontinuities within the excavation damaged zone of the adjacent rock mass. Groundwater flow can result in the erosion and transport of bentonite colloids, resulting in bentonite mass loss which can negatively impact the long-term integrity and safety of the overall engineered barrier system. The hydro -mechanicalchemical interactions between the buffer, surrounding host rock and groundwater influence the erosion characteristics of the bentonite buffer. Hence, assessing the critical hydro-mechanicalchemical factors that negatively affect bentonite erosion is crucial for the safety design of the deep geological repository. In this study, the effects of initial bentonite density, aperture, discontinuity angle and groundwater chemistry on the erosion characteristics of Bentonil WRK are investigated via bentonite extrusion and artificial fracture experiments. Both experiments examine bentonite swelling and intrusion into simulated rock discontinuities; cylindrical holes for bentonite extrusion experiments and plane surfaces for artificial fracture experiments. Compacted bentonite blocks and bentonite pellets are manufactured using a compaction press and granulation compactor respectively and installed in the transparent extrusion cells and artificial fracture cells. The reference test condition is set to be 1.6 g/cm3 dry density and saturation using distilled water. After distilled water or solution injection, the axial and radial expansion of the bentonite specimens into the simulated rock discontinuities are monitored for one month under free swelling conditions with no groundwater flow. Subsequent flow tests are conducted using the artificial fracture cell to determine the critical flow rate for bentonite erosion. The intrusion and erosion characteristics are modelled using a modified hydro-mechanicalchemical coupled dynamic bentonite diffusion model and a fluid-based hydro-mechanical penetration model.

Buffer materials play an important role in preventing the leakage of radionuclides from the residue. The mineralogical properties of these buffer materials are critical in repository design. This study presents the fundamental properties of Na-type MX80 and a novel Ca-type Bentonil- WRK. The CaO to MgO ratio in Bentonil-WRK was approximately 1:1, and the CaO to Na2O ratio was approximately 2.8:1. These results suggest that Bentonil-WRK demonstrates a lower swelling index compared to Gyeongju bentonite due to its CaO-to-MgO ratio’s proximity to 1:1, despite having a higher montmorillonite content than Gyeongju bentonite. The results of this research can provide useful foundational data for the evaluation of the thermal-hydraulic-mechanical-chemical behavior of buffer materials.

Buffer materials are one of the engineering barrier components in high-level radioactive waste disposal facilities. Compacted bentonite has been known as the most suitable buffer material so far, and research is being conducted worldwide to determine the characteristics of suitable bentonite blocks in each country. Therefore, this study aims to compare and analyze various properties of different buffer material candidates, including thermal, hydraulic, and mechanical properties. Buffer material candidates for Korea disposal system, Kyungju Bentonite (KJ-I, KJ-II), and Bentonil- WRK were compared. The properties were compared and analyzed based on experimental and literature data. The data obtained from this report can be used to understand the behavior of buffer materials and assess whether they meet various criteria, such as temperature and saturation, and ultimately serve as an essential input variable database for safety evaluations of disposal systems.

The presence of technological voids in deep geological repositories for high-level radioactive nuclear waste can have negative effects on the hydro-mechanical properties of the engineered barrier system when groundwater infiltrates from the surrounding rock. This study conducted hydration tests along with image acquisition and X-ray CT analysis on compacted Korean bentonite samples, which simulated technological voids filling to investigate the behavior of fracturing (piping erosion) and cracking deterioration. We utilized a dual syringe pump to inject water into a cell consisting of a bentonite block and technological voids at a consistent flow rate. The results showed that water inflow to fill technological voids led to partial hydration and self-sealing, followed by the formation of an erosional piping channel along the wetting front. After the piping channel generated, the cyclic filling-piping stage is characterized by the repetitive accumulation and drop of water pressure, accompanied by the opening and closing of piping channels. The stoppage of water inflow leads to the formation of macro- and micro cracks in bentonite due to moisture migration caused by high suction pressure. These cracks create preferential flow paths that promote longterm groundwater infiltration. The experimental test and analysis are currently ongoing. Further experiments will be conducted to investigate the effects of different dry density in bentonite, flow rate, and chemical composition of injected water.

In this study, a third metal layer with a higher corrosion potential than copper was introduced between the copper and cast iron layer to strengthen the corrosion resistance of the copper layer which is considered as a corrosion resistant barrier in the disposal container for spent nuclear fuel. Three types of corrosion-resistant metals, silver, nickel, and titanium, were selected as the intermediate insertion layer, and the galvanic specimens of two bonded metals were exposed to KURT (KAERI Underground Research Tunnel) groundwater and a high voltage of 1.0 V was applied to corrode the specimens at electrochemically accelerated condition. Corrosion of copper part was confirmed in Cu-Ti, Cu-Ni, and Cu-Ag galvanic specimens, but copper part was not corroded in Cu-Fe galvanic specimen. If the corrosion-resistant intermediate layer proposed in this study works properly, the local corrosion problem of copper disposal canister is expected to be some degree solved, which can apply to a welding part or a stress concentrated part.

Korea Atomic Energy Research Institute (KAERI) has been operating the Post Irradiation Examination Facility (PIEF). The facility has many PIE equipment and one of them is a hydrogen analyzer for measuring hydrogen contents in Zr cladding of spent fuel. The cladding tube of fuel is oxidized in the core environment of high temperature and pressure and absorbs some of the hydrogen generated during the oxidation. The hydrogen content increases with the increase of burn-up, and causes hydriding of the material, which degrades the mechanical properties. Therefore, hydrogen content analysis of the cladding tube is required for the performance and integrity evaluation of spent fuel. In PIEF, the hydrogen analyzer extracts hydrogen gas from Zr cladding by the hot extraction method. The hydrogen gas flows with inert gas and oxidizes to H2O through a CuO reagent. Finally, an IR detector measures the hydrogen amount from the absorbed IR intensity at a specific wavelength. Because the equipment is in the glove box and has some consumable parts, the maintenance work was performed as a radiation work.

Nuclear fuel assemblies are exposed to high temperature and high pressure environments underwater for long periods of time in a reactor, leading to deterioration of the assembly structure. These assembly consists of fuel rods, grids, a top nozzle, a bottom nozzle and guide tubes. In particular, the integrity of the guide tube made of Zircaloy-4 is a very important part in handling the assembly. In the Post Irradiation Examination Facility (PIEF), there are 14×14 Westinghouse STD assemblies that have lost their handleability due to the top nozzle being removed for damaged fuel rod test. To handle these assemblies, it is reasonable to use cut guide tubes whenever possible. Therefore, it is necessary to determine the irradiation embrittlement state of the guide tube before designing or manufacturing parts that can connect the top nozzle and the guide tubes. Therefore, in this paper, the location for installing the top nozzle-guide tube connection parts was selected in the height range of 3,460 to 3,713 mm, and guide tube specimens were made within that range. Offset strain was derived from the load-displacement curve obtained through compression testing to confirm whether the ductility of guide tubes was maintained. As a result, there was no significant difference in strength and ductility of the guide tube within the above length range. In addition, it was confirmed that the ductility was maintained enough to install the top nozzle-guide tube connection parts. Therefore, it is judged that there will be no problem even if the top nozzle-guide tube connection parts are installed in the guide tube to restore the handleability of the assemblies.

For efficient design and manufacture of PWR spent fuel burnup detector, data simulated with various condition of spent fuel in the NPP storage pool is required. In this paper, to derive performance requirements of spent fuel burnup detector for neutron flux and dose rates were evaluated at various distances from CE16 and WH17 types of fuel, representatively. The evaluation was performed by the following steps. First, the specifications of the spent fuel, such as enrichment, burnup, cooling time, and fuel type, were analyzed to find the conditions that emit maximum radioactivity. Second, gamma and neutron source terms of spent fuel were analyzed. The gamma source terms by actinides and fission products and neutron source terms by spontaneous and (α, n) reactions were calculated by SCALE6 ORIGAMI module. Third, simulation input data and model were applied to the evaluation. The material composition and dose conversion factor were referred as PNNL-15870 and ICRP-74 data, respectively and dose rates were displayed with the MCNP output data. It was assumed that there was only one fuel modeled by MCNP 6.2 code in pool. The evaluation positions for each distance were selected as 5 cm, 10 cm, 25 cm, 50 cm, and 1 m apart from the side of fuel, respectively. Fourth, neutron flux and dose rates were evaluated at distance from each fuel type by MCNP 6.2 code. For WH 17 types with a 50 GWd/MTU burnup from 5 cm distance close to fuel, the maximum neutron flux, gamma dose rates and neutron dose rates are evaluated as 1.01×105 neutrons/sec, 1.41×105 mSv/hr and 1.61×101 mSv/hr, respectively. The flux and dose rate of WH type were evaluated to be larger than those of CE type by difference in number of fuel rods. The relative error for result was less than 3~7% on average secured the reliability. It is expected that the simulated data in this paper could contribute to accumulate the basic data required to derive performance requirements of spent fuel burnup detector.