• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.11-19
• /
• 2005

First, methods of numerical analysis of gas-particle flows is classified into micro, meso and macro scale approaches based on the concept of multi-scale mechanics. Next, the explanation moves on to discrete particle simulation where motion of individual particles is calculated numerically using the Newtonian equations of motion. The author focuses on the cases where particle-to-particle interaction has significant effects on the phenomena. Concerning the particle-to-particle interaction, two cases are considered: the one is collision-dominated flows and the other is the contact-dominated flows. To treat this interaction mathematically, techniques named DEM(Distinct Element Method) or DSMC (Direct Simulation Monte Carlo) have been developed DEM, which has been developed in the field of soil mechanics, is useful for the contact -dominated flows and DSMC method, developed in molecular gas flows, is for the collision-dominated flows. Combining DEM or DSMC with CFD (computer fluid dynamics), the discrete particle simulation becomes a more practical tool for industrial flows because not only the particle-particle interaction but particle-fluid interaction can be handled. As examples of simulations, various results are shown, such as hopper flows, particle segregation phenomena, particle mixing in a rotating drum, dense phase pneumatic conveying, spouted bed, dense phase fluidized bed, fast circulating fluidized bed and so on.

• Computers and Concrete
• /
• v.5 no.4
• /
• pp.329-342
• /
• 2008

A local behavior law, which includes elasticity, plasticity and damage, is developed in a three dimensional numerical model for concrete. The model is based on the Discrete Element Method (DEM)and the computational implementation has been carried out in the numerical Code YADE. This model was used to study the response of a concrete slab impacted by a rigid missile, and focuses on the extension of the compacted zone. To do so, the model was first used to simulate compression and hydrostatic tests. Once the local constitutive law parameters of the discrete element model were calibrated, the numerical model simulated the impact of a rigid missile used as a reference case to be compared to an experimental data set. From this reference case, simulations were carried out to show the importance of compaction during an impact and how it expands depending on the different impact conditions. Moreover, the numerical results were compared to empirical predictive formulae for penetration and perforation cases, demonstrating the importance of taking into account the local compaction process in the local interaction law between discrete elements.

• Journal of the Korean Geosynthetics Society
• /
• v.10 no.2
• /
• pp.55-66
• /
• 2011

This paper presents the results of an investigation on the failure mechanism of geosynthetic reinforced soil walls in back-to-back configuration using 1-g reduced-scale model tests as well as discrete element method-based numerical investigation. In the 1-g reduced scale model tests, 1/10 scale back-to-back walls were constructed so that the wall can be brought to failure by its own weight and the effect of reinforcement length on the failure mechanism was investigated. In addition, a validated discrete element method-based numerical model was used to further investigate the failure mechanism of back-to-back walls with different boundary conditions. The results were then compared with the failure mechanisms defined in the FHWA design guideline.

• Tunnel and Underground Space
• /
• v.30 no.5
• /
• pp.484-495
• /
• 2020

To numerically simulate the advance of EPB TBM, various type of numerical analysis methods have been adopted including discrete element method (DEM), finite element method (FEM), and finite difference method (FDM). In this paper, an EPB TBM driving model was proposed by using coupled DEM-FDM. In the numerical model, DEM was applied in the TBM excavation area, and contact properties of particles were calibrated by a series of triaxial tests. Since the ground around the excavation area was coupled with FDM, the horizontal stress considering the coefficient of earth pressure at rest could be applied. Also, the number of required particles was reduced and the efficiency of the analysis was increased. The proposed model can control the advance rate and rotational speed of the cutter head and screw conveyor, and derive the torque, thrust force, chamber pressure, and discharging during TBM tunnelling.

• Geomechanics and Engineering
• /
• v.7 no.6
• /
• pp.665-678
• /
• 2014

In discrete element modeling, 2D software has been widely used in order to gain further insights into the fundamental mechanisms with less computational time. The porosities used in 2D DEM studies should be determined with appropriate approaches based on 3D laboratory porosities. This paper summarizes the main approaches for converting porosities from 3D to 2D for DEM studies and theoretical evaluations show that none of the current approaches can be widely used in dealing with soil mechanical problems. Therefore, a parabolic equation and a criterion have been suggested for the determination of 2D porosities in this paper. Moreover, a case study has been used to validate that the 2D porosity obtained from the above suggestion to be rational with both the realistic contact force distribution in the specimen and the good agreement of the DEM simulation results of direct shear tests with the corresponding experimental data. Therefore, the parabolic equation and the criterion are suggested for the determination of 2D porosities in a wide range of polydisperse particle systems, especially in dealing with soil mechanical problems.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2009.09a
• /
• pp.591-604
• /
• 2009

