• Tunnel and Underground Space
• /
• 제26권2호
• /
• pp.59-67
• /
• 2016

In this paper, we perform a numerical analysis to evaluate the potential of fault reactivation during gas production from hydrate bearing sediments and the moment magnitude of induced seismicity. For the numerical analysis, sequential coupling of TOUGH+Hydrate and FLAC3D was used and the change in effective stress and consequent geomechanical deformation including fault reactivation was simulated by assuming that Mohr-Coulomb shear resistance criterion is valid. From the test production simulation of 30 days, we showed that pore pressure reduction as well as effective stress change hardly induces the fault reactivation in the vicinity of a production well. We also investigated the influence of stress state conditions to a fault reactivation, and showed that normal fault stress regime, where vertical stress is relatively greater than horizontal, may have the largest potential for the reactivation. We tested one simulation that earthquake can be induced during gas production and calculated the moment magnitude of the seismicity. Our calculation presented that all the magnitudes from the calculation were negative values, which indicates that induced earthquakes can be grouped into micro-seismic and as small as hardly perceived by human beings. However, it should be noted that the current simulation was carried out using the highly simplified geometric model and assumptions such that the further simulations for a scheduled test production and commercial scale production considering complex geometric conditions may produce different results.

• Smart Structures and Systems
• /
• 제25권4호
• /
• pp.487-499
• /
• 2020

This paper proposes a novel approach to model updating for a large-scale cable-stayed bridge based on ambient vibration tests coupled with a hybrid metaheuristic search algorithm. Vibration measurements are carried out under excitation sources of passing vehicles and wind. Based on the measured structural dynamic characteristics, a finite element (FE) model is updated. For long-span bridges, ambient vibration test (AVT) is the most effective vibration testing technique because ambient excitation is freely available, whereas a forced vibration test (FVT) requires considerable efforts to install actuators such as shakers to produce measurable responses. Particle swarm optimization (PSO) is a famous metaheuristic algorithm applied successfully in numerous fields over the last decades. However, PSO has big drawbacks that may decrease its efficiency in tackling the optimization problems. A possible drawback of PSO is premature convergence leading to low convergence level, particularly in complicated multi-peak search issues. On the other hand, PSO not only depends crucially on the quality of initial populations, but also it is impossible to improve the quality of new generations. If the positions of initial particles are far from the global best, it may be difficult to seek the best solution. To overcome the drawbacks of PSO, we propose a hybrid algorithm combining GA with an improved PSO (HGAIPSO). Two striking characteristics of HGAIPSO are briefly described as follows: (1) because of possessing crossover and mutation operators, GA is applied to generate the initial elite populations and (2) those populations are then employed to seek the best solution based on the global search capacity of IPSO that can tackle the problem of premature convergence of PSO. The results show that HGAIPSO not only identifies uncertain parameters of the considered bridge accurately, but also outperforms than PSO, improved PSO (IPSO), and a combination of GA and PSO (HGAPSO) in terms of convergence level and accuracy.

• Membrane Journal
• /
• 제10권3호
• /
• pp.139-147
• /
• 2000

Non-Fickian (or anomalous) diffusion was observed in transient sorption of aromatic solvents(such as benzene, toluene, and chlorobenzene) in fluoropolymers (such as ETFE, ECTFE and PVDF). In this study, five other transient sorption models (Crank, Long ＆ Richman, Berens ＆ Hopfenberg, Neogi, Li) based on Fick's law were employed to fit the anomalous sorption data for aromatic solvents. The adjustable parameters were determined by least square analysis of the measured and predicted fractional uptake. For ETFE sorption data slightly deviating from Fickian behavior, all the models exhibited satisfactory results in fitting the anomalous sorption data. In particular, Neogj model predicted intrinsic diffusivity (0.4～0.8$\times$10$^{-5}$ $\textrm{cm}^2$/day) and equilibrium diffusivity (0.13～0.31$\times$10$^{-4}$ $\textrm{cm}^2$/day) as well as relaxation kinetics related to non-Fickain diffusion. For a typical sigmoidal sorption behavior in PVDF, only Crank's model could give the reasonable evaluation on transport properties. The ratio of intial diffusivity (D$_{i}$) to final equilibrium diffusivity (D$_{\infty}$) was ranged from 80 to 200. For the final stage of uptake In ECTFE with drastic acceleration, all the models exhibited significant deviations from the sorption data. New diffusion models based on thermodynamics and continuum mechanics should be employed to get valuable information on transport properties as well as relaxation kinetics coupled with non-Fickian diffusion.

• Tunnel and Underground Space
• /
• 제23권5호
• /
• pp.372-382
• /
• 2013

Thermal-hydro-mechanical (THM) modeling is a critical R&D issue in the performance and safety assessment of a high-level waste repository. With an $\ddot{A}$sp$\ddot{o}$ prototype repository, its thermal behavior was analyzed and then compared with in-situ experimental data for its validation. A model simulation was used to calculate the temperature distributions in the deposition holes, deposition tunnel, and surrounding host rock. A comparison of the simulation results with the experimental data was made for deposition hole DH-6, which showed that there was a temperature difference of $2{\sim}5^{\circ}C$ depending on the location of the measuring points, but there was a similar trend in the evolution curves of temperature as a function of time. It was expected that the coupled modeling of the thermal behavior with the hydro-mechanical behavior in the buffer and backfill of the $\ddot{A}$sp$\ddot{o}$ prototype repository would give a better agreement between the experimental and model calculation results.

• Tunnel and Underground Space
• /
• 제23권6호
• /
• pp.493-505
• /
• 2013

