Statistical Properties of Granular Materials near Jamming ESMC

  • Slides: 49
Download presentation
Statistical Properties of Granular Materials near Jamming ESMC 2009, Lisbon September 8, 2009 R.

Statistical Properties of Granular Materials near Jamming ESMC 2009, Lisbon September 8, 2009 R. P. Behringer Duke University Collaborators: Max Bi, Chris Bizon, Karen Daniels, Julien Dervaux, Somayeh Farhadi, Junfei Geng, Bob Hartley, Silke Henkes, Dan Howell, Trush Majmudar, Guillaume Reydellet, Trevor Shannon, Matthias Sperl, Junyao Tang, Sarath Tennakoon, Brian Tighe, John Wambaugh, Brian Utter, Peidong Yu, Jie Ren, Jie Zhang, Bulbul Chakraborty, Eric Clément, Isaac Goldhirsch, Lou Kondic, Stefan Luding, Guy Metcalfe, Corey O’Hern, David Schaeffer, Josh Socolar, Antoinette Tordesillas Support: NSF, ARO, NASA, IFPRI, BSF

Roadmap • What/Why granular materials? • Behavior of disordered solids—possible universal behavior? • Where

Roadmap • What/Why granular materials? • Behavior of disordered solids—possible universal behavior? • Where granular materials and molecular matter part company—open questions of relevant scales Experiments at Duke explore: • Forces, force fluctuations • Isotropic jamming • Effect of shear—anisotropic stresses

What are Granular Materials? • Collections of macroscopic ‘hard’ particles: interactions are dissipative/frictional –

What are Granular Materials? • Collections of macroscopic ‘hard’ particles: interactions are dissipative/frictional – – – – Classical h 0 Highly dissipative Draw energy for fluctuations from macroscopic flow Large collective systems, but outside normal statistical physics Although many-body, a-thermal in the usual sense Exist in phases: granular gases, fluids and solids Analogues to other disordered solids: glasses, colloids. . May be dry or wet

Broader context: Disordered N-body systems— far from equilibrium • There exist a number of

Broader context: Disordered N-body systems— far from equilibrium • There exist a number of such systems: glasses, foams, colloidal suspensions, granular materials, … • For various reasons, these systems are not in ordinary thermal equilibrium (although temperature may still play a role) • Energy may not be conserved • Other conservation rules—e. g. stress, may apply

Standard picture of jamming • Jamming—how disordered N-body systems becomes solid-like as particles are

Standard picture of jamming • Jamming—how disordered N-body systems becomes solid-like as particles are brought into contact, or fluid-like when grains are separated • Density is implicated as a key parameter, expressed as packing (solid fraction) φ • Marginal stability (isostaticity) for spherical particles (disks in 2 D) contact number, Z, attains a critical value, Ziso at φiso

Jamming How do disordered solids lose/gain their solidity? Bouchaud et al. Liu and Nagel

Jamming How do disordered solids lose/gain their solidity? Bouchaud et al. Liu and Nagel

Theoretical tools: Statistical ensembles: what to do when energy is not conserved? • Edwards

Theoretical tools: Statistical ensembles: what to do when energy is not conserved? • Edwards ensemble for rigid particles: V replaces E • Real particles are deformable (elastic): forces, stresses, and torques matter, and are ‘conserved’ for static systems—hence force/stress ensembles emerge • Snoeijer, van Saarloos et al. Tighe, Socolar et al. , Henkes and Chakraborty and O’Hern, Makse et al.

Experimental tools: what to measure, and how to look inside complex systems • Confocal

Experimental tools: what to measure, and how to look inside complex systems • Confocal techniques in 3 D—with fluidsuspended particles—for colloids, emulsions, fluidized granular systems • Bulk measurements— 2 D and 3 D • Measurements at boundaries— 3 D • 2 D measurements: particle tracking, Photoelastic techniques (this talk) • Numerical experiments—MD/DEM

What happens when shear is applied? Plasticity—irreversible deformation when a material is sheared •

What happens when shear is applied? Plasticity—irreversible deformation when a material is sheared • System becomes anisotropic—e. g. long force chains form • Shear causes irreversible (plastic) deformation. Particles move ‘around’ each other • What is the microscopic nature of this process for granular materials?

Different types methods of applying shear • Example 1: pure shear • Example 2:

Different types methods of applying shear • Example 1: pure shear • Example 2: simple shear • Example 3: steady shear

Dense Granular Material Phases-Some simple observations Forces are carried preferentially on force chains multiscale

Dense Granular Material Phases-Some simple observations Forces are carried preferentially on force chains multiscale phenomena Friction matters Howell et al. PRL 82, 5241 (1999)

