Using the DEM-CFD method to predict Brownian particle deposition in a constricted tube

Modelling of the agglomeration and deposition on a constricted tube collector of colloidal size particles immersed in a liquid is investigated using the discrete element method (DEM). The ability of this method to represent surface interactions allows the simulation of agglomeration and deposition at the particle scale. The numerical model adopts a mechanistic approach to represent the forces involved in colloidal suspensions by including near-wall drag retardation, surface interaction and Brownian forces. The model is implemented using the commercially available DEM package EDEM 2.3^®, so that results can be replicated in a standard and user-friendly framework. The effects of various particle-to-collector size ratios, inlet fluid flow-rates and particle concentrations are examined and it is found that deposition efficiency is strongly dependent on the inter-relation of these parameters. Particle deposition and re-suspension mechanisms have been identified and analyzed thanks to EDEM's post processing capability. One-way coupling with computational fluid dynamics (CFD) is considered and results are compared with a two-way coupling between EDEM 2.3^® and FLUENT 12.1^®. It is found that two-way coupling requires circa 500% more time than one-way coupling for similar results.     

Mixing and segregation assessment of bi-disperse solid particles in a double paddle mixer

A double paddle blender's flow patterns and mixing mechanisms were analyzed using discrete element method (DEM) and experiments. The mixing performance of this type of the blender containing bi-disperse particles has been rarely studied in the literature. Plackett-Burman design of experiments (DoE) methodology was used to calibrate the DEM input parameters. Subsequently, the impact of the particle number ratio, vessel fill level, and paddle rotational speed on mixing performance was investigated using the calibrated DEM model. The mixing performance was assessed using relative standard deviation and segregation intensity. Mixing performance was significantly affected by the paddle rotational speed and particle number ratio. Moreover, the Peclet number and diffusivity coefficient were used to evaluate the mixing mechanism in the blender. Results revealed that the diffusion was the predominant mixing mechanism, and the best mixing performance was observed when the diffusivity coefficients of 3 mm and 5 mm particles were almost equal.     

Numerical simulation of tetrahedral particle mixing and motion in rotating drums

A regular tetrahedron is the simplest three-dimensional structure and has the largest non-sphericity. Mixing of tetrahedral particles in a thin drum mixer was studied by the soft-sphere-imbedded pseudo-hard particle model and compared with that of spherical particles. The two particle types were simulated with different rotation speeds and drum filling levels. The Lacey mixing index and Shannon information entropy were used to explore the effects of sphericity on the mixing and motion of particles. Moreover, the probability density functions and mean values and variances of motion velocities, including translational and rotational, were computed to quantify the differences between the motion features of tetrahedra and spheres. We found that the flow regime depended on the particle shape in addition to the rotation speed and filling level of the drum. The mixing of tetrahedral particles was better than that of spherical particles in the rolling and cascading regimes at a high filling level, whereas it may be poorer when the filling level was low. The Shannon information entropy is better than the Lacey mixing index to evaluate mixing because it can reflect the real change of flow regime from the cataracting to the centrifugal regime, whereas the mixing index cannot.     

Using the discrete element method to assess the mixing of polydisperse solid particles in a rotary drum

Despite the wide applications of powder and solid mixing in industry, knowledge on the mixing of polydisperse solid particles in rotary drum blenders is lacking. This study investigates the mixing of monodisperse, bidisperse, tridisperse, and polydisperse solid particles in a rotary drum using the discrete element method. To validate the model developed in this study, experimental and simulation results were compared. The validated model was then employed to investigate the effects of the drum rotational speed, particle size, and initial loading method on the mixing quality. The degree of mixing of polydisperse particles was smaller than that for monodisperse particles owing to the segregation phenomenon. The mixing index increased from an initial value to a maximum and decreased slightly before reaching a plateau for bidisperse, tridisperse, and polydisperse particles as a direct result of the segregation of particles of different sizes. Final mixing indices were higher for polydisperse particles than for tridisperse and bidisperse particles. Additionally, segregation was weakened by introducing additional particles of intermediate size. The best mixing of bidisperse and tridisperse particles was achieved for top–bottom smaller-to-larger initial loading, while that of polydisperse systems was achieved using top–bottom smaller-to-larger and top–bottom larger-to-smaller initial loading methods.     

Discrete element method study of shear-driven granular segregation in a slowly rotating horizontal drum

Segregation and mixing of granular materials are complex processes and are not fully understood. Motivated by industrial need, we performed a simulation using the discrete element method to study size segregation of a binary mixture of granular particles in a horizontal rotating drum. Particles of two different sizes were poured into the drum until it was 50% full. Shear-driven segregation was induced by rotating the side-plates of the drum in the opposite direction to that of the cylindrical wall. We found that radial segregation diminished in these systems but did not completely vanish. In an ordinary rotating drum, a radial core of smaller particles is formed in the center of the drum, surrounded by larger revolving particles. In our system, however, the smaller particles were found to migrate toward the side-plates. The shear from anti-spinning side-plates reduces the voidage and increases the bulk density. As such, smaller particles in the mixer tend to move to denser regions. We varied the shear by changing the coefficient of friction on the side-plates to study the influence of shear rate on this migration. We also compared the extent of radial segregation with stationary side-plates and with side-plates moving in different angular directions.     

