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

The seven-year research project entitled “Development of workflow for integrated 3D geological site descriptive modeling” is being carried out from 2023. This research is funded by Ministry of Trade, Industry, and Energy (MOTIE). Progress of the research is discussed here. The integrated 3D geological SDM (site descriptive model; GSDM hereafter) consists of three part; 1) three dimensional representation of geologic elements, 2) database for material properties and modeling results from SDMs of other disciplines (e.g., rock mechanics), and 3) a visualization tool for geology, material properties and modeling results. The GSDM is comparable to the GDSMs of SKB and POSIVA in its representation of geology by volume of geologic elements. However, our GSDM is different in that extra information of material properties and an extra tool for visualization is included in the GDSM. The rationale for incorporating material properties and a visualization tool into the GSDM is to expedite the development of the GSDM and SDMs of other disciplines by allowing single institution to integrate database and visualization with the GSDM. SKUA-GOCAD is used for representation of geologic surfaces for ductile and brittle shear zones, and also for surfaces for delineation of volumes of rock units. We have adopted SKUAGOCAD because the program offers powerful functions of interpolation including borehole data and geophysical prospecting. So far, we have tested the program for five different geologies, including sedimentary, high-grade metamorphic, and intrusive igneous geology. The test results are promising. Incorporation of data and modeling results for the SDMs of other disciplines is at conceptual stage. The working conceptual model involves the following steps, 1) to provide the modeler of other disciplines with surface information representing geologic elements, 2) the modeler returns not only material properties but the results of numerical analysis, and 3) incorporation of material properties and modeling results into database. Since the numerical codes in other disciplines adopt different types of formats for 3D geology, we plan to adopt the widely used FEM format prepared by Gmsh. The visualization tool will also adopt Gmsh for graphical representation of 3D geology as well as database for material properties and modeling results. When the working model of GSDM becomes available, rapid and significant progress is expected in the SDMs of other disciplines and related areas, for example, geotechnical investigation for deep geological repository.

Since the first operation of the Gori No. 1 nuclear power plant in Korea was started to operate in 1978, currently 24 nuclear power plants have been being operated, out of which 21 plants are PWR types and the rest are CANDU types. About 30% of total electricity consumed in Korea is from all these nuclear power plants. The accumulated spent nuclear fuels (SFs) generated from each site are temporarily being stored as wet or dry storage type at each plant site. These SFs with their high radiotoxicity, heat generating, and long-lived radioactivity are actually the only type of high-level radioactive waste (HLW) in Korea, which urgently requires to be disposed of in deep geological repository. Studies on disposal of HLW in various kind of geological repositories have been carried out in such countries as Sweden, Finland, United States, and etc. with their own methodologies and management policies in consideration of their situations. In Korea long-term R&D research program for safe management of SF has also been conducted during last couple of decades since around 1997, during which several various alternative type of disposal concepts for disposal of SNFs in deep geological formations have been investigated and developed. The first concept developed was KAERI Reference Disposal System (KRS) which is actually very much similar to Swedish KBS-3, a famous concept of direct disposal of SF in stable crystalline rock at a depth of around 500 m which has been regarded as one of the most plausible method worldwide. The world first Finnish repository which is expected to begin to operate sooner or later will be also this type. Since the characteristics of SF discharged from domestic nuclear reactors have been changed and improved, and burnup has sometimes increased, a more advanced deep geological repository system has been needed, KRS-HB (KRS with High Burnup SF) has been developed and in consideration of the dimensions of SNFs and the cooling period at the time point of the disposal time, KRS+, a rather improved disposal concept has also been subsequently developed which is especially focused on the efficient disposal area. Recently research has concentrated on rather advanced disposal technology focused on a safer and more economical repository system in recent view of the rapidly growing amount of accumulated SF. Especially in Korea the rock mass and the footprint area for the repository extremely limited for disposal site. Some preliminary studies to achieve rather higher efficiency repository concept for disposal of SF recently have already been emphasized. Among many possible ones for consideration of design for high-efficiency repository system, a double-layered system has been focused which is expected to maximize disposal capacity within the minimum footprint disposal area. Based on such disposal strategy a rather newly designed performance assessment methodology might be required to show long-term safety of the repository. Through the study some prerequisites for such methodological development has been being roughly checked and investigated, which covers FEP identification and pathway and scenario analyses as well as preliminary conceptual modeling for the nuclide release and transport in nearfield, far-field, and even biosphere in and around the conceptual repository system. Through the study such scenarios and models has been implemented to development of a safety assessment by utilizing GoldSim development tool for a rough quantitative comparison with existing disposal options and simple illustration purpose as well as for showing how to develop and implementation of the model to GoldSim templet.

In the design of HLW repositories, it is important to confirm the performance and safety of buffer materials at high temperatures. Most existing models for predicting hydraulic conductivity of bentonite buffer materials have been derived using the results of tests conducted below 100°C. However, they cannot be applied to temperatures above 100°C. This study suggests a prediction model for the hydraulic conductivity of bentonite buffer materials, valid at temperatures between 100°C and 125°C, based on different test results and values reported in literature. Among several factors, dry density and temperature were the most relevant to hydraulic conductivity and were used as important independent variables for the prediction model. The effect of temperature, which positively correlates with hydraulic conductivity, was greater than that of dry density, which negatively correlates with hydraulic conductivity. Finally, to enhance the prediction accuracy, a new parameter reflecting the effect of dry density and temperature was proposed and included in the final prediction model. Compared to the existing model, the predicted result of the final suggested model was closer to the measured values.

