Force Distributions near Jamming and Glass Transitions

We calculate the distribution of interparticle normal forces P(F) near the glass and jamming transitions in model supercooled liquids and foams, respectively. P(F) develops a peak that appears near the glass or jamming transitions, whose height increases with decreasing temperature, decreasing shear stress and increasing packing density. A similar shape of P(F) was observed in experiments on static granular packings. We propose that the appearance of this peak signals the development of a yield stress. The sensitivity of the peak to temperature, shear stress, and density lends credence to the recently proposed generalized jamming phase diagram.     

Jamming, two-fluid behavior, and "self-filtration" in concentrated particulate suspensions

We study the flow of model hard-sphere colloidal suspensions at high volume fraction Phi driven through a constriction by a pressure gradient. Above a particle-size dependent limit Phi(0), direct microscopic observations demonstrate jamming and unjamming-conversion of fluid to solid and vice versa-during flow. We show that such a jamming flow produces a reduction in colloid concentration Phi(x) downstream of the constriction. We propose that this "self-filtration" effect is due to a combination of jamming of the particulate part of the system and continuing flow of the liquid part, i.e., the solvent, through the pores of the jammed solid. Thus we link jamming in colloidal and granular media with a "two-fluid-like" picture of the flow of concentrated suspensions. Results are also discussed in the light of the original experiments of Reynolds on dilation in granular materials.     

Correlation time and diffusion coefficient imaging: application to a granular flow system

A parametric method for spatially resolved measurements for velocity autocorrelation functions, R(u)(tau) = , expressed as a sum of exponentials, is presented. The method is applied to a granular flow system of 2-mm oil-filled spheres rotated in a half-filled horizontal cylinder, which is an Ornstein-Uhlenbeck process with velocity autocorrelation function R(u)(tau) = e(- ||tau ||/tau(c)), where tau(c) is the correlation time and D = tau(c) is the diffusion coefficient. The pulsed-field-gradient NMR method consists of applying three different gradient pulse sequences of varying motion sensitivity to distinguish the range of correlation times present for particle motion. Time-dependent apparent diffusion coefficients are measured for these three sequences and tau(c) and D are then calculated from the apparent diffusion coefficient images. For the cylinder rotation rate of 2.3 rad/s, the axial diffusion coefficient at the top center of the free surface was 5.5 x 10(-6) m(2)/s, the correlation time was 3 ms, and the velocity fluctuation or granular temperature was 1.8 x 10(-3) m(2)/s(2). This method is also applicable to study transport in systems involving turbulence and porous media flows.     

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

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

Sheared force networks: anisotropies, yielding, and geometry

A scenario for the yielding of granular matter is presented by considering the ensemble of force networks for a given contact network and applied shear stress tau. As tau is increased, the probability distribution of contact forces becomes highly anisotropic, the difference between average contact forces along minor and major axes grows, and the allowed networks span a shrinking subspace of all force networks. Eventually, contacts start to break, and at the maximal shear stress the packing becomes effectively isostatic. The size of the allowed subspace exhibits simple scaling properties, which lead to a prediction for the yield stress for packings of an arbitrary contact number.     

Jamming of Granular Flow in a Two-Dimensional Hopper

We study experimentally the jamming phenomenon of granular flow of monodisperse disks of D = 5 mm diameter in a two-dimensional hopper with opening R. The jamming probability J(d) is measured where d identical withR/D. We found that J(d) decreases from 1 to zero when d increases from 2 to 5. From observing the disk configurations of the arch in the jamming events, the jamming probability can be explained quantitatively by treating the arch as the trajectory of a restricted random walker.     

Phase transition and bistable phenomenon of granular flows down a chute with successive turnings

This paper studies the granular flow down a chute with two successive turnings, which play the role of bottlenecks for the granular flow system and determine the granular flow state in main section between them. With the increase of main section width D, phase transition from dilute to dense granular flow is observed: When the main section width D is small (large), the granular flow at upper (lower) bottleneck is dense and the granular flow is dilute (dense) in the main section. More interestingly, a bistable region is exhibited, in which either dilute flow or dense flow may occur and continue for the entire run. In this region, the packing in the reservoir will affect initial flow rate and then affect the flow pattern. This study can be viewed as a paradigm for the jamming and unjamming transitions under shear due to gravity.     

