DEM study of the effect of impeller design on mixing performance in a U-shape ribbon mixer

The mixing of powders in a U-shape mixer is significantly influenced by the mixer design, especially impellers, but the studies on the mixing processes are still insufficient. In this study, the effect of impeller designs on mixing performance in an industrial-scale U-shaped ribbon mixer is studied using DEM simulations. Three impeller designs are studied: 2-bladed impeller spiralling in the same direction (i.e., Design I) and the opposite direction (i.e., Design II), and 4-bladed impeller (i.e., Design III). Different particle mixing behaviours in three different impeller designs are studied in aspects of mixing status, particle path line, velocity distribution, and forces. The radial direction has the highest dispersion coefficient while the axial direction has the lowest dispersion coefficient. Most particles in the mixers are imposed a weak force. Design III shows the best mixing performance among the three with the front-by-back and top-by-bottom loading used. Design II shows a better mixing performance used than Design I and III with the side-by-side loading but takes a longer time to reach the stable status. This work evaluates the effect of different impeller designs on the mixing performance in an industrial-scale U-shaped ribbon mixer and provides an effective way to assist industrial design in an economical and safe manner.     

Discrete element simulation for mixing performances and power consumption in a twin-blade planetary mixer with non-cohesive particles

The present study aims to characterize the mixing performances and power consumption of a twin-blade planetary mixer with non-cohesive particles through the discrete element method (DEM). A DEM model used for simulating the particle flow and mixing kinetics of the mixer was experimentally verified. The particle velocity and mixing mechanism are elaborated quantitatively, indicating that particle mixing is realized under the combined actions of radial, circumferential and vertical circulations, and some local collisions and mergers. Increasing the absolute speed N and the speed ratio i promotes the radial circulation, while the tangential and vertical circulations are strengthened with the increase of N and the decrease of i. The mixing time required for the homogeneous state decreases, and the power consumption increases as N increases and i decreases. Thus, increasing N and decreasing i can improve the mixing performance but require more energy to reach the homogeneous state. Also, the mixing performance shows a strong correlation with the swept volume of blades, which proves that the dominant mixing mechanism of the mixer is convection.     

Mixing assessment of non-cohesive particles in a paddle mixer through experiments and discrete element method (DEM)

In this study the mixing kinetics and flow patterns of non-cohesive, monodisperse, spherical particles in a horizontal paddle blender were investigated using experiments, statistical analysis and discrete element method (DEM). EDEM 2.7 commercial software was used as the DEM solver. The experiment and simulation results were found to be in a good agreement. The calibrated DEM model was then utilized to examine the effects of the impeller rotational speed, vessel fill level and particle loading arrangement on the overall mixing quality quantified by the relative standard deviation (RSD) mixing index. The simulation results revealed as the impeller rotational speed was increased from 10?RPM to 40?RPM, generally a better degree of mixing was reached for all particle loading arrangements and vessel fill levels. As the impeller rotational speed was increased further from 40?RPM to 70?RPM the mixing quality was affected, for a vessel fill level of 60% and irrespective of the particle loading arrangement. Increasing the vessel fill level from 40% to 60% enhanced the mixing performance when impeller rotational speed of 40?RPM and 70?RPM were used. However, the mixing quality was independent of vessel fill level for almost all simulation cases when 10?RPM was applied, regardless of the particle loading arrangement. Furthermore, it was concluded that the particle loading arrangement did not have a considerable effect on the mixing index. ANOVA showed that impeller rotational speed had the strongest influence on the mixing quality, followed by the quadratic effect of impeller rotational speed, and lastly the vessel fill level. The granular temperature data indicated that increasing the impeller rotational speed from 10?RPM to 70?RPM resulted in higher granular temperature values. By evaluating the diffusivity coefficient and Peclet number, it was concluded that the dominant mixing mechanism in the current mixing system was diffusion.     

Dry mixing of cathode composite powder for all-solid-state batteries using a high-shear mixer

