Effect of particle polydispersity on micromechanical properties and energy dissipation in granular mixtures

A series of numerical tests was conducted to study the micromechanical properties and energy dissipation in polydisperse assemblies of spherical particles subjected to uniaxial compression. In general, distributed particle size assemblies with standard deviations ranging from 0% to 80% of the particle mean diameter were examined. The microscale analyses included the trace of the fabric tensor, magnitude and orien- tation of the contact forces, trace of stress, number of contacts and degree of mobilization of friction in contacts between particles. In polydisperse samples, the average coordination numbers were lower than in monodisperse assemblies, and the mobilization of friction was higher than in monodisperse assemblies due to the non-uniform spatial rearrangement of spheres in the samples and the smaller displacements of the particles. The effect of particle size heterogeneity on both the energy density and energy dissipation in systems was also investigated.     

Effect of particle polydispersity on micromechanical properties and energy dissipation in granular mixtures

A series of numerical tests was conducted to study the micromechanical properties and energy dissipation in polydisperse assemblies of spherical particles subjected to uniaxial compression. In general, distributed particle size assemblies with standard deviations ranging from 0% to 80% of the particle mean diameter were examined. The microscale analyses included the trace of the fabric tensor, magnitude and orientation of the contact forces, trace of stress, number of contacts and degree of mobilization of friction in contacts between particles. In polydisperse samples, the average coordination numbers were lower than in monodisperse assemblies, and the mobilization of friction was higher than in monodisperse assemblies due to the non-uniform spatial rearrangement of spheres in the samples and the smaller displacements of the particles. The effect of particle size heterogeneity on both the energy density and energy dissipation in systems was also investigated.     

Spectrum of density fluctuations in a particlefluid system—II. Polydisperse spheres

This paper presents an extension of the analysis shown in Part I to a polydisperse particle-fluid system. The density autocorrelation is shown to be a function of two quantities, a generalized Overlap function for which an analytical expression is derived, and the radial distribution function (RDF). In Fourier transform space, the density spectrum again appears to be a strong function of the mean particle size, and secondarily the mean particle separation distance. One unusual result is previously observed oscillations in the density spectrum of a monodisperse system of particles are severely dampened or even eliminated in the polydisperse case, depending on the width of the particle size distribution. Apparently contributions from different particle correlations interfere with each other, thereby reducing the coherent oscillations seen in the monodisperse particle-fluid system. Furthermore at large wavenumbers, the spectrum decays with a −2 power-law, independent of the shape of the particle size distribution. This behavior can be traced to the Overlap function which controls the behavior of the spectrum beyond the first peak. Remarkably the −2 power-law spectrum is determined by the shape of the particles (i.e. spheres) rather than their spatial distribution (RDF).

The effect of an asymptotically large pressure gradient on the correlation of several important higher-order moments is revisited for the polydisperse system. The relatively simple relationships developed for the monodisperse system are lost in the polydisperse case because particles of different sizes will be influenced differently by an applied pressure gradient. The result is moments that are of different order in velocity can no longer be related to each other (as they were in the monodisperse system), even in this idealized flow. A more comprehensive understanding of this phenomenon can only be achieved through direct numerical simulation or experiment.     

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.     

Modeling of the random packing of a loose layer of polydisperse spherical particles

A method for calculating the loose packing structure of polydisperse spherical particles with a predetermined size distribution function is proposed. The coordinates of the particle centers in the loose layer are determined as the result of random fall of single spheres on a substrate under the action of gravity, assuming the inelastic collision of the spheres and considering the force of their adhesive interaction, and also assuming that the motion of one sphere on the surface of the other is pure slip. Numerical simulation is used to obtain the pattern of arrangement of polydisperse spherical particles in the loose powder layer, whose porosity depends on the particle size distribution function. The results are compared with experimental data.     

Numerical solution of an inverse problem of nozzle theory for a two-phase gas-particle mixture

In the present paper gas flows with monodisperse and polydisperse particles in plane and axisymmetric nozzles are calculated by the inverse method [1, 2]. The gas velocity distribution is specified on the axis of symmetry of the nozzle, while the gas and particle parameters are specified in the entrance section. As a result of the numerical integration of a system of equations describing a flow of gas with condensate particles in it we determine the gas and particle parameters, the gas streamlines, and the particle trajectories with allowance for the mutual influence of the gas and particles. One of the gas streamlines is taken as the nozzle contour and the limiting trajectories and pure gas zone are found. A difference method is described which makes it possible to calculate the subsonic, transonic, and supersonic flow regions using a single algorithm, its features are noted, and the results of the calculation for monodisperse mixtures with particle diameters 1 and 5 m and fractions by weight 0.3 are given. A comparison is made with the results of calculations by other methods.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 106–114, July–August, 1986.The authors express their gratitude to N. B. Ponomarev and G. E. Dumnov for their useful discussions and help in carrying out the calculations.     

