An Innovative Coupled Common-Node Discrete Element Method-Smoothed Particle Hydrodynamics Model Developed with LS-DYNA and Its Applications

In this study,a common-node DEM-SPH coupling model based on the shared node method is proposed,and a fluid–structure coupling method using the common-node discrete element method-smoothed particle hydrodynamics(DS-SPH)method is developed using LS-DYNA software.The DEM and SPH are established on the same node to create common-node DEM-SPH particles,allowing for fluid–structure interactions.Numerical simulations of various scenarios,including water entry of a rigid sphere,dam-break propagation over wet beds,impact on an ice plate floating on water and ice accumulation on offshore structures,are conducted.The interaction between DS particles and SPH fluid and the crack generation mechanism and expansion characteristics of the ice plate under the interaction of structure and fluid are also studied.The results are compared with available data to verify the proposed coupling method.Notably,the simulation results demonstrated that controlling the cutoff pressure of internal SPH particles could effectively control particle splashing during ice crushing failure.

Study on the particle breakage of ballast based on a GPU accelerated discrete element method

Breakage of particles will have greatly influence on mechanical behavior of granular material(GM)under external loads,such as ballast,rockfill and sand.The discrete element method(DEM)is one of the most popular methods for simulating GM as each particle is represented on its own.To study breakage mechanism of particle breakage,a cohesive contact mode is developed based on the GPU accelerated DEM code-Blaze-DEM.A database of the 3D geometry model of rock blocks is established based on the 3D scanning method.And an agglomerate describing the rock block with a series of non-overlapping spherical particles is used to build the DEM numerical model of a railway ballast sample,which is used to the DEM oedometric test to study the particles’breakage characteristics of the sample under external load.Furthermore,to obtain the meso-mechanical parameters used in DEM,a black-analysis method is used based on the laboratory tests of the rock sample.Based on the DEM numerical tests,the particle breakage process and mechanisms of the railway ballast are studied.All results show that the developed code can better used for large scale simulation of the particle breakage analysis of granular material.

Predicting the Dynamic Behavior of Asphalt Concrete Using Three-dimensional Discrete Element Method

A user-defined three-dimensional （3D） discrete element model was presented to predict the dynamic modulus and phase angle of asphalt concrete （AC）. The 3D discrete element method （DEM） model of AC was constructed employing a user-defined computer program developed using the ＂Fish＂ language in PFC3D. Important microstructural features of AC were modeled, including aggregate gradation, air voids and mastic. The irregular shape of aggregate particle was modeled using a clump of spheres. The developed model was validated through comparing with experimental measurements and then used to simulate the cyclic uniaxial compression test, based on which the dynamic modulus and phase angle were calculated from the output stress- strain relationship. The effects of air void content, aggregate stiffness and volumetric fraction on AC modulus were further investigated. The experimental results show that the 3D DEM model is able to accurately predict both dynamic modulus and phase angle of AC across a range of temperature and loading frequencies. The user- defined 3D model also demonstrated significant improvement over the general existing two-dimensional models.

Micromechanical modeling of asphalt concrete fracture using a user-defined three-dimensional discrete element method

A user-defined micromechanical model was developed to investigate the fracture mechanism of asphalt concrete （AC） using the discrete element method （DEM）. A three-dimensional （3D） AC beam was built using the ＂Fish＂ language provided by PFC3D and was employed to simulate the three-point bending beam test at two temperature levels： -10 ℃ and 15℃. The AC beam was modeled with the consideration of the microstructural features of asphalt mixtures. Uniaxial complex modulus test and indirect tensile strength test were conducted to obtain material input parameters for numerical modeling. The 3D predictions were validated using laboratory experimental measurements of AC beams prepared by the same mixture design. Effects of mastic stiffness, cohesive and adhesive strength on AC fracture behavior were investigated using the DEM model. The results show that the 3D DEM fracture model can accurately predict the fracture patterns of asphalt concrete. The ratio of stress at interfaces to the stress in mastics increases as the mastic stiffness decreases; however, the increase in the cohesive strength or adhesive strength shows no significant influence on the tensile strength.

Microstructure modeling and virtual test of asphalt mixture based on three-dimensional discrete element method

