Determination of contact parameters for discrete element method simulations of granular systems

Both linear-spring-dashpot (LSD) and non-linear Hertzian-spring-dnshpot (HSD) contact models are commonly used for the calculation of contact forces in Discrete Element Method (DEM) simulations of granular systems.Despite the popularity of these models, determination of suitable values for the contact parameters of the simulated particles such as stiffness, damping coefficient, coefficient of restitution, and simulation time step,is not altogether obvious.In this work the relationships between these contact parameters for a model system where a particle impacts on a flat base are examined.Recommendations are made concerning the determination of these contact parameters for use in DEM simulations.     

基于OpenGL的DEM地形可视性分析算法研究

针对规则网格模型的DEM可视性问题,提出了一种基于仰角坡度比较法的LOS方法,以LOS视线为基准,根据DEM地形实际的高程和坡度值来判断其是否可视,并给出了可视化区域的确定方式.利用OpenGL图形库提供的网格化方式创建可视区域,给出了VC++开发环境下的可视化区域绘制程序.     

Modelling of massive particulates for breakwater engineering using coupled FEMDEM and CFD

The seaward slope of many breakwaters consists of thousands of interlocking units of rock or concrete comprising a massive granular system of large elements each weighing tens of tonnes. The dumped quarry materials in the core are protected by progressively coarser particulates. The outer armour layer of freely placed units is intended to both dissipate wave energy and remain structurally stable as strong flows are drawn in and out of the particulate core. Design guidance on the mass and shape of these units is based on empirical equations derived from scaled physical model tests. The main failure mode for armour layers exposed to severe storms is hydraulic instability where the armour units of concrete or rock are subjected to uplift and drag forces which can in turn lead to rocking, displacement and collisions sufficient to cause breakage of units. Recently invented armour unit designs making up such granular layered system owe much of their success to the desirable emergent properties of interlock and porosity and how these combine with individual unit structural strength and inertial mass. Fundamental understanding of the forces governing such wave-structure interaction remains poor. We use discrete element and combined finite-discrete element methods to model the granular solid skeleton of randomly packed units coupled to a CFD code which resolves the wave dynamics through an interface tracking technique. The CFD code exploits several methods including a compressive advection scheme, node movement, and general mesh optimization. We provide the engineering context and report progress towards the numerical modelling of instability in these massive granular systems.     

DEM simulation of particle mixing in a sheared granular flow

Mixing behaviors of particles are simulated in a sheared granular flow using differently colored but otherwise identical glass spheres, with five different bottom wall velocities. By DEM simulation, the solid fractions, velocities, velocity fluctuations and granular temperatures are measured. The mixing layer thicknesses are compared with the calculations from a simple diffusion equation using the data of apparent self-diffusion coefficients obtained from the current simulation measurements. The calculations and simulation results showed good agreements, demonstrating that the mixing process of granular materials occurred through the diffusion mechanism.     

DEM simulations of die filling during pharmaceutical tabletting

The flow behaviour of powders from a stationary shoe into a moving die, which mimics the die filling process in a rotary tablet press, was analysed using a discrete element method （DEM）, in which 2D irregular shaped particles were considered. The influence of the particle shape, size and size distribution, the number of particles used in the simulation, the initial height of powder bed in the shoe, and the filling speed on the average mass flow rate and the critical filling speed （the highest speed at which the die can be completely filled） were explored. It has been found that a maximum flow rate is obtained at the critical filling speed for all systems investigated and poly-disperse systems have higher mass flow rates and higher critical filling speeds than mono-disperse systems. In addition, the powder with particles which can tessellate generally has a lower filling rate and a lower critical titling speed.     

Modelling of massive particulates for breakwater engineering using coupled FEMDEM and CFD

The seaward slope of many breakwaters consists of thousands of interlocking units of rock or concrete comprising a massive granular system of large elements each weighing tens of tonnes.The dumped quarry materials in the core are protected by progressively coarser particulates.The outer armour layer of freely placed units is intended to both dissipate wave energy and remain structurally stable as strong flows are drawn in and out of the particulate core.Design guidance on the mass and shape of these units is based on empirical equations derived from sealed physical model tests.The main failure mode for armour layers exposed to severe storms is hydraulic instability where the armour units of concrete or rock are subjected to uplift and drag forces which can in turn lead to rocking,displacement and collisions sufficient to cause breakage of units.Recently invented armour unit designs making up such granular layered system owe much of their success to the desirable emergent properties of interlock and porosity and how these combine with individual unit structural strength and inertial mass. Fundamental understanding of the forces governing such wave-structure interaction remains poor.We use discrete element and combined finite-discrete element methods to model the granular solid skeleton of randomly packed units coupled to a CFD code which resolves the wave dynamics through an interface tracking technique.The CFD code exploits several methods including a compressive advection scheme, node movement, and general mesh optimization.We provide the engineering context and report progress towards the numerical modelling of instability in these massive granular systems.     

