• Tunnel and Underground Space
• /
• v.30 no.5
• /
• pp.446-461
• /
• 2020

Numerical simulations of CIEMAT column test in Spain are performed to investigate the coupled thermo-hydro-mechanical (THM) behavior of MX80 bentonite pellets using TOUGH2-FLAC3D. The heater power and injection pressure of water in the numerical simulations are identical to those in the laboratory test. To investigate the applicability of the thermo-hydraulic (TH) model used in TOUGH2 code to prediction of the coupled TH behavior, the simulation results are compared with the observations of temperature and relative humidity with time. The tendencies of the coupled behavior observed in the test are well represented by the numerical models and the simulator in terms of temperature and relative humidity evolutions. Moreover, the performance of the models for the reproduction and prediction of the coupled TH behavior is globally satisfactory compared with the observations. However, the calculated stress change is relatively small and slow due to the limitations of the simple elastic and swelling pressure model used in numerical simulations. It seems that the two models are insufficient to realistically reproduce the complex coupled THM behavior in the bentonite pellets.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.17 no.2
• /
• pp.183-202
• /
• 2019

For design and performance assessment of a high-level radioactive waste (HLW) disposal system, it is necessary to understand the characteristics of coupled thermo-hydro-mechanical (THM) behavior. However, in previous studies for the Korean Reference HLW Disposal System (KRS), thermal analysis was performed to determine the spacing of disposal tunnels and interval of disposition holes without consideration of the coupled THM behavior. Therefore, in this study, TOUGH2-MP/FLAC3D is used to conduct THM modeling for performance assessment of the Korean Reference HLW Disposal System (KRS). The peak temperature remains below the temperature limit of $100^{\circ}C$ for the whole period. A rapid rise of temperature caused by decay heat occurs in the early years, and then temperature begins to decrease as decay heat from the waste decreases. The peak temperature at the bentonite buffer is around $96.2^{\circ}C$ after about 3 years, and peak temperature at the rockmass is $68.2^{\circ}C$ after about 17 years. Saturation of the bentonite block near the canister decreases in the early stage, because water evaporation occurs owing to temperature increase. Then, saturation of the bentonite buffer and backfill increases because of water intake from the rockmass, and bentonite buffer and backfill are fully saturated after about 266 years. The stress is calculated to investigate the effect of thermal stress and swelling pressure on the mechanical behavior of the rockmass. The calculated stress is compared to a spalling criterion and the Mohr-Coulumb criterion for investigation of potential failure. The stress at the rockmass remains below the spalling strength and Mohr-Coulumb criterion for the whole period. The methodology of using the TOUGH2-MP/FLAC3D simulator can be applied to predict the long-term behavior of the KRS under various conditions; these methods will be useful for the design and performance assessment of alternative concepts such as multi-layer and multi-canister concepts for geological spent fuel repositories.

• Tunnel and Underground Space
• /
• v.21 no.5
• /
• pp.394-405
• /
• 2011

In this paper, we applied coupled non-isothermal, multiphase fluid flow and geomechanical numerical modeling using TOUGH-FLAC coupled analysis to study the complex thermodynamic and geomechanical performance of underground lined rock caverns (LRC) for compressed air energy storage (CAES). Mechanical stress in concrete linings as well as pressure and temperature within a storage cavern were examined during initial and long-term operation of the storage cavern for CAES. Our geomechanical analysis showed that effective stresses could decrease due to air penetration pressure, and tangential tensile stress could develop in the linings as a result of the air pressure exerted on the inner surface of the lining, which would result in tensile fracturing. According to the simulation in which the tensile tangential stresses resulted in radial cracks, increment of linings' permeability and air leakage though the linings, tensile fracturing occurred at the top and at the side wall of the cavern, and the permeability could increase to $5.0{\times}10^{-13}m^2$ from initially prescribed $10{\times}10^{-20}m^2$. However, this air leakage was minor (about 0.02% of the daily air injection rate) and did not significantly impact the overall storage pressure that was kept constant thanks to sufficiently air tight surrounding rocks, which supports the validity of the concrete-lined underground caverns for CAES.

• Tunnel and Underground Space
• /
• v.31 no.4
• /
• pp.221-242
• /
• 2021