The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete element method.A randomly generating algorithm was proposed to capture the three-dimensional irregular shape of coarse aggregate.And then,modeling algorithm and method for graded aggregates were built.Based on the combination of modeling of coarse aggregates,asphalt mastic and air voids,three-dimensional virtual sample of asphalt mixture was modeled by using PFC^(3D).Virtual tests for penetration test of aggregate and uniaxial creep test of asphalt mixture were built and conducted by using PFC^(3D).By comparison of the testing results between virtual tests and actual laboratory tests,the validity of the microstructure modeling and virtual test built in this study was verified.Additionally,compared with laboratory test,the virtual test is easier to conduct and has less variability.It is proved that microstructure modeling and virtual test based on three-dimensional discrete element method is a promising way to conduct research of asphalt mixture.

Applicability of discrete element method with spherical and clumped particles for constitutive study of granular materials

Discrete element method(DEM)has been intensively used to study the constitutive behaviour of granular materials.However,to what extent a real granular material can be reproduced by virtual DEM simulations remains unclear.This study attempts to answer this question by comparing DEM simulations with typical features of experimental granular materials.Three groups of models with spherical and clumped particles are investigated from four perspectives:(i)deviatoric stress and volumetric behaviour;(ii)critical state behaviour;(iii)stress-dilatancy relationship;and(iv)the evolution of principal stress ratio against axial strain.The results demonstrate that DEM with spherical or clumped particles is capable of qualitatively describing macroscopic deviatoric stress responses,volumetric behaviour,and critical state behaviour observed in experiments for granular materials.On the other hand,some qualitative deviations between experiments and the investigated DEM simulations are also observed,in terms of the stress-dilatancy behaviour and principal stress ratio against axial strain,which are proven to be critical for constitutive modelling.The results demonstrate that DEM with spherical or clumped particles may not necessarily fully capture experimental features of granular materials even from a qualitative perspective.It is thus encouraged to thoroughly validate DEM with experiments when developing constitutive models based on DEM observations.

A new method to simulate irregular particles by discrete element method

Granular materials are ubiquitous in nature and important in various applications such as road and railway engineering. Granular materials exhibit complicated mechanical behaviors, which are affected significantly by the irregular shape of particles. Currently, the discrete element method (DEM) has been accepted as an effective approach to investigate the mechanical behaviors of granular materials. However, there are scarce simulations based on DEM in literatures considering the irregularity of particle shape. A new method is proposed to simulate individual real particle with irregular shape using clump constituted by overlapping spheres. First, the geometric model of real particle with surface nodes and inner nodes is established through digitally processing the computerized tomography (CT) scanning data. Second, a clump consisting of spheres is generated to simulate the real particle using a minimum distance criterion. The criterion is implemented by tree optimization algorithm. Influential factors are also introduced to balance the model accuracy and computing cost. Effects of the influential factors, including the density of geometric grid and the minimum distance, on simulations are discussed. Results show that this new method is simpler and more efficient than the previous methods in terms of the model accuracy and computing cost.

Coupling of discrete-element method and smoothed particle hydrodynamics for liquid-solid flows

Particle based methods can be used for both the simulations of solid and fluid phases in multiphase medium, such as the discrete-element method for solid phase and the smoothed particle hydrodynamics for fluid phase. This paper presents a computational method combining these two methods for solid-liquid medium. The two phases are coupled by using an improved model from a reported Lagrangian-Eulerian method. The technique is verified by simulating liquid-solid flows in a two-dimensional lid-driven cavity.

Discrete Element Modeling of Asphalt Concrete Cracking Using a User-def ined Three-dimensional Micromechanical Approach

We established a user-defined micromechanical model using discrete element method （DEM） to investigate the cracking behavior of asphalt concrete （AC）. Using the ＂Fish＂ language provided in the particle flow code in 3-Demensions （PFC3D）, the air voids and mastics in asphalt concrete were realistically built as two distinct phases. With the irregular shape of individual aggregate particles modeled using a clump of spheres of different sizes, the three-dimensional （3D） discrete element model was able to account for aggregate gradation and fraction. Laboratory uniaxial complex modulus test and indirect tensile strength test were performed to obtain input material parameters for the numerical simulation. A set of the indirect tensile test were simulated to study the cracking behavior of AC at two levels of temperature, i e, -10 ℃ and 15 ℃. The predicted results of the numerical simulation were compared with laboratory experimental measurements. Results show that the 3D DEM model is able to predict accurately the fracture pattern of different asphalt mixtures. Based on the DEM model, the effects of air void content and aggregate volumetric fraction on the cracking behavior of asphalt concrete were evaluated.

