Interparticle Collisions Research Articles

Performance Evaluation of Polygonal and Cylindrical Ball Mills Using Discrete Element Method

ABSTRACTBall mills are frequently utilized in the chemical and mineral processing industries to reduce the size of particles. In mineral processing, cylindrical rotating mills with wear‐resistant materials are used at various operation scales. Polygon‐shaped rotating mills are used in artisanal and small‐scale mining due to technological and economic reasons. It is imperative to thoroughly evaluate polygon‐shaped mills to inform decisions regarding their improvement or replacement based on technological data. Even a minor enhancement in grinding efficiency and mill durability can yield substantial economic and environmental advantages. In this study, the performance of polygonal‐shaped mills was evaluated and compared with cylindrical profiles, taking into account power draw, collision energy dissipation, relative wear, and operational stability. The effect of installing lifters in the mills was also assessed. Simulations using the discrete element method (DEM) were carried out to describe the mechanical behavior of grinding media and its interaction with the mill walls. Notably, polygonal‐shaped mills without lifters facilitated interparticle interaction with limited centrifuging of particles compared to cylindrical mills without lifters. The installation of lifters in polygonal mills increased interparticle collisions and reduced power draw and wear rate. Cylindrical mills with lifters had high interparticle collisions similar to polygonal profiles but with much lower wear rate and high operational stability.

A multiscale investigation on protein addition toward steering agglomeration and yield in spray drying

The positive effect of protein drying aids that reduce the stickiness of sugar-rich products on the spray drying yield is known. However, as agglomeration also depends on stickiness, one can expect an effect of these drying aids on interparticle collisions. Therefore, the effect of protein addition to maltodextrin systems on agglomeration and yield in spray drying was investigated using single droplet drying and pilot-scale spray drying experiments. Single droplet drying was used to compare the sticking regimes of whey and pea proteins at different concentrations. Adding protein increased the chance of creating agglomerates upon forced collisions between a drying droplet and a glass bead. This was reflected in an increased fraction of agglomerated particles in pilot-scale experiments with injection of fines. Pea protein decreased the fraction of non-agglomerated primary particles from 38% to 22% when the protein content in the liquid feed was 50 g/kg dry basis. In pilot-scale spray drying, protein addition improved the yield. Industrially, these results can be used to steer formulations regarding agglomeration and yield.

Establishment of particle motion model and study of particle volume fraction distribution in the turbo air classifier based on DSMC

In order to study the particles' motion law and collision behavior in the flow field of the turbo air classifier, A particle motion simulation platform is developed using Compute Unified Device Architecture (CUDA) based on the Visual Studio. The particle-eddy interaction model, the particle-wall collision model, and the inter-particle collision Direct Simulation Monte Carlo (DSMC) solver are implemented. The results show that the deviation between calculated cut size d50 and experimental d50 using DSMC is smallest compared to DSMC, DDPM and DPM. DSMC has better acceleration performance compared to DDPM. The more particles need to be calculated, the more significant this advantage is. The influences of feeding rate on particle volume fractions and inter-particle collision number are analyzed. The increase of feeding rate leads to the increase of the particle volume fractions and inter-particle collision number.

An improved DSMC method for acoustic agglomeration of solid particles assisted by spray droplets

The direct simulation Monte Carlo (DSMC) method was improved for describing the process of acoustic agglomeration assisted by spray droplets. The agglomeration and rebound as consequences of inter-particle collisions were modeled by considering all possible particle types, in particular, mixed-phase particles associated with the immersion and distribution mechanisms. Numerical predictions by the present method were validated by both analytical solutions and experimental data. The dynamic process of the acoustic agglomeration in terms of the size and type evolution as well as the evolution of number concentration of different particle types were examined. Results suggest that the strong interaction between the solid particles and the spray droplets causes the significant enhancement of acoustic agglomeration. Furthermore, the effect of frequency varying from audible to ultrasonic range on the acoustic agglomeration was evaluated. Overall, a better performance of acoustic agglomeration can be achieved at a lower frequency, regardless of addition of spray droplets. Agglomeration efficiency of 64.8 % can be achieve at acoustic frequency of 1000 Hz in case with spray droplets. This study provides not only a numerical model for describing the agglomeration of complex particles in particle-laden gas flows, but also important insights into the acoustic agglomeration assisted by spray droplets.

