• Tunnel and Underground Space
• /
• v.30 no.6
• /
• pp.573-590
• /
• 2020

This study presents the current status of DECOVALEX-2023 project Task G and our research results so far. Task G, named 'Safety ImplicAtions of Fluid Flow, Shear, Thermal and Reaction Processes within Crystalline Rock Fracture NETworks (SAFENET)' aims at developing a numerical method to simulate the fracture creation and propagation, and the coupled thermohydro-mechanical processes in fracture in crystalline rocks. The first research step of Task G is a benchmark simulation, which is designed for research teams to make their modelling codes more robust and verify whether the models can represent an analytical solution for displacements of a single rock fracture. We reproduced the mechanical behavior of rock and embedded single fracture using a three-dimensional grain-based distinct element model for the simulations. In this method, the structure of the rock was represented by an assembly of rigid tetrahedral grains moving independently of each other, and the mechanical interactions at the grains and their contacts were calculated using 3DEC. The simulation results revealed that the stresses induced along the embedded fracture in the model were relatively low compared to those calculated by stress analysis due to stress redistribution and constrained fracture displacements. The fracture normal and shear displacements of the numerical model showed good agreement with the analytical solutions. The numerical model will be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated using various experiments in a further study.

• Tunnel and Underground Space
• /
• v.31 no.1
• /
• pp.52-65
• /
• 2021

In this study, authors analyzed the vapor cloud explosion induced by propane leak at the PEMIX Terminal, which is the propane storage facility outside of Mexico City. TNT equivalence mass for the leaked 4750 kg propane was estimated to be 9398 kg. Blast parameters such as peak overpressure, positive phase duration, and impact at 40-400 (m) away from the center of the explosion were calculated by applying TNT Equivalency Method and Multi-Energy Method. The probability of damage due to lung damage, eardrum rupture, head impact, and whole-body displacement impact by applying the probit function obtained using blast parameters was evaluated. The peak overpressure obtained using Multi-Energy Method was found to be greater than the peak overpressure obtained by applying the TNT Equivalency Method at all distances considered, but it was evaluated that there was no significant difference from the points above 200 m. The peak overpressure obtained by Multi-Energy Method was computed to assess the extent of damage to the structure, and it was shown that structures within 100 m of the explosion center would collapse completely, and that the glasses of the structures 400 m away would be almost broken. The probability of death due to lung damage was shown to vary depending on a human body's position located in the propagating direction of shock wave, and if there is a reflecting surface in the immediate surroundings of a human body, the probability of death was estimated to be the greatest. The impact of shock wave on lung damage, eardrum rupture, head impact, and whole-body displacement impact was evaluated and found to affect whole-body impact < lung damage < eardrum rupture

• Tunnel and Underground Space
• /
• v.30 no.6
• /
• pp.509-518
• /
• 2020

This study constructed a new simulation platform, including mesh generation process, numerical simulation, and post-processing for results analysis based on exploration data to perform real-scale numerical analysis considering the actual geological structure efficiently. To build the simulation platform, we applied for open-source programs. The source code is open to be available for code modification according to the researcher's needs and compatibility with various numerical simulation programs. First, a three-dimensional model(3D) is acquired based on the exploration data obtained using a drone. Then, the domain's mesh density was adjusted to an interpretable level using Blender, the free and open-source 3D creation suite. The next step is to create a 3D numerical model by creating a tetrahedral volume mesh inside the domain using Gmsh, a finite element mesh generation program. To use the mesh information obtained through Gmsh in a numerical simulation program, a converting process to conform to the program's mesh creation protocol is required. We applied a Python code for the procedure. After we completed the stability analysis, we have created various visualization of the study using ParaView, another open-source visualization and data analysis program. We successfully performed a preliminary stability analysis on the full-scale Dokdo model based on drone-acquired data to confirm the usefulness of the proposed platform. The proposed simulation platform in this study can be of various analysis processes in future research.