Discrete particle methods for engineering simulation:Reproducing mesoscale structures in multiphase systems

Most natural resources are processed as particle-fluid multiphase systems in chemical,mineral and material indus-tries,therefore,discrete particles methods(DPM)are reasonable choices of simulation method for engineering the relevant processes and equipments.However,direct application of these methods is challenged by the complex multiscale behavior of such systems,which leads to enormous computational cost or otherwise qualitatively inac-curate description of the mesoscale structures.The coarse-grained DPM based on the energy-minimization multi-scale(EMMS)model,or EMMS-DPM,was proposed to reduce the computational cost by several orders while main-taining an accurate description of the mesoscale structures,which paves the way for its engineering applications.Further empowered by the high-efficiency multi-scale DEM software DEMms and the corresponding customized heterogeneous supercomputing facilities with graphics processing units(GPUs),it may even approach realtime simulation of industrial reactors.This short review will introduce the principle of DPM,in particular,EMMS-DPM,and the recent developments in modeling,numerical implementation and application of large-scale DPM which aims to reach industrial scale on one hand and resolves mesoscale structures critical to reaction-transport coupling on the other hand.This review finally prospects on the future developments of DPM in this direction.

THE APPLICATION OF DISCRETE ELEMENT METHOD IN SOLVING THREE-DIMENTIONAL IMPACT DYNAMICS PROBLEMS

A three-dimensional discrete element model of the connective type is presented. Moreover,a three-dimensional numerical analysis code,which can carry out the transitional pro- cess from connective model(for continuum)to contact model(for non-continuum),is developed for simulating the mechanical process from continuum to non-continuum.The wave propagation process in a concrete block(as continuum)made of cement grout under impact loading is numer- ically simulated with this code.By comparing its numerical results with those by LS-DYNA,the calculation accuracy of the model and algorithm is proved.Furthermore,the failure process of the concrete block under quasi-static loading is demonstrated,showing the basic dynamic tran- sitional process from continuum to non-continuum.The results of calculation can be displayed by animation.The damage modes are similar to the experimental results.The two numerical examples above prove that our model and its code are powerful and efficient in simulating the dynamic failure problems accompanying the transition from continuum to non-continuum.It also shows that the discrete element method(DEM)will have broad prospects for development and application.

Numerical analysis of behavior of active layer in rotary kilns by discrete element method

The behavior of the active layer of material bed within rotary kilns plays a key role in industrial applications.To obtain its influences on industrial process,different regimes of particle motion have been simulated by discrete element method(DEM) in three dimensions under variant rotation speeds,filling degree,based on the background of induration process of iron ore pellets.The influences of the mentioned factors on the maximum thickness of the active layer and the average velocity of particles have been investigated.The average velocity of particles increases with Froude number following the power function over a wide range,and the maximum thickness rises with increasing rotation speed in a way of logarithm.The influence of the filling degree f on the maximum thickness exhibits a good linearity under two classic regimes,but the increasing of the average velocity of the active layer is limited at f=0.4.This basic research highlights the impact of the active layer within rotary kilns,and lays a good foundation for the further investigation in mixing and heat transfer within the particle bed inside rotary kilns.

Parallel Computing of Discrete Element Method on GPU

General purpose computing on GPU for scientific computing has been rapidly growing in recent years. We investigate the applicability of GPU to discrete element method (DEM) often used in particle motion simulation. NVIDIA provides a sample code for this type of simulation, which obtained superior performance than CPU in computational time. A computational model of the contact force in NVIDIA’s sample code is, however, too simple to use in practice. This paper modifies the NVIDIA’s simple model by replacing it with the practical model. The computing speed of the practical model on GPU is compared with the simple one on GPU and with the practical one on CPU in numerical experiments. The result shows that the practical model on GPU obtains the computing speed 6 times faster than the practical one on CPU while 7 times slower than that of the simple one on GPU. The effects of the GPU architectures on the computing speed are analyzed.

