• Tunnel and Underground Space
• /
• v.29 no.6
• /
• pp.439-455
• /
• 2019

A buffer is the major component of a high level radioactive waste repository. Due to their thermal conductivity and low permeability, bentonites have been considered as a key component of a buffer system in most countries. The deep geological condition generates ground water inflow and results in swelling pressure in the buffer and backfill. Investigation of swelling pressure of bentonite buffer is an important task for the safe disposal system. The swelling pressure that can be critical is affected by mechanical and hydro properties of the system. Therefore, in this study, a sensitivity analysis was conducted to examine the effect of hydro-mechanical (HM) behaviors in the MX-80 bentonite. Based on the results of the swelling pressure generation with HM model parameters, a coupled HM analysis of an unsaturated buffer and backfill in a deep geological repository was also carried out to investigate the major factor of the swelling pressure generation.

• 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.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.4 no.2
• /
• pp.117-131
• /
• 2006

A coupled T-H-M(Thermo-Hydro-Mechanical) analysis was carried out for KENTEX (KAERI Engineering-scale T-H-M Experiment for Engineered Barrier System), which is a facility for validating the coupled T-H-M behavior in the engineered barrier system of the Korean reference HLW(high-level waste) disposal system. The changes of temperature, water saturation, and stress were estimated based on the coupled T-H-M analysis, and the influence of the types of mechanical constitutive material laws was investigated by using elastic model, poroelastic model, and poroelastic-plastic model. The analysis was done using ABAQUS, which is a commercial finite element code for general purposes. From the analysis, it was observed that the temperature in the bentonite increased sharply for a couple of days after heating the heater and then slowly increased to a constant value. The temperatures at all locations were nearly at a steady state after about 37.5 days. In the steady state, the temperature was maintained at $90^{\circ}C$ at the interface between the heater and the bentonite and at about $70^{\circ}C$ at the interface between the bentonite and the confining cylinder. The variation of the water saturation with time in bentonite was almost same independent of the material laws used in the coupled T-H-M processes. By comparing the saturation change of T-H-M and that of H-M(Hydro-Mechanical) processes using elastic and poroelastic material mod31 respectively, it was found that the degree of saturation near the heater from T-H-M calculation was higher than that from the coupled H-M calculation mainly because of the thermal flux, which seemed to speed up the saturation. The stresses in three cases with different material laws were increased with time. By comparing the stress change in H-M calculation using poroelasetic and poroelasetic-plastic model, it was possible to conclude that the influence of saturation on the stress change is higher than the influence of temperature. It is, therefore, recommended to use a material law, which can model the elastic-plastic behavior of buffer, since the coupled T-H-M processes in buffer is affected by the variation of void ratio, thermal expansion, as well as swelling pressure.

• Coupled systems mechanics
• /
• v.6 no.4
• /
• pp.439-464
• /
• 2017

The hydro-elastic system consisting of a pre-stretched highly elastic plate, compressible Newtonian viscous fluid, and the rigid wall is considered and it is assumed that on the plate a lineal-located time-harmonic force acts. It is required to investigate the dynamic behavior of this system and determine how the problem parameters and especially the pre-straining of the plate acts on this behavior. The elasticity relations of the plate are described through the harmonic potential and linearized (with respect to perturbations caused by external time-harmonic force) form of these relations is used in the present investigation. The plane-strain state in the plate is considered and the motion of that is described within the scope of the three-dimensional linearized equations of elastic waves in elastic bodies with initial stresses. The motion of the fluid is described by the linearized Navier-Stokes equations and it is considered the plane-parallel flow of this fluid. The Fourier transform with respect to the space coordinate is applied for a solution to the corresponding boundary-value problem. Numerical results on the frequency response of the interface normal stress and normal velocity and the influence of the initial stretching of the plate on this response are presented and discussed. In particular, it is established that the initial stretching of the plate can decrease significantly the absolute values of the aforementioned quantities.

• Journal of the Korean Geotechnical Society
• /
• v.31 no.9
• /
• pp.5-15
• /
• 2015

In this study, rainfall-induced slope stability and coupling effect are investigated using hydro-mechanical finite element model. This model is developed by formulating constitutive and coupled balance equations and is verified by comparing the numerical results with field matric suction. The homogeneous soil layer (soil column) and soil slope are modeled by this model, and the results of variation in matric suction, mean effective stress, porosity, displacement, factor of safety are compared with those of staggered analysis. It is found that the vertical and horizontal displacement from coupling analysis considering change in porosity is larger than that of staggered analysis. The displacement and matric suction from coupling analysis by rainfall infiltration can affect slope instability, which shows a progressive failure behavior. The lowest factor of safety is observed under short-term rainfall. This results confirm the fact that coupling analysis is needed to design soil slope under severe rain condition.

• Journal of the Korean Geotechnical Society
• /
• v.27 no.11
• /
• pp.55-69
• /
• 2011

The importance of unsaturated state in various geo-engineering problems has led to the advance of mechanical constitutive model emulating behavior of unsaturated soils in response to thermo-hydro-mechanical loading. Elasto-plastic mechanical constitutive model for unsaturated soil is formulated based on Bishop's effective stress. Effective stress and temperature are main variables in constitutive equation, and incremental formulation of constitutive relationship is derived to compute stress update and stiffness tensor. Numerical simulations involving coupled THM processes are conducted to discuss numerical stability and applicability of developed constitutive model: one-dimensional test, tri-axial compression test, and clay-buffering at high level radioactive waste disposal. Numerical results demonstrated that developed model can predict very complex behavior of coupled THM phenomena and is applicable to geo-engineering problems under various environmental conditions, as well as interpret typical behavior of unsaturated soils.