With the recent concern regarding cellulose enhancing radionuclide mobility upon its degradation to ISA, disposal of cellulosic wastes is being held off until the disposal safety is vindicated. Thus, a rational assessment should be conducted, applying an appropriate cellulose degradation model considering the disposal environment and cellulose degradation mechanisms. In this paper cellulose degradation mechanisms and the disposal environment are studied to propose the best-suitable cellulose degradation model for the domestic 1st phase repository. For the cellulose to readily degrade, the pH should be greater than 12.5. As in the case of SKB, 1BLA is excluded from the safety assessment because the pH of 1BLA remains below 12.5. Furthermore, despite cellulose degradation occurring, it does not always produce ISA. At low Ca2+ concentration, the ISA yield rate is around 25%, but at high Ca2+ concentration, the ISA yield rate increases up to 90%. Thus, for the cellulose to be a major concern, both pH and Ca2+ concentration conditions must be satisfied. To satisfy both conditions, the cement hydration must be in 2nd phase, when the porewater pH remains around 12.5 and a significant amount of Ca2+ ion is leaching out from the cement. However, according to the safety evaluation and domestic research, 2nd phase of cement hydration for silo concrete would achieve a pH of around 12.4, dissatisfying cellulose degradation condition like in 1BLA. Thus, cellulose degradation would be unlikely to occur in the domestic 1st phase repository. To derive waste acceptance criteria, a quantitative evaluation should be conducted, conservatively assuming cellulose is degraded. To conduct a safety evaluation, an appropriate degradation model should be applied to determine the degradation rate of cellulose. According to overseas research, despite the mid-chain scission being yet to be seen in the experiments, the degradation model considering mid-chain scission is applied, resulting in an almost 100% degradation rate. The model is selected because the repositories are backfilled with cement, achieving a pH greater than 13, so extensive degradation is reasonably conservative. However, under the domestic disposal condition, where cellulose degradation is unlikely to occur, applying such model would be excessively conservative. Thus, the peeling and stopping model derived by Van Loon and Haas, which suggests 10~25% degradation rate, is reasonably conservative. Based on this model, cellulose would not be a major concern in the domestic 1st phase repository. In the future, this study could be used as fundamental data for planning waste acceptance criteria.

For the performance analysis of deep geological repository systems, numerical simulation with multi-physics is required, which specifically covers Thermal (T), Hydraulic (H), and Mechanical (M) behaviors in the disposal environment. Numerous simulation models have been developed so far, each of which varies in the approach and methodology for solving THM problems. Fully-coupled THM simulation codes such as ROCMAS, THAMES, and CODE_BRIGHT were mainly developed in the initial stage of DEvelopment of COupled models and their VALidation against EXperiments (DECOVALEX), with the advantage of thorough calculations consisting of correlated several variables on different physics. Due to the difficulty of solving the complex Jacobian Matrix and the following burden for the computational calculation, weakly-coupled THM models have been suggested in recent researches: TOUGH2-MP with FLAC3D, TOUGH2 with UDEC and OpenGeoSys with FLAC3D. This methodology of loose coupling allows the practical use of computational code optimized for each physics, thereby increasing the efficiency in simulation. However, these suggested models require two different numerical codes to calculate THM behaviors, which leads to several inherent issues: compatibility during maintenance, updating and dependency between two codes. In this study, therefore, the authors build a unified code for simulating THM behaviors in the deep geological repository. The concept involves the iterative sequential coupling between TH and M for calculation efficiency. As having developed the simulation code, High-level rAdiowaste Disposal Evaluation System (HADES), to describe TH behavior based on Multi-physics Object-Oriented Simulation Environment (MOOSE) software, the authors make a milestone to develop and couple the MOOSE-based new code for M behavior as Sub-app, with the previous HADES set to be Main-app. New model for M behavior will be verified with the benchmark case of DECOVALEX-THMC Task D, comparing the mechanical simulation results: stress evolution over time, profiles of stress and vertical displacement. The existing simulation results from HADES will also be updated with the coupled calculations, with regard to temperature and saturation. Additionally, the effective stress evolution can be assessed in terms of repository’s stability with Spalling Strength and Mohr-Coulomb failure criterion. This concept for new simulation model has its meaning in that it aims to demonstrate the specific methodology of loosely coupling multi-physics in unified simulation code and analyze THM complex interactions with considering mutual influence on various physics. It is expected that HADES can be renewed as an integral simulation model for deep geological repository systems by possessing the capacity for analyzing and assessing mechanical behavior.

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.

There was researches on methodologies that provides procedures of developing corporate models and patterns of step, tasks, tactics and technologies in Information Strategy Planning(ISP). However these could not provide data input and output relations in ISP procedures, effective data and knowledge management tactics which enhance knowledge reusability. Therefore, this study aims at developing Repository model supporting Information Strategic Planning process and showing its usage and effectiveness.