Complex force network in marginally and deeply jammed solids

This paper studies the force network properties of marginally and deeply jammed packings of frictionless soft particlesfrom the perspective of complex network theory. We generate zero-temperature granular packings at different pressures by minimizing the inter-particle potential energy. The force networks are constructed as nodes representing particles and links representing normal forces between the particles. Deeply jammed solids show remarkably different behavior from marginally jammed solids in their degree distribution, strength distribution, degree correlation, and clustering coefficient. Bimodal and multi-modal distributions emerge when the system enters the deep jamming region. The results also show that small and large particles can show different correlation behavior in this simple system.     

Flow, ordering, and jamming of sheared granular suspensions

We study the rheological properties of a granular suspension subject to constant shear stress by constant volume molecular dynamics simulations. We derive the system "flow diagram" in the volume fraction or stress plane (phi, F): at low phi the flow is disordered, with the viscosity obeying a Bagnold-like scaling only at small F and diverging as the jamming point is approached; if the shear stress is strong enough, at higher phi an ordered flow regime is found, the order-disorder transition being marked by a sharp drop of the viscosity. A broad jamming region is also observed where, in analogy with the glassy region of thermal systems, slow dynamics followed by kinetic arrest occurs when the ordering transition is prevented.     

Jamming: A simple introduction

Amorphous materials as diverse as foams, emulsions, colloidal suspensions and granular media can jam into a rigid, disordered state where they withstand finite shear stresses before yielding. Here we give a simple introduction to the surprising physics displayed by a very simple model system for the jamming transition: frictionless, soft spheres at zero temperature and zero shear that act through purely repulsive contact forces. This system starts to become rigid, i.e. goes through the jamming transition, whenever the confining pressure becomes positive. We highlight some of the remarkable geometrical features of the zero pressure jamming point and discuss the peculiar mechanical properties of these systems for small pressures.     

Finite frequency quantum noise in an interacting mesoscopic conductor

We present a quantum calculation based on scattering theory of the frequency-dependent noise of current in an interacting chaotic cavity. We include interactions of the electron system via long range Coulomb forces between the conductor and a gate with capacitance C. We obtain explicit results exhibiting the two time scales of the problem, the cavity's dwell time tau(D) and the RC time tau(C) of the cavity in relation to the gate. The noise shows peculiarities at frequencies of the order and exceeding the inverse charge relaxation time tau(-1) = tau(D)(-1) + tau(C)(-1).     

Probing weak forces in granular media through nonlinear dynamic dilatancy: clapping contacts and polarization anisotropy

Rectification (demodulation) of high-frequency shear acoustic bursts is applied to probe the distribution of contact forces in 3D granular media. Symmetry principles allow for rectification of the shear waves only with their conversion into longitudinal mode. The rectification is due to nonlinear dynamic dilatancy, which is found to follow a quadratic or Hertzian power law in the shear wave amplitude. Evidence is given that a significant portion of weak contact forces is localized below 10(-2) of the mean force-a range previously inaccessible by experiment. Strong anisotropy of nonlinearity for shear waves with different polarization is observed.     

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.     

A column of grains in the jamming limit: glassy dynamics in the compaction process

We investigate a stochastic model describing a column of grains in the jamming limit, in the presence of a low vibrational intensity. The key control parameter of the model, , is a representation of granular shape, related to the reduced void space. Regularity and irregularity in grain shapes, respectively corresponding to rational and irrational values of , are shown to be centrally important in determining the statics and dynamics of the compaction process.Received: 12 May 2003, Published online: 15 October 2003PACS: 45.70.-n Granular systems - 45.70.Cc Static sandpiles; granular compaction - 45.70.Mg Granular flow: mixing, segregation and stratification     

Jamming and flow of random-close-packed spherical bubbles: an analogy with granular materials

Wet foams are random-close-packed assemblies of approximately spherical gas bubbles in a liquid. We report rheological experiments with this material, showing that even though the stiffness and frictional interactions of bubbles strongly distinguish them from solid spherical grains, jamming and flow of wet foams and granular materials are governed by closely analogous laws.     

Granular flow from a silo: Discrete-particle simulations in three dimensions

