• Tunnel and Underground Space
• /
• v.30 no.3
• /
• pp.242-255
• /
• 2020

The numerical simulations of Heater Experiment-D (HE-D) at the Mont Terri rock laboratory in Switzerland were performed to investigate an applicability of FLAC3D to reproduce the coupled thermo-hydro-mechanical (THM) behaviour in Opalinus Clay, as part of the DECOVLEX-2015 project Task B. To investigate the reliability of numerical simulations of the coupled behaviour using FLAC3D code, the simulation results were compared with the observations from the in-situ experiment, such as temperature at 16 sensors, pore pressure at 6 sensors, and strain at 22 measurement points. An anisotropic heat conduction model, fluid flow model, and transversely isotropic elastic model in FLAC3D successfully represented the coupled thermo-hydraulic behaviour in terms of evolution for temperature and pore pressure, however, performance of the models for mechanical behavior is not satisfactory compared with the measured strain.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2011.06a
• /
• pp.1481-1482
• /
• 2011

• Tunnel and Underground Space
• /
• v.30 no.4
• /
• pp.335-358
• /
• 2020

We present a numerical model to simulate coupled hydro-mechanical behavior of fault using TOUGH-FLAC simulator. This study aims to develop a numerical method to estimate fluid injection-induced fault reactivation in low permeability rock and to access the relevant hydro-mechanical stability in rock as part of DECOVALEX-2019 Task B. A coupled fluid flow and mechanical interface model to explicitly represent a fault was suggested and validated from the applications to benchmark simulations and the field experiment at Mont Terri underground laboratory in Switzerland. The pressure build-up, hydraulic aperture evolution, displacement, and stress responses matched those obtained at the site, which indicates the capability of the model to appropriately capture the hydro-mechanical processes in rock fault.

• Tunnel and Underground Space
• /
• v.26 no.3
• /
• pp.166-180
• /
• 2016

The present study numerically simulated the CO2 injection into the saline aquifer of CO2CRC Otway pilot project and the resulting hydrological-mechanical coupled process in the storage site by TOUGH-FLAC simulator. A three-dimensional numerical model was generated using the stochastic geological model which was established based on well log and core data. It was estimated that the CO2 injection of 30,000t over a period of 200 days increased the pressure near the injection point by 0.5 MPa at the most. The pressure increased rapidly and tended to approach a certain value at an early stage of the injection. The hydrological and mechanical behavior observed from the CO2 flow, effective stress change and stress-strength ratio revealed that the CO2 injection into the saline aquifer under the given condition would not have significant effects on the mechanical safety of the storage site and the hydrological state around the adjacent fault.

• Tunnel and Underground Space
• /
• v.31 no.6
• /
• pp.593-609
• /
• 2021

We proposed a numerical method for the thermo-mechanical behavior of rock fracture using a grain-based distinct element model (GBDEM) and simulated thermally induced fracture slip. The present study is the benchmark simulation performed as part of DECOVALEX-2023 Task G, which aims to develop a numerical method to estimate the coupled thermo-hydro-mechanical processes within the crystalline rock fracture network. We represented the rock sample as an assembly of Voronoi grains and calculated the interaction of the grains (blocks) and their interfaces (contacts) using a distinct element code, 3DEC. Based on an equivalent continuum approach, the micro-parameters of grains and contacts were determined to reproduce rock as an elastic material. Then, the behavior of the fracture embedded in the rock was characterized by the contacts with Coulomb shear strength and tensile strength. In the benchmark simulation, we quantitatively examined the effects of the boundary stress and thermal stress due to heat conduction on fracture behavior, focusing on the mechanism of thermally induced fracture slip. The simulation results showed that the developed numerical model reasonably reproduced the thermal expansion and thermal stress increment, the fracture stress and displacement and the effect of boundary condition. We expect the numerical model to be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated in further study experiments.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.5 no.1
• /
• pp.43-54
• /
• 2003

In the deign of cut and cover tunnel, the structural analysis such as rigid frame analysis has been used for its simplicity and convenience. The structural analysis, however, can not consider the geological and geotechnical factors such as soil arching effect. In this study, the dominant factors influencing the behavior of cut and cover tunnel such as interface element, slope of excavation plane, distance between slope and tunnel lining, and location of slope of covered soil, were investigated by the numerical analysis to develop the analysis technique and design technology. Based on the results, the variation of bending moment, shear stress, axial force and displacements were evaluated and analyzed for each factor.

• The KIPS Transactions:PartB
• /
• v.9B no.3
• /
• pp.277-286
• /
• 2002

Genetic programming (GP) has been combined with quantum mechanical perturbation theory to make a new algorithm to construct mathematical models and perform predictions for chaotic time series from real world. Procedural similarities between time series modeling and perturbation theory to solve quantum mechanical wave equations are discussed, and the exemplary GP approach for implementing them is proposed. The approach is based on multiple populations and uses orthogonal functions for GP function set. GP is applied to original time series to get the first mathematical model. Numerical values of the model are subtracted from the original time series data to form a residual time series which is again subject to GP modeling procedure. The process is repeated until predetermined terminating conditions are met. The algorithm has been successfully applied to construct highly effective mathematical models for many real world chaotic time series. Comparisons with other methodologies and topics for further study are also introduced.