Capability of discrete element method to investigate the macro-micro mechanical behaviours of granular soils considering different stress conditions and morphological gene mutation

Discrete element method(DEM)has been widely utilised to model the mechanical behaviours of granular materials.However,with simplified particle morphology or rheology-based rolling resistance models,DEM failed to describe some responses,such as the particle kinematics at the grain-scale and the principal stress ratio against axial strain at the macro-scale.This paper adopts a computed tomography(CT)-based DEM technique,including particle morphology data acquisition from micro-CT(mCT),spherical harmonic-based principal component analysis(SH-PCA)-based particle morphology reconstruction and DEM simulations,to investigate the capability of DEM with realistic particle morphology for modelling granular soils’micro-macro mechanical responses with a consideration of the initial packing state,the morphological gene mutation degree,and the confining stress condition.It is found that DEM with realistic particle morphology can reasonably reproduce granular materials’micro-macro mechanical behaviours,including the deviatoric stressevolumetric straineaxial strain response,critical state behaviour,particle kinematics,and shear band evolution.Meanwhile,the role of multiscale particle morphology in granular soils depends on the initial packing state and the confining stress condition.For the same granular soils,rougher particle surfaces with a denser initial packing state and a higher confining stress condition result in a higher degree of shear strain localisation.

Computer vision-aided DEM study on the compaction characteristics of graded subgrade filler considering realistic coarse particle shapes

The compaction quality of subgrade filler strongly affects subgrade settlement.The main objective of this research is to analyze the macro-and micro-mechanical compaction characteristics of subgrade filler based on the real shape of coarse particles.First,an improved Viola-Jones algorithm is employed to establish a digitalized 2D particle database for coarse particle shape evaluation and discrete modeling purposes of subgrade filler.Shape indexes of 2D subgrade filler are then computed and statistically analyzed.Finally,numerical simulations are performed to quantitatively investigate the effects of the aspect ratio(AR)and interparticle friction coefficient(μ)on the macro-and micro-mechanical compaction characteristics of subgrade filler based on the discrete element method(DEM).The results show that with the increasing AR,the coarse particles are narrower,leading to the increasing movement of fine particles during compaction,which indicates that it is difficult for slender coarse particles to inhibit the migration of fine particles.Moreover,the average displacement of particles is strongly influenced by the AR,indicating that their occlusion under power relies on particle shapes.The dis-placement and velocity of fine particles are much greater than those of the coarse particles,which shows that compaction is primarily a migration of fine particles.Under the cyclic load,the interparticle friction coefficientμhas little effect on the internal structure of the sample;under the quasi-static loads,however,the increase inμwill lead to a significant increase in the porosity of the sample.This study could not only provide a novel approach to investigate the compaction mechanism but also establish a new theoretical basis for the evaluation of intelligent subgrade compaction.

Comparison of Construction Method for DEM Simulation of Ellipsoidal Particles

Discrete element model was developed to simulate the ellipsoidal particles moving in the moving bed. Multi-element model was used to describe a ellipsoidal particle, the contact detection algorithm of ellipsoidal particle was developed, and both contact force and gravity force were considered in the models. The simulation results were validated by our experiment. Three algorithms for representing an ellipsoidal particle were compared in macro and micro aspects. The results show that there exists big difference in the microscopic parameters such as kinetic energy, rotational kinetic energy, deformation, contact force and collision number which leads to the difference of macroscopic parameters. The relative error in the discharge rate and tracer particle position is the largest between 3-tangent-element representation and experimental results. The flow pattern is similar for the 5-element and 3-intersection representations. The only difference is the discharge rate of 5-element representation is larger than the experimental value and that of the 3-intersection representation has the contrary result. Finally the 3-intersection- element reoresentation is chosen in the simulation due to less comouting time than that of the 5-element renresentation.

Multi-level coarse-grained discrete element method modeling of cylinder particle flow in a rotating drum

