As part of the third ATM Challenge, we performed a series of atmospheric dispersion simulations for routine releases of Xe-133 from ordinarily operating nuclear facilities such as Medical Isotope Production Facilities (MIPFs), Nuclear Power Plants (NPPs), and Research Reactors (RRs) in the Northern Hemisphere using our ATM, Lagrangian Atmospheric Dose Assessment System (LADAS), with Numerical Weather Prediction (NWP) data produced by the Korea Meteorological Administration (KMA). The simulation time period is 6 months, from June to November in 2014, and we used the stack emission data except for CNL (Canada) and IRE (Belgium) in accordance with the scenario of the third ATM Challenge 2019. In addition, the simulations were done individually for all MIPFs, NPPs, and RRs. We utilized 3-hourly KMA’s Unified Model Global Data Assimilation and Prediction System (UM-GDAPS) data with 0.35°×0.23° horizontal resolution as input meteorological fields and extracted hourly time series results for Xe-133 activity concentrations with few different resolutions such as 0.5°×0.5°, 0.35°×0.23°, and 0.1°×0.1° at several IMS stations in the Northern Hemisphere which were in normal operation in 2014. Considering previously reported values of daily Xe-133 release amounts for CNL and IRE, measured signals at some IMS stations (such as CAX17, DEX33, SEX63, and USX75) were well reproduced from the simulation results.

Considering the Fukushima nuclear accident and the marine discharge plan of contaminated (or treated) water, it is necessary a seafood monitoring system for radioactive material screening. Currently, radioactivity tests in seafood are conducting in Korea. Although current method using a HPGe detector can provide very low uncertainty in determining radioactivity, there is a limitation in that rapid inspection cannot be performed because of a time-consuming pretreatment process as well as long measurement time (typically 10,000 s). To overcome this limitation, we are developing an insitu inspection device, a kind of screening system, which can monitor the radioactivity in seafood in a near real-time basis. In this study, the actual seafood with a check source was measured to verify the reliability of the Monte Carlo simulation model. The detector used in the experiment was a 5-cm-thick polyvinyl toluene (PVT) plastic scintillator with a 0.5-cm-thick lead shield for background reduction. A Cs-137 check source was placed within seafood. The seafood used in the experiment was fishcake, raw oyster, and dried laver, which is the representative seafood of fish, shellfish, and seaweed. These three seafood products of the same size and shape as the manufacturing process were used to predict the performance realistically. We compared the energy spectrum of the Cs-137 check source obtained from measurements and simulations. The region of interest (ROI) around the Compton edge was set to reduce the influence of the background and electronic noise. The results showed that the measured and simulated spectrum were in good agreement.

A simulation model was developed for heavy water pre-enrichment and detritiation by the Combined Electrolysis and Catalytic Exchange (CECE) process. In the CECE process, heavy water enrichment and detritiation are based on the principle that concentrated in to water phase through an isotopic exchange reaction between water vapor and hydrogen gas produced by a water electrolysis. An operational analysis for a liquid phase catalytic exchange column was carried out by the model equations, composed of a material balance and combined equilibrium relationships for a scrubbing and catalyst bed, respectively. As a result of simulation, the optimum flow ratio of water to the rising hydrogen gas in contact with the down-coming water was predicted as the key variables in the separation performance analysis at a given feed flow rate and isotopic composition. From a graphical approach based on this model, the operating conditions can be determined within the range where the operating line does not meet the combined equilibrium line before reaching the specified target concentration.

Molten Salt Reactor, which employs molten salt mixture as fuel, has many advantages in reactor size and operation compared to conventional nuclear reactor. In developing Molten Salt Reactor, Offgas system should be properly designed since the fission products in off-gas accelerates the corrosion in reactor structure materials and deteriorates the purity of liquid fuel. The design of off-gas system therefore requires the preliminary study of the behavior of evolved fission products in off-gas units and the development of off-gas model is crucial in developing such system. In this study, we corrected the off-gas illustrative model proposed by ORNL (Nuclear Engineering and Design, vol 385(15) 111529, 2021) by employing physically consistent concept of capture rate of fission product and holdup. For the application of the corrected off-gas model to Chloride-based 6 MW Molten Salt Reactor, major fission products were firstly determined from OpenMC based neutronics calculation and chain reaction related to the major fission products were defined. Based on these data, the holdup behavior of fission products in off-gas units (decay tank, caustic scrubber, Halide trap, H2O trap and charcoal bad) were investigated.

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.

