Temperature Properties in Polydisperse Granular Mixtures

Numerical simulations are employed to consider the problem of determining the granular temperatures of the species of a homogeneous heated granular mixture with a power-law size distribution. The partial granular temperature ratios are studied as functions of the fractal dimension D, the restitution coefficient e, the rescaled viscosity time, the average occupied area fraction φ, the total particle number N and the number fraction. Different species of particles in a power-law system typically do not have the same mean kinetic energy, namely the granular temperature. It is found that the extent of nonequipartition of kinetic energy is determined by the fractal dimension D, the restitution coefficient e and the rescaled viscosity time, while is insensitive to the total particle number N , the area fraction φ and the number fraction.     

Dynamic behavior of non-uniform granular gases system with the fractal characteristic in one dimension

A one-dimensional dynamic model of polydisperse granular mixture with a power-law size distribution is presented, in which the particles are subject to inelastic mutual collisions and driven by Gaussian white noise. The particle size distribution of the mixture has the fractal characteristic, and a fractal dimension D as a measurement of the inhomogeneity of the particle size distribution is introduced. We define the global granular temperature and the kinetic pressure of the mixture, and obtain their expressions. By molecular dynamics simulations, we have mainly investigated how the inhomogeneity of the particle size distribution and the inelasticity of collisions influence the steady-state dynamic properties of the system, focusing on the global granular temperature, kinetic pressure, velocity distribution and distribution of interparticle spacing. Some novel results are found that, with the increase of the fractal dimension D, the global granular temperature and the kinetic pressure decrease, the velocity distribution deviates more obviously from the Gaussian one and the particles cluster more pronouncedly at the same value of the restitution coefficient e (0<e<1). On the other hand, as the restitution coefficient e decreases, the dynamic behavior has the similar evolution as above at the fixed fractal dimension D. The dynamic behavior changing with e and D is, respectively, presented.     

Dynamic Properties of Two-Dimensional Polydisperse Granular Gases

We propose a two-dimensional model of polydisperse granular mixtures with a power-law size distribution in the presence of stochastic driving. A fractal dimension D is introduced as a measurement of the inhomogeneity of the size distribution of particles. We define the global and partial granular temperatures of the multi-component mixture. By direct simulation Monte Carlo, we investigate how the inhomogeneity of the size distribution influences the dynamic properties of the mixture, focusing on the granular temperature, dissipated energy, velocity distribution, spatial clusterization, and collision time. We get the following results： a single granular temperature does not characterize a multi-component mixture and each species attains its own ＂granular temperature＂; The velocity deviation from Gaussian distribution becomes more and more pronounced and the partial density of the assembly is more inhomogeneous with the increasing value of the fractal dimension D; The global granular temperature decreases and average dissipated energy per particle increases as the value olD augments.     

Dynamic Properties of Two-Dimensional Polydisperse Granular Gases

We propose a two-dimensional model of polydisperse granular mixtures with a power-law size distribution in the presence of stochastic driving. A fractal dimension D is introduced as a measurement of the inhomogeneity of the size distribution of particles. We define the global and partial granular temperatures of the multi-component mixture. By direct simulation Monte Carlo, we investigate how the inhomogeneity of the size distribution influences the dynamic properties of the mixture, focusing on the granular temperature, dissipated energy, velocity distribution, spatial clusterization, and collision time. We get the following results: a single granular temperature does not characterize a multi-component mixture and each species attains its own "granular temperature"; The velocity deviation from Gaussian distribution becomes more and more pronounced and the partial density of the assembly is more inhomogeneous with the increasing value of the fractal dimension D; The global granular temperature decreases and average dissipated energy per particle increases as the value of D augments.     

Steady-State Properties of Driven Polydisperse Granular Mixtures

We represent a two-dimensional model of polydisperse granular mixtures with a power-law size distribution. The model consists of smooth hard disks in a rectangular box with inelastic collisions, driven by a homogeneous heat bath at zero gravity. The width of particle size distribution is characterized by the only
parameter, namely, the fractal dimension D. The energy dissipation of the mixture is increased as D increases or as e decreases. Furthermore, it is found that the steady-state properties of the mixture such as the collision rate, granular temperature, kinetic pressure and velocity distribution depend sensitively on size distribution parameter D.     

Numerical simulation of particle size effects on contact electrification in granular systems