• Tunnel and Underground Space
• /
• v.31 no.3
• /
• pp.145-166
• /
• 2021

Korea Stress Map database is built by integrating actual data of 1,400 in-situ stress measurements using hydraulic fracturing and overcoring method in South Korea. Korea Stress Map 2020 is presented based on the guideline proposed by World Stress Map Project. As detailed data, stress ratio and maximum horizontal stress direction distribution for each region are also presented. The dominant maximum horizontal stress direction in the Korean Peninsula is from northeast to southeast, and the magnitude of the in-situ stress is relatively distributed. There is some stress heterogeneity caused by local characteristics such as topographical and geological properties. We investigated case studies in which the in-situ stress was affected by mountainous topography, difference in rock quality of fracture zone, presence of mine or underground cavities, and geological structure of fault zone.

• Wind and Structures
• /
• v.34 no.1
• /
• pp.127-136
• /
• 2022

In this work a multi-fidelity non-intrusive polynomial chaos (MF-NIPC) has been applied to a structural wind engineering problem in architectural design for the first time. In architectural design it is important to design structures that are safe in a range of wind directions and speeds. For this reason, the computational models used to design buildings and bridges must account for the uncertainties associated with the interaction between the structure and wind. In order to use the numerical simulations for the design, the numerical models must be validated by experi-mental data, and uncertainties contained in the experiments should also be taken into account. Uncertainty Quantifi-cation has been increasingly used for CFD simulations to consider such uncertainties. Typically, CFD simulations are computationally expensive, motivating the increased interest in multi-fidelity methods due to their ability to lev-erage limited data sets of high-fidelity data with evaluations of more computationally inexpensive models. Previous-ly, the multi-fidelity framework has been applied to CFD simulations for the purposes of optimization, rather than for the statistical assessment of candidate design. In this paper MF-NIPC method is applied to flow around a rectan-gular 5:1 cylinder, which has been thoroughly investigated for architectural design. The purpose of UQ is validation of numerical simulation results with experimental data, therefore the radius of curvature of the rectangular cylinder corners and the angle of attack are considered to be random variables, which are known to contain uncertainties when wind tunnel tests are carried out. Computational Fluid Dynamics (CFD) simulations are solved by a solver that employs the Finite Element Method (FEM) for two turbulence modeling approaches of the incompressible Navier-Stokes equations: Unsteady Reynolds Averaged Navier Stokes (URANS) and the Large Eddy simulation (LES). The results of the uncertainty analysis with CFD are compared to experimental data in terms of time-averaged pressure coefficients and bulk parameters. In addition, the accuracy and efficiency of the multi-fidelity framework is demonstrated through a comparison with the results of the high-fidelity model.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.25 no.2
• /
• pp.141-155
• /
• 2023

Monitoring technology based on electrical resistivity is widely used for non-destructive data collection and health analysis of underground structures and tunnels. Vulnerable sections such as fault zone generates many problems during construction of the tunnel. These problems cause displacement and stress changes of the ground. Therefore, it is necessary to predict the state of the fault zone section to ensure the mechanical stability of the underground structure. Monitoring the size of joints and the porosity of the fillings is essential for rocks. Previous studies have not considered the variety of fillings in rock joints. In this study, electrical resistivity tests were conducted according to the particle mixing state of the sandy fillings. When the size of fillings is decreased at the constant porosity, the electrical resistivity tends to increase. The results of this study are expected to be useful as basic electrical resistivity data for predicting the ground conditions and evaluation of the ground behavior that is containing sandy fillings in the rock joint for tunnels.

• Korean Journal of Mineralogy and Petrology
• /
• v.36 no.1
• /
• pp.33-40
• /
• 2023