For safety assessment of a high-level radioactive waste disposal system, it is important to predict and analyze the coupled thermo-hydro-mechanical (THM) behaviors of bentonite, which is a buffer candidate material in the engineered barrier system. The Barcelona Basic Model (BBM) is a constitutive model to describe the geomechanical behaviors of partially saturated soils. Complicated tests are required to directly measure BBM parameters of bentonite. In this study, we demonstrate that probable BBM parameters can be sought by calibrating the BBM parameters to match simulation results to observed ones for two kinds of simple tests (swelling pressure test and free swelling test) instead of the complicated direct tests. In the swelling pressure test and free swelling test that were conducted by Japan Atomic Energy Agency (JAEA), water was injected into constrained and unconstrained bentonite core samples, and then swelling pressure and displacements were measured, respectively. We find optimal BBM parameters using a quasi-Newton optimization method that reproduce the observed swelling pressures and displacements in hydro-mechanical simulations. The optimal BBM parameters that are sought in the inversion process can be used to predict the THM behaviors of bentonite barriers in a high-level radioactive waste disposal system.

Radiation dose rates for spent fuel storage casks and storage facilities of them are typically calculated using Monte Carlo calculation codes. In particular, Monte Carlo computer code has the advantage of being able to analyze radiation transport very similar to the actual situation and accurately simulate complex structures. However, to evaluate the radiation dose rate for models such as ISFSI (Independent Spent Fuel Storage Installation) with a lot of spent fuel storage casks using Monte Carlo computational techniques has a disadvantage that it takes considerable computational time. This is because the radiation dose rate from the cask located at the outermost part of the storage facility to hundreds of meters must be calculated. In addition, if a building is considered in addition to many storage casks, more analysis time is required. Therefore, it is necessary to improve the efficiency of the computational techniques in order to evaluate the radiation dose rate for the ISFSI using Monte Carlo computational codes. The radiation dose rate evaluation of storage facilities using evaluation techniques for improving calculation efficiency is performed in the following steps. (1) simplified change in detailed analysis model for single storage cask, (2) create source term for the outermost side and top surface of the storage cask, (3) full modeling for storage facilities using casks with surface sources, (4) evaluation of radiation dose rate by distance corresponding to the dose rate limit. Using this calculation method, the dose rate according to the distance was evaluated by assuming that the concrete storage cask (KORAD21C) and the horizontal storage module (NUHOMS-HSM) were stored in the storage facility. As a result of calculation, the distance to boundary of the radiation control area and restricted area of the storage facility is respectively 75 m / 530 m (KORAD21C case), and 20 m / 350 m (NUHOMS-HSM case).

Considering the domestic condition with small land area and high population density, it is necessary to develop technology that can reduce the disposal area than the deep geological disposal method. For this, KAERI is developing a nuclide management process that can reduce the environmental burden of spent fuel, and establishing an evaluation model that can evaluate the performance of various process options. It is expected that an optimal option of the nuclide management process can be derived from disposal perspective by applying the evaluation model. The mass flow between processing steps of the radionuclide management process is the basic quantity required to quantify the evaluation criteria. Therefore, we built a generalized block model on GoldSim, which can simulate mass flow of various radionuclide management process options. In addition to the mass flow, this model was established to derive the amount of wastes generated by each processing step, the composition of nuclides, and radiological properties (decay heat, radioactivity, etc.). The mass flow and waste property derived from the models are closely related to the factors that determine the area of disposal concepts. Based on this, a disposal area calculation model was established as a model to evaluate the effectiveness of the radionuclide management process on environmental burden reduction. For verification, three process options, which can manage radionuclides having high decay heat (Cs, Sr) or large volume (U), were selected and evaluated as reference processes. And two disposal options, deep geological disposal and deep borehole disposal concepts were considered to be linked with the processes. As a result, it was confirmed that the disposal area could be reduced in the process separating radionuclides having high decay heat. In the future, other evaluation models for economic viability and safety will be added in the GoldSim model.