DEM simulation of polydisperse systems of particles in a fluidized bed

Numerical simulations based on three-dimensional discrete element model (DEM) are conducted for mono-disperse, binary and ternary systems of particles in a fluidized bed. Fluid drag force acting on each particle depending on its size and relative velocity is assigned. The drag coefficient corresponding to Ergun’s correlation is applied to the system of fluidized bed with particle size ratios of 1:1 for the mono-disperse system, 1:1.2, 1:1.4 and 1:2 for the binary system and 1:1.33:2 for the ternary system b...     

Euler-Lagrange/DEM simulation of wood gasification in a bubbling fluidized bed reactor

We present an Euler–Lagrange method for the simulation of wood gasification in a bubbling fluidized bed. The gas phase is modeled as a continuum using the 2D Navier–Stokes equations and the solid phase is modeled by a Discrete Element Method(DEM)using a soft-sphere approach for the particle collision dynamic. Turbulence is included via a Large-Eddy approach using the Smagorinsky sub-grid model.The model takes into account detailed gas phase chemistry,zero-dimensional modeling of the pyrolysis and gasification of each individual particle,particle shrinkage,and heat and mass transfer between the gas phase and the particulate phase.We investigate the influence of wood feeding rate and compare exhaust gas compositions and temperature results obtained with the model against experimental data of a laboratory scale bubbling fluidized bed reactor.     

Eulerian-Lagrangian simulation of distinct clustering phenomena and RTDs in riser and downer

Numerical simulation of fully developed hydrodynamics of a riser and a downer was carried out using an Eulerian-Lagrangian model, where the particles are modeled by the discrete element method (DEM) and the gas by the Navier-Stokes equations. Periodic flow domain with two side walls was adopted to simulate the fully developed dynamics in a 2D channel of 10 cm in width. All the simulations were carried out under the same superficial gas velocity and solids holdup in the domain, starting with a homogenous sta...     

Particle convective heat transfer near the wall in a supercritical water fluidized bed by single particle model coupled with CFD-DEM

Supercritical water fluidized bed (SCWFB) is a promising reactor to gasify biomass or coal. Its optimization design is closely related to wall-to-bed heat transfer, where particle convective heat transfer plays an important role. This paper evaluates the particle convective heat transfer coefficient (h_pc) at the wall in SCWFB using the single particle model. The critical parameters in the single particle model which is difficult to get experimentally are obtained by the computational fluid dynamics-discrete element method (CFD-DEM). The contact statistics related to particle-to-wall heat transfer, such as contact number and contact distance, are also presented. The results show that particle residence time (τ), as the key parameter to evaluate h_pc, is found to decrease with rising velocity, while increase with larger thermal boundary layer thickness. τ follows a gamma function initially adopted in the gas–solid fluidized bed, making it possible to evaluate h_pc in SCWFB by a simplified single particle model. The theoretical predicted h_pc tends to increase with rising thermal gradient thickness at a lower velocity (1.5 U_mf), while first decreases and then increases at higher velocity (1.75 and 2 U_mf). h_pc occupies 30%–57% of the overall wall-to-bed heat transfer coefficient for a particle diameter of 0.25 mm. The results are helpful to predict the overall wall-to-bed heat transfer coefficient in SCWFB combined with a reasonable fluid convective heat transfer model from a theoretical perspective.     

Influences of the fluidizing and spouting pulsation on particle motion in spout-fluid beds

Effects of variable airflow on particle motion in spout-fluid beds are studied. Computational fluid dynamics using Navier–Stokes equations for the gas phase coupled with the discrete element method using Newton’s laws for the solid phase have been employed. Results indicate that increasing the fluidizing velocity diminishes dead zones and increases both the total height of the bed and the traversed distance by particles in the steady spout-fluid bed. In pulsed airflows, two configurations are investigated, namely, the spouted pulsed-fluidized bed with pulsed flow of the fluidizing velocity, and the pulsed-spouted fluidized bed with pulsed flow of the spouting velocity. The positive effect of pulsation on particle motion is shown and the effects of parameters, such as amplitude and frequency, on the dynamics of the bed are investigated in each configuration. An increase of up to 19% in traversed distance is found for the range studied, which suggests flow pulsation as a promising technique for increasing particle mixing in spout-fluid beds.     

Euler–Lagrange/DEM simulation of wood gasification in a bubbling fluidized bed reactor

We present an Euler–Lagrange method for the simulation of wood gasification in a bubbling fluidized bed. The gas phase is modeled as a continuum using the 2D Navier–Stokes equations and the solid phase is modeled by a Discrete Element Method (DEM) using a soft-sphere approach for the particle collision dynamic. Turbulence is included via a Large-Eddy approach using the Smagorinsky sub-grid model. The model takes into account detailed gas phase chemistry, zero-dimensional modeling of the pyrolysis and gasification of each individual particle, particle shrinkage, and heat and mass transfer between the gas phase and the particulate phase. We investigate the influence of wood feeding rate and compare exhaust gas compositions and temperature results obtained with the model against experimental data of a laboratory scale bubbling fluidized bed reactor.