Dynamic Simulation of Random Packing of Polydispersive Fine Particles

In this paper, we perform molecular dynamic (MD) simulations to study the two-dimensional packing process of both monosized and random size particles with radii ranging from 1.0 to 7.0 μm. The initial positions as well as the radii of five thousand fine particles were defined inside a rectangular box by using a random number generator. Both the translational and rotational movements of each particle were considered in the simulations. In order to deal with interacting fine particles, we take into account both the contact forces and the long-range dispersive forces. We account for normal and static/sliding tangential friction forces between particles and between particle and wall by means of a linear model approach, while the long-range dispersive forces are computed by using a Lennard-Jones-like potential. The packing processes were studied assuming different long-range interaction strengths. We carry out statistical calculations of the different quantities studied such as packing density, mean coordination number, kinetic energy, and radial distribution function as the system evolves over time. We find that the long-range dispersive forces can strongly influence the packing process dynamics as they might form large particle clusters, depending on the intensity of the long-range interaction strength.     

高压电场内细颗粒堆积机理研究

细颗粒的堆积既是电厂尾部除尘系统的核心问题,也是航天领域空间环境内限制动力或光学元件性能的关键问题。电场力因具有长程有效性和强可操作性成为控制细颗粒堆积的主要手段。本文通过微观实验研究了高压平行板电场间细颗粒堆积的机理,观察到颗粒在预极化、预荷电、变外电场电压等工况下的堆积形貌,并发展了图像处理的方法统计堆积颗粒数。研究表明,偶极力是外电场下颗粒成链的主因,而颗粒的倒伏则是来自曳力的作用,颗粒链的极限高度主要受外电场场强和颗粒堆积结构的影响。     

基于超二次曲面的颗粒材料缓冲性能离散元分析

自然界或工业中普遍是由非球形颗粒组成的复杂体系,与球形颗粒相比,非球形颗粒间的高离散和咬合互锁可使冲击载荷引起的能量有效衰减实现缓冲作用.基于连续函数包络的超二次曲面单元能准确地描述非球形颗粒的几何形态,并可精确地计算单元间的接触碰撞作用.本文采用离散元方法对冲击载荷作用下非球形颗粒物质的缓冲性能进行数值分析,并与圆柱体冲击的理论结果和球体冲击的实验结果进行对比验证.在此基础之上,进一步研究了筒底作用力在不同颗粒层厚度和形状等因素影响下的变化规律.计算结果表明:不同颗粒形状都存在一个临界厚度H_c.当HH_c时,缓冲率随H的增加而增加;当HH_c时,缓冲率的变化不再显著并趋于稳定值.此外,减小颗粒表面尖锐度和增加或减小圆柱形和长方形颗粒的长宽比都会提高颗粒材料的缓冲效果.     

Role of interparticle forces in the formation of random loose packing

We present a physical and numerical study of the settling of uniform spheres in liquids and show that interparticle forces play a critical role in forming the so-called random loose packing (RLP). Different packing conditions give different interparticle forces and, hence, different RLP. Two types of interparticle forces are identified: process dependent and process independent. The van der Waals force, as the major cohesive force in the present study, plays a critical role in effecting the process-dependent forces such as drag and lift forces. An equation is formulated to describe the relationship between the macroscopic packing fraction and microscopic interparticle forces in a packing. We argue there is no lowest packing fraction for a mechanically stable RLP; hence, the packing fractions of RLP can range from 0 to 0.64 depending on the cohesive and frictional conditions between particles.     

湿颗粒堆力学特性的离散元法模拟研究

利用基于线性黏聚接触模型的离散元法对不同颗粒系统的堆积过程进行了数值模拟研究,分析了颗粒形状和湿颗粒间液桥力对颗粒堆积形态的影响机理,获得了球形和块状湿颗粒堆基底表面所受的法向力以及堆中颗粒间的法向力和切向力"中心凹陷"式的分布规律,讨论了颗粒形状和黏聚能量密度对基底表面作用力和颗粒间作用力的影响.研究结果表明,颗粒形状和液桥力对颗粒堆的堆积形态具有显著的影响.堆积角随着黏聚能量密度的增加而增大,并且相同条件下的块状颗粒堆积角大于球形颗粒.颗粒形状和黏聚能量密度对基底表面所受作用力和堆中颗粒间的作用力变化及最大幅值均有影响作用.当黏聚能量密度值逐渐增大时,颗粒堆的作用力最大幅值均逐渐增大,并且块状颗粒堆的作用力最大幅值大于球形颗粒堆.当黏聚能量密度值过大时,颗粒堆力学特性更加复杂,液桥力对颗粒堆积特性的影响作用大于颗粒形状的影响.     

