Jamming as a critical phenomenon: a field theory of zero-temperature grain packings

A field theory of frictionless grain packings in two dimensions is shown to exhibit a zero-temperature critical point at a nonzero value of the packing fraction. The zero-temperature constraint of force balance plays a crucial role in determining the nature of the transition. Two order parameters, , the deviation of the average number of contacts from the isostatic value, and , the average magnitude of the force per contact, characterize the transition from the jammed (high packing fraction) to the unjammed (low packing fraction state). The critical point has a mixed character with the order parameters showing a jump discontinuity but with fluctuations of the contact force diverging. At the critical point, the distribution of phi shows the characteristic plateau observed in static granular piles. The theory makes falsifiable predictions about the spatial fluctuations of the contact forces.     

Influence of the long-range forces in non-Gaussian random-packing dynamics

In this paper, we perform molecular dynamics (MD) simulations to study the random packing of spheres with different particle size distributions. In particular, we deal with non-Gaussian distributions by means of the Lévy distributions. The initial positions as well as the radii of five thousand non-overlapping particles are assigned inside a confining rectangular box. After that, the system is allowed to settle under gravity towards the bottom of the box. Both the translational and rotational movements of each particle are considered in the simulations. In order to deal with interacting particles, we take into account both the contact and long-range cohesive forces. The normal viscoelastic force is calculated according to the nonlinear Hertz model, whereas the tangential force is calculated through an accurate nonlinear-spring model. Assuming a molecular approach, we account for the long-range cohesive forces using a Lennard-Jones (LJ)-like potential. The packing processes are studied assuming different long-range interaction strengths.     

Random packings of frictionless particles

We conduct numerical simulations of random packings of frictionless particles at T = 0. The packing fraction where the pressure becomes nonzero is the same as the jamming threshold, where the static shear modulus becomes nonzero. The distribution of threshold packing fractions narrows, and its peak approaches random close packing as the system size increases. For packing fractions within the peak, there is no self-averaging, leading to exponential decay of the interparticle force distribution.     

Jamming III: Characterizing randomness via the entropy of jammed matter

The nature of randomness in disordered packings of frictional and frictionless spheres is investigated using theory and simulations of identical spherical grains. The entropy of the packings is defined through the force and volume ensemble of jammed matter and this is shown to be difficult to calculate analytically. A mesoscopic ensemble of isostatic states is then utilized in an effort to predict the entropy through the definition of a volume function that is dependent on the coordination number. Equations of state are obtained relating entropy, volume fraction and compactivity characterizing the different states of jammed matter, and elucidating the phase diagram for jammed granular matter. Analytical calculations are compared to numerical simulations using volume fluctuation analysis and graph theoretical methods, with reasonable agreement. The entropy of the jammed system reveals that random loose packings are more disordered than random close packings, allowing for an unambiguous interpretation of both limits. Ensemble calculations show that the entropy vanishes at random close packing (RCP), while numerical simulations show that a finite entropy remains in the microscopic states at RCP. The notion of a negative compactivity, which explores states with volume fractions below those achievable by existing simulation protocols, is also explored, expanding the equations of state. The mesoscopic theory reproduces the simulations results in shape well, though a difference in magnitude implies that the entire entropy of the packing may not be captured by the methods presented herein. We discuss possible extensions to the present mesoscopic approach describing packings from random loose packing (RLP) to RCP to the ordered branch of the equation of state in an effort to understand the entropy of jammed matter in the full range of densities from RLP to face-centered cubic (FCC) packing.     

Formation of superconductivity through interparticle interactions in ferrimagnetic-like Sn nanoparticle assemblies

We report on the observation of the development of superconductivity through interparticle interactions in 3, 5, 7, and 23 nm ferrimagnetic-like Sn nanoparticle assemblies. The Sn nanoparticles are fabricated using the gas condensation method. Each sample consists of a macroscopic amount of individual Sn nanoparticles without a capping molecule. Ferrimagnetism is found but no sign of superconductivity can be detected when the 3 nm particles are very loosely assembled. A reduction in the mean particle moment results when the packing fraction of the assembly is increased. Superconductivity occurs when a critical packing fraction is reached. Beyond this, the superconducting transition temperature T _C continues to increase and noticeably exceeds that of the bulk T _C. The enhancement of superconductivity by interparticle interactions has also been observed in 5, 7, and 23 nm particle assemblies, with the effect becoming less significant in larger particles. We attribute these observations to the transfer of electrons between the surface and the core regions of the nanoparticles triggered by finite size effects and interparticle interactions.     

