Discrete particle simulation of particulate systems: A review of major applications and findings

Modelling the Multiphase Flow in Dense Medium Cyclones

Dense medium cyclone (DMC) is widely used in mineral industry to separate solids by density. It is simple in design but the flow pattern within it is complex due to the size and density distributions of the feed and process medium solids, and the turbulent vortex formed. Recently, the so-called combined computational fluid dynamics (CFD) and discrete element method (DEM) (CFD-DEM) was extended from two-phase flow to model the flow in DMCs at the University of New South Wales (UNSW). In the CFD-DEM model, the flow of coal particles is modelled by DEM and that of medium flow by CFD, allowing consideration of medium-coal mutual interaction and particle-particle collisions. In the DEM model, Newton's laws of motion are applied to individual particles, and in the CFD model the local-averaged Navier-Stokes equations combined with the volume of fluid (VOF) and mixture multiphase flow models are solved. The application to the DMC studies requires intensive computational effort. Therefore, various simplified versions have been proposed, corresponding to the approaches such as Lagrangian particle tracking (LPT) method where dilute phase flow is assumed so that the interaction between particles can be ignored, one-way coupling where the effect of particle flow on fluid flow is ignored, and the use of the concept of parcel particles whose properties are empirically determined. In this paper, the previous works on the modelling of DMCs at UNSW are summarized and the features and applicability of the models used are discussed.

An integrated Quality by Design (QbD) approach towards design space definition of a blending unit operation by Discrete Element Method (DEM) simulation

A combined Quality by Design (QbD) and Discrete Element Model (DEM) simulation-approach is presented to characterize a blending unit operation by evaluating the impact of formulation parameters and process variables on the blending quality and blending end point. Understanding the variability of both the API and the excipients, as well as their impact on the blending process are critical elements for blending QbD. In a first step, the QbD-methodology is systematically used to (1) establish the critical quality attribute content uniformity and to link this CQA to its surrogate blend homogeneity, (2) identify potentially critical input factors that may affect blending operation quality and (3) risk-rank these factors to define activities for process characterization. Subsequently, a DEM-simulation-based characterization of the blending process is performed. A statistical evaluation is finally presented, relating blend homogeneity of systems with low particle number to the regulatory requirements. Data are then used to map out a three-dimensional knowledge space, providing parameters to define a design space and set up an appropriate control strategy.

Microstructure evolution in density relaxation by tapping

The density relaxation phenomenon is modeled using both Monte Carlo and discrete element simulations to investigate the effects of regular taps applied to a vessel having a planar floor filled with monodisperse spheres. Results suggest the existence of a critical tap intensity which produces a maximum bulk solids fraction. We find that the mechanism responsible for the relaxation phenomenon is an evolving ordered packing structure propagating upwards from the plane floor.

Numerical study of the effects of particle shape and polydispersity on permeability

We study through numerical simulations the dependence of the hydraulic permeability of granular materials on the particle shape and the grain size distribution. Several models of sand are constructed by simulating the settling under gravity of the grains; the friction coefficient is varied to construct packs of different porosity. The size distribution and shapes of the grains mimic real sands. Fluid flow is simulated in the resulting packs using a finite element method and the permeability of the packs is successfully compared with available experimental data. Packs of nonspherical particles are less permeable than sphere packs of the same porosity. Our results indicate that the details of grain shape and size distribution have only a small effect on the permeabilty of the systems studied.

Three-dimensional simulations of a vertically vibrated granular bed including interstitial air

We present a numerical study of the effect of interstitial air on a vertically vibrated granular bed within one period of oscillation. We use a three-dimensional molecular-dynamics simulation including air phenomenologically. The simulations are validated with experiments made with spherical glass beads in a rectangular container. After validation, results are reported for a granular column of 9000 grains and approximately 50 layers deep (at rest), agitated with a sinusoidal excitation with maximal acceleration 4.7g at 11.7 Hz. We report the evolution of density, granular temperature, and coordination number within a vibration cycle, and the effect of interstitial air on those parameters. In three-dimensional computer simulations we found that the presence of interstitial air can promote the collective motion of the granular material as a whole.

Vibrated granular bed on a bumpy surface

We investigate effects of physical characteristics of the vibrating base in a vibrated granular bed system on the bulk granular behavior of the bed using molecular dynamics simulations. With a vibrating base that exhibits low coefficient of restitution for particle-base collisions, a monolayer of granular materials was observed to form dynamically on the surface of the base. Such a system was shown to be both qualitatively and quantitatively equivalent to one with a bumpy base composed of discrete solid particles undergoing the same type of imposed oscillatory motion. The "bumpiness" of such a base, as defined by the size of particles constituting the vibrating base, was also found to have insignificant effects on the behavior of the bulk granular materials. The observations made in this study may point toward a possible methodology to model vibrating granular bed systems with inelastic bases using continuum theories.