Collective motion of polar active particles on a sphere

Collective motion of active particles with polar alignment is investigated on a sphere. We discussed the factors that affect particle swarm motion and define an order parameter that can show the degree of particle swarm motion. In the model, we added a polar alignment strength, along with Gaussian curvature, affecting particles swarm motion. We find that when the force exceeds a certain limit, the order parameter will decrease with the increase of the force. Combined with our definition of order parameter and observation of the model, the reason is that particles begin to move side by side under the influence of polar forces. In addition, the effects of velocity, rotational diffusion coefficient, and packing fraction on particle swarm motion are discussed. It is found that the rotational diffusion coefficient and the packing fraction have a great influence on the clustering motion of particles, while the velocity has little influence on the clustering motion of particles.     

The dynamics of multiparticle collisions in motion of a granular material

We consider the complex problem of how to calculate particle motions taking into account multiparticle collisions. Multiparticle contacts occur when a particle collides with neighbouring particles, so that those contacts have a direct influence on each other. We will focus on the molecular dynamics method. Particularly, we will analyse what happens in cohesive materials during multiparticle contacts. We investigated the expression of repulsive force formulated under fractional calculus which is able to control dynamically the transfer and dissipation of energy in granular media. Such approach allows to perform simulations of arbitrary multiparticle collisions and also granular cohesion dynamics.     

Coupling of ferromagnetic nanoparticles through dipolar interactions

We consider two ferromagnetic nanoparticles coupled via long-range dipolar interactions. We model each particle by a three-dimensional array of classical spin vectors, with a central spin surrounded by a variable number of shells. Within each particle only ferromagnetic coupling between nearest neighbor spins is considered. The interaction between particles is of the dipolar type and the magnetic properties of the system is studied as a function of temperature and distance between the centers of the particles. We perform Monte Carlo simulations for particles with different number of shells, and the magnetic properties are calculated via two routes concerning the dipolar contribution: one assuming a mean-field like coupling between effective magnetic moments at the center of the particles, and other one, where we take into account interactions among all the pairs of spins, one in each particle. We show that the dipolar coupling between the particles enhances the critical temperature of the system relative to the case in which the particles are very far apart. The dipolar energy between the particles is smaller when the assumption of effective magnetic moment of the particles is used in the calculations.     

On hydraulic permeability of random packs of monodisperse spheres: Direct flow simulations versus correlations

Hydraulic permeability is studied in porous media consisting of randomly distributed monodisperse spheres by means of computational fluid dynamics (CFD) simulations. The packing of spheres is generated by inserting a certain number of nonoverlapping spherical particles inside a cubic box at both low and high packing fractions using proper algorithms. Fluid flow simulations are performed within the interparticulate porous space by solving Navier-Stokes equations in a low-Reynolds laminar flow regime. The hydraulic permeability is calculated from the Darcy equation once the mean values of velocity and pressure gradient are calculated across the particle packing. The simulation results for the pressure drop across the packing are verified by the Ergun equation for the lower range of porosities (<0.75), and the Stokes equation for higher porosities (∼1). Using the results of simulations, the effects of porosity and particle diameters on the hydraulic permeability are investigated. Simulations precisely specified the range of applicability of empirical or semi-empirical correlations for hydraulic permeability, namely the Carman-Kozeny, Rumpf-Gupte, and Howells-Hinch formulas. The number of spheres in the model is gradually decreased from 2000 to 20 to discover the finite-size effect of pores on the hydraulic permeability of spherical packing, which has not been clearly addressed in the literature. In addition, the scale dependence of hydraulic permeability is studied via simulations of the packing of spheres shrunk to lower scales. The results of this work not only reveal the validity range of the aforementioned correlations, but also show the finite-size effect of pores and the scale-independence of direct CFD simulations for hydraulic permeability.     

粉末材料堆积的物理模型与仿真系统