Three-dimensional measurement method of binary particle collisions under dry and wet conditions

The description, calculation and design of solids processes, such as fluidization, pneumatic conveying or mixing, require knowledge of collision dynamics between particles and apparatus walls or interparticle collisions without, but often also with the presence of an additional liquid. For the study of such dry and wet interparticle collisions, an experimental setup is presented that allows the performance of dry or wet free binary collisions of spherical particles down to a minimum diameter of 1 mm. The recording of the particle motions in all three spatial directions is provided by two synchronized highspeed cameras. Methods and algorithms are presented for analyzing the translational and rotational motion in three-dimensional space, as well as for measuring the amount of liquid on the particles in the case of wet collisions. Motion data from dry collisions between equal sized and unequal sized spherical particles are compared with a three-parameter collision model. In addition, a first series of measurements of wetted collisions is used to compare the collision dynamics with dry collisions.

Structure evolution analysis of asphalt mixtures during compaction based on particle tracking and granular properties

Pavement compaction is a process of structure evolution of materials accompanied by energy dissipation, which is closely related to the inherent granular properties of asphalt mixture. This research achieved dynamic characterization of particle migration energy and contact force development through particle tracking technique and discrete element simulation. Analytical theories for multi-stage structure evolution and compaction mechanisms were proposed based on granular properties. Results indicate that kinetic, dissipative, and potential energies dynamically represent stages of structure evolution, determined by particle collisions, self-organization into arches, and interlocking reinforcement. The initial density state is established after interparticle collisions, while self-organization occurs as particles adaptively rearrange spatially by rotation in response to unbalanced forces, accompanied by a cyclical strengthening trend of arch formation and collapse. Arches facilitate stress redistribution toward isotropy, with strong contact forces exhibiting a top-down and side-center trend, progressively supported by particles larger than 4.75 mm. Theoretical measures for compaction control were proposed, including pre-calculation of IC, selection of compaction temperature corresponding to the optimum energy efficiency, and control of principal stress deflection within the friction cone. This research aims to provide insight into granular properties and compaction dynamics, thereby contributing to intelligent compaction.

Instabilities of a circular moderately dense particle cloud impacted by an incident shock

Shock-driven gas-particle flows exist in a wide variety of physical systems, covering the range from dilute to dense regime. There are a lot of fundamental questions remaining to be answered, like the evolutions of flow instability, two-phase interactions, and inter-particle collision effects. In this work, the Eulerian-Lagrangian approach with gas-particle four-way coupling is applied to model two-dimensional shock-driven gas-particle flow under the moderately dense regime. For the gas flow, the formation mechanism of primary and secondary vortexes is discussed. The primary vortex is caused by the roll-up of slip line, while the secondary vortexes are induced by the velocity and density difference between the penetrating flow of particle-laden zone and the outside mainstream. Parametric studies of different Mach numbers, particle diameters and initial volume fractions have been performed to investigate the effects on shock behaviours and slip line instabilities. For the particle phase, a high-volume-fraction compression region forms at the upstream cloud zone, which blocks the penetrating gas, and results in the formation of slip line. The cloud tail is mainly composed of the upstream edge particles. The inter-particle collisions have impact on the particle distributions, and greatly affect the slip line stability at the upstream cloud edge. This work provides valuable insights into the instability mechanism of compressible gas-solid flow and the particle dispersion behaviours. However, it is important to note that our analysis is based on a two-dimensional configuration. In a real three-dimensional context, the generation of vortices is likely to differ, and may affect the particle dynamics to a certain degree.

Viscosity-modulated clustering of heated bidispersed particles in a turbulent gas