Based on the contact charge transfer model between two particles due to a single collision proposed by Apodaca, the contact charges carried on a particle is derived due to multiple collisions, including the repeat collisions between two particles and the collisions with different particles, in mixed-size granular system of identical material. The effect of the particle size on the charges carried on the particle is simulated. The results indicate that for a mixed-size granular system, due to multiple collisions among particles, there exists a threshold particle radius, the particles with radius higher than which and the particles with radius lower than which carry opposite charges. The threshold particle radius is equal to mean value of particle size in the mixed-size granular system. Basically, the polarity of the charges carried on the largest particle is same as the polarity of the transfer charge carrier, and in case of the positive charge transferred, the largest particle will be positively charged and the smallest particle will be negatively charged, and vice versa. In the same size region, the more dispersive the particle size is, the more the net charges can be produced. In normal-distributed granular system, the magnitude of contact charge is determined mainly by the particle size distribution, size region, total particle number and the relative impact velocity.     

Spatial Density Distributions and Correlations in a Quasi-one-Dimensional Polydisperse Granular Gas

By Monte Carlo simulations, the effect of the dispersion of particle size distribution on the spatial density distributions and correlations of a quasi one-dimensional polydisperse granular gas with fractal size distribution is investigated in the same inelasticity. The dispersive degree of the particle size distribution can be measured by a fractal dimension dr, and the smooth particles are constrained to move along a circle of length L, colliding inelastically with each other and thermalized by a viscosity heat bath. When the typical relaxation time τ of the driving Brownian process is longer than the mean collision time To, the system can reach a nonequilibrium steady state. The average energy of the system decays exponentially with time towards a stable asymptotic value, and the energy relaxation time τB to the steady state becomes shorter with increasing values of df. In the steady state, the spatial density distribution becomes more clusterized as df increases, which can be quantitatively characterized by statistical entropy of the system. Furthermore, the spatial correlation functions of density and velocities are found to be a power-law form for small separation distance of particles, and both of the correlations become stronger with the increase of df. Also, tile density clusterization is explained from the correlations.     

Particle dynamics simulations of triboelectric charging in granular insulator systems

Particle dynamics simulations are carried out to study triboelectric charging in granular systems composed of a single insulating material. The simulations implement a model in which electrons trapped in localized high energy states can be transferred during collisions to low energy states in the other particle. It is shown that this effect alone can generate electrostatic charging in the system, and cause net electron transfer from larger particles to smaller particles. The magnitude of charging is small for systems of a single particle size but becomes much greater for a system with polydispersal particle sizes, due to the net electron transfer from larger to smaller particles. The negative charge of smaller particles, and positive charge of larger particles has been observed in field studies and laboratory experiments of granular systems.     

Soot particles in piston-top pool fires and exhaust at 5 and 15 MPa injection pressure in a gasoline direct-injection engine

This study shows the structure of soot particles sampled directly from wall wetting-induced pool fires formed on the piston top in a spark-ignition direct-injection (SIDI) engine. Of particular interest is its variation with injection pressure considering the current trend of high-pressure DI system development to reduce engine-out particulate emissions. Thermophoretic particle sampling was performed for transmission electron microscope (TEM) imaging, which was post-processed for statistical analysis of key morphology and internal structure parameters. These include the size distribution of soot aggregates and primary particles as well as carbon-layer fringe-to-fringe gap and concentricity. With the fixed engine speed and load conditions, in-flame soot particles are compared to the exhaust particles sampled simultaneously at selected 5 and 15 MPa injection pressures corresponding to low-speed/low-load and high-speed/high-load in practical engine operation. From the TEM images and statistical analysis, it was found that more concentrated and taller pool fire developed for 5-MPa injection leads to smaller soot aggregates composed of smaller soot primary particles due to suppressed soot growth in fuel-rich flames. However, the soot particles formed in 15-MPa injection-induced pool fires are at a more reactive status evidenced by less defined core-shell boundaries and higher fringe separation. The higher soot reactivity results in enhanced soot oxidation, which explains smaller soot aggregates and primary particles found for the 15-MPa injection in the exhaust sample.     

Spray flow-network flow transition of binary Lennard-Jones particle system

We simulate gas-liquid flows caused by rapid depressurization using a molecular dynamics model. The model consists of two types of Lennard-Jones particles, which we call liquid particles and gas particles. These two types of particles are distinguished by their mass and strength of interaction: a liquid particle has heavier mass and stronger interaction than a gas particle. By simulations with various initial number densities of these particles, we found that there is a transition from a spray flow to a network flow with an increase of the number density of the liquid particles. At the transition point, the size of the liquid droplets follows a power-law distribution, while it follows an exponential distribution when the number density of the liquid particles is lower than the critical value. The comparison between the transition of the model and that of models of percolation is discussed. The change of the average droplet size with the initial number density of the gas particles is also presented.     

