• International Journal of Highway Engineering
• /
• v.15 no.5
• /
• pp.21-29
• /
• 2013

PURPOSES: Simulation of aggregate slump test using equivalent sphere particle in DEM and its validity evaluation against lab aggregate slump test METHODS : In this research, aggregate slump tests are performed and compared with DEM simulation. To utilize spheric particles in YADE, equivalent sphere diameter concept is applied. As verification measures, the volume in slump cone filled with aggregate is used and it is compared with volume in slump cone filled with equivalent sphere particle. Slump height and diameter are also used to evaluate the suggested numerical method with equivalent concept RESULTS : Simulation test results show good agrement with lab test results in terms of loose packing volume, height and diameter of slumped particle clump. CONCLUSIONS : It is concluded that numerical simulation using DEM is applicable to evaluate the effect of aggregate morphological property in loose packing and optimum gradation determination based on the aggregate slump test simulation result.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.18 no.3
• /
• pp.59-65
• /
• 2019

In many industries, particle packing is adopted quite frequently. In the particle packing process, the Discrete Element Method (DEM) can analyze the multi-collision of particles efficiently. Two types of contact models are frequently used for the DEM. One is the linear spring model, which has the fastest calculation time, and the other is the Hertz-Mindlin model, which is the most frequently used contact model employing the DEM. Meanwhile, very tiny particles in the micrometer order are used in modern industries. In the micro length order, surface force is important to decreased particle size. To consider the effect of surface force in this study, we performed a simulation with the Hertz-Mindlin model and added the Johnson-Kendall-Roberts (JKR) theory depicting surface force with surface energy. In addition, three contact models were compared with several parameters. As a result, it was found that the JKR model has larger residual stress than the general contact models because of the pull-off force. We also validated that surface force can influence particle behavior if the particles are small.

• Journal of the Korean Society for Railway
• /
• v.6 no.3
• /
• pp.174-178
• /
• 2003

This paper presents the development of 2-dimensional discrete element algorithm to generate circle and line elements for the simulation of the ballast and sleeper in railway. An example of randomly distributed circle elements show a good applicability of this algorithm for the modeling of the behaviors of ballast. The output about unbalaned force, particle velocity, and total energy conservation from the code is evaluated to check if the calculation is conducted properly.

• Korean Journal of Materials Research
• /
• v.26 no.11
• /
• pp.611-622
• /
• 2016

Particle morphology change and different experimental condition analysis during composite fabrication process by traditional ball milling with discrete element method (DEM) simulation were investigated. A simulation of the three dimensional motion of balls in a traditional ball mill for research on the grinding mechanism was carried out by DEM simulation. We studied the motion of the balls, the ball behavior energy and velocity; the forces acting on the balls were calculated using traditional ball milling as simulated by DEM. The effect of the operational variables such as the rotational speed, ball material and size on the flow velocity, collision force and total impact energy were analyzed. The results showed that increased rotation speed with interaction impact energy between balls and balls, balls and pots and walls and balls. The rotation speed increases with an increase of the impact energy. Experiments were conducted to quantify the grinding performance under the same conditions. Furthermore, the results showed that ball motion affects the particle morphology, which changed from irregular type to plate type with increasing rotation speed. The evolution was also found to depend on the impact energy increase of the grinding media. These findings are useful to understand and optimize the particle motion and grinding behavior of traditional ball mills.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.3197-3202
• /
• 2007

An analysis program for pedestrian flow has been developed to investigate the flow patterns of passenger in railway stations. Analysis algorithms for pedestrian flow based on DEM(Discrete Element Method) are newly developed. There are lots of similarity between particle-laden two phase flow and passenger flow. The velocity component of 1st phase corresponds to the unit vector of calculation cell, each particle to passenger, volume fraction to population density and the particle velocity to the walking velocity, etc. And, the walking velocity of passenger is also represented by the function of population density. Key algorithms are developed to determine the position of passenger, population density and numbering to each passenger. By using the developed program, we compared the simulation results of the effects of the location and size of exit and elapsed time.

• International Journal of Highway Engineering
• /
• v.19 no.4
• /
• pp.37-44
• /
• 2017