Todorokite is a manganese oxide mineral containing Mg²⁺ in a tunnel structure in which MnO₆ octahedra share corners. In order to investigate the suitability and efficiency of high temperature-treated todorokite as a material for adsorption and fixation of Cs, Cs was ion exchanged and the amount of leached Cs from todorokite was measured. The todorokite used in this study was synthesized by transforming Na-birnessite to Mg-buserite and used as a precursor. After high temperature treatment, Cs exchanged todorokite changed to birnessite and hausmannite as the temperature increased. The amount of leached Cs was investigated for Cs exchanged todorokite which was reacted with distilled water and 1 M NaCl solution at different reaction times. In general, for the samples reacted with 1 M NaCl solution, the fixation of Cs was quite effective, although the amount of leached Cs was greater due to the ion exchange reaction with Na. As the treatment temperature increased, the amount of leached Cs increased and then decreased again, which was related to the mineral phases formed at each temperature. As birnessite was formed, the amount of leached Cs increased, but as birnessite decreased, that decreased again. As the mineral phase changed to hausemanite, the amount of Cs decreased rapidly. The results of our study show that Cs exchanged todorokite can be used as a material that effectively fixes Cs and prevents its diffusion by high temperature treatment.

• Journal of the Society of Disaster Information
• /
• v.20 no.1
• /
• pp.156-168
• /
• 2024

Purpose: The trend of highway traffic accidents shows a repeating pattern of increase and decrease, with the fatality rate being highest on highways among all road types. Therefore, there is a need to establish improvement measures that reflect the situation within the country. Method: We conducted accident severity analysis using Random Forest on data from accidents occurring on 10 specific routes with high accident rates among national highways from 2019 to 2021. Factors influencing accident severity were identified. Result: The analysis, conducted using the SHAP package to determine the top 10 variable importance, revealed that among highway traffic accidents, the variables with a significant impact on accident severity are the age of the perpetrator being between 20 and less than 39 years, the time period being daytime (06:00-18:00), occurrence on weekends (Sat-Sun), seasons being summer and winter, violation of traffic regulations (failure to comply with safe driving), road type being a tunnel, geometric structure having a high number of lanes and a high speed limit. We identified a total of 10 independent variables that showed a positive correlation with highway traffic accident severity. Conclusion: As accidents on highways occur due to the complex interaction of various factors, predicting accidents poses significant challenges. However, utilizing the results obtained from this study, there is a need for in-depth analysis of the factors influencing the severity of highway traffic accidents. Efforts should be made to establish efficient and rational response measures based on the findings of this research.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.20 no.3
• /
• pp.543-560
• /
• 2018

With the increased development in downtown underground space facilities that vertically cross under a railway at a shallow depth, the demand for non-open cut method is increasing. However, most construction sites still adopt the pipe roof method, where medium and large diameter steel pipes are pressed in to form a roof, enabling excavation of the inside space. Among the many factors that influence the loosening region and loads that occur while pressing in steel pipes, the size of the pipe has the largest impact, and this factor may correspond to the magnitude of load applied to the underground structure inside the steel pipe roof. The super equilibrium method (SEM) has been developed to minimize ground disturbance and loosening load, and uses small diameter pipes of approximately 114 mm instead of conventional medium and large diameter pipes. This small diameter steel pipe is called an SEM pile. After SEM piles are pressed in and the grouting reinforcement is constructed, a crossing structure is pressed in by using a hydraulic jack without ground subsidence or heaving. The SEM pile, which plays the role of timbering, is a fore-poling pile of approximately 5 m length that prevents ground collapse and supports surface load during excavation of toe part. The loosening region should be adequately calculated to estimate the spacing and construction length of the piles and stiffness of members. In this paper, we conducted a comparative analysis of calculations of loosening load that occurs during the press-in of SEM pile to obtain an optimal design of SEM. We analyzed the influence of factors in main theoretical and empirical formulas applied for calculating loosening regions, and carried out FEM analysis to see an appropriate loosening load to the SEM pile. In order to estimate the soil loosening caused by actual SEM-pile indentation and excavation, a steel pipe indentation reduction model test was conducted. Soil subsidence and soil loosening were investigated quantitatively according to soil/steel pipe (H/D).

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