Molecular (or granular) dynamics methods are used to study the gravity-driven flow of granular material through a horizontal aperture in three dimensions. The grains are spherical and modeled using a short-range repulsive interaction, together with normal and tangential frictional damping forces. The material is contained in a rough-walled cylindrical container with a circular hole in its base, and to permit flow measurements under steady-state conditions a continuous feed approach is employed in which exiting grains are replaced at the upper surface of the material. The dependence of flow velocity and discharge rate on aperture diameter is found to agree with experiment; other quantities such as the kinetic energy and pressure distributions are also examined. Received 5 June 2000 and Received in final form 21 September 2000     

Scaling laws for melting ice avalanches

This Letter describes an investigation of interfacial melting in ice-bearing granular flows. It is proposed that energy associated with granular collisions causes melting at an ice particle's surface, which can thus occur at temperatures well below freezing. A laboratory experiment has been designed that allows quantification of this process and its effect on the dynamics of a granular shear flow of ice spheres. This experiment employs a rotating drum, half filled with ice particles, situated in a temperature controlled laboratory. Capillary forces between the wetted melted particle surfaces lead to the clumping of particles and enhanced flow speeds, in turn leading to further melting. Dimensional analysis defines a parameter space for further experimentation.     

离散颗粒层被横向推移过程中的力学行为研究

研究发现,颗粒物质层被匀速推移挤压过程中,所需推移力先以线性规律增加,在某一确定点后,则会以指数规律增加.而颗粒物质是由众多离散颗粒组成的软凝聚态物质,其宏观上反映的是离散颗粒的个体性质和凝聚态物质的集体效应.颗粒与颗粒之间以及颗粒与边界之间的细观尺度接触力链的构成以及演变规律将会直接影响各种宏观受力情况,其摩擦力与挤压力便是力链的主要构成形式.围绕着定量描述细观力链特征,从而揭示力的变化规律这一目标,采用计算机模拟的方法,依照球形颗粒Hertz法向接触理论和Mindlin-Deresiewicz切向接触理论,对重力作用下不同数目的三维等径球体颗粒层的推移情况进行了离散元仿真模拟,量化分析了推移力变化规律、各摩擦力变化规律以及力链分布规律,发现摩擦力与挤压力在颗粒堆积的不同阶段对力链的构成起到了不同的主次作用,使得力链发生强弱演变,从而发现了推移颗粒物质层时推移力的变化规律及原因.这些结果有利于从力链角度揭示颗粒内部和颗粒与各边界之间的受力情况.     

Global nature of dilute-to-dense transition of granular flows in a 2D channel

The dilute-to-dense transition of granular flow of particle size d(0) is studied experimentally in a two-dimensional channel (width D) with confined exit (width d). Our results show that with fixed d and D there is a maximum inflow rate Q(c) above which the flow changes from dilute to dense and the outflow rate drops abruptly from Q(c) to a dense rate Q(d). A rescaled critical rate q(c) is found to be a function of a scaling variable lambda only: q(c) approximately F(lambda), where lambda identical with d/d(0) d/D-d. This form of lambda suggests that the dilute-to-dense transition is a global property of the flow, unlike the jamming transition which depends only on d/d(0). Furthermore, the transition is found to occur when the area fraction of particles near the exit exceeds a critical value which is close to 0.65+/-0.03.     

Repulsive contact interactions make jammed particulate systems inherently nonharmonic

Many jammed particulate systems, such as granular and colloidal materials, interact via repulsive contact forces. We find that these systems possess no harmonic regime in the large system limit (N→∞) for all compressions Δ? studied, and at jamming onset Δ?→0 for all N. We perform fixed energy simulations following perturbations with amplitude δ along eigendirections of the dynamical matrix. The fluctuations abruptly spread to all modes for δ≈δ(c) (where a single contact breaks) in contrast to linear and weakly nonlinear behavior. For δ > δ(c), all discrete modes disappear into a continuous frequency band. <δ(c)> scales with 1/N and Δ?, which limits harmonic behavior to only overcompressed systems. The density of vibrational modes deviates strongly from that predicted from the dynamical matrix when the system enters the nonharmonic regime, which significantly affects its mechanical and transport properties.