• Journal of the Korean Institute of Gas
• /
• v.20 no.4
• /
• pp.65-77
• /
• 2016

Modeling for complex reacting flow in Fischer-Tropsch reactor is one of the challenges in the field of Computational Fluid Dynamics (CFD). It is hard to derive each and every reaction rate for all chemical species because Fisher-Tropsch reaction produces many kinds of hydrocarbons which include lots of isomers. To overcome this problem, after analyzing the existing methodologies for reaction rate modeling, non-Anderson-Schulz-Flory methodology is selected to model the detailed reaction rates. In addition, the inside flow has feature of multi-phase flow, and the methodologies for modeling multi-phase flow depend on the interference between the phases, distribution of the dispersed phase, flow pattern, etc. However, existing studies have used a variety of inside flow modeling methodologies with no basis or rationale for the feasibility. Modeling inside flow based on the experimental observation of the flow would be the best way, however, with limited resources we infer the probable regime of inside flow based on conventional fluid dynamics theory; select the appropriate methodology of Mixture model; and perform systematic CFD modeling. The model presented in this study is validated through comparisons between experimental data and simulation results for 10 experimental conditions.

• Journal of the KSME
• /
• v.33 no.1
• /
• pp.22-32
• /
• 1993

본드선도는 물리적 등가가 적용되는 모든 상사 시스템을 대단히 조직이며 일관성 있게 모델링할 수있을 뿐만 아니라 시스템의 동적 방정식을 제어동역학 관점에 유리한 상태방정식으로 직접 유도할 수 있다. 또한 본드선도는 모델링 단계에서 무시된 각 요소들의 동적 특성 즉, 기계시스 템에서 회전축의 탄성처짐이나 운동 물체들 사이의 마찰 효과 전기 . 자기시스템에서 자속유출량, 그리고 유체 . 유압시스템에서 밀봉부분의 유량손실이나 압축성 유체 특성 등을 기곤에 구성된 본드선도의 변경 없이 각 특성들이 나타나는 본드상에 적절한 접합요소와 함께 단순히 첨가하여 고려할 수 있기 때문에 모델링 과정을 다시 반복하지 않고도 무시된 동적 특성들이 전체시스템에 미치는 영향을 파악할 수 있다. 특히 본드선도가 에너지 변환 장치나 에너지 유동 메카니즘이 복잡한 다에너지역 시스템 등에 거동이 전체시스템에 미치는 효과를 시각적으로 보다 세밀히 파악할 수 있을 뿐만 아니라 구성된 본드선도에 인과관계를 할당하여 본드선도의 접합요소를 통하여 유도된 최종적인 동적 방정식이 현실적으로 실현 가능한 물리시스템인지 파악할 수 있어 본드선도 모델링 단계에서 시스템을 묘사하는 동적 방정식의 옳고 그름을 평가할 수 있다. 이와 같이 본드선도는 기계, 전기 . 자기, 유체 . 유압, 열 시스템 및 이들이 조합된 복잡한 다에너지역 시스템 등을 효과적으로 모데링할 수 있는 매우 유용한 동적 시스템 모델링 방법이다.

• Proceedings of the Korean Geotechical Society Conference
• /
• 1998.05a
• /
• pp.35-81
• /
• 1998

Distinct Element Method(DEM) has a great advantage to model the discontinuous behaviour of jointed rock masses such as rotation, sliding, and separation of rock blocks. Geometrical data of joints by a field monitoring is not enough to model the jointed rock mass though the results of DE analysis for the jointed rock mass is most sensitive to the distributional properties of joints. Also, it is important to use a properly joint law in evaluating the stability of a jointed rock mass because the joint is considered as the contact between blocks in DEM. In this study, a stochastic modelling technique is developed and the dilatant rock joint is numerically modelled in order to consider th geometrical and mechanical properties of joints in DE analysis. The stochastic modelling technique provides a assemblage of rock blocks by reproducing the joint distribution from insufficient joint data. Numerical Modelling of joint dilatancy in a edge-edge contact of DEM enable to consider not only mechanical properties but also various boundary conditions of joint. Preprocess Procedure for a stochastic DE model is composed of a statistical process of raw data of joints, a joint generation, and a block boundary generation. This stochastic DE model is used to analyze the effect of deviations of geometrical joint parameters on .the behaviour of jointed rock masses. This modelling method may be one tool for the consistency of DE analysis because it keeps the objectivity of the numerical model. In the joint constitutive law with a dilatancy, the normal and shear behaviour of a joint are fully coupled due to dilatation. It is easy to quantify the input Parameters used in the joint law from laboratory tests. The boundary effect on the behaviour of a joint is verified from shear tests under CNL and CNS using the numerical model of a single joint. The numerical model developed is applied to jointed rock masses to evaluate the effect of joint dilation on tunnel stability.