• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.1
• /
• pp.15-23
• /
• 1998

Characteristics of turbulence beneath the free surface around a blunt bow are numerically-investigated. Three dimensional wavier-Stokes and continuity equations are solved for the simulations.. The Large Eddy Simulation(LES) with the external disturbance is performed to simulate the turbulent free surface flow called sub-breaking wave. The result shows that the free surface fluctuates beyond a certain critical condition and the characteristics of the fluctuation are similar to the turbulent boundary layer flow around a solid body.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.43 no.5
• /
• pp.621-630
• /
• 2023

In this study, to investigate an energy reduction effect by field application of cylindrical baffle arrays, the 3D Debris flow numerical analysis was conducted with various baffle configurations for the simulation of a real-scale valley, where the cylindrical baffle arrays were installed. For this, the valley of the watershed was modeled using terrestrial LiDAR data from the real-scale experiment site. Numerical analysis simulated the flow behavior of debris flow and the structures using Smooth Particle Hydrodynamics (SPH) technique of ABAQUS (Ver. 2021). The numerical analysis results that the case without cylindrical baffle arrays had a similar velocity change to that of the real-scale experiment. Also, the installation of baffles significantly reduced the frontal velocity of debris flow. Furthermore, increasing the baffle height increased the downstream energy reduction because of the higher flow impedance of taller baffles.

• Tunnel and Underground Space
• /
• v.19 no.2
• /
• pp.86-94
• /
• 2009

A three dimensional rock joint element was developed considering fracture mechanics and subcritical crack growth to simulate non-linear behavior and the progressive failure of rock joints. Using this 3-D joint element, joint shear tests of rock discontinuities were simulated by a numerical method. The asperities on the joint surface began to fail at stress levels lower than the rock fracture toughness and continued progressively due to subcritical crack growth. As a result of progressive failing in each and every asperity, the joint showed non-linear stress-time behavior including stress hardening/softening and the reaching of a residual stress.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.15 no.3
• /
• pp.345-356
• /
• 2013

Numerical simulation on rock cutting by a TBM disc cutter was carried out using SPH (Smoothed Particle Hydrodynamics) code. AUTODYN3D, a commercial software program based on finite element method, was used in this study. The three-dimensional geometry of a disc cutter and a rock specimen were modeled by Lagrange and SPH code respectively. The numerical simulation was carried out for Hwangdeung granite for 10 different cutting conditions. The results of the numerical simulation, i.e. the relation between cutter force and failure behavior, had a good agreement with those from LCM test. The cutter forces measured in the numerical simulation had 10% deviation from the LCM test results. Moreover, the optimum cutter spacing was almost identical with the experimental results. These results indicate that SPH code can be successfully used had applicability for simulation on rock cutting by a TBM disc cutter. However, further study on Lagrange-SPH coupled modelling would be necessary to reduce the computation time.

• Journal of the Korean GEO-environmental Society
• /
• v.15 no.3
• /
• pp.41-48
• /
• 2014

Experimental study of sedimentation and self-weight consolidation has been primary research area in dredged soil. However, good quality of the dredged soil and minimum water pollution caused by the pumping of reclaimed soil require intensive study of the flow characteristics of dredged material due to dumping. In this study, continuity and the equilibrium equations for mass flow assuming single phase was derived to simulate mass flow in dredged containment area. To optimize computation and modeling time for three dimensional geometry and boundary conditions, depth integration is applied to governing equations to consider three dimensional topography of the site. Petrov-Galerkin formulation is applied in spatial discretization of governing equations. Generalized trapezoidal rule is used for time integration, and Newton iteration process approximated the solution. DG and CDG technique were used for weighting matrix in discontinuous test function in dredged flow analysis, and numerical stability was evaluated by performed a square slump simulation. A comparative analysis for numerical methods showed that DG method applied to SU / PG formulation gives minimal pseudo oscillation and reliable numerical results.

• Tunnel and Underground Space
• /
• v.25 no.3
• /
• pp.275-283
• /
• 2015

