• Structural Engineering and Mechanics
• /
• v.29 no.1
• /
• pp.91-111
• /
• 2008

This paper outlines the development of a computational model in order to analyze the dynamic behaviour of coupled fluid-structure systems such as a) liquid containers, b) a set of parallel or radial plates. In this work a hybrid fluid-solid element is developed, capable of simulating both membrane and bending effects of the plate. The structural mass and stiffness matrices are determined using exact integration of governing equations which are derived using a combination of classical plate theory and a finite element approach. The Bernoulli equation and velocity potential function are used to describe the liquid pressure applied on the solid-fluid element. An impermeability condition assures a permanent contact at the fluid-structure interface. Applications of this model are presented for both parallel and radial plates as well as fluid-filled rectangular reservoir. The effect of physical parameters on the dynamic behaviour of a coupled fluid-structure system is investigated. The results obtained using the presented approach for dynamic characteristics such as natural frequency are in agreement to those calculated using other theories and experiments.

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

• Tunnel and Underground Space
• /
• v.28 no.5
• /
• pp.400-425
• /
• 2018

This study presents the research results and current status of the DECOVALEX-2019 project Task B. Task B named 'Fault slip modelling' is aiming at developing a numerical method to simulate the coupled hydro-mechanical behavior of fault, including slip or reactivation, induced by water injection. The first research step of Task B is a benchmark simulation which is designed for the modelling teams to familiarize themselves with the problem and to set up their own codes to reproduce the hydro-mechanical coupling between the fault hydraulic transmissivity and the mechanically-induced displacement. We reproduced the coupled hydro-mechanical process of fault slip using TOUGH-FLAC simulator. The fluid flow along a fault was modelled with solid elements and governed by Darcy's law with the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by creating interface elements between two separating rock blocks in FLAC3D. A methodology to formulate the hydro-mechanical coupling relations of two different hydraulic aperture models and link the solid element of TOUGH2 and the interface element of FLAC3D was suggested. In addition, we developed a coupling module to update the changes in geometric features (mesh) and hydrological properties of fault caused by water injection at every calculation step for TOUGH-FLAC simulator. Then, the transient responses of the fault, including elastic deformation, reactivation, progressive evolutions of pathway, pressure distribution and water injection rate, to stepwise pressurization were examined during the simulations. The results of the simulations suggest that the developed model can provide a reasonable prediction of the hydro-mechanical behavior related to fault reactivation. The numerical model will be enhanced by continuing collaboration and interaction with other research teams of DECOLVAEX-2019 Task B and validated using the field data from fault activation experiments in a further study.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.21 no.1
• /
• pp.23-41
• /
• 2023

Coupled thermo-hydraulic-mechanical (THM) processes are essential for the long-term performance of deep geological disposal of high-level radioactive waste. In this study, a numerical sensitivity analysis was performed to analyze the effect of rock properties on THM responses after the execution of the heater test at the Kamaishi mine in Japan. The TOUGH-FLAC simulator was applied for the numerical simulation assuming a continuum model for coupled THM analysis. The rock properties included in the sensitivity study were the Young's modulus, permeability, thermal conductivity, and thermal expansion coefficients of crystalline rock, rock salt, and clay. The responses, i.e., temperature, water content, displacement, and stress, were measured at monitoring points in the buffer and near-field rock mass during the simulations. The thermal conductivity had an overarching impact on THM responses. The influence of Young's modulus was evident in the mechanical behavior, whereas that of permeability was noticed through the change in the temperature and water content. The difference in the THM responses of the three rock type models implies the importance of the appropriate characterization of rock mass properties with regard to the performance assessment of the deep geological disposal of high-level radioactive waste.

• Tunnel and Underground Space
• /
• v.31 no.5
• /
• pp.309-332
• /
• 2021

Natural barrier systems surrounding the geological repository for the high-level radioactive waste should guarantee the hydraulic performance for preventing or delaying the leakage of radionuclide. In the case of the behavior of a crystalline rock, the hydraulic performance tends to be decided by the existence of discontinuities, so the coupled hydro-mechanical(HM) processes on the discontinuities should be characterized. The discontinuum modelling can describe the complicated behavior of discontinuities including creation, propagation, deformation and slip, so it is appropriate to model the behavior of a crystalline rock. This paper investigated the coupled HM processes in discontinuum modelling such as UDEC, 3DEC, PFC, DDA, FRACOD and TOUGH-UDEC. Block-based discontinuum methods tend to describe the HM processes based on the fluid flow through the discontinuities, and some methods are combined with another numerical tool specialized in hydraulic analysis. Particle-based discontinuum modelling describes the overall HM processes based on the fluid flow among the particles. The discontinuum methods that are currently available have limitations: exclusive simulations for two-dimension, low hydraulic simulation efficiency, fracture-dominated fluid flow and simplified hydraulic analysis, so it could be improper to the modelling the geological repository. Based on the concepts of various discontinuum modelling compiled in this paper, the advanced numerical tools for describing the accurate coupled HM processes of the deep geological repository should be developed.

• Tunnel and Underground Space
• /
• v.9 no.1
• /
• pp.1-11
• /
• 1999

A finite element code was developed to analyze coupled thermo-hydro-mechanical phenomena. This code is based on the finite element formulation provided by Noorishad et al. (1984) and Joint behavior was simulated Goodman's joint constitutive model. The developed code was applied for T-H-M coupling analysis for two kinds of shaft models, with a joint or without a joint respectively. For a model without a joint, temperature increased from the shaft wall to outward evidently. The radial displacement showed opposite directions of outward and inward at some distance from shaft wall. For a model with a joint, closure of joint was found due to thermal expansion. The temperature distribution along a joint showed relatively lower than that of rock matrix because of low thermal conductivity and high specific heat of water. And it could be concluded that effects of thermal flow to joint were more than that of hydraulic flow in a rock mass.

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

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.22 no.4
• /
• pp.419-434
• /
• 2020

In this study, the damage parameters of Mazars model for KURT (KAERI Underground Research Tunnel) granite are measured form uniaxial compressive and Brazilian tests under the simulated coupled condition of a deep geological disposal. The tests are conducted in three different temperatures (15℃, 45℃, and 75℃) and dry/saturated conditions. Major model parameters such as maximum effective tensile strain (𝜖_d0), A_t, B_t, A_c, and B_c differ from the typical reference values of concrete specimens. This is likely due to the difference in elastic modulus between rock and concrete. It is found that the saturation of specimens causes an increase in value of B_t and B_c while, the rise in temperature increases 𝜖_d0 and B_t and decreases B_c. The damage model obtained from this study will be used as the primary input parameters in the development of coupled Thermo-Hydro-Mechanical Damage numerical model in KAERI.

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

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