• Tunnel and Underground Space
• /
• v.30 no.4
• /
• pp.320-334
• /
• 2020

This study presents an application of hydro-mechanical coupled Particle Flow Code 3D (PFC3D) to simulation of fluid injection induced fault slip experiment conducted in Mont Terri Switzerland as a part of a task in an international research project DECOVALEX-2019. We also aimed as identifying the current limitations of the modelling method and issues for further development. A fluid flow algorithm was developed and implemented in a 3D pore-pipe network model in a 3D bonded particle assembly using PFC3D v5, and was applied to Mont Terri Step 2 minor fault activation experiment. The simulated results showed that the injected fluid migrates through the permeable fault zone and induces fault deformation, demonstrating a full hydro-mechanical coupled behavior. The simulated results were, however, partially matching with the field measurement. The simulated pressure build-up at the monitoring location showed linear and progressive increase, whereas the field measurement showed an abrupt increase associated with the fault slip We conclude that such difference between the modelling and the field test is due to the structure of the fault in the model which was represented as a combination of damage zone and core fractures. The modelled fault is likely larger in size than the real fault in Mont Terri site. Therefore, the modelled fault allows several path ways of fluid flow from the injection location to the pressure monitoring location, leading to smooth pressure build-up at the monitoring location while the injection pressure increases, and an early start of pressure decay even before the injection pressure reaches the maximum. We also conclude that the clay filling in the real fault could have acted as a fluid barrier which may have resulted in formation of fluid over-pressurization locally in the fault. Unlike the pressure result, the simulated fault deformations were matching with the field measurements. A better way of modelling a heterogeneous clay-filled fault structure with a narrow zone should be studied further to improve the applicability of the modelling method to fluid injection induced fault activation.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.11 no.2
• /
• pp.141-150
• /
• 2009

When constructing a subsea tunnel in discontinuous rock mass, fluid flow in joints has a great influence on the behavior of the tunnel so that hydro-mechanical coupled analysis should be performed for the stability estimation. In practice, relaxed rock load is generally used for the design of tunnel concrete lining. In a continuum analysis, a method based on the distribution of local safety factor around a tunnel was proposed for the estimation of a potential relaxed zone. However, in the case of discontinuous rock mass in which joints are developed, the whole stability of tunnels depends on the behavior of the joints. In this study, therefore, a method is proposed for the estimation of a potential relaxed zone occurred by the excavation of a tunnel in discontinuous rock mass. The suggested method is validated by sensitivity analysis and the comparison with the results of continuum analysis.

• Tunnel and Underground Space
• /
• v.30 no.1
• /
• pp.39-62
• /
• 2020

In this study, Barcelona Basic Model (BBM) was implemented into TOUGH2-MP/FLAC3D for the numerical analysis of coupled thermo-hydro-mechanical (THM) behavior of unsaturated soils and the prediction of long-term behaviors. Similar to the methodology described in a previous study for the implementation of BBM into TOUGH-FLAC, the User Defined Model (UDM) of FLAC based on the Modified Cam Clay Model (MCCM) and the FISH function of FLAC3D were used to extend the existing MCCM module in FLAC3D for the implementation of BBM into TOUGH2-MP/FLAC3D. In the developed BBM module in TOUGH2-MP/FLAC3D, the plastic strains due to change in suction increase (SI) in addition to mean effective stress are calculated. In addition to loading-collapse (LC) yield surface, suction increase (SI) yield surface is changed by hardening rules in the developed BBM module. Several numerical simulations were conducted to verify and validate the implementation of BBM: using an example presented in the FLAC3D manual for the standard MCCM, simulation results using COMSOL, and experimental data presented in SKB Reports. In addition, the developed BBM analysis module was validated by simultaneously performing a series of modeling tests that were performed for the validation of the Quick tools developed for the purpose of effectively deriving BBM parameters, and by comparing the Quick tools and Code_Bright results reported in a previous study.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.33 no.4
• /
• pp.1489-1498
• /
• 2013

The monolithically coupled finite element analysis for a deformable unsaturated soil slope is performed to investigate the effect of antecedent rainfall which is assumed by initial conditions varying degree of saturation (36, 51, 77%) in finite element analysis. The distributions of matric suction and deformation on slope surface obtained from numerical simulation show the instability of antecedent rainfall-induced unsaturated soil slope. Moreover, the numerical analysis using Drucker-Prager model can be checked if a soil slope has reached failure (trial failure criterion $f^{tr}$ >0, plastic behavior) or not (trial failure criterion $f^{tr}$ < 0, elastic behavior). It is found that displacement of slope surface layer increases and the matric suction on soil slope decreases with an increase of initial degree of saturation by antecedent rainfall. Especially, the matric suction of the soil slope in dry condition (S=36%) rapidly decreases rather than that in wet condition (S=51%) at the same rainfall duration. The results of the trial failure criterion ($f^{tr}$ > 0) show slope instability in the toe region and surface of the slopes.