Effects of Fractal Size Distributions on Velocity Distributions and Correlations of a Polydisperse Granular Gas

By the Monte Carlo method, the effect of dispersion of disc size distribution on the velocity distributions and correlations of a polydisperse granular gas with fractal size distribution is investigated in the same inelasticity. The dispersion can be described by a fractal dimension D, and the smooth hard discs are engaged in a two- dimensional horizontal rectangular box, colliding inelastically with each other and driven by a homogeneous heat bath. In the steady state, the tails of the velocity distribution functions rise more significantly above a Gaussian as D increases, but the non-Gaussian velocity distribution functions do not demonstrate any apparent universal form for any value of D. The spatial velocity correlations are apparently stronger with the increase of D. The perpendicular correlations are about half the parallel correlations, and the two correlations are a power-law decay function of dimensionless distance and are of a long range. Moreover, the parallel velocity correlations of postcollisional state at contact are more than twice as large as the precollisional correlations, and both of them show almost linear behaviour of the fractal dimension D.     

Strength and failure of cemented granular matter

Cemented granular materials (CGMs) consist of densely packed solid particles and a pore-filling solid matrix sticking to the particles. We use a sub-particle lattice discretization method to investigate the particle-scale origins of strength and failure properties of CGMs. We show that jamming of the particles leads to highly inhomogeneous stress fields. The stress probability density functions are increasingly wider for a decreasing matrix volume fraction, the stresses being more and more concentrated in the interparticle contact zones with an exponential distribution as in cohesionless granular media. Under uniaxial loading, pronounced asymmetry can occur between tension and compression both in strength and in the initial stiffness as a result of the presence of bare contacts (with no matrix interposed) between the particles. Damage growth is analyzed by considering the evolution of stiffness degradation and the number of broken bonds in the particle phase. A brutal degradation appears in tension as a consequence of brittle fracture in contrast to the more progressive nature of damage growth in compression. We also carry out a detailed parametric study in order to assess the combined influence of the matrix volume fraction and particle-matrix adherence. Three regimes of crack propagation can be distinguished corresponding to no particle damage, particle abrasion and particle fragmentation, respectively. We find that particle damage scales well with the relative toughness of the particle-matrix interface with respect to the particle toughness. This relative toughness is a function of both matrix volume fraction and particle-matrix adherence and it appears therefore to be the unique control parameter governing transition from soft to hard behavior.     

梯度界面层对复合材料有效光学性质的影响

 研究了球形颗粒被无规地浸没在基质中,具有梯度分界面和光滑分界面的球形颗粒复合物界面层对介电响应和光学性质的影响;利用第一原理近似,计算了在稀释极限下具有梯度分布界面层复合物的有效介电常数,讨论了其实部和吸收光谱随外加场频率的变化关系。对结果进行比较,发现梯度分布的界面层使复合材料在等离子体共振频率处的有效介电常数的涨落明显减小,吸收谱线明显宽化并伴有“蓝移”现象,梯度界面层能很好地控制等离子体共振和减小有效介电响应的涨落。     

Coexistence of two granular temperatures in binary vibrofluidized beds

An investigation into the granular temperature distributions of a binary vibrofluidized granular bed has been conducted using positron emission particle tracking. By repeating each experiment with the tracer selected in turn from the two size components, the granular temperature and packing fraction distributions for each phase were determined. It was found that, for a range of size fractions, the granular temperature of the larger particles was higher than that of the smaller diameter grains, a result which was supported by a simple theoretical analysis based on the steady state energy equation.     

Size and Shape Prediction of Colloidal Metal Oxide MgBaFeO Particles from Light Scattering Measurements

The size and structure of colloidal metal oxide (MgBaFeO) particles are determined using an Elliptically Polarized Light Scattering (EPLS) technique. The approach is based on a hybrid experimental/theoretical study where the experimental data are compared against predictions obtained using a T-Matrix model that accounts for particle shape irregularities. A power-law distribution function with two parameters is employed to account for the particle size distribution. The refractive index of the particles in ethyl alcohol is calculated based on the Maxwell-Garnet formula. The experiments are conducted using a second-generation nephelometer. It is shown that the current EPLS measurements can effectively be used for identification of both the shape and the size of the colloids.     

Effect of particle size distribution on the polarity of triboelectric charging in granular insulator systems