New discrete element models for elastoplastic problems

The discrete element method （DEM） has attractive features for problems with severe damages, but lack of theoretical basis for continua behavior especially for nonlinear behavior has seriously restricted its application. The present study proposes a new approach to developing the DEM as a general and robust technique for modeling the elastoplastic behavior of solid materials. New types of connective links between elements are proposed, the interelement parameters are theoretically determined based on the principle of energy equivalence and a yield criterion and a flow rule for DEM are given for describing nonlinear behavior of materials. Moreover, a numerical scheme, which can be applied to modeling the behavior of a continuum as well as the transformation from a continuum to a discontinuum, is obtained by introducing a fracture criterion and a contact model into the DEM. The elastoplastic stress wave propagations and the tensile failure process of a steel plate are simulated, and the numerical results agree well with those obtained from the finite element method （FEM） and corresponding experiment, and thus the accuracy and efficiency of the DEM scheme are demonstrated.     

Experimental and DEM analyses on wheel-soil interaction

In this paper, the wheel-soil interaction for a future lunar exploration mission is investigated by physical model tests and numerical simulations. Firstly, a series of physical model tests was conducted using the TJ-1 lunar soil simulant with various driving conditions, wheel configurations and ground void ratios. Then the corresponding numerical simulations were performed in a terrestrial environment using the Distinct Element Method (DEM) with a new contact model for lunar soil, where the rolling resistance and van der Waals force were implemented. In addition, DEM simulations in an extraterrestrial (lunar) environment were performed. The results indicate that tractive efficiency does not depend on wheel rotational velocity, but decreases with increasing extra vertical load on the wheel and ground void ratio. Rover performance improves when wheels are equipped with lugs. The DEM simulations in terrestrial environment can qualitatively reproduce the soil deformation pattern as observed in the physical model tests. The variations of traction efficiency against the driving condition, wheel configuration and ground void ratio attained in the DEM simulations match the experimental observations qualitatively. Moreover, the wheel track is found to be less evident and the tractive efficiency is higher in the extraterrestrial environment compared to the performance on Earth.     

No more conventional reference electrode: Transition time for determining chloride ion concentration

Ion selective electrodes (ISE) are used extensively for the potentiometric determination of ion concentrations in electrolytes. However, the inherent drift in these measurements and the requirement of a stable reference electrode restrict the feasibility of this method for long-term in-situ applications. This work presents a chronopotentiometric approach to minimize drift and avoid the use of a conventional reference electrode for measuring chloride ion concentration. An anodic current pulse is applied to a Ag/AgCl working electrode which initiates a faradaic reaction that depletes the chloride ions near the electrode surface. The rate of change in potential at the Ag/AgCl electrode, due to chloride ion depletion, reaches an inflection point once the chloride ions deplete completely near the electrode surface. The moment of the inflection point, also known as the transition time, is a function of the chloride ion concentration and is described by the Sand equation. It is shown that the square root of the transition time is linearly proportional to the chloride ion concentration. Drift in the response over two weeks is negligible: 59 μM/day when measuring 1 mM of Cl⁻ ions using a 10 A m⁻² current pulse. The transition time at a specific ion concentration can be tuned by the applied current pulse, e.g., in a solution containing 5 mM chloride ions, the transition times with current pulses of 10 and 20 A m⁻² are 1.56 and 0.25 s, respectively. The moment of inflection determines the response, and thus is independent of the absolute potential of reference electrode. Therefore, any metal wire can act as a pseudo-reference electrode, enabling this approach for long-term and integrated-sensor applications such as measurement inside concrete structures.     

Development and validation of SuperDEM for non-spherical particulate systems using a superquadric particle method

This article presents the development and validation of the Superquadric Discrete Element Method (SuperDEM) for non-spherical particle simulation using a superquadric particle method in open-source CFD suite MFiX. A superquadric particle–particle contact algorithm with accelerating and stabilizing strategy was developed. A superquadric particle–arbitrary wall contact algorithm was developed, which enables the simulation in complex geometry. The solver was validated by comparing with experimental data generated in this study or available in the literature. Tests include cylinder contacting with a wall, static packing of M&M chocolate candies in a cylindrical container, static packing of cylinders in a cylindrical container, dynamic angle of repose of cylinders in a rotating drum, and discharging of chocolate candies from a hopper. Besides, MPI parallelization of the solver was implemented and the parallel performance of the solver using MPI was assessed through large-scale simulations of 1 million, 10 million, and 100 million particles on up to 6800 cores, which demonstrates that the SuperDEM solver has great potential for industrial-scale systems simulation.