TRANSMISSION AND DIFFRACTION PATTERN BY A MONOLAYER OF OPAQUE SPHERES

Transmission and diffraction of radiation through a monolayer of monodisperse opaque spherical particles at different particle concentrations are studied theoretically and by numerical calculation. It is found that the transmission and the diffraction pattern are strongly dependent on the receiving aperture and the monolayer density.     

Differentially weighted direct simulation Monte Carlo method for particle collision in gas-solid flows

In gas-solid flows,particle-particle interaction(typical,particle collision) is highly significant,despite the small particles fractional volume.Widely distributed polydisperse particle population is a typical characteristic during dynamic evolution of particles(e.g.,agglomeration and fragmentation) in spite of their initial monodisperse particle distribution.The conventional direct simulation Monte Carlo(DSMC)method for particle collision tracks equally weighted simulation particles,which results in high statistical noise for particle fields if there are insufficient simulation particles in less-populated regions.In this study,a new differentially weighted DSMC(DW-DSMC) method for collisions of particles with different number weight is proposed within the framework of the general Eulerian-Lagrangian models for hydrodynamics.Three schemes(mass,momentum and energy conservation) were developed to restore the numbers of simulation particle while keeping total mass,momentum or energy of the whole system unchanged respectively.A limiting case of high-inertia particle flow was numerically simulated to validate the DW-DSMC method in terms of computational precision and efficiency.The momentum conservation scheme which leads to little fluctuation around the mass and energy of the whole system performed best.Improved resolution in particle fields and dynamic behavior could be attained simultaneously using DW-DSMC,compared with the equally weighted DSMC.Meanwhile,computational cost can be largely reduced in contrast with direct numerical simulation.     

Differentially weighted direct simulation Monte Carlo method for particle collision in gas–solid flows

In gas–solid flows, particle–particle interaction (typical, particle collision) is highly significant, despite the small particles fractional volume. Widely distributed polydisperse particle population is a typical characteristic during dynamic evolution of particles (e.g., agglomeration and fragmentation) in spite of their initial monodisperse particle distribution. The conventional direct simulation Monte Carlo (DSMC) method for particle collision tracks equally weighted simulation particles, which results in high statistical noise for particle fields if there are insufficient simulation particles in less-populated regions. In this study, a new differentially weighted DSMC (DW-DSMC) method for collisions of particles with different number weight is proposed within the framework of the general Eulerian–Lagrangian models for hydrodynamics. Three schemes (mass, momentum and energy conservation) were developed to restore the numbers of simulation particle while keeping total mass, momentum or energy of the whole system unchanged respectively. A limiting case of high-inertia particle flow was numerically simulated to validate the DW-DSMC method in terms of computational precision and efficiency. The momentum conservation scheme which leads to little fluctuation around the mass and energy of the whole system performed best. Improved resolution in particle fields and dynamic behavior could be attained simultaneously using DW-DSMC, compared with the equally weighted DSMC. Meanwhile, computational cost can be largely reduced in contrast with direct numerical simulation.     

Numerical simulation of oxide nanoparticle growth characteristics under the gas detonation chemical reaction by space-time conservation element–solution element method

Under harsh conditions (such as high temperature, high pressure, and millisecond lifetime chemical reaction), a long-standing challenge remains to accurately predict the growth characteristics of nanosize spherical particles and to determine the rapid chemical reaction flow field characteristics. The growth characteristics of similar spherical oxide nanoparticles are further studied by successfully introducing the space-time conservation element–solution element (CE/SE) algorithm with the monodisperse Kruis model. This approach overcomes the nanosize particle rapid growth limit set and successfully captures the characteristics of the rapid gaseous chemical reaction process. The results show that this approach quantitatively captures the characteristics of the rapid chemical reaction, nanosize particle growth and size distribution. To reveal the growth mechanism for numerous types of oxide nanoparticles, it is very important to choose a rational numerical method and particle physics model.     

Fluid Permeability in Stratified Unconsolidated Particulate Bed