Triboelectric charging occurs in granular insulating systems even when all particles are composed of identical material. A simple model is used here to address triboelectric charging in such systems. The basis of the model is the existence of electrons trapped in high-energy states, which can be released during collisions with another particle and transferred to the other particle. This model shows that triboelectric charging in insulator systems composed of particles of identical material can be attributed to a distribution of particle sizes, such that smaller particles tend to charge negatively and larger particles tend to charge positively. This polarity of charging has been observed in field studies of sand storms, dust devils and volcanic plumes, and most laboratory experiments on triboelectric charging in granular systems.     

双分散颗粒体系在临界堵塞态的结构特征

堵塞行为是颗粒体系中一种常见的现象,其力学性质与堆积结构的关联非常复杂.本文采用离散元法研究了由两种不同半径颗粒组成的二维双分散无摩擦球形颗粒体系在临界堵塞态所呈现的结构特征,讨论了大小颗粒粒径比与大颗粒百分比对临界堵塞态的影响.数值模拟结果表明,当粒径比小于1.4时,临界平均接触数与大颗粒百分比关系不大,当粒径比大于1.4时随着大颗粒百分比的增大临界平均接触数先减小再增大.而临界体积分数在粒径比小于1.8时随着大颗粒百分比的增加先减小后增大,大于1.8时又基本不随大颗粒百分比而变化.大颗粒百分比在接近0或1时,系统近似为单分散体系,临界平均接触数与体积分数基本不随半径比的增大而变化;在接近0.5时,临界平均接触数随着半径比的增大逐渐减小,而临界体积分数则是先减小后增大.文中对大-小颗粒这一接触类型的百分比也进行了探讨,其值随着大颗粒百分比的增大呈二次函数的变化趋势,粒径比对这一变化趋势只有较小的影响.     

Bidisperse granular avalanches on inclined planes: A rich variety of behaviors

Experiments were performed to provide insight into the flow behavior and structure of bimodal mixtures of grains in gravity-driven, free-surface flows. Unsteady unconfined flows were produced by releasing instantaneously a dry granular mass, composed of two particle sizes, over a rough inclined plane. As a result of size segregation, the small particles are found at the bottom of the flow and final deposit, the large particles are found at the free surface, but also on the lateral borders and at the front of the flow. The lateral and vertical inhomogeneous repartitions of particles lead to two main effects that are completely absent in monodispersed flows. The outline effect results from the accumulation of large beads on the periphery of the flow depending on the value of the relative friction of each particle species on the plane. This effect in turn causes a narrowing of the flow and/or an increase of length of the final deposit. The interface effect results of the interaction between layers of different size particles and causes the modification of the thickness of the deposit. These effects occur simultaneously and their combination leads to a great variety of behaviors. In this investigation, evidence of the diversity of behaviors is presented as the size ratio, relative friction and concentration of each particle species are varied.     

非球形椭球粒子参数变化对光偏振特性的影响

针对自然界中多数沙尘、烟煤粒子的非球形问题, 在球形粒子偏振特性的基础上, 进一步研究非球形椭球粒子的折射率、有效半径、粒子形状等参数变化对光偏振特性的影响, 采用基于T矩阵的非球形粒子仿真方法, 模拟非偏振光经椭球粒子传输后光的偏振特性及其与球形粒子间的差异, 并以实际沙尘、海洋、烟煤三种气溶胶粒子为例说明结果的正确性. 结果表明: 当折射率实部越小, 虚部越大时, 球形粒子与非球形粒子的偏振差异越不明显; 当粒子有效半径增加时, 球形粒子偏振度的变化比非球形粒子更为明显, 且最大值分别出现在散射角为150°和120°的位置; 当粒子形状不同时, 不同形状椭球及球形粒子的差异在散射角小于60° 时并不明显, 且当椭球粒子纵横比互为倒数时, 两种粒子的偏振特性近似相同. 通过以上分析可知, 在光传输过程中, 椭球粒子多数情况下无法被近似为球形粒子进行计算.     

The Fractal Characteristic of Particles in Granular Material Flows and Its Effect on Effective Thermal Conductivity

According to the fact that many pulverized particles possess fractal characteristic, a fractal model for studying fine particles in granular material flows is first proposed. An expression of particles' fractal distribution is derived to describe the relationship between the particle fractal dimensions and particle velocity distribution function. In accordance with this model, the theoretical particle effective thermal conductivity is derived. The analytical results show that for the small Biot-Fourier number, the effective thermal conductivity increases with the square root of the granular temperature. For very large Biot-Fourier number, the effective thermal conductivity linearly increases with the granular temperature. Numerically calculated results show that the thermal conductivity increases with the particle size fractal dimensions and decreases with the particle surface fractal dimensions.