A numerical study of the performance assesment of coupled thermo-hydro-mechanical (THM) processes in improved Korean reference disposal system (KRS⁺) for high-level radioactive waste is conducted using TOUGH2-MP/FLAC3D simulator. Decay heat from high-level radioactive waste increases the temperature of the repository, and it decreases as decay heat is reduced. The maximum temperature of the repository is below a maximum temperature criterion of 100℃. Saturation of bentonite buffer adjacent to the canister is initially reduced due to pore water evaporation induced by temperature increase. Bentonite buffer is saturated 250 years after the disposal of high-level radioactive waste by inflow of groundwater from the surrounding rock mass. Initial saturation of rock mass decreases as groundwater in rock mass is moved to bentnonite buffer by suction, but rock mass is saturated after inflow of groundwater from the far-field area. Stress changes at rock mass are compared to the Mohr-Coulomb failure criterion and the spalling strength in order to investigate the potential rock failure by thermal stress and swelling pressure. Additional simulations are conducted with the reduced spacing of deposition holes. The maximum temperature of bentonite buffer exceeds 100℃ as deposition hole spacing is smaller than 5.5 m. However, temperature of about 56.1% volume of bentonite buffer is below 90℃. The methodology of numerical modeling used in this study can be applied to the performance assessment of coupled THM processes for high-level radioactive waste repositories with various input parameters and geological conditions such as site-specific stress models and geothermal gradients.

• Tunnel and Underground Space
• /
• v.30 no.4
• /
• pp.359-381
• /
• 2020

Within the framework of DECOVALEX-2019 Task D, full-scale engineered barriers experiment (FEBEX) at Grimsel Test Site was numerically simulated to investigate an applicability of implemented Barcelona basic model (BBM) into TOUGH2-MP/FLAC3D simulator, which was developed for the prediction of the coupled thermo-hydro-mechanical behavior of bentonite buffer. And the calculated heater power, temperature, relative humidity, total stress, saturation, water content and dry density were compared with in situ data monitored in the various sections. In general, the calculated heater power and temperature provided a fairly good agreement with experimental observations, however, the difference between power of heater #1 and that of heater #2 could not captured in the numerical analysis. It is necessary to consider lamprophyre with low thermal conductivity around heater #1 and non-simplified installation progresses of bentonite blocks in the tunnel for better modeling results. The evolutions and distributions of relative humidity were well reproduced, but hydraulic model needs to be modified because the re-saturation process was relatively fast near the heaters. In case of stress evolutions due to the thermal and hydraulic expansions, the computed stress was in good agreement with the data. But, the stress is slightly higher than the measured in situ data at the early stage of the operation, because gap between rock mass and bentonite blocks have not been considered in the numerical simulations. The calculated distribution of saturation, water content, and dry density along the radial distance showed good agreement with the observations after the first and final dismantling. The calculated dry density near the center of the FEBEX tunnel and heaters were overestimated compared with the observations. As a result, the saturation and water content were underestimated with the measurements. Therefore, numerical model of permeability is needed to modify for the production of better numerical results. It will be possible to produce the better analysis results and more realistically predict the coupled THM behavior in the bentonite blocks by performing the additional studies and modifying the numerical model based on the results of this study.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.12 no.4
• /
• pp.267-274
• /
• 2014

The first Low- and Intermediate-Level Waste (LILW) disposal facility with 6 silos has been constructed in granite host rock saturated with groundwater in Korea. A two-dimensional numerical modeling on gas migration was carried out using TOUGH2 with EOS5 module in the disposal facility. Laboratory-scale experiments were also performed to measure the important properties of silo concrete related with gas migration. The gas entry pressure and relative gas permeability of the concrete was determined to be $0.97{\pm}0.15bar$ and $2.44{\times}10^{-17}m^2$, respectively. The results of the numerical modeling showed that hydrogen gas generated from radioactive wastes was dissolved in groundwater and migrated to biosphere as an aqueous phase. Only a small portion of hydrogen appeared as a gas phase after 1,000 years of gas generation. The results strongly suggested that hydrogen gas does not accumulate inside the disposal facility as a gas phase. Therefore, it is expected that there would be no harmful effects on the integrity of the silo concrete due to gas generation.

• Geomechanics and Engineering
• /
• v.13 no.1
• /
• pp.1-23
• /
• 2017

A coupled liquid-gas-solid three-phase model, linking two numerical codes (TOUGH2/EOS3 and $FLAC^{3D}$), was firstly established and validated by simulating an in-situ air flow test in Essen. Then the coupled model was employed to investigate responses of multiphase flow and soil skeleton deformation to compressed air or freshwater injection using the same simulation conditions in an aquifer of Tianjin, China. The simulation results show that with injecting pressurized fluids, the vertical effective stress in some area decreases owing to the pore pressure increasing, an expansion of soil skeleton appears, and land uplift occurs due to support actions from lower deformed soils. After fluids injection stops, soil deformation decreases overall due to injecting fluids dissipating. With the same applied pressure, changes in multiphase flow and geo-mechanical deformation caused by compressed air injection are relatively greater than those by freshwater injection. Furthermore, the expansion of soil skeleton induced by compressed air injection transfers upward and laterally continuously with time, while during and after freshwater injection, this expansion reaches rapidly a quasi-steady state. These differences induced by two fluids injection are mainly because air could spread upward and laterally easily for its lower density and phase state transition appears for compressed air injection.