PURPOSES : Conductive and convective heat transfer simulations for an asphalt mixture were made by using discrete element method (DEM) and similarity principle. METHODS : In this research, virtual specimens composed of discrete element method particles were generated according to four different predetermined particle size distribution curves. Temperature variations of the four different particles for a given condition were estimated and were compared with measurements and analytical solutions. RESULTS : The virtual specimen with mixed particles and with the smallest particle show very good agreement with laboratory test results and analytical solutions. As particle size decreases, better heat transfer simulation can be performed due to smaller void ratio and more contact points and areas. In addition, by utilizing the similarity principle of thermal properties and corresponding time unit, analytical time can be drastically reduced. CONCLUSIONS : It is concluded that the DEM asphalt mixture specimens with similarity principle could be used to predict the temperature variation for a given condition. It is observed that the void ratio has critical effect on prediction of temperature variation. Comparing the prediction for a 4 mm particle specimen with a mixed particle specimen, it is also concluded that predicting the mixed particle specimen temperature is much more efficient considering the number of particles that are directly associated with computational time in DEM analysis.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.1-1
• /
• 2009

In this presentation, I talk about various fluid simulation methods that have been developed for computer graphics special effects since 1996. They are all based on CFD but sacrifice physical reality for visual plausability and time. But as the speed of computer increases rapidly and the capability of GPU (graphics processing unit) improves, methods for more physical realism have been tried. In this talk, I will focus on four aspects of fluid simulation methods for computer graphics: (1) particle level-set methods, (2) particle-based simulation, (3) methods for exact satisfaction of incompressibility constraint, and (4) GPU-based simulation. (1) Particle level-set methods evolve the surface of fluid by means of the zero-level set and a band of massless marker particles on both sides of it. The evolution of the zero-level set captures the surface in an approximate manner and the evolution of marker particles captures the fine details of the surface, and the zero-level set is modified based on the particle positions in each step of evolution. (2) Recently the particle-based Lagrangian approach to fluid simulation gains some popularity, because it automatically respects mass conservation and the difficulty of tracking the surface geometry has been somewhat addressed. (3) Until recently fluid simulation algorithm was dominated by approximate fractional step methods. They split the Navier-Stoke equation into two, so that the first one solves the equation without considering the incompressibility constraint and the second finds the pressure which satisfies the constraint. In this approach, the first step introduces error inevitably, producing numerical diffusion in solution. But recently exact fractional step methods without error have been developed by fluid mechanics scholars), and another method was introduced which satisfies the incompressibility constraint by formulating fluid in terms of vorticity field rather than velocity field (by computer graphics scholars). (4) Finally, I want to mention GPU implementation of fluid simulation, which takes advantage of the fact that discrete fluid equations can be solved in parallel.

• Proceedings of the KSME Conference
• /
• 2008.11a
• /
• pp.1922-1927
• /
• 2008

In the present study, deposition of discrete and small particles, which diameter is less than $1{\mu}m$, on a filter element was simulated by stochastic method. Trajectory of each particle was numerically solved by Langevin equation and Brownian random motion was treated by Brownian dynamics. Lattice Boltzmann method (LBM) was used to solve flow field around the filter collector and deposit layer. Interaction between flow field and deposit layer was obtained from a converged solution from an inner-loop calculation. Simulation method is properly validated and collection efficiency due to different filtration parameters are examined and discussed. Morphology of deposit layer and its evolution was visualized in terms of the particle size. The particle loaded effect on collection efficiency was also discussed.

• International Journal of Highway Engineering
• /
• v.15 no.5
• /
• pp.31-45
• /
• 2013

PURPOSES : Determination of particle packing model variables that can be used for formulation of new DEM based particle packing model by examining existing particle packing models METHODS : Existing particle packing models are thoroughly examined by analytical reformulation and sensitivity analysis in order to set up DEM based new particle packing model and to determine its variables. All model equations considered in this examination are represented with consistent expressions and are compared to each others to find mathematical and conceptual similarity in expressions. RESULTS : From the examination of existing models, it is observed that the models are very similar in their shapes although the derivation of the models may be different. As well, it is observed that variables used in some existing models are comprehensive enough to estimate particle packing but not applicable to DEM simulation. CONCLUSIONS : A set of variables that can be used in DEM based particle packing model is determined.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.38 no.7
• /
• pp.751-756
• /
• 2014

The dynamic behavior of numerous grains interacting with each other can be easily observed. In this study, this dynamic behavior was analyzed based on the contact between numerous grains. The discrete element method was used for analyzing the dynamic behavior of each particle and the neighboring-cell algorithm was employed for detecting their contact. The Hertzian and tangential sliding friction contact models were used for calculating the contact force acting between the particles. A GPU-based parallel program was developed for conducting the computer simulation and calculating the numerous contacts. The dam break experiment was performed to verify the simulation results. The reliability of the program was verified by comparing the results of the simulation with those of the experiment.