Clustering of externally and evenly heated particles is enhanced by the increased viscosity of heated fluid in the vicinity of these clusters – a phenomenon known as viscous capturing (VC). Herein we study, via direct numerical simulations of decaying turbulence, the effect of temperature-driven viscosity on clustering with different particle loading densities. We employ a two-way momentum and energy coupling, and gas viscosity is modelled by a power law to understand the role of the increased drag and particle back-reaction force on the clustering intensity. For the continuum and dispersed phases, Eulerian and Lagrangian point particle schemes have been used, neglecting inter-particle collisions. We found that the enhanced viscosity-driven clustering is a strong function of particle loading density, as the increase in particle number density enables the formation of large uneven clusters before heating, which is the main condition for VC to take effect. Higher number density should result in greater turbulence modulation and negate local temperature-based viscous effects leading to VC. However, due to higher local particle number density in the clusters and interphase heat transfer, increased drag force prevails in such cases and delivers excessive clustering. By sampling conditionally the particle velocity and temperature inside the clusters, it is found that the thermodynamic and kinematic properties of the particles in the clusters are highly correlated, and this correlation increases with the particle loading density. Therefore, based on the particle number density, temperature-based viscosity can enhance considerably the clustering of heated particles and alter the effect of particles on the underlying turbulence.

Mediation of collisionless turbulent dissipation through cyclotron resonance

The dissipation of turbulence in astrophysical systems is fundamental to energy transfer and heating in environments ranging from the solar wind and corona to accretion disks and the intracluster medium. Although turbulent dissipation is relatively well understood in fluid dynamics, astrophysical plasmas often exhibit exotic behaviour, arising from the lack of interparticle collisions, which complicates turbulent dissipation and heating in these systems. Recent observations by NASA’s Parker Solar Probe mission in the inner heliosphere have shed new light on the role of ion cyclotron resonance as a potential candidate for turbulent dissipation and plasma heating. Here, using in situ observations of turbulence and wave populations, we show that ion cyclotron waves provide a major pathway for dissipation and plasma heating in the solar wind. Our results support recent theoretical predictions of turbulence in the inner heliosphere, known as the helicity barrier, that suggest a role of cyclotron resonance in ion-scale dissipation. Taken together, these results provide important constraints for turbulent dissipation and acceleration efficiency in astrophysical plasmas.

Modelling dense particle streams during free-fall electrostatic separation

Electrostatic particle separation is an attractive option for dry processing in the resource and recycling industries. Free-fall electrostatic separators are particularly simple and versatile member of this family of technologies. This paper presents numerical simulations of the movement of particles with different charges and masses within a stream of other particles in an electric field. The modelled system comprises much larger numbers and higher spatial densities of particles than previous models, and allows a realistic assessment of the effect of interparticle electrostatic forces and collisions. The effect of an angled feed with both co- and counter-directional fields is investigated, and the results translated to grade-recovery curves for different configurations. Among the important results were the degradation of separation performance with particle feed rate, the potential importance of Coulombic interparticle reactions in dense particle streams, and the benefits of an angled feed with a counter-directional electric field.

Effect of two-way coupling on clustering and settling of heavy particles in homogeneous turbulence

When inertial particles are dispersed in a turbulent flow at sufficiently high concentrations, the continuous and dispersed phases are two-way coupled. Here, we show via laboratory measurements how, as the suspended particles modify the turbulence, their behaviour is also profoundly changed. In particular, we investigate the spatial distribution and motion of sub-Kolmogorov particles falling in homogeneous air turbulence. We focus on the regime considered in Hassaini & Coletti (J. Fluid Mech., vol. 949, 2022, A30), where the turbulent kinetic energy and dissipation rate were found to increase as the particle volume fraction increases from $10^{-6}$ to $5\times 10^{-5}$ . This leads to strong intensification of the clustering, encompassing a larger fraction of the particles and over a wider range of scales. The settling rate is approximately doubled over the considered range of concentrations, with particles in large clusters falling even faster. The settling enhancement is due in comparable measure to the predominantly downward fluid velocity at the particle location (attributed to the collective drag effect) and to the larger slip velocity between the particles and the fluid. With increasing loading, the particles become less able to respond to the fluid fluctuations, and the random uncorrelated component of their motion grows. Taken together, the results indicate that the concentrated particles possess an effectively higher Stokes number, which is a consequence of the amplified dissipation induced by two-way coupling. The larger relative velocities and accelerations due to the increased fall speed may have far-reaching consequences for the inter-particle collision probability.

