S. Wirtz

Selection of an appropriate time integration scheme for the discrete element method (DEM)

Abstract With increasing computer power simulation methods addressing discrete problems in a broad range of scientific fields become more and more available. The discrete element method is one of these discontinuous approaches used for modeling granular assemblies. Within this method the dynamics of a system of particles is modeled by tracking the motion of individual particles and their interaction with their adjacencies over time. For the interaction of particles, force models need to be specified. The resulting equations of motion are of coupled ordinary differential configuration, which are usually solved by explicit numerical schemes. In large-scale systems like avalanches, planetary rings, hoppers or chemical reactors vast numbers of particles need to be addressed. Therefore, integration schemes need to be accurate on the one hand, but also numerically efficient on the other hand. This numerical efficiency is characterized by the methods demand for memory and CPU-time. In this paper a number of mostly explicit numerical integration schemes are reviewed and applied to the benchmark problem of a particle impacting a fixed wall as investigated experimentally by Gorham and Kharaz [Gorham, D. A., & Kharaz, A. H. (2000). The measurement of particle rebound characteristics. Powder Technology , 112 (3), 193–202]. The accurate modeling which includes the correct integration of the equations of motion is essential. In discrete simulation methods the accuracy of properties on the single particle level directly influence the global properties of the granular assembly like velocity distributions, porosities or flow rates, whereas their correct knowledge is often of key interest in engineering applications. The impact experiment is modeled with simple force displacement approaches which allow an analytical solution of the problem. Aspects discussed are the dependency of the step size on the accuracy of certain collision properties and the related computing time. The effect of a fixed time step is analyzed. Guidelines for the efficient selection of an integration scheme considering the additional computational cost by contact detection and force calculation are presented.

Modeling of Granular Flow and Combined Heat Transfer in Hoppers by the Discrete Element Method (DEM)

The discrete element method can be used for modeling moving granular media in which heat and mass transport takes place. In this paper the concept of discrete element modeling with special emphasis on applicable force laws is introduced and the necessary equations for heat transport within particle assemblies are derived. Possible flow regimes in moving granular media are discussed. The developed discrete element model is applied to a new staged reforming process for biomass and waste utilization which employs a solid heat carrier. Results are presented for the flow regime and heat transport in substantial vessels of the process.

A Comparative Numerical Study of Particle Mixing on Different Grate Designs Through the Discrete Element Method

Based on LEAT’s discrete element codes, granular flow and mixing on conveying equipment are studied in two and three dimensions. Discrete element simulations, which are briefly introduced, provide detailed information on particle positions and velocities over time. This information is used to derive quantities characterizing the dynamic process of mixing. The main focus of the study presented is the mixing process of inhomogeneous particle ensembles on different grate types. For this purpose, the introduced mixing parameters are used to compare the mixing in a 3D situation with the corresponding 2D approximation on identical grates and to compare different grate designs in two dimensions.

An Experimental and Numerical Study of Transversal Dispersion of Granular Material on a Vibrating Conveyor

The mixing of thin granular layers transported on the surface of an oscillating trough is experimentally and numerically examined. The particle dispersion was experimentally quantified by an image processing system recording the growth of the mixing layer thickness of two differently colored but otherwise identical sand particle streams along the longitudinal position within the transporting channel. Granular flow and dispersion on the vibrating conveyor were studied numerically based on a three-dimensional discrete element code. Both experiments and simulations were used to derive quantities characterizing the transversal dispersion. The mixing was found to be directly proportional to the vertical acceleration of the conveyor and inversely proportional to the mass flow of the transported material. Keeping the above-mentioned parameters constant, the dispersion increases with increasing mean particle diameter. When performing the experiments with materials of different mean particle diameters and tuning the mass flow to achieve the same level of dimensionless bed height, the magnitude of the dispersion coefficient remains constant, as was also confirmed by the numerical simulation.

Applicable Contact Force Models for the Discrete Element Method: The Single Particle Perspective