Fluid permeability of polydisperse particulate bed with finite thickness has been examined. On the assumption of creeping flow, the permeability of monodisperse particles with arbitrary arrangement is calculated by means of Stokesian dynamics approach in which the interaction between individual particles and interstitial fluid is described by multipole expansion of the Oseen tensor. We have extended such calculation method to polydisperse particulate systems which have not so dense structures (up to particle volume fraction ${\phi \sim}$ 0.2). The particles are located infinitely in space and their interaction has been taken into account by Ewald summation technique. For the spatial distribution of polydisperse particles, we consider locally stratified particulate beds and define stratification degree as a parameter which apparently and mathematically represents the thickness of the mixing region of different-sized particles. The permeability profiles in the particulate beds with different stratification degree show the dependence of local permeability on the spatial and size distribution of particles. Consequently, the calculation results indicate that the permeability of non-uniform polydisperse particulate bed can be predicted by integrating the local permeation resistance which is determined by the local specific surface area.     

Effects of Particle Size Non-Uniformity on Transport and Retention in Saturated Porous Media

Short-pulse injection experiments are investigated to study the effects of particle size non-uniformity on the transport and retention in saturated porous media. Monodisperse particles (3, 10, and 16 \(\upmu \hbox {m}\) latex microspheres) and polydisperse particles (containing 3, 10, and 16 latex microspheres) were explored. The obtained results suggest considering not only the particle sizes but also their polydispersivity (particle size non-uniformity) in transport and retention. Although, the density of the suspended particles is close to that of water, results reveal a slow transport of particles compared to the dissolved tracer whatever their size and flow velocity. The recovered particles in the mixture experiments show that the retention of large particles (10 and 16 \(\upmu \hbox {m}\)) enhances the retention of small ones (3 \(\upmu \hbox {m}\)). However, the straining of 10 and 16 \(\upmu \hbox {m}\) particles in “mixture experiments” is smaller than their straining in “monodisperse experiments”. A linear relationship summarizing the simultaneous effect of particle sizes and flow velocity on deposition kinetics coefficient is proposed.     

Regimes of sedimentation of a consolidated system of solid spherical particles

The results of the experimental investigation of gravitational sedimentation of a consolidated system of solid monodisperse spherical particles in a viscous liquid are represented over wide ranges of the particle number density and the Reynolds and Stokes numbers. Empirical dependences of the velocity of sedimentation of a particle aggregate and the drag coefficient of a system of particles as functions of the initial volume number density are obtained. The boundary values of the particle number density separating the sedimentation regimes are determined.     

Study of fluid cell coarsening for CFD-DEM simulations of polydisperse gas–solid flows

Particle polydispersity is ubiquitous in industrial fluidized beds, which possesses a significant impact on hydrodynamics of gas–solid flow. Computational fluid dynamics-discrete element method (CFD-DEM) is promising to adequately simulate gas–solid flows with continuous particle size distribution (PSD) while it still suffers from high computational cost. Corresponding coarsening models are thereby desired. This work extends the coarse-grid model to polydisperse systems. Well-resolved simulations with different PSDs are processed through a filtering procedure to modify the gas–particle drag force in coarse-grid simulations. We reveal that the drag correction of individual particle exhibits a dependence on filtered solid volume fraction and filtered slip velocity for both monodisperse and polydisperse systems. Subsequently, the effect of particle size and surrounding PSD is quantified by the ratio of particle size to Sauter mean diameter. Drag correction models for systems with monodisperse and continuous PSD are developed. A priori analysis demonstrates that the developed models exhibit reliable prediction accuracy.     

Particle size distribution in CPFD modeling of gas–solid flows in a CFB riser

A computational particle fluid dynamics (CPFD) numerical method to model gas–solid flows in a circulating fluidized bed (CFB) riser was used to assess the effects of particle size distribution (PSD) on solids distribution and flow. We investigated a binary PSD and a polydisperse PSD case. Our simulations were compared with measured solids concentrations and velocity profiles from experiments, as well as with a published Eulerian-Eulerian simulation. Overall flow patterns were similar for both simulation cases, as confirmed by experimental measurements. However, our fine-mesh CPFD simulations failed to predict a dense bottom region in the riser, as seen in other numerical studies. Above this bottom region, distributions of particle volume fraction and particle vertical velocity were consistent with our experiments, and the simulated average particle diameter decreased as a power function with riser height. Interactions between particles and walls also were successfully modeled, with accurate predictions for the lateral profiles of particle vertical velocity. It was easy to implement PSD into the CPFD numerical model, and it required fewer computational resources compared with other models, especially when particles with a polydisperse PSD were present in the heterogeneous flow.     