All-solid-state lithium-ion batteries (ASSLIBs) are promising candidates for next-generation batteries because of their various attractive properties. The uniform mixing of active materials (AMs) and solid electrolytes (SEs) is important for high-performance ASSLIBs. However, most AMs and SEs have poor flowability owing to their small particle size, which makes it difficult to uniformly mix the AM and SE particles. This study is focused on a high-shear mixer (HSM) as a scalable method to uniformly mix the AM and SE particles. The objective of this study is to determine the optimal operating conditions for HSM and its effectiveness in AM-SE mixing. The higher rotating speed of the chopper caused uniform SE dispersion by deagglomerating the SE particles and improving the adhesion of SE onto the AM particles, affording an electrode with well-balanced electrical/ionic conductivity and lower internal resistance. The ASSLIB with this electrode exhibited lower electrode polarization and excellent rate and cycle performance. Additionally, it has been demonstrated that the HSM could lead to a more uniform SE dispersion than conventional lab-scale mixing methods, resulting in significantly improved battery performance. Moreover, insights into the process-homogeneity-performance relations have been obtained.     

DEM study on identification of mixing mechanisms in a pot blender

A pot blender with both blending and storage capabilities offers an advantage over a conventional rotating drum. However, the mixing mechanism of the pot blender is extremely complicated because the pot blender rotates and swings simultaneously. Owing to the lack of systematic investigations, the mixing mechanism of the pot blender has not been fully elucidated. In this study, we clarify the mixing mechanism of the pot blender by using the discrete element method. Simulation results reveal that the main mixing mechanism is convective mixing in the rotational direction and shear mixing in the axial direction. Moreover, the mixing performance is unaffected by particle density, whereas the velocity gradient in the axial direction, which mainly determines the axial mixing performance, is affected by the particle filling ratio. Considering the relationship between the variance of axial particle velocity and granular temperature, the filling ratio is shown to significantly influence the mixing efficiency in the pot blender. In addition, the dependency of shear and diffusive mixing on Lacey’s mixing index in the pot blender is newly clarified. Consequently, this study demonstrates essential insights into the mixing mechanism of the pot blender and the pot blender as an effective industrial mixer.     

Numerical study of the mixing process of binary-density particles in a bladed mixer

Discrete element method (DEM) simulations of binary mixing of particles with different densities were conducted to study the influence of density ratio, blade speed, and filling level on the particle dynamics and mixing performance in a bladed mixer. Four particles with different densities at different locations were tagged to discuss the influence of three factors on the particle trajectory and velocity field in the mixer. A method based on cubic polynomial fitting of relative standard deviation was used to determine the critical revolution during the mixing process. It was found that the non-dimensional tangential velocity decreases with the increase of the blade speed and filling level, the fluctuation of vertical velocity increases with the radial location, blade speed, and filling level, and it is more pronounced than the fluctuation of tangential and radial velocity during the mixing process. Results obtained indicate that the mixing performance of particles with different density increases with the decrease of density ratio and filling level, while it increases with the increase of blade speed.     

Scale-up and flow behavior of cohesive granular material in a four-bladed mixer: effect of system and particle size

Flow of cohesive granular materials with different moisture contents was examined in a four-bladed mixer via the discrete element method (DEM). Firstly, the mixer diameter (D) was increased while keeping the particle diameter (d) constant. It was observed that when the mixer diameter to the particle diameter ratio (D/d) was larger than a certain critical size (D/d ≥ 75), granular flow behaviors and mixing kinetics followed simple scaling relations. For D/d ≥ 75, flow patterns and mixing kinetics were found to be independent of system size, and velocities of particles scaled linearly with the tip speed of the impeller blades and particle diffusivities scaled with the tip speed of the blades and mixer diameter. These results suggest that past a certain system size the flow and mixing of cohesive particles in large-scale units can be predicted from smaller systems. Secondly, system size was kept constant and particle diameter was changed and it was observed that by keeping the Bond number constant (by changing the level of cohesion) the flow behavior and mixing patterns did not change, showing that larger particles can be used to simulate flow of smaller cohesive particles in a bladed mixer by matching the Bond numbers.     

A review of current techniques for the evaluation of powder mixing

Blending a mixture of powders to a homogeneous system is a crucial step in many manufacturing processes. To achieve a high quality of the end product, powder mixtures should be made with high content uniformity. For instance, producing uniform tablets depends on the homogeneous dispersion of active pharmaceutical ingredient (API), often in low level quantities, into excipients. To control the uniformity of a powder mixture, the first required step is to estimate the powder content information during blending. There are several powder homogeneity evaluation techniques which differ in accuracy, fundamental basis, cost and operating conditions. In this article, emerging techniques for the analysis of powder content and powder blend uniformity, are explained and compared. The advantages and drawbacks of all the techniques are reviewed to help the readers to select the appropriate equipment for the powder mixing evaluation. In addition, the paper highlights the recent innovative on-line measurement techniques used for the non-invasive evaluation of the mixing performance.     