Magnetic field effect on the sedimentation process of two non-magnetic particles inside a ferrofluid

The behaviours of non-magnetic particles inside a ferrofluid have a significant effect on the performance of ferrofluid-based biosensors, such as antibody detection. The present work numerically investigates the sedimentation process of two non-magnetic spherical particles inside a ferrofluid under a uniform external magnetic field. A three-dimensional lattice Boltzmann model coupled with the bounce-back schemes is employed to simulate the corresponding flows. The present work particularly focuses on the effects of the magnetic field orientation and intensity on the sedimentation behaviours. The results show that the two particles are accelerated in a vertical magnetic field due to the reduction of drag, but decelerated in the magnetic field with the orientation of 30∘≤θ≤90∘. The inter-particle collision is encouraged by the attractive dipole force for 0∘≤θ≤45∘, resulting in an earlier occurrence of the ‘kissing’ stage, but is suppressed by the repulsive dipole force for 60∘≤θ≤90∘, leading to a postponement or even the absence of the ‘kissing’ stage. Meanwhile, it is found that the time duration of the ‘drafting’ stage decreases with the increasing of the magnetic field intensity. When the magnetic field intensity is beyond a critical value, the two particles self-assemble into a ‘dumbbell-like’ structure and the ‘tumbling’ stage will disappear during the sedimentation process. The present work demonstrates that the sedimentation behaviour of non-magnetic particles immersed in a ferrofluid can be controlled by varying the orientation and intensity of the external magnetic field, which can provide a valuable technical guidance for the industrial applications of ferrofluids.

Considering solid shear stress in MP-PIC simulation of a CFB riser

To enable successful multi-phase particle-in-cell (MP-PIC) simulations of circulating fluidized beds (CFBs), both the interphase drag force and interparticle collisions need to be carefully modeled. Particle collisions are usually represented by the solid stress consisting of the normal and shear components, in which the normal stress was found to have a leading role in the numerical stability of MP-PIC simulations, whereas the impact of the shear stress has seldom been reported. In this work, the effects of the solid shear stress are investigated by using two-dimensional simulations of the moderately dense laboratory-scale CFB riser with in-house MP-PIC code implemented on MFIX open-source platform. Cases including only normal solid stress and those considering both normal and shear solid stresses are simulated. The results of solids flux, axial and radial solids profiles are compared with available experimental data. The results show that the solid shear stress plays a minor role on the accuracy of simulation, and the increase of PPP (number particles per parcel) leads to a lower accuracy of simulation.

CFD-DEM model of plugging in flow with cohesive particles

Plugging in flows with cohesive particles is crucial in many industrial and real-life applications such as hemodynamics, water distribution, and petroleum flow assurance. Although probabilistic models for plugging risk estimation are presented in the literature, multiple details of the process remain unclear. In this paper, we present a CFD-DEM model of plugging validated against several experimental benchmarks. Using the simulations, we consider the process of plugging in a slurry of ice in decane, focusing on inter-particle collisions and plugging dynamics. We conduct a parametric study altering the Reynolds number (3000...9000), particle concentration (1.6...7.3%), and surface energy (21...541 mJ/m^2). We note the process possesses complex non-linear behaviour for the cases where particle-wall adhesion reduces by more than 20% relative to inter-particle cohesion. Finally, we demonstrate how the simulation results match the flow maps based on the third-party experiments.

Experimental signatures of the transition from acoustic plasmon to electronic sound in graphene.