Particle size distribution in CPFD modeling of gas-solid flows in a CFB riser

A computational particle fluid dynamics(CPFD) numerical method to model gas-solid flows in a circulating fluidized bed(CFB) riser was used to assess the effects of particle size distribution(PSD) on solids distribution and flow.We investigated a binary PSD and a polydisperse PSD case.Our simulations were compared with measured solids concentrations and velocity profiles from experiments,as well as with a published Eulerian-Eulerian simulation.Overall flow patterns were similar for both simulation cases,as confirmed by experimental measurements.However,our fine-mesh CPFD simulations failed to predict a dense bottom region in the riser,as seen in other numerical studies.Above this bottom region,distributions of particle volume fraction and particle vertical velocity were consistent with our experiments,and the simulated average particle diameter decreased as a power function with riser height.Interactions between particles and walls also were successfully modeled,with accurate predictions for the lateral profiles of particle vertical velocity.It was easy to implement PSD into the CPFD numerical model,and it required fewer computational resources compared with other models,especially when particles with a polydisperse PSD were present in the heterogeneous flow.     

蒙特卡洛方法数值研究大气颗粒物动力学效应和辐射传输性质

爆发性增强的雾天,空气污染严重能见度低,这与大气边界层湍流性质、悬浮颗粒的动力学及散射性质密切相关.文中基于颗粒群平衡方程和Mie理论,采取加权蒙特卡洛方法,自行开发了Fortran程序.文中计算所得的颗粒尺度分布函数、颗粒散射性质与实验值、理论解一致,验证了数值模型和方法的正确性.此外,数值研究了雾爆发性增强阶段雾滴谱拓宽、能见度降低的机理,讨论湍流输运和颗粒局部聚集效应下颗粒间的碰并过程,并耦合颗粒散射性质,数值分析雾发展中湍流耗散率对颗粒对径向相对速度、系统透过率的影响;以及颗粒对径向相对速度与系统透过率、颗粒尺度的关系.研究结果表明:随着湍流耗散率的增大,颗粒的径向相对速度呈现先缓慢而后快速增大的变化趋势.1000s时刻,湍流的耗散率为1.0×10-2m2/s3,颗粒径向相对速度(无量纲)为0.0969;对于0.6μm的可见光,雾环境颗粒系统的透过率为0.47.此外,雾发展中雾滴易与气溶胶碰并,系统的散射性质与水组成的雾滴系统不同,天气的能见度明显降低.      

Transport of Neutrally Buoyant and Dense Variably Sized Colloids in a Two-Dimensional Fracture with Anisotropic Aperture

The transport of monodisperse as well as polydisperse colloid suspensions in a two-dimensional, water saturated fracture with spatially variable and anisotropic aperture is investigated with a particle tracking model. Both neutrally buoyant and dense colloid suspensions are considered. Although flow and transport in fractured subsurface formations have been studied extensively by numerous investigators, the transport of dense, polydisperse colloid suspensions in a fracture with spatially variable and anisotropic aperture has not been previously explored. Simulated snapshots and breakthrough curves of ensemble averages of several realizations of a log-normally distributed aperture field show that polydisperse colloids exhibit greater spreading than monodisperse colloids, and dense colloids show greater retardation than neutrally buoyant colloids. Moreover, it is demonstrated that aperture anisotropy oriented along the flow direction substantially increases colloid spreading; whereas, aperture anisotropy oriented transverse to the flow direction retards colloid movement.     

Evolution of shock waves in polydisperse gas suspensions with a nonuniform particle concentration distribution

The results of a numerical investigation of the laws of shock wave propagation in polydisperse (two-fraction) gas suspensions with a non-uniform initial particle concentration distribution are presented. Examples of shock wave propagation in extended layers of a gas suspension with linearly increasing, linearly decreasing and sinusoidal laws of variation of the particle concentration are considered. It is shown that when shock waves pass through layers of a gas suspension with increasing and decreasing laws of variation of the particle concentration, respectively, amplification and attenuation of the waves are observed; when shock waves travel through gas suspensions with a periodic law of variation of the particle concentration the pressure distribution behind the wave fronts is nonmonotonic. The solutions corresponding to polydisperse and monodisperse gas suspensions with an effective particle size are examined. The nonequilibrium and thermodynamic-equilibrium solutions are compared.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 183–190, September–October, 1991.