Development of a method for estimating particles mixing curves in short DEM simulation time

A new method for estimating particles mixing curves by simulating the particles behavior for a short period of time using discrete element method (DEM) was developed. The mixing curve is the time variation of the mixing degree, and can be divided into the following two stages: one is the mixing degree increasing stage, and the other is the mixing degree stagnation stage. The mixing curves could be estimated from time variations of mixing and de-mixing rates by assuming that the stagnation occurs when the mixing and de-mixing are balanced. Assuming the mixing and de-mixing rates are represented by the first-order exponential functions, each rate was estimated by the simulations of particles behavior for a short period of time. The estimated mixing curves agree with experimental ones and the proposed estimation method reduced the calculation time to approximately one-fifth of the conventional one that does not use a method to reduce the mixing time. Therefore, the proposed estimation method could estimate the mixing curve in a short period of time.     

DEM simulation of the packing structure and wall load in a 2-dimensional silo

The filling and discharge of a two-dimensional wedged-bottom silo holding circular objects was modelled using DEM technique to examine the influence of method of silo filling on distribution of orientations of unit vectors normal to contact points (contact normals) and normal contact forces. It was found that packing structure determined through method of generation of grain bedding significantly influenced distribution of contact normals. Nearly hexagonal network of contact normals was obtained for central filling of silo while sprinkle filling provided higher anisotropy of contact normals. The significance of frictional conditions and number of particles in system on distribution of contact normals was analysed. Increase in number of grains reduced disturbance from boundaries on behaviour of assembly. Distribution of loads on silo bottom obtained in simulation for different wall roughness was found in qualitative agreement with experimental data.     

Influence of particle shape on mixing rate in rotating drums based on super-quadric DEM simulations

Fundamental research on the flow and mixing of non-spherical particles is critical for industrial production and design. In this paper, the Discrete Element Method (DEM) is used to study the flow and mixing of granular materials in the horizontal rotating drum, and the periodic boundary condition is employed to eliminate end wall effect. Super-quadric elements are adopted to describe spherical and non-spherical particles. The influences of rotating speed, blockiness, and aspect ratio on the mixing rate are investigated by the Lacey mixing index. The results show that the rotating speed has a primary effect on the mixing rate, whereas the effect of the particle shape on the mixing rate is a secondary factor for non-spherical granular systems. Moreover, the mixing rate of spherical and non-spherical particle systems is significantly different. The mixing rate of spheres is the lowest, and the cubes have a higher mixing rate than the cylinders. As the blockiness decreases or aspect ratio deviates from 1.0, the mixing rate decreases. Ordered face-to-face contacts and dense packing structures result in a higher mixing rate. The analysis of kinetic energy shows that particle shape affects the transfer efficiency of external energy to the granular systems. The translational kinetic energy of non-spherical particles is higher than that of spherical particles, and their rotational kinetic energy is lower than that of spheres. Meanwhile, the blockiness enhances the transfer efficiency of external energy to the non-spherical systems; in contrast, the aspect ratio reduces the energy conversion efficiency.     

橡胶连续混炼粉料混合均匀性仿真及实验研究

为了提高橡胶连续混炼中混炼胶质量稳定性,实现炭黑等粉体物料精确配比和均匀性混合问题,针对粉体物料在混合和输送过程存在复杂的物理性质,建立了炭黑等混合粉料的球体颗粒模型和Hertz接触力-位移模型,采用EDEM对典型粉体物料混合均匀性进行模拟仿真和粉体物料混合实验,对炭黑等粉体物料在橡胶连续混炼工艺中的混合均匀性进行分析,探求橡胶粉料连续混合和输送机理.研究发现:粉体物料混合仿真与实验测试结果具有较高的拟合性,表明在橡胶连续混炼工艺中可以在保证混合均匀性的前提下实现多粉体混合物的连续称量和输送,同时也验证了运用EDEM数据模拟仿真粉体物料混合的可行性.     

DEM simulation and experimental validation for mechanical response of ellipsoidal particles under confined compression