Excavation of underground space in soft ground implicate to the structure, such as subsidence. As a result, it has been acting as a serious risk to the stability of the roads and facilities. Therefore, in order to stabilize the soil stabilization and reinforcement of the structure, we have been using a number of methods and injecting material. In this study, we compared and analyzed the amount of subsidence regarding the ground reinforcement during underground excavation in soft ground by performing model test. And three-dimensional numerical analysis was performed using FLAC 3D. The subsidence was simulated numerically according to the tunnel excavation. The subsidence results of the model tests and numerical analyzes were relatively consistent. Thus comparing the ground subsidence by varying the reinforcement area on the numerical analysis was analyzed. As a results, three-dimensional numerical simulation could be regarded to simulate better on the ground subsidence by various kinds of underground excavation and it can be used as a material of subsidence prevention methods.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.24 no.3
• /
• pp.231-245
• /
• 2022

When vertical shafts are constructed in soft clay with low strength, there is a risk of basal heave, which causes the excavation surface to heave due to the low bearing capacity of the ground against the imbalance of earth pressure at the excavation surface. Methods of deriving a safety factor have been proposed to evaluate the stability against the basal heave. However, there are limitations in that it is difficult to accurately evaluate the heave stability because many assumptions are included in the theoretical derivation. In this study, assuming that a circular vertical shaft is constructed in soft clay, the existing safety factor equation proposed through a theoretical approach was supplemented. Bearing capacity according to the shaft geometry, inhomogeneity of the soil, and the effect of soil plug were considered theoretically and applied in a previous safety factor equation. A three-dimensional numerical analysis was conducted to simulate the occurrence of basal heave and review the supplemented equation through various case studies. Several series of case studies were conducted targeting various factors affecting heave stability. It was verified that the additionally considered characteristics were properly reflected in the supplemented equation. Furthermore, the effects of each factor constituting the safety factor equation were examined using the results of the numerical analysis performed by simulating various cases. It was confirmed that considering the undrained shear strength increment according to depth had the most significant effect on the calculation of the safety factor.

• Journal of the Korean Geotechnical Society
• /
• v.31 no.8
• /
• pp.29-38
• /
• 2015

This paper presents the application of the Coupled Eulerian-Lagrangian (CEL) numerical technique to simulate the driving of open-ended piles into sandy soil. The main objective of this study was to investigate the applicability of CEL technique to the behavior of the driven pile penetration. Comprehensive studies to verify the behavior of driven pile penetration are presented in this paper. Through comparison with results of field load tests, the CEL methodology was found to be in good agreement with the general trend observed by in situ measurement, and the CEL approach accurately simulated the behavior of driven pipe piles.

• Journal of the Korean Geotechnical Society
• /
• v.27 no.5
• /
• pp.85-92
• /
• 2011

The dynamic behavior of piles becomes very complex due to soil-pile dynamic interaction, soil non-linearity, resonance phenomena of soil-pile system and so on. Therefore, the proper numerical simulation of the pile behavior needs much effort and calculation time. In this research, a new modeling method, which can be applied to the conventional finite difference analysis program FLAC 3D, was developed to reduce the calculation time. The soil domain in this method is divided into a near-field region and a far-field region, which is not influenced by the soil-pile dynamic interaction. Then, the ground motion of the far-field is applied to the boundaries of the near-field instead of modeling the far-field region as finite meshes. In addition, the soil non-linearity behavior is modeled by using the hysteretic damping model, which determines the soil tangent modulus as a function of shear strain and the interface element was applied to simulate the separation and slip between the soil and pile. The proposed method reduced the calculation time by as much as one third compared with a usual modeling method and maintained the accuracy of the calculated results. The calculated results by the proposed method showed a good agreement with the prototype pile behavior, which was obtained by applying a similitude law to the 1-g shaking table test results.

• Tunnel and Underground Space
• /
• v.33 no.1
• /
• pp.57-69
• /
• 2023

In this study, we simulated induced seismicity in the Groningen natural gas reservoir using 2D hydro-mechanical coupled discrete element modelling (DEM). The code used is PFC2D (Particle Flow Code 2D), a commercial software developed by Itasca, and in order to apply to this study we further developed 1)initialization of inhomogeneous reservoir pressure distribution, 2)a non-linear pressure-time history boundary condition, 3)local stress field monitoring logic. We generated a 2D reservoir model with a size of 40 × 50 km² and a complex fault system, and simulated years of pressure depletion with a time range between 1960 and 2020. We simulated fault system failure induced by pressure depletion and reproduced the spatiotemporal distribution of induced seismicity and assessed its failure mechanism. Also, we estimated the ground subsidence distribution and confirmed its similarity to the field measurements in the Groningen region. Through this study, we confirm the feasibility of the presented 2D hydro-mechanical coupled DEM in simulating the deformation of a complex fault system by hydro-mechanical coupled processes.