研究了粉末材料堆积过程仿真的物理模型和系统,并探讨了适合多种不同粒径颗粒混合堆积过程仿真的高性能计算方法.在该仿真系统中,考虑了重力、接触力、阻尼力、摩擦力和范德瓦耳斯力等多种作用力的影响,集成了多种接触力模型和阻尼模型,使其适用于三维大规模粉末材料堆积过程的计算机仿真.利用该系统对粉末材料领域中的两个典型应用进行了模拟研究.模拟了两种相同密度不同粒径颗粒（粒径比为10）的混合堆积过程.当小颗粒数为大颗粒数的300倍时,得到最大的堆积密度（体积分数）为0.82.另外,还模拟了两种不同密度相同粒径颗粒的混合堆积过程.当堆积结束时,出现了明显的分离（segregation）现象和团聚现象.所研究的物理模型和仿真系统既可用于粉末材料堆积过程研究,亦可用于普通的球形物体堆积过程的模拟研究. 关键词： 粉末堆积 物理模型 仿真系统 离散元法     

Relaxation of a test particle in systems with long-range interactions: diffusion coefficient and dynamical friction

We study the relaxation of a test particle immersed in a bath of field particles interacting via weak long-range forces. To order 1/N in the N→+∞ limit, the velocity distribution of the test particle satisfies a Fokker-Planck equation whose form is related to the Landau and Lenard-Balescu equations in plasma physics. We provide explicit expressions for the diffusion coefficient and friction force in the case where the velocity distribution of the field particles is isotropic. We consider (i) various dimensions of space d=3,2 and 1; (ii) a discret spectrum of masses among the particles; (iii) different distributions of the bath including the Maxwell distribution of statistical equilibrium (thermal bath) and the step function (water bag). Specific applications are given for self-gravitating systems in three dimensions, Coulombian systems in two dimensions and for the HMF model in one dimension.     

Numerical Simulation of Random Close Packings in Particle Deformation from Spheres to Cubes

Variation of packing density in particle deforming from spheres to cubes is studied. A new model is presented to describe particle deformation between different particle shapes. Deformation is simulated by relative motion of component spheres in the sphere assembly model of a particle. Random close packings of particles in deformation form spheres to cubes are simulated with an improved relaxation algorithm. Packings in both 2D and 3D cases are simulated. With the simulations, we find that the packing density increases while the particle sphericity decreases in the deformation. Spheres and cubes give the minimum (0.6404) and maximum (0.7755) of packing density in the deformation respectively. In each deforming step, packings starting from a random configuration and from the final packing of last deforming step are both simulated. The packing density in the latter case is larger than the former in two dimensions, but is smaller in three dimensions. The deformation model can be applied to other particle shapes as well.     

Fluctuations in granular media

Dense slowly evolving or static granular materials exhibit strong force fluctuations even though the spatial disorder of the grains is relatively weak. Typically, forces are carried preferentially along a network of "force chains." These consist of linearly aligned grains with larger-than-average force. A growing body of work has explored the nature of these fluctuations. We first briefly review recent work concerning stress fluctuations. We then focus on a series of experiments in both two- and three-dimension [(2D) and (3D)] to characterize force fluctuations in slowly sheared systems. Both sets of experiments show strong temporal fluctuations in the local stress/force; the length scales of these fluctuations extend up to 10(2) grains. In 2D, we use photoelastic disks that permit visualization of the internal force structure. From this we can make comparisons to recent models and calculations that predict the distributions of forces. Typically, these models indicate that the distributions should fall off exponentially at large force. We find in the experiments that the force distributions change systematically as we change the mean packing fraction, gamma. For gamma's typical of dense packings of nondeformable grains, we see distributions that are consistent with an exponential decrease at large forces. For both lower and higher gamma, the observed force distributions appear to differ from this prediction, with a more Gaussian distribution at larger gamma and perhaps a power law at lower gamma. For high gamma, the distributions differ from this prediction because the grains begin to deform, allowing more grains to carry the applied force, and causing the distributions to have a local maximum at nonzero force. It is less clear why the distributions differ from the models at lower gamma. An exploration in gamma has led to the discovery of an interesting continuous or "critical" transition (the strengthening/softening transition) in which the mean stress is the order parameter, and the mean packing fraction, gamma, must be adjusted to a value gamma(c) to reach the "critical point." We also follow the motion of individual disks and obtain detailed statistical information on the kinematics, including velocities and particle rotations or spin. Distributions for the azimuthal velocity, V(theta), and spin, S, of the particles are nearly rate invariant, which is consistent with conventional wisdom. Near gamma(c), the grain motion becomes intermittent causing the mean velocity of grains to slow down. Also, the length of stress chains grows as gamma-->gamma(c). The 3D experiments show statistical rate invariance for the stress in the sense that when the power spectra and spectral frequencies of the stress time series are appropriately scaled by the shear rate, Omega, all spectra collapse onto a single curve for given particle and sample sizes. The frequency dependence of the spectra can be characterized by two different power laws, P proportional, variant omega(-alpha), in the high and low frequency regimes: alpha approximately 2 at high omega; alpha<2 at low omega. The force distributions computed from the 3D stress time series are at least qualitatively consistent with exponential fall-off at large stresses. (c) 1999 American Institute of Physics.     