The coarse-grained discrete element method(DEM)is probably a feasible option for simulating an actual drum-type biomass boiler,which contains over 10 million cylinder particles.A multi-level study was conducted based on particle and coarse-grained level data to evaluate the adequacy of the coarse-grained approach in terms of geometrical characteristics,kinematic features,and dynamic properties.Two scaling laws for contact parameters were used and compared during the simulations.The results show that the coarse-grained approach can accurately predict the positions of the free surface and active-passive interface,the mixing index,and the orientation properties.Deviations in the velocity fields may occur due to the worse flowability of coarse-grained particles near the free surface.The efficiency is significantly improved by the coarse-grained model compared with the corresponding original case(the same DEM code without a coarse-grained model was used for the original simulations).

Particle size spatial distribution in landslide dams

The particle composition and spatial distribution of landslide-induced dam bodies are critical geotechnical parameters for studying the hazards of dam-break floods.However,current research often neglects the influence of the initial particle composition and spatial distribution of the landslide on the particle composition and spatial distribution of the landslide dam.This study investigated the impact of initial particle size distribution,volume,and sliding length on the energy and velocity changes of characteristic particles during the sliding process and the spatial distribution of particle sizes in the landslide dam body.Numerical simulations and physical models were employed to examine the effects of sequential gradient arrangements(where particle sizes decrease from top to bottom)and four other different initial particle arrangements on the energy and velocity changes of particles and the spatial distribution of particle sizes in the dam body.The study reveals the characteristics of translational and rotational energy of different particles and the laws of mechanical energy conversion,obtaining the spatial distribution patterns of particle sizes in landslide-induced dams.The results show that under the sequential gradient arrangement,the energy dissipation of the landslide movement is lower,with larger particles mainly distributed at the distal end and smaller particles at the proximal end of the landslide dam.In contrast,under the reverse gradient arrangement,the energy dissipation of the landslide movement is higher,and the distribution pattern of the dam particles is opposite to that of the sequential gradient arrangement.For the other arrangement modes,the spatial distribution of dam particles falls between the aforementioned two.There is a positive correlation between particle size and translational kinetic energy within the particle flow during the landslide process,and rotational motion increases energy dissipation.Under constant slope conditions,sliding length does not affect the movement pattern of the particle flow or the spatial distribution of particles in the dam body.The findings of this study provide a scientific basis for the accurate simulation and prediction of dam-break flood processes.

Oblique impact breakage unification of nonspherical particles using discrete element method

Particle breakage commonly occurs during processing of particulate materials,but a mechanistic model of particle impact breakage is not fully established.This article presents oblique impact breakage characteristics of nonspherical particles using discrete element method(DEM)simulations.Three different particle shapes,i.e.spherical,cuboidal and cylindrical,are investigated.Constituent spheres are agglomerated with bridging bonds to model the breakage characteristics under impact conditions.The effect of agglomerate shapes on the breakage pattern,damage ratio,and fragment size distribution is fully investigated.By using a newly proposed oblique impact model,unified breakage master surfaces are theoretically constructed for all the particle shapes under oblique impact conditions.The developed approach can be applied to modelling particulate processes where nonspherical particles and oblique impact breakage are prevailing.

Discrete Element Modelling of Dynamic Behaviour of Rockfills forResisting High Speed Projectile Penetration

This paper presents a convex polyhedral based discrete element method for modelling the dynamic behaviour ofrockfills for resisting high speed projectile penetration. The contact between two convex polyhedra is defined by theMinkowski overlap and determined by the GJK and EPA algorithm. The contact force is calculated by a Minkowskioverlap based normal model. The rotational motion of polyhedral particles is solved by employing a quaternionbased orientation representation scheme. The energy-conserving nature of the polyhedral DEM method ensures arobust and effective modelling of convex particle systems. The method is applied to simulate the dynamic behaviourof a rockfill system under impact of a high speed projectile. The rockfill sample is generated by a three-dimensionalVoronoi meso method with a specific particle size distribution. The penetrating process of the projectile strikingthe rockfill target is simulated. Some physical quantities associated with the projectile such as the residual velocity,penetration resistance, and deflection angle are monitored which can reflect the influence of the characteristics ofthe rockfill target on its anti-penetration performance. It can be concluded that the developed polyhedral DEMmethod is a very promising numerical approach in analysing the dynamic behaviour of rockfill systems subject tohigh speed projectile impact.