Effect of interparticle forces on the coagulation of ultrafine SiO2 particles

The mechanism of ultrafine (<100 nm) SiO₂ particle production under thermal arc plasma conditions is studied by modelling. Two cases of the process are considered in the model: (i) when it is determined as a pure free-molecular coagulation; (ii) when the coagulation is influenced by the interparticle forces. The Hamaker formula is used to present the van der Waals forces between the particles. Particle size distribution functions (PSDF) are calculated for both cases. It is shown that inclusion of the interparticle interactions does not affect the self-preservation of the PSDF. The mean particle sizes are obtained from the PSDF and compared. Higher values are observed in the case that includes van der Waals forces. Comparison between experimental and calculated PSDF shows better agreement in the case considering interparticle forces.     

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

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

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.     

Physics of compaction of fine cohesive particles

Fluidized fractal clusters of fine particles display critical-like dynamics at the jamming transition, characterized by a power law relating consolidation stress with volume fraction increment [sigma--(c) proportional, variant(Deltaphi)(beta)]. At a critical stress clusters are disrupted and there is a crossover to a logarithmic law (Deltaphi = nu logsigma--(c)) resembling the phenomenology of soils. We measure lambda identical with- partial differentialDelta(1/phi)/ partial log(sigma--(c) proportional, variant Bo(0.2)(g), where Bo(g) is the ratio of interparticle attractive force (in the fluidlike regime) to particle weight. This law suggests that compaction is ruled by the internal packing structure of the jammed clusters at nearly zero consolidation.     

Jamming II: Edwards’ statistical mechanics of random packings of hard spheres

The problem of finding the most efficient way to pack spheres has an illustrious history, dating back to the crystalline arrays conjectured by Kepler and the random geometries explored by Bernal in the 1960s. This problem finds applications spanning from the mathematician’s pencil, the processing of granular materials, the jamming and glass transitions, all the way to fruit packing in every grocery. There are presently numerous experiments showing that the loosest way to pack spheres gives a density of ∼55% (named random loose packing, RLP) while filling all the loose voids results in a maximum density of ∼63%-64% (named random close packing, RCP). While those values seem robustly true, to this date there is no well-accepted physical explanation or theoretical prediction for them. Here we develop a common framework for understanding the random packings of monodisperse hard spheres whose limits can be interpreted as the experimentally observed RLP and RCP. The reason for these limits arises from a statistical picture of jammed states in which the RCP can be interpreted as the ground state of the ensemble of jammed matter with zero compactivity, while the RLP arises in the infinite compactivity limit. We combine an extended statistical mechanics approach ‘a la Edwards’ (where the role traditionally played by the energy and temperature in thermal systems is substituted by the volume and compactivity) with a constraint on mechanical stability imposed by the isostatic condition. We show how such approaches can bring results that can be compared to experiments and allow for an exploitation of the statistical mechanics framework. The key result is the use of a relation between the local Voronoi volumes of the constituent grains (denoted the volume function) and the number of neighbors in contact that permits us to simply combine the two approaches to develop a theory of volume fluctuations in jammed matter. Ultimately, our results lead to a phase diagram that provides a unifying view of the disordered hard sphere packing problem and further sheds light on a diverse spectrum of data, including the RLP state. Theoretical results are well reproduced by numerical simulations that confirm the essential role played by friction in determining both the RLP and RCP limits. The RLP values depend on friction, explaining why varied experimental results can be obtained.     

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.     

Onset of mechanical stability in random packings of frictional spheres

Using sedimentation to obtain precisely controlled packings of noncohesive spheres, we find that the volume fraction phiRLP of the loosest mechanically stable packing is in an operational sense well defined by a limit process. This random loose packing volume fraction decreases with decreasing pressure p and increasing interparticle friction coefficient mu. Using x-ray tomography to correct for a container boundary effect that depends on particle size, we find for rough particles in the limit p-->0 a new lower bound, phiRLP=0.550+/-0.001.     

Jamming of packings of frictionless particles with and without shear

By minimizing the enthalpy of packings of frictionless particles, we obtain jammed solids at desired pressures and hence investigate the jamming transition with and without shear. Typical scaling relations of the jamming transition are recovered in both cases. In contrast to systems without shear, shear-driven jamming transition occurs at a higher packing fraction and the jammed solids are more rigid with an anisotropic force network. Furthermore, by introducing the macrofriction coefficient, we propose an explanation of the packing fraction gap between sheared and non-sheared systems at fixed pressure.     

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.     

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

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

剪切振动下湿颗粒的力学谱

在恒温25 ℃剪切振动条件下,测量不同水分含量的NaCl湿颗粒体系的力学谱（能量耗散tanφ和剪切模量G）.研究发现,随着剪切振幅增大,NaCl湿颗粒体系的剪切模量G和能量耗散tanφ都表现出类似于干颗粒体系的阻塞（Jamming）转变行为.随着体系中水含量的增大,湿颗粒体系的剪切模量G和能量耗散tanφ在质量分数约等于11%的临界水浓度下均出现一个峰值,且峰位与应变振幅无关,表明此时颗粒之间主要的作用力发生了变化.     

A Dugdale–Barenblatt model for a strip with a semi-infinite crack embedded in decagonal quasicrystals

The present study is to determine the solution of a strip with a semi-infinite crack embedded in decagonal quasicrystals, which transforms a physically and mathematically daunting problem. Then cohesive forces are incorporated into a plastic strip in the elastic body for nonlinear deformation. By superposing the two linear elastic fields, one is evaluated with internal loadings and the other with cohesive forces, the problem is treated in Dugdale-Barenblatt manner. A simple but yet rigorous version of the complex analysis theory is employed here, which involves conformal mapping technique. The analytical approach leads to the establishment of a few equations, which allows the exact calculation of the size of cohesive force zone and the most important physical quantity in crack theory: stress intensity factor. The analytical results of the present study may be used as the basis of fracture theory of decagonal quasicrystals.     

二维颗粒体系单轴压缩形成的力链结构

从接触力、能量分布和接触网络结构特点出发,提出了强力链的力大小(Fc)判据(Fc大于平均接触力〈F〉)和角度θc判据(θc=180/〈Z〉,其中〈Z〉是平均配位数),指出强力链和弱力链是本质不同的两类结构存在于颗粒体系中,其中强力链网络与体系的宏观性质直接相关.以二维颗粒体系的单轴压缩为例,计算发现了强力链长度的幂率分布规律,分析了侧向压力系数与相应强力链网络结构的关系:当内部强力链网络充分发育而不再变化时,侧向压力系数趋于稳定数值.     

The relationship between attractive interparticle forces and bulk behaviour in dry and uncharged fine powders

Memento, homo, qui pulvis est et pulverem reverteris.

Genesis 3

Polvos serán, mas polvo enamorado.

Francisco de Quevedo

The physics of granular materials in ambient gases is governed by interparticle forces, gas–particle interaction, geometry of particle positions and geometry of particle contacts. At low consolidations these are strongly dependent on the external forces, boundary conditions and on the assembling procedure. For dry fine powders of micron and sub-micron particle size interparticle attractive forces are typically much higher than particle weight, and particles tend to aggregate. Because of this, cohesive powders fracture before breaking, flow and avalanche in coherent blocks much larger than the particle size. Similarly the drag force for micron sized particles is large compared to their weight for velocities as low as 1?mm/s. Due to this extreme sensitivity to interstitial gas flow, powders transit directly from plastic dense flows to fluidization without passing through collisional regimes with negligible gas interaction. These two features, strong attractive forces and strong gas interaction make powder behaviour differ qualitatively from the behaviour of large, noncohesive grains. In this paper we investigate the implications of these two features on the bulk powder behaviour. More in particular, the aim of this paper is to examine the relationship between attractive interparticle forces at grain level, with solid bulk properties at low consolidations (solid fraction, stresses), fluidization (aggregation, settling) and flow regime boundaries (plastic flow, inertial flow, fluidization and suspension). Many of the experimental results reported here are for dry and uncharged fine powders made of polymer particles of the order of 10 microns in diameter. However, the basic concepts and methodology are of general applicability.

The relationship between attractive interparticle forces and bulk behaviour in dry and uncharged fine powders

All authors

A. Castellanos

https://doi.org/10.1080/17461390500402657

Published online:

19 February 2007

Table     

Bose-einstein condensation in trapped dipolar gases

We discuss Bose-Einstein condensation in a trapped gas of bosonic particles interacting dominantly via dipole-dipole forces. We find that in this case the mean-field interparticle interaction and, hence, the stability diagram are governed by the trapping geometry. Possible physical realizations include ultracold heteronuclear molecules, or atoms with laser induced electric dipole moments.