Roadmap • What/Why granular materials? • Where granular materials and molecular matter part company—open

Roadmap • What/Why granular materials? • Where granular materials and molecular matter part company—open questions of relevant scales • Dense granular materials: need statistical approach Use experiments to explore: • Forces, force fluctuations ◄ • Jamming—isotropic • Shear and anisotropic stresses

Experiments to determine vector contact forces, distribution P 1(F) is example of particle-scale statistical

Experiments to determine vector contact forces, distribution P 1(F) is example of particle-scale statistical measure Experiments use biaxial tester and photoelastic particles (Trush Majmudar and RPB, Nature, June 23, 2005)

Overview of Experiments Biax schematic Compression Shear Image of Single disk ~2500 particles, bi-disperse,

Overview of Experiments Biax schematic Compression Shear Image of Single disk ~2500 particles, bi-disperse, d. L=0. 9 cm, d. S= 0. 8 cm, NS /NL = 4

Measuring forces by photoelasticity

Measuring forces by photoelasticity

Basic principles of technique • Process images to obtain particle centers and contacts •

Basic principles of technique • Process images to obtain particle centers and contacts • Interparticle contact forces determine stresses within each particle, including principal stresses, σi • Stresses determine photoelastic response: I = Iosin 2[(σ2 - σ1)CT/λ] • Now go backwards, using nonlinear inverse technique to obtain contact forces • In the previous step, invoke force and torque balance • Newton’s 3 d law provides error checking

Examples of Experimental and ‘Fitted’ Images Experiment—color Filtered Experiment--original Fitted

Examples of Experimental and ‘Fitted’ Images Experiment—color Filtered Experiment--original Fitted

Force distributionsfor shear and compression Shear εxx = -εyy =0. 04; Zavg = 3.

Force distributionsfor shear and compression Shear εxx = -εyy =0. 04; Zavg = 3. 1 Compression εxx = -εyy =0. 016; Zavg = 3. 7 From T. Majmudar and RPB, Nature, 2005

Stress ensemble models for P(f) • Consider all possible states consistent with applied external

Stress ensemble models for P(f) • Consider all possible states consistent with applied external forces, or other boundary conditions— assume all possible states occur with equal probability • Compute Fraction where at least one contact force has value f P(f) • E. g. Snoeier et al. PRL 92, 054302 (2004) • Tighe et al. Phys. Rev. E, 72, 031306 (2005)

Some Typical Cases—isotropic compression and shear Snoeijer et al. ↓ Tigue et al ↓.

Some Typical Cases—isotropic compression and shear Snoeijer et al. ↓ Tigue et al ↓. Compression Shear Latest update—B. Tigue, this session

What about force correlations? Compression Shear

What about force correlations? Compression Shear

Correlation functions determine important scales • C(r) = <Q(r + r’) Q(r’)> • <>

Correlation functions determine important scales • C(r) = <Q(r + r’) Q(r’)> • <> average over all vector displacements r’ • For isotropic cases, average over all directions in r. • Angular averages should not be done for anisotropic systems

Spatial correlations of forces—angle dependent Shear Chain direction Compression Both directions equivalent Direction normal

Spatial correlations of forces—angle dependent Shear Chain direction Compression Both directions equivalent Direction normal To chains New work: S. Henkes, B. Chakraborty, G. Lois, C. O’Hern, J. Zhang, RPB

Roadmap • What/Why granular materials? • Where granular materials and molecular matter part company—open

Roadmap • What/Why granular materials? • Where granular materials and molecular matter part company—open questions of relevant scales • Dense granular materials: need statistical approach Use experiments to explore: • Forces, force fluctuations ◄ • Jamming--isotropic ◄ • Shear—anisotropic stresses

Jamming How do disordered solids lose/gain their solidity? Bouchaud et al. Liu and Nagel

Jamming How do disordered solids lose/gain their solidity? Bouchaud et al. Liu and Nagel

The Isotropic Jamming Transition—Point J • Simple question: What happens to key properties such