Fermi liquids respond differently to perturbations depending on whether their frequency is higher (collisionless regime) or lower (hydrodynamic regime) than the interparticle collision rate. This results in a different phase velocity between the collisionless zero sound and the hydrodynamic first sound. We performed terahertz photocurrent nanoscopy measurements on graphene devices, with a metallic gate close to the graphene layer, to probe the dispersion of propagating acoustic plasmons, the counterpart of sound modes in electronic Fermi liquids. We report the observation of a change in the plasmon phase velocity when the excitation frequency approaches the electron-electron collision rate that is compatible with the transition between the zero and the first sound mode.

Fluctuation analysis and spatial distribution model of particle concentration in a lab-scale scrubbing-cooling chamber

The motion and distribution of ash from coal gasification during the scrubbing and cooling process are crucial for revealing the fluid flow pattern and improving the efficiency of multiphase separation. Experimental measurements of the local solid holdup (εs) in the liquid pool of a laboratory-scale gas-liquid-solid three-phase scrubbing-cooling chamber were carried out using an optical fiber probe. The fluctuations of εs in time and frequency domains under different operating conditions were discussed. The axial and radial distribution of εs was quantitatively analyzed. The results showed that inter-particle collision and particle-fluid interaction are the key factors leading to a fluctuation in particle concentration. Axial uniformity of particle concentration can be quantified by solid-phase volume fraction, particle Reynolds number, Archimedes number, and Morton number; εs in different regions is expressed according to the degree of influence of bubble wake as a correlation equation consisting only of cross-sectional average solid holdup and dimensionless radial position; a dimensionless expression for εs concerning operating conditions was established.

Numerical investigation on the influence of immersed tube bundles on biomass gasification in industrial-scale dual fluidized bed gasifier

Immersed tubes play a pivotal role in enhancing fluidization quality within fluidized beds by effectively disrupting bubbles. However, the effect of such internals within a large-scale gasification apparatus remains unclear. Consequently, this study focuses on numerically simulating biomass steam gasification within an industrial-scale dual fluidized bed gasifier that features immersed tube bundles using the multiphase particle-in-cell method. In this method, particles with identical properties are grouped, and the interparticle collisions are model via the solid stress to improve computational efficiency. The numerical model has been well verified with experiment in terms of bed pressure drop, solid recirculation and chemical reactions. Some main conclusions are shown as follows. Immersed tubes significantly alter the motion of gas and particles in the gasifier. The presence of immersed tubes can ameliorate solid entrainment and avoid excess particles escaping from the gasifier. The asymmetrical structure of the gasifier results in asymmetrical distribution of gas temperature and gas species. Increasing the number of immersed tubes slightly enhances particle temperature near the loop seal outlet. The radial biomass dispersion coefficient and axial sand dispersion coefficient decrease about 40.9 % and 28.9 % when increasing the immersed tube number from three to nine. The results provide enhanced understanding of the role of immersed tubes in altering gas-particle interactions and fluidization dynamics in industrial-scale fluidized bed gasifier.

Theoretical model and numerical simulation of high-gravity backwashing process for deep bed filter in oilfield

With the large-scale application of polymer flooding technology, the polymer concentration in produced water is increasing, which leads to incomplete regeneration of filter media by gravity field hydraulic backwashing technology only. In this study, the high-gravity backwashing technology is constructed by coupling gravity field and swirl field, which promotes the development of the oilfield. Based on the cascade theory, this paper establishes a mathematical expression of radial collision force in the high-gravity backwashing process in the vortex tube and solves a mathematical equation of high-gravity backwashing velocity gradient G based on the efficient zone of collision force action. Meanwhile, based on the discrete element method and computational fluid dynamics theory, the inter-particle collision regulation of the high-gravity backwashing process is simulated by the coupled FLUENT-EDEM method, and the high-efficiency zone of the backwashing process is verified. In this high-efficiency zone, the theoretical and simulated values of velocity gradient G for different high-gravity values match well, with differences ranging from –5.8 % to +7.1 %, indicating that the theoretical model of high-gravity backwashing is correct. When the high-gravity value is 210–240 g, the total impulse of the filter bed per unit time reaches a great value, and the backwashing velocity gradient G is 3047.1 s−1, which is much larger than the G value of combined air-water backwashing. The high-gravity backwashing technology breaks through the technical bottleneck of low efficiency in the action of the gravity field. It enriches and develops the gravity field hydraulic backwashing technology of the filter bed.