Several processes in nature as well as many industrial applications involve static or dynamic granular materials. Granulates can adopt solid-, liquid-, or gaslike states and thereby reveal intriguing physical phenomena not observable in its versatility for any other form of matter. The frequent occurrence of phase transitions and the related characteristics thereby strongly affect their processing quality and economics. This situation demands for prediction methods for the behavior of granulates. In this context simulations provide a feasible alternative to experimental investigations. Several different simulation approaches are applicable to granular materials. The time-driven discrete element method turns out to be not only the most complex but also the most general simulation approach. Discrete element simulations have been used in a wide variety of scientific fields for more than 30 years. With the tremendous increase in available computer power, especially in the past years, the method is more and more developing to the state of the art simulation technique for granular materials not only in science but also in industrial applications. Several commercial software packages utilizing the time-driven discrete element method have emerged and are becoming more and more popular within the engineering community. Despite the long time of usage of the time-driven discrete element method, model advances derived and theoretical and experimental studies performed in the different branches of application lack harmonization. They thereby provide potential for improvements. Therefore, the scope of this paper is a review of methods and models for contact forces based on theoretical considerations and experimental data from literature. Particles considered are of spherical shape. Through model advances it is intended to contribute to a general enhancement of simulation techniques, which help improve products and the design of the related equipment.

Convective Drying of Agitated Silica Gel and Beech Wood Particle Beds—Experiments and Transient DEM-CFD Simulations

With coupled discrete element (DEM)–computational fluid dynamics (CFD) simulations, drying processes can be simultaneously described on the system scale while resolving detailed subprocesses on the particle scale. In this contribution, DEM-CFD simulations are used to analyze the transient heat and mass transfer in mechanically agitated particle beds during drying. Results are compared to convective batch-drying experiments with silica gel and beech wood spheres and mixing effects are studied in detail. A good agreement with the measurements of both single-particle and particle bed drying is achieved by resolving heat and moisture transport three-dimensionally inside each particle.

Fast multitarget tracking via strategy switching for sensor-based sorting

State-of-the-art sensor-based sorting systems provide solutions to sort various products according to quality aspects. Such systems face the challenge of an existing delay between perception and separation of the material. To reliably predict an objects position when reaching the separation stage, information regarding its movement needs to be derived. Multitarget tracking offers approaches through which this can be achieved. However, processing time is typically limited since the sorting decision for each object needs to be derived sufficiently early before it reaches the separation stage. In this paper, an approach for multitarget tracking in sensor-based sorting is proposed which supports establishing an upper bound regarding processing time required for solving the measurement to track association problem. To demonstrate the success of the proposed method, experiments are conducted for data-sets obtained via simulation of a sorting system. This way, it is possible to not only demonstrate the impact on required runtime but also on the quality of the association.

Simulation-based evaluation of predictive tracking for sorting bulk materials

Multitarget tracking problems arise in many real-world applications. The performance of the utilized algorithm strongly depends both on how the data association problem is handled and on the suitability of the motion models employed. Especially the motion models can be hard to validate. Previously, we have proposed to use multitarget tracking to improve optical belt sorters. In this paper, we evaluate both the suitability of our model and the tracking and then of our entire system incorporating the image processing component via the use of highly realistic numerical simulations. We first assess the model using noise-free measurements generated by the simulation and then evaluate the entire system by using synthetically generated image data.

HEAT TRANSFER EXPERIMENTS IN A ROTARY DRUM FOR A VARIETY OF GRANULAR MATERIALS

The influence of material properties on the contact heat transfer coefficient between the covered wall surface and the solid bed was investigated. The contact heat transfer coefficients were calculated from the measured radial and circumferential temperature profiles. Experiments were carried out with six different materials, including steel spheres, animal powder, cement clinker, quartz sand, glass beads, and expanded clay. The rotational speeds were varied from 1 to 6 rpm to evaluate the influence of rotational speed on the contact heat transfer coefficient. The measured contact heat transfer coefficients were compared with four models from the literature.

Simulation of Heat Transfer in Moving Granular Material by the Discrete Element Method With Special Emphasis on Inner Particle Heat Transfer

Physical processes involving static or dynamic granular assemblies are best modeled on the particle scale by Discrete Element Methods (DEM) rather than continuum approaches. Due to the high computational effort of DEM simulations, present studies assume the inner particle temperature to be spatially uniform and neglect the inner particle heat transfer. For this reason the Radial Temperature Model was introduced [1, 2] It assumes radial temperature distributions within the particles and is based on an analytical solution of the heat conduction equation in a spherical particle. The scope of this paper is to present the further development of the Radial Temperature Model that allows to simulate granular systems of particles of different sizes and materials, enabling the use of DEM in various applications. The contact heat transfer is modeled making additional material-specific data unnecessary. It is shown that a very good accuracy for the contact heat transfer between different spherical particles is achieved for binary contacts. DEM simulations were performed using the Radial Temperature Model and uniform particle temperatures, respectively. The results demonstrate that the Radial Temperature Model that has been developed and incorporated in the Discrete Element Method allows for an improved calculation of the transient thermal behavior of granular assemblies even for large numbers of particles.Copyright