The Isotropic Jamming Transition—Point J • Simple question: What happens to key properties such as pressure, contact number as a sample is isotropically compressed/dilated through the point of mechanical stability? Z = contacts/particle; Φ = packing fraction Predictions for spherical frictionless particles (e. g. O’Hern et al. Torquato et al. , Schwarz et al. , Wyart et al. Z ~ ZI +(φ – φc)ά (discontinuity) Exponent ά ≈ 1/2 P ~(φ – φc)β β depends on force law (= 1 for ideal disks) S. Henkes and B. Chakraborty: entropy-based model gives P and Z in terms of a field conjugate to entropy. Can eliminate to get P(z)

Experiment: Characterizing the Jamming Transition —Isotropic compression ßIsotropic compression Majmudar et al. PRL 98,

Experiment: Characterizing the Jamming Transition —Isotropic compression ßIsotropic compression Majmudar et al. PRL 98, 058001 (2007)

How do we obtain stresses and Z? Fabric tensor Rij = Sk, c ncik

How do we obtain stresses and Z? Fabric tensor Rij = Sk, c ncik ncjk Z = trace[R] Stress tensor : sij = (1/A) Sk, c rcik fcjk A is system area, trace of stress tensor gives P

LSQ Fits for Z give an exponent of 0. 5 to 0. 6

LSQ Fits for Z give an exponent of 0. 5 to 0. 6

LSQ Fits for P give β ≈ 1. 0 to 1. 1

LSQ Fits for P give β ≈ 1. 0 to 1. 1

Comparison to Henkes and Chakraborty prediction

Comparison to Henkes and Chakraborty prediction

Roadmap • What/Why granular materials? • Where granular materials and molecular matter part company—open

Roadmap • What/Why granular materials? • Where granular materials and molecular matter part company—open questions of relevant scales • Dense granular materials: need statistical approach Use experiments to explore: • Forces, force fluctuations ◄ • Isotropic jamming ◄ • Shear—anisotropic stresses◄

A different part of the jamming diagram σ2 σ1 What happens here or here,

A different part of the jamming diagram σ2 σ1 What happens here or here, when shear is applied to a granular material? Note: P = (s 2 + s 1)/2 Coulomb failure: : t = (s 2 – s 1)/2 |t|/P = m

Shear near jamming • Example 1: pure shear • Example 2: simple shear •

Shear near jamming • Example 1: pure shear • Example 2: simple shear • Example 3: steady shear

Starting point is anisotropic, unjammed state Isotroic jamming point at φ = 0. 82

Starting point is anisotropic, unjammed state Isotroic jamming point at φ = 0. 82 Start here, with φ = 0. 76

What happens for granular materials subject to pure shear? Use biax and photoelastic particles

What happens for granular materials subject to pure shear? Use biax and photoelastic particles Mark particles with UV-sensitive Dye for tracking J. Zhang et al. to appear, Granular Matter

Apply Cyclic Pure Shear—starting from an unjammed state Resulting state with polarizer And without

Apply Cyclic Pure Shear—starting from an unjammed state Resulting state with polarizer And without polarizer

Consider one cycle of shear

Consider one cycle of shear

Particle Displacements and Rotations Forward shear—under UV

Particle Displacements and Rotations Forward shear—under UV

Deformation Field—Shear band forms At strain = 0. 085 At strain = 0. 111

Deformation Field—Shear band forms At strain = 0. 085 At strain = 0. 111 At strain = 0. 105—largest plastic event

Particle displacement and rotation (forward shear) • Green arrows are displace- ment of particle

Particle displacement and rotation (forward shear) • Green arrows are displace- ment of particle center • Blob size stands for rotation magnitude • Blue color—clockwise rotation • Brown color— counterclockwise rotation • Mean displacement subtracted

Stresses and Z

Stresses and Z

Hysteresis in stress-strain and Z-strain curves

Hysteresis in stress-strain and Z-strain curves

But!! Apparent scaling for stresses vs. Z

But!! Apparent scaling for stresses vs. Z

Force Distributions Tangential forces Normal forces

Force Distributions Tangential forces Normal forces

Corrections for missed contacts Estimate of missed contacts: ∫ 0 F-min P(F)d. F/∫ 0∞

Corrections for missed contacts Estimate of missed contacts: ∫ 0 F-min P(F)d. F/∫ 0∞ P(F)d. F Similar approach corrects P, but Effect is much smaller, since P ~ ∫ F P(F) d. F

How do we contemplate jamming in frictional granular materials? σ2 Sheared granular materials fail

How do we contemplate jamming in frictional granular materials? σ2 Sheared granular materials fail to other stable states (σ2 - σ1) J ? ? 1/φ ? ? σ1 Note that Reynolds Dilatancy weakly confined samples dilate under shear— Hence, rigidly confined materials show an iσncrease of P under shear

Conclusions • Granular materials show important features due to friction • Distributions of forces

Conclusions • Granular materials show important features due to friction • Distributions of forces show sensitivity to stress state but agree reasonably with stress-ensemble approach • Correlations forces in sheared systems—thus, force chains can be mesoscopic • Predictions for jamming (mostly) verified • Granular states near jamming show jamming under shear at densities below isotropic jamming • Friction (and also preparation history) implicated • Z may be key variable for shear failure/jamming • Generalized SGR explains rate-dependence

What is actual force law for our disks?

What is actual force law for our disks?