The representation of non-spherical particles in discrete element method (DEM) has not been addressed adequately. Although the multiple sphere method (MSM) is the most popular approach to describe non-spherical particle shape, the validity of the MSM has not been established yet. The purpose of this study is to examine the validity and adequacy of the MSM. A uni-axial confined compression test was designed and set up to study the mechanical behaviour of an ellipsoidal granular assembly under vertical loading and the load transfer to the contacting boundary. Four levels of multi-sphere approximation for an axi-symmetric ellipsoidal particle were employed in DEM simulation to investigate the adequacy of multi-sphere approximation. A comparison on compression characteristics between the numerical and experimental results was made and discussed in this paper. Most of the compared physical properties showed reasonable agreement, indicating that capturing the key linear dimensions of a non-spherical particle may be sufficient to predict reasonable results. A small number of sub-spheres (say, N?≥?5) for representing an axi-symmetric ellipsoidal particle can give plausible results. However, the DEM simulations also produced a certain extent of discrepancy in loading stiffness with experiments. Plausible explanations are provided and require further investigation.     

Recent progress on the discrete element method simulations for powder transport systems: A review

Powder transport systems are ubiquitous in various industries, where they can encounter single powder flow, two-phase flow with solids carried by gas or liquid, and gas–solid–liquid three-phase flow. System geometry, operating conditions, and particle properties have significant impacts on the flow behavior, making it difficult to achieve good transportation of granular materials. Compared to experimental trials and theoretical studies, the numerical approach provides unparalleled advantages over the investigation and prediction of detailed flow behavior, of which the discrete element method (DEM) can precisely capture complex particle-scale information and attract a plethora of research interests. This is the first study to review recent progress in the DEM and coupled DEM with computational fluid dynamics for extensive powder transport systems, including single-particle, gas–solid/solid–liquid, and gas–solid–liquid flows. Some important aspects (i.e., powder electrification during pneumatic conveying, pipe bend erosion, non-spherical particle transport) that have not been well summarized previously are given special attention, as is the application in some new-rising fields (ocean mining, hydraulic fracturing, and gas/oil production). Studies involving important large-scale computation methods, such as the coarse grained DEM, graphical processing unit-based technique, and periodic boundary condition, are also introduced to provide insight for industrial application. This review study conducts a comprehensive survey of the DEM studies in powder transport systems.     

Investigating the upper limit for applying the coarse grain model in a discrete element method examining mixing processes in a rolling drum

Mixing is an essential manufacturing process in various industries. The processing procedure and final product quality depend on the homogeneity of mixing. Because it is difficult to evaluate mixing systems experimentally, the discrete element method is commonly employed. However, as the number of particles increases, this approach incurs huge computational costs. The coarse grain model offers a potential solution, but its applicability has not been widely demonstrated; this study aimed to elucidate the upper limit for applying the coarse grain model. To determine the appropriate simulation parameters, calibrations were performed by comparing the powder bed in experiments versus simulations. Various mixing processes were numerically evaluated, and the mixing characteristics were qualitatively consistent among all coarse-grained ratios. These mixing systems were also evaluated quantitatively based on Lacey’s mixing index, which indicated that the upper limit of the coarse-grained ratio was five times. It is therefore important to secure a sufficient number of particles in each cell and to use an appropriate number of cells. This study clarified the upper application limit and criteria for the coarse grain model and verified the maximum coarse-grained ratio (five times). This approach can be used to determine the coarse-grained ratio and reduce computational costs.     

DEM numerical studies on the design and efficiency of the continuous operating triboelectric separator

Triboelectric separation has gained importance as a dry separator for mineral and coal beneficiation. The different grade components of coal powder can be classified, beneficiated, and segregated with the help of the triboelectric separator. A batch-mode triboelectric separator is being used in CSIR-Institute of Minerals & Materials Technology for separating pure coal particles from gangue quartz and kaolinite particles. The sticking of particles on the electrode plates forces the separator to be paused for cleaning the electrode plates. A moving belt over an electrode plate can prevent its direct contact with particles and also carries the particle down towards collection bins. The feasibility study of installing moving belts along with the optimum conditions for design and operating parameters were studied with the help of DEM simulations. The bi-modal frequency distributions of coal and quartz in collection bins ware observed. The rationale for observation of bi-modal distribution is attributed to the varying initial angles and velocities of the particles entering the active electric field region. The non-uniform probabilistic opportunities of the particles to attach with the respective belts resulted in two groups of particles; (a) particles attached to belt followed by their collection in the expected bins (b) particles with initial inclination opposite to the direction of external force and their collection in the bins on the other side. The second group of particles followed trajectories governed by external forces and got collected into the bins on the other side. The phenomena resulted in the bi-modal distribution of mass distribution.     