Comprehensive Euler/Lagrange modelling including particle erosion for confined gas-solid flows

The present research aims to assess the capability of a comprehensive Euler/Lagrange approach for predicting gas-solid flows and the associated solid particle erosion. The open-source code OpenFOAM® 4.1 was used to carry out the numerical simulations, where the standard Lagrangian libraries were substantially extended to account for all necessary models. Particles are tracked considering both translational and rotational motion as well as all relevant forces, such as gravity/buoyancy, drag and transverse lift due to shear and particle rotation. The tracking time step was dynamically adapted according to the locally relevant time scales, which drastically reduces computational times. Stochastic approaches are adopted to model particle turbulent dispersion, particle collisions with rough walls and particle-particle interactions. Five solid particle erosion models, available in the literature, were considered to estimate pipe bend erosion. Three study cases are provided to validate the adopted numerical approach and erosion models extensively. The first case intends to evaluate the ability of the extended CFD code to predict the behaviour of gas-solid flows in pneumatic conveying systems. This goal is achieved by comparing the numerical results with the experimental data obtained by Huber (1997) and Huber and Sommerfeld (1994, 1998) in a pneumatic conveying system. Here, the importance of considering inter-particle collisions and surface roughness for predicting particle velocity, mass flux and mean diameter distributions in gas-solid flows is highlighted. The second and third case intend to evaluate the ability of the erosion models in estimating bend erosion in diluted gas-solid flows. The erosion data obtained experimentally by Mazumder et al. (2008) and Solnordal et al. (2015) in very dilut pneumatic conveying systems is used for validating the numerical results, neglecting now inter-particle collisions and two-way coupling. Besides a comprehensive analysis of the different influential properties on erosion, the innovation of the present study is as follows. For the first time also a temporal modification of the surface roughness due to the erosion was considered in the simulations obtained from previous measurements (Novelletto Ricardo & Sommerfeld, 2020). As the surface roughness is increased due to erosion, eventually erosion rate becomes lower. This is the result of diminishing wall collision frequency. Simulations for several degrees of surface roughness showed that larger roughness is coupled with a drastic reduction of erosion. Hence, numerical simulations neglecting wall surface roughness are not realistic. The consideration of a particle size distribution instead of mono-sized computations showed a possible reduction of erosion rate. The detailed analysis of the different single-particle erosion models revealed that the model proposed by Oka et al. (2005) and Oka and Yoshida (2005) yields the best agreement with the measurements, however particle and wall properties are needed.

CFD-DEM modeling of turbidity current propagation in channels with two different topographic configurations

Submarine turbidity currents are a special type of sediment gravity flow responsible for turbidite deposits, attracting great interests from scientists and engineers in marine and petroleum geology. This paper presents a fully coupled computational fluid dynamics (CFD) and discrete element method (DEM) model to quantitatively analyze the turbidity current propagation in channels with two different topographic configurations. An appropriate drag force model is first incorporated in the CFD-DEM scheme, and two benchmark cases, including a single-particle sedimentation case and an immersed granular collapse case, are conducted to verify the accuracy of the developed CFD-DEM model. The model is then employed to investigate the fluid and particle dynamics of turbidity currents flowing over a flat bed (FB), and three obstacle-placed beds with different heights (OPB, OPB_1 and OPB_2). The CFD-DEM results indicate that the front position of turbidity current in the FB case is well consistent with the classic lock-exchange experiment. Results also show that the presence of the obstacle can clearly diminish the inter-particle collisions and the particle kinetic energy, weaken the particle-fluid interactions, and further make more sediment particles settle in front of the obstacle. Increase of obstacle height can result in diverse flow morphology of particles and fluids, and intensify the influences of obstacle on particle dynamics of turbidity currents. We show that our models enable reproducing the typical process of turbidity current propagation, and further can provide more valuable insights in understanding the turbidite-related geological phenomena from the point of view of particulate flow.