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.
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.
In this study, TbDyFe thin films with the thickness of 1000 Å are fabricated by DC magnetron sputtering. TbDyFe thin films are prepared by DC magnetron sputtering method. The pressure of Ar gas below 1.33 kPa and DC input power of 200 W are used for the sputtering conditions. During sputtering process the substrate holder is heated up to 150 ℃. The thin films are deposited to a thickness of 1000 Å on polyimide substrate with a thickness of 2 μm. The fabricated microstructures are observed by X-ray diffraction (XRD) and the film thickness is measured. Magnetostrictions are determined from the curvature of the thin films which are measured by the optical cantilever method. The experimental results are discussed with numerical data.