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.

Geologic disposal at deep depth is an acceptable way to dispose of high-level radioactive waste and isolate it from the biosphere. The geological repository system comprises an engineered barrier system (EBS) and the host rock. The system aims to delay radionuclide migration through groundwater flow, and also, the flow affects the saturation of the bentonite in the EBS. The thermal conductivity of bentonite is a function of saturation, so the temperature in the EBS is directly related to the flow system. High-temperature results in the two-phase flow, and the two-phase flow system also affects the flow system. Therefore, comprehending the influencing parameters on the flow system is critical to ensure the safety of the disposal system. Various studies have been performed to figure out the complex two-phase flow characteristics, and numerical simulation is considered an effective way to predict the coupled behavior. DECOVALEX (DEvelopment of COupled models and their VALidation against EXperiments) is one of the most famous international cooperating projects to develop numerical methods for thermo-hydro-mechanicalchemical interaction, and Task C in the DECOVALEX-2023 has the purpose of simulating the Fullscale Emplacement (FE) experiment at the Mont-Terri underground research laboratory. We used OGS-FLAC, a self-developed numerical simulator combining OpenGeoSys and FLAC3D, for the simulation and targeted to analyze the effecting parameters on the two-phase flow system. We focused on the parameters of bentonite, a key component of the disposal system, and analyzed the effect of compressibility and air entry pressure on the flow system. Compressibility is a parameter included in the storage term, defining the fluid storage capacity of the medium. While air entry pressure is a crucial value of the water retention curve, defining the relation between saturation and capillary pressure. From a series of sensitivity analyses, low compressibility resulted in faster flow due to low storage term, while low air entry pressure slowed flow inflow into the bentonite. Low air entry pressure means the air easily enters the medium; hence the flow rate becomes lower based on the relativity permeability definition. Based on the sensitivity analysis, we further investigate the effect of shotcrete around the tunnel and excavation damaged zone. Also, long-term analysis considering heat decay of the radioactive waste will be considered in future studies.

In high-level radioactive waste disposal, a high temperature is generated from the canister containing the waste in the engineered barrier, while groundwater flows into the buffer system from the host rock. The temperature increase and groundwater inflow result in the water phase change and saturation variation. Saturation change is related to the thermal conductivity of buffer material; hence the phase change and saturation strongly interact with the temperature evolution. The complex coupled behavior affects the stability of the whole disposal system, and the security of the repository is critical to human-being life. However, it is difficult to predict the long-term coupled behavior in the disposal system due to the considerable field-test scale, and therefore a numerical simulation is a suitable method having repeatability and cost-effectiveness. DECOVALEX is an international cooperating project for developing numerical methods and models for thermo-hydro-mechanical-chemical (THMC) interaction. DECOVALEX has a four-year cycle with various topics. At the current phase, Task C aims to simulate the full-scale emplacement (FE) experiment performed at Mont Terri underground rock laboratory. Nine research groups are participating in the task, and among them, KAERI simulates the experiment using OGS-FLAC. The simulator combines OpenGeoSys for TH simulation and FLAC3D for M simulation. Through the benchmark simulation, we verified OGS-FLAC for the two-phase flow analysis in the disposal system and finally modeled the FE experiment with a three-dimensional grid. We performed a simple sensitivity analysis to investigate the effect of input parameters on the two-phase flow system and confirmed that the compressibility and permeability affected the flow behavior. We also compared the simulation results to the field data and obtained well-matched results from a series of simulation.

Geologic disposal of high-level radioactive waste is considered the most effective method to isolate high-level radioactive waste from the biosphere. A high-level radioactive waste repository is designed to be placed at a deep depth and generally consists of canisters, buffer material, and host rock. In the disposal system, the heat from the canister occurs for millions of years due to the long half-life of the high-level radioactive waste, and the heat induces vaporization of groundwater in the buffer material. The resaturation process also occurs due to groundwater inflow from the host rock by the hydraulic head and capillarity. The saturation variation leads to the heat transfer and multi-phase flow in the buffer material, and thermal pressurization of groundwater due to the heat affects the effective stress change in the host rock. The stress change can make the porosity and permeability change in the flow system of the host rock, and the flow system affects the nuclide migration to the biosphere. Therefore, it is crucial to understand the complex thermo-hydro-mechanical-chemical (THMC) coupled behavior to secure the repository’s long-term safety. DECOVALEX is an international cooperating project to develop numerical methods and models for predicting the THMC interactions in the disposal systems through validation and comparison with test results. In Task C of DECOVALEX-2023, nine participating groups (BGR, BGE, CAS, ENSI, GRS, KAERI, LBNL, NWMO, Sandia) models the full-scale emplacement (FE) experiments at the Mont Terri underground rock laboratory and focus on understanding pore pressure development, heat transfer, thermal pressurization, vaporization and resaturation process in the disposal system. In the FE experiment, three heaters generated heat with constant power for five years at a 1:1 scale in the emplacement tunnel based on Nagra’s reference repository design. KAERI used OGS-FLAC3D for the numerical simulation, combining OpenGeoSys for TH simulation and FLAC3D for M simulation. We generated a full-scale three-dimensional numerical model with a dimension of 100 by 100 by 60 meters. The pressure and temperature distribution were well simulated with the host rock's anisotropy. Based on the capillarity, we observed vaporization and resaturation in the bentonite under the twophase flow system. We plan to compare the simulation results with the field data and investigate the effect of input parameters, including thermal conductivity and pore compressibility affecting the thermal and flow system.

We present an optical-infrared photometric study of galaxies in six nearby clusters of galaxies at z=0.041∼0.098 (A1436, A1773, A1809, A2048, A2142, and A2152). Using BV I photometry obtained at the Bohyunsan Optical Astronomical Observatory and JHKS photometry extracted from the 2-Micron All-Sky Survey catalog, we investigate the colors of galaxies in the clusters. Using the (B − V ) versus (I −KS) color-color diagrams in comparison with the simple stellar population model, we estimate the ages and metallicities of bright early-type member galaxies. Early-type galaxies in each cluster show the color-magnitude relation. Ages and metallicities of early-type members show little dependence on their velocity dispersions. Mean ages of early-types in the clusters range from 3 Gyr to 20 Gyr, showing a large dispersion, and mean metallicities range from Z = 0.03 to 0.05 above the solar value, showing a negligible dispersion.

We present a photometric study of galaxies in the central regions of six nearby galaxy clusters at redshift z=0.0231~0.0951. We have derived BVI photometry of the galaxies from the CCD images obtained at the Bohyunsan Optical Astronomical Observatory (BOAO) in Korea, and JHKs photometry of the bright galaxies from the 2MASS extended source catalog. Comparing the galaxy photometry results with the simple stellar population model of Bruzual & Charlot (2003) in the optical & NIR color-color diagrams, we have estimated the ages and metallicities of early type galaxies. We have found that the observed galaxies had recent star-formation mostly 5 ~ 7 Gyrs ago but the spread in age estimation is rather large. The average metallicities are [Fe/H]=0.l~0.5 dex. These results support the hypothesis that large early type galaxies in clusters are formed via hierarchical merging of smaller galaxies.