Discrete element modeling and experimental investigation of hot pressing of intermetallic NiAl powder

This paper presents the numerical and experimental analysis of hot pressing of NiAl powder with an emphasis on the best possible representation of its main stages: initial powder compaction and pressure-assisted sintering. The numerical study has been performed within the discrete element framework. In the paper, an original viscoelastic model of hot pressing has been used. In order to ensure that the applied values of material parameters in numerical simulations are appropriate, the reference literature has been reviewed. It produced the relations and equations to estimate the values of all required sintering material parameters of the considered viscoelastic model. Numerical simulations have employed the geometrical model of the initial dense specimen generated by a special algorithm which uses the real grain distribution of powder. The numerical model has been calibrated and validated through simulations of the real process of hot pressing of intermetallic NiAl material. The kinetics of compaction, sintering and cooling stage indicated by the evolution of density, shrinkage and densification rate have been studied. The comparison of numerical and experimental results has shown a good performance of the developed numerical model.     

A study of principle stress rotation on granular soils using DEM simulation of hollow cylinder test

This study presents the numerical modeling of a hollow cylinder test (HCT) on granular soils by modifying the TRUBAL code. The discrete element method (DEM) was employed for this purpose. Owing to the significant expenditure of the HCT, a verified numerical modeling for this test was developed. The introduced numerical model (HCTBALL) defines plane and cylindrical walls for the boundary conditions to be applied. In addition, the article presents an efficient method to apply the torque. The displacements of the inner and outer walls are interdependent while torsion was applied to control the intermediate principal stress parameter (b). To verify the model, the paper employs results from experimental HCTs on Firoozkooh sand under both monotonic loading and drained conditions. A comparison of the presented model and the experimental results shows that both are closely concordant.It was shown that the deviatoric stress decreases as the principal stress angle (α) increases. In addition, it was observed that by increasing the confining pressure, the internal friction angle (φ′) decreases; however, at higher confining pressures, this reduction is insignificant. Furthermore, this study investigates the coordination number and its relationship to volumetric strain variations.     

Numerical study of particle mixing in a tilted three-dimensional tumbler and a new particle-size mixing index

A new mixing index is proposed, which is an improved Lacey index based on coordination number fraction. The differences and similarities among many mixing indices are compared, including the new mixing index, the information entropy based on coordination number fraction, the Lacey index based on local concentration, and the information entropy based on local concentration. The first two indices are microscopic since the coordination number fraction is on particle-scale, whereas the latter two are mesoscopic as the local concentration is mesoscopic scale. The newly proposed mixing evaluation indices does not include inauthentic temporal oscillations. Moreover, using mixing index, the mixing characteristics of particles in a tilted tumbler are studied by discrete element method (DEM). The tumbler’s angle of tilt α = 0°, 10°, 20°, 30°, 40°, 50°, 60° and 70°, at five rotating velocities ω = 0.175, 0.35, 0.5, 0.6, 0.7 and 1.4 rad/s corresponding to Froude number F_r = 0.0025, 0.001, 0.002, 0.003, 0.004, 0.016 respectively are simulated. It is found that both increasing the tilt angle and the rotating speed have negative effects on the particle mixing within the scope of this study.     

Mathematical modeling of coating time in dry particulate coating using mild vibration field with bead media described by DEM simulation

To theoretically understand the previously reported dry particulate coating process using a mild vibration field with a bead media, a mathematical analysis model of the dry coating system was developed. In this coating process, an ordered mixture with coarse host particles (drug-loaded ion exchange resin, diameter approximately 100 µm) and fine guest particles (acrylic polymer particle, primary particle size of approximately 100 nm) is formed using a vibrating a vessel. Second, the guest particles on the host particulate surface are firmly fixed using the collision of coated particles zirconia beads (diameter 1.5 mm). Our model assumes that the unfixed guest particles are fixed by particle-to-particle collisions (C_c) provided by the apparatus, thereby increasing the coating ratio. C_c was estimated using the discrete element method and some experimental results. The model includes parameters such as the number of C_c, host particles and unfixed guest particles. The coating time simulated by the established model equation in this study fits well with the experimental results of the dry process. It depends on the ratio of the number of collisions contributing to the increased coating ratio to the number of unfixed guest particles on the surface of host particles.