For the decommissioning or continuous long-term power generation of nuclear power plants, it is necessary to transfer the spent nuclear fuel from the wet storage pool to the dry storage. Spent nuclear fuel should go through the drying process, which is the first step of dry storage. The most important part in the drying process is the removal of the residual water. The spent fuel might be stored in a dry storage system for a long time. The integrity of internal components and spent fuel cladding should be maintained during the storage period. If residual water is present, problems such as aging of metal materials, oxidation of cladding, and the hydride-reorientation could occur. The presence or absence of residual water after vacuum drying is evaluated by pressure. If there is residual water in the vacuum drying process, it evaporates easily at low pressure to form water vapor pressure and the internal pressure rises. In the recent EPRI High burn up demonstration test, the gas inside the canister that satisfied the dryness criteria was extracted and analyzed. It showed that the water content was higher than the expected value. We are conducting verification studies on the pressure evaluation method, which is an indirect evaluation method of vacuum drying. The vacuum drying test was performed on small specimens at Sandia National Laboratory, and quantitative residual water evaluation was also performed. The report did not mention a detailed method for the assessment of residual water. Based on the test results of SNL, direct residual water evaluation was performed using energy balance. If the dryness criteria were satisfied, the quantitative amount of residual water was also evaluated. As a result, almost the same result as the evaluation result of SNL was derived, and it was confirmed that most of the water was removed when the dryness criteria was satisfied.

For the spent fuel modeling, the plastic model of the cladding used in FRAPCON uses the σ􀷥 = K􀟝̃􀯡 􁉂 􀰌􁈶 􀬵􀬴􀰷􀰯􁉃 􀯠 format. Strength coefficient (K), strain hardening exponent (n), strain rate sensitivity constant (m) are derived as the function of temperature. The coefficient m related to the strain rate shows dependence on the strain rate only in the α-β phase transition section, 1,172.5~1,255 K. But this is the analysis range of the FRAPTRAN code, which is an accident condition nuclear fuel behavior evaluation code. It does not apply to evaluate spent fuel. This coefficient in FRAPCON is used as a constant value (0.015) below 750 K (476.85°C), and at a temperature above 750 K, it is assumed that it is linearly proportional to the temperature without considering the strain rate dependence, also. In order to confirm the effect of strain rate, multiple test data performed under various conditions are required. However, since the strain rate dependence is not critical and test specimen limitation in the case of spent fuel, it is needed to replace with a new plastic model that does not include the strain rate term. In the new plastic model, the basic form of the Ramberg-Osgood equation (RO equation) is the same as ε􀷤 = 􀰙􀷥 􀮾 + 􀜭􀯥 􁉀􀰙􀷥 􀮾􁉁 􀯡􀳝. If the new variable α is defined as α = 􀜭􀯥􁈺􀟪􀯢/􀜧􁈻􀯡􀳝􀬿􀬵, this equation can be transformed into ε􀷤 = 􀰙􀷥 􀮾 + 􀟙 􀰙􀷥 􀮾 􁉀 􀰙􀷥 􀰙􀰬 􁉁 􀯡􀳝􀬿􀬵 . The procedure for expressing the stress-strain curve of the cladding with the RO equation is as follows. First, convert the engineering stress-strain into true stress-strain. Second, using a data analysis program such as EXCEL or ORIGIN, obtain the slope of the linear trend-line on the linear part and use it as the elastic modulus. Third, using the 0.2% offset method, find the yield point and the yield stress. Finally, using the solver function of EXCEL, find the optimal values of α and 􀝊􀯥 that minimize the sum of errors. The applicability of the suggested RO equation was evaluated using the results of the Zircaloy-4 plate room temperature tensile test performed by the KAERI and the Zircaloy cladding uniaxial tensile test results presented in the PNNL report. Through this, the RO equation was able to express the tensile test results within the uncertainty range of ±0.005. In this paper, the RO equation is suggested as a new plastic model with limited test data due to the test specimen limitation of spent fuel and its applicability is confirmed.