This paper presents some important issues in modeling rock behaviour around an underground opening at depth which characterized as stress-induced fractural failure of rock. Unlike other conventional modeling approaches, stress-induced rock failure is highly complex process due to its own heterogeneous and discrete natures. Because of this complexity, many researchers has been struggled to mimic such processes as close as possible to reality with various approaches in both analytical, and numerical approaches for past few decades. Such approaches which are based on continuum mechanics, analytical fracture mechanics, and DEM(Discrete Element Method) were explored in this paper, and fundamental shortcomings for each approaches were illustrated here. In addition, DEM approach using $PFC^{2D}$(Particle Flow Code) was also implemented and illuminated in this paper and discuss the improvement and considerations for the future research.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.19 no.6
• /
• pp.1029-1044
• /
• 2017

The Discrete Element Method (DEM) is one of the useful numerical methods to analyze the behavior of the ground formation by computing the motion and interaction using particles. The DEM has not been applied in civil engineering but also a wide range of industrial fields, such as chemical engineering, pharmacy, material science, food engineering, etc. In this study, to review a performance of the spoke-type earth pressure balance (EPB) shield TBM (Tunnel Boring Machine), the commercial software based on the DEM technology was used. An analysis of the TBM during excavation was conducted according to two pre-defined excavation conditions with the different rotation speed of a cutterhead. During the analysis, the resistant torque at the face of the cutterhead, the compressive force at the cutterhead and shield surface, the muck discharge at the screw auger were measured and compared. Upon the two kinds of excavation conditions, the applicability of the DEM analysis was reviewed as a modelling method for the TBM.

• Ocean Systems Engineering
• /
• v.11 no.4
• /
• pp.313-330
• /
• 2021

In this study, the dynamic interaction between an Arctic Spar and drifting level ice is examined in time domain using the newly developed ice-hull-mooring coupled dynamics program. The in-house program, CHARM3D, which is the hull-riser-mooring coupled dynamic simulator is extended by coupling with the open-source discrete element method (DEM) simulator, LIGGGHTS. In the LIGGGHTS module, the parallel-bonding method is implemented to model the level ice using an assembly of multiple bonded spherical particles. As a case study, a spread-moored Artic Spar platform, whose hull surface near waterline is the inverted conical shape, is chosen. To determine the breaking-related DEM parameter (the critical bonding strength), the four-point numerical bending test is used. A series of numerical simulations is systematically performed under the various ice conditions including ice drift velocity, flexural strength, and thickness. Then, the effects of these parameters on the ice force, platform motions, and mooring tensions are discussed. The simulations reveal various features of dynamic interactions between the drifting ice and moored platform for various ice conditions including the novel synchronous resonance at low ice speed. The newly developed simulator is promising and can repeatedly be used for the future design and analysis including ice-floater-mooring coupled dynamics.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.33 no.3
• /
• pp.187-192
• /
• 2020

To reduce fuel consumption by reducing the weight of the deck of a dump truck and to design an eco-friendly deck, accurate structural analysis is required. To date, the load on the deck has been calculated based on the hydrostatic pressure or by applying the earth pressure theory. However, these methods cannot be used to determine the non-uniformity of the load on the deck. Load distribution varies depending on the size distribution and interaction of aggregate particles. Compared with the finite element method, the discrete element method can simulate the behavior of aggregate particles more effectively. In this study, major properties were obtained by measuring bulk density and repose. The deck of a 15 ton dump truck was simulated using the obtained properties and bumping, breaking, and turning load conditions were applied. EDEM, which is a discrete element analysis software, was employed. The stress and strain distribution of the deck were calculated by NASTRAN and compared with the measured values. The study revealed that the results derived from a DEM simulation were more accurate than those based on mathematical assumption.

• Geomechanics and Engineering
• /
• v.35 no.5
• /
• pp.475-486
• /
• 2023

Ballast particles have an irregular shape and are discrete in nature. Due to the discrete nature of ballast, it exhibits complex mechanical behaviour under loading conditions. The discrete element method (DEM) can model the behaviour of discrete particles under a multitude of loading conditions. DEM is used in this paper to simulate a series of three-dimensional direct shear tests in order to investigate the shear behaviour of railway ballast and its interaction at the microscopic level. Particle flow code in three dimension (PFC3D) models the irregular shape of ballast particles as clump particles. To investigate the influence of particle size distribution (PSD), real PSD of Indian railway ballast specification IRS:GE:1:2004, China high-speed rail (HSR) and French rail specifications are generated. PFC3D built-in linear contact model is used to simulate the interaction of ballast particles under various normal stresses, shearing rate and shear box sizes. The results indicate how shear resistance and volumetric changes in ballast assembly are affected by normal stress, shearing rate, PSD and shear box size. In addition to macroscopic behaviour, DEM represents the microscopic behaviour of ballast particles in the form of particle displacement at different stages of the shearing process.