This numerical study investigates seismicity and fault slip induced by fluid injection in deep geothermal reservoir with pre-existing fractures and fault. Particle Flow Code 2D is used with additionally implemented hydro-mechanical coupled fluid flow algorithm and acoustic emission moment tensor inversion algorithm. The output of the model includes spatio-temporal evolution of induced seismicity (hypocenter locations and magnitudes) and fault deformation (failure and slip) in relation to fluid pressure distribution. The model is applied to a case of fluid injection with constant rates changing in three steps using different fluid characters, i.e. the viscosity, and different injection locations. In fractured reservoir, spatio-temporal distribution of the induced seismicity differs significantly depending on the viscosity of the fracturing fluid. In a fractured reservoir, injection of low viscosity fluid results in larger volume of induced seismicity cloud as the fluid can migrate easily to the reservoir and cause large number and magnitude of induced seismicity in the post-shut-in period. In a faulted reservoir, fault deformation (co-seismic failure and aseismic slip) can occur by a small perturbation of fracturing fluid (<0.1 MPa) can be induced when the injection location is set close to the fault. The presented numerical model technique can practically be used in geothermal industry to predict the induced seismicity pattern and magnitude distribution resulting from hydraulic stimulation of geothermal reservoirs prior to actual injection operation.

• Journal of the Microelectronics and Packaging Society
• /
• 제16권2호
• /
• pp.43-50
• /
• 2009

Although thin PCBs(Printed Circuit Boards) have recently been required for high density interconnection, high electrical performance, and low manufacturing cost, the utilization of thin PCBs is severely limited by warpage and reliability issues. Warpage of the thin PCB leads to failure in solder-joints and chip. The $B^2it$(Buried Bump Interconnection Technology) for PCB has been developed to achieve a competitive manufacturing price. In this study, chip temperature, package warpage, chip stress and solder-joints stress characteristics of the PCB prepared with $B^2it$ process have been calculated using thermo-mechanical coupled analysis by the FEM(Finite Element Method). FEM computation was carried out with the variations in bump shapes and kinds of materials under 1.5W power of chip and constant convection heat transfer. The results show that chip temperature distribution reached more quickly steady-state status with PCB prepared with $B^2it$ process than PCB prepared with conventional via interconnection structure. Although $B^2it$ structures are effective on low package warpage and chip stress, with high strength bump materials arc disadvantage for low stress of solder-joints. Therefore, it is recommended that optimized bump shapes and materials in PCB design should be considered in terms of reliability characteristics in the packaging level.

• Journal of the Korean Society of Safety
• /
• 제11권4호
• /
• pp.143-150
• /
• 1996

This is an explicit-Implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE ) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for and existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The Ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory Darcy's law Is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full newton-Raphson method is used for solving the nonlinear governing equations.

• Tunnel and Underground Space
• /
• 제29권3호
• /
• pp.197-213
• /
• 2019

We simulated the fault reactivation experiment conducted at 'Main Fault' intersecting the low permeability clay formations of Mont Terri Underground Research Laboratory in Switzerland 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. We formulate the hydro-mechanical coupling relation of hydraulic aperture to consider the elastic fracture opening and failure-induced dilation for reproducing the abrupt changes in injection flow rate and monitoring pressure at fracture opening pressure. A parametric study was conducted to examine the effects of in-situ stress condition and fault deformation and strength parameters and to find the optimal parameter set to reproduce the field observations. In the best matching simulation, the fracture opening pressure and variations of injection flow rate and monitoring pressure showed good agreement with field experiment results, which suggests the capability of the numerical model to reasonably capture the fracture opening and propagation process. The model overestimated the fault displacement in shear direction and the range of reactivated zone, which was attributed to the progressive shear failures along the fault at high injection pressure. In the field experiment results, however, fracture tensile opening seems the dominant mechanism affecting the hydraulic aperture increase.

• Tunnel and Underground Space
• /
• 제29권4호
• /
• pp.262-279
• /
• 2019

In the engineering barriers of high-level radioactive waste disposal, gases could be generated through a number of processes. If the gas production rate exceeds the gas diffusion rate, the pressure of the gas increases and gases could migrate through the bentonite buffer. Because people and the environment can be exposed to radioactivity, it is very important to clarify gas migration in terms of long-term integrity of the engineered barrier system. In particular, it is necessary to identify the hydro-mechanical mechanism for the dilation flow, which is a very important gas flow phenomenon only in medium containing large amounts of clay materials such as bentonite buffer, and to develop and validate new numerical approach for the quantitative evaluation of the gas migration phenomenon. Therefore, in this study, we developed a two-phase flow model considering the mechanical damage model in order to simulate the gas migration in the engineered barrier system, and validated with 1D gas flow modelling through saturated bentonite under constant volume boundary conditions. As a result of numerical analysis, the rapid increase in pore water pressure, stress, and gas outflow could be simulated when the dilation flow was occurred.

• Tunnel and Underground Space
• /
• 제31권5호
• /
• pp.333-359
• /
• 2021

Gases such as hydrogen and radon can be generated around the canister in high-level radioactive waste disposal systems due to several reasons including the corrosion of metal materials. When the gas generation rate exceeds the gas diffusion rate in the low-permeability bentonite buffer, the gas phase will form and accumulate in the engineered barrier system. If the gas pressure exceeds the gas entry pressure, gas can migrate into the bentonite buffer, resulting in pathway dilation flow and advective flow. Because a sudden occurrence of dilation flow can cause radionuclide leakage out of the engineered barrier of the radioactive waste disposal system, it is necessary to understand the gas migration behavior in the bentonite buffer to quantitatively evaluate the long-term safety of the engineered barrier. Experimental research investigating the characteristics of gas migration in saturated bentonite and research developing numerical models capable of simulating such behaviors are being actively conducted worldwide. In this technical note, previous gas injection experiments and the numerical models proposed to verify such behaviors are introduced, and the future challenges necessary for the investigation of gas migration are summarized.