Simulation of random packing of spherical particles with different size distributions

A numerical model for a loose packing process of spherical particles is presented. The simulation model starts with randomly choosing a sphere according to a pregenerated continuous particle-size distribution, and then dropping the sphere into a dimension-specified box, and obtaining its final position by using dropping and rolling rules which are derived from a similar physical process of spheres dropping in the gravitational field to minimize its gravity potential. Effects of three different particle-size distributions on the packing structure were investigated. Analysis on the physical background of the powder-based manufacturing process is additionally applied to produce optimal packing parameters of bimodal and Gaussian distributions to improve the quality of the fabricated parts. The results showed that higher packing density can be obtained using bimodal size distribution with a particle-size ratio from 1.5 to 2.0 and the mixture composition around n ₂:n ₁=6:4. For particle size with a Gaussian distribution, the particle radii should be limited in a narrow range around 0.67 to 1.5.     

Contact forces in a granular packing

We present the results of a systematic numerical investigation of force distributions in granular packings. We find that all the main features of force transmission previously established for two-dimensional systems of hard particles hold in three-dimensional systems and for soft particles, too. In particular, the probability distribution of normal forces falls off exponentially for forces above the mean force. For forces below the mean, this distribution is either a decreasing power law when the system is far from static equilibrium, or nearly uniform at static equilibrium, in agreement with recent experiments. Moreover, we show that the forces below the mean do not contribute to the shear stress. The subnetwork of the contacts carrying a force below the mean thus plays a role similar to a fluid surrounding the solid backbone composed of the contacts carrying a force above the mean. We address the issue of the computation of contact forces in a packing at static equilibrium. We introduce a model with no local simplifying force rules, that allows for an exact computation of contact forces for given granular texture and boundary conditions. (c) 1999 American Institute of Physics.     

Electrohydrodynamic controlled assembly and fracturing of thin colloidal particle films confined at drop interfaces

Particles can adsorb strongly at liquid interfaces due to capillary forces, which in practice can confine the particles to the interface. Here we investigate the electrohydrodynamic flow driven packing and deformation of colloidal particle layers confined at the surface of liquid drops. The electrohydrodynamic flow has a stagnation point at the drop equator, leading to assembly of particles in a ribbon shaped film. The flow is entirely controlled by the electric field, and we demonstrate that AC fields can be used to induce hydrodynamic “shaking” of the colloidal particle film. We find that the mechanical properties of the film is highly dependent on the particles: monodisperse polystyrene beads form packed granular monolayers which “liquefies” upon shaking, whereas clay mineral particles form cohesive films that fracture upon shaking. The results are expected to be relevant for understanding the mechanics and rheology of particle stabilized emulsions.     

Direct observation of anisotropic interparticle forces in nematic colloids with optical tweezers

Interparticle forces in a nematic liquid-crystal colloid have been directly observed by the dual beam laser trapping method with pN sensitivity. We introduce two different types of spatial distributions of forces, detected between the particles accompanied by hyperbolic hedgehog defects. These force distributions lead to specific particle arrangements, which are both stabilized by the balance of the orientational stress field of nematics. On the basis of these results, we propose novel artificial construction for multiparticle regular arrangements.     

Short-time dynamics of a packing of polyhedral grains under horizontal vibrations

We analyze the dynamics of a 3D granular packing composed of particles of irregular polyhedral shape confined inside a rectangular box with a retaining wall subjected to horizontal harmonic forcing. The simulations are performed by means of the contact dynamics method for a broad set of loading parameters. We explore the vibrational dynamics of the packing, the evolution of solid fraction and the scaling of dynamics with the loading parameters. We show that the motion of the retaining wall is strongly anharmonic as a result of jamming and grain rearrangements. It is found that the mean particle displacement scales with inverse square of frequency, the inverse of the force amplitude and the square of gravity. The short-time compaction rate grows in proportion to frequency up to a characteristic frequency, corresponding to collective particle rearrangements between equilibrium states, and then it declines in inverse proportion to frequency.