Numerous spent nuclear fuels are generated every year in Korea. To solve the spent nuclear fuel problem within saturated temporary storage, the authorities are readying to build an interim storage and a permanent disposal facility in the country. At the same time, the authorities are readying to establish a management procedure for spent nuclear fuel. In the future, the authorities need to make and apply the Database of spent nuclear fuel to practice the management procedure. However, the structure of a traditional database is not reasonable for information management because it has a problem with listing data and identifying data features due to its structure. In addition, the traditional database always exists human error from working in Excel program by a human. Therefore, this research proposes a new standard information management model based on Semantic Web technique. Semantic Web uses a data structure named ontology. By using the ontology in the information database of the spent nuclear fuel, users, such as institutions related to management, could more easily recognize and understand the Database. Furthermore, since this task proceeds in the ontology construction program, the human error in the new model reduces rather than an environment of the traditional database.

According to the ‘Basic Plan for High-Level Radioactive Waste Management (draft)’, the total amount of CANDU spent nuclear fuel is expected to be approximately 660,000 bundles. To safely and efficiently transport this amount to interim storage facilities, it is essential to develop a large-capacity transport cask. Therefore, we have been developing a large-capacity PHWR spent nuclear fuel transport cask, called the KTC-360 transport cask. According to the transport-cask related regulations, the KTC-360 transport cask was classified as a Type B package, and such packages must be able to withstand a temperature of 800°C for a period of 30 min. It is desirable to conduct a test using a fullscale model of a shipping package when performing tests to evaluate its integrity. However, it is costly to perform a test using a full-scale model. Therefore, to evaluate the thermal integrity of the KTC-360 transport cask, the fire test was conducted using a slice model. For comparison purposes, the fire test was also carried out using a 1/4 scale model. In the fire test using a slice model and in the fire test using a 1/4 scale model, the maximum temperature of the cask body was lower than the permitted maximum temperature limit. Therefore, the thermal integrity of the KTC-360 transport cask could be considered to be maintained. The temperature results from the fire test using a slice model were higher than those of the fire test using a 1/4 scale model. Therefore, the effect of flame on a transport cask without combustible materials, such as the KTC-360 transport cask, seems to be affected by the reduction in the time rather than the size reduction.

As the zircaloy cladding absorbs an excessive amount of hydrogen and cooled down under hoop stress, radial hydride may be precipitated by hydride reorientation phenomenon. There have been many previous studies about the threshold stress of the reorientation, but it is known that the quantitative degree of hydride reorientation rather than the threshold is important for the prediction of mechanical properties. A thermodynamic model for Radial Hydride Fraction (RHF) prediction has been developed in this study. The model calculates RHF with respect to temperature, cooling rate, hydrogen content, and applied stresses. Once the cooling rate is given, the solid solution concentration at each temperature is determined by Hydrogen-Nucleation-Growth-Dissolution model. Subsequently, the increment of radial hydride is derived by nucleation and growth theory. The code based on the thermodynamic theory can provide the prediction of RHF under hoop stress, as well as a change in precipitation behavior over time. RHF of the zircaloy cladding in long-term dry storage can be obtained by the implementation of the code and the degradation of the cladding is directly estimated according to the correlation between RHF and mechanical properties. Ongoing experimental validation of the developed model is discussed.

As the saturation rate of temporary storage facilities for spent nuclear fuel increases, regulatory demands such as interim storage and permanent disposal of spent nuclear fuel are expected to begin in earnest. Considering the domestic situation where all nuclear power plants are located on the waterfront site, the interim storage site is also likely to be located on the waterfront site, and maritime transportation is one of the essential management stages. Currently, there are no independently developed maritime transportation risk assessment code in Korea, and no research has been conducted to evaluate the release of radioactive waste due to the sinking of transport container. Therefore, it is necessary to secure technology to properly reflect the domestic maritime transportation environment and to assess the impact of the sinking accident and to carry out safety regulations. To accurately calculate the releaser rate of radionuclides contained in a cask with breached containment boundary, the flow rate through the gap generated in the containment boundary should be calculated. The fluid flow through this gap which is probably in micro scale in most situations should be evaluated combining the fluid flow inside and outside the cask. In this study, a detailed computational fluid dynamics model to evaluate the internal fluid flow in the cask and a simplified model to capture the fluid flow and the heat transfer around the cask in the sea are constructed. The results for the large scale model are compared with the analytic formula for verification of heat transfer coefficient and they showed good agreements. The heat transfer coefficient thus found can be used in the detailed model to provide more realistic data than those obtained from assumed heat transfer coefficient around the surface of the cask. In the future, fluid flow through the gap between the lid and the body of the cask will be evaluated coupling the models developed in this work.