Skip to content
PaperFren

Granular matter

How do force chains in flowing sand become rigid?

Dashti H, Saberi AA, Rahbari SHE, et al. · Science advances · 2023

Open access · cc by · source: Europe PMC

In slowly sheared granular material, the network of strong contact forces becomes connected in exactly the way predicted by rigidity percolation theory, but faster flow changes the critical behaviour.

Study at a glance

Design
Computational / modelling — Molecular dynamics (LAMMPS) of two-dimensional frictionless bidisperse disks under simple shear, with percolation analysis of force networks and finite-size scaling across system sizes and shear rates.
N
No participant N; system sizes range from 2048 to 65,536 disks, each averaged over 10 independent simulations at each of several shear rates.
Population
Simulated dense (above jamming) packings of frictionless bidisperse disks in two dimensions
Outcome
Critical exponents (nu, beta, gamma, eta) and fractal dimension of the percolating interparticle-force network as a function of shear rate

Structured fields used in claim comparison tables when every cited study has a complete layer.

Key findings

At the slowest shear rates the exponents (for example nu of about 1.21, beta of about 0.20 and gamma of about 2.15) matched standard rigidity percolation, which the authors describe as the first verification of this class in an off-lattice molecular dynamics simulation. Above a crossover shear rate the correlation-length exponent rose steadily, so flow rate acts as a relevant perturbation producing a continuous line of exponents, passing through random-percolation values at intermediate rates. The force-force correlation exponent and the fractal dimension of the spanning cluster, however, stayed essentially the same across all flow rates.

Methodology

The authors simulated a dense two-dimensional packing of frictionless disks of two sizes being continuously sheared, using system sizes from 2048 to 65,536 particles and shear rates spanning three orders of magnitude. For each configuration they kept only contacts whose force exceeded a threshold and asked at which threshold a system-spanning cluster of force chains appears. Using finite-size scaling, they extracted the critical exponents of this percolation transition and compared them with known universality classes, and interpreted the shear-rate dependence with the extended Harris criterion for correlated disorder.

Limitations

The system is idealised: two-dimensional, frictionless, athermal disks with simple linear contact forces, so the results may not carry over directly to frictional sand, three-dimensional flows or thermal glasses, which the authors flag as open questions. Percolation is defined by an arbitrary force threshold on contacts rather than by a direct rigidity (pebble-game) test of each cluster. The explanation of changing exponents via long-range flow-induced correlations is a theoretical argument fitted to the data, and the beta exponent had error bars too large for a conclusive test.

How this study connects

Role on claims

Each row is a claim on a concept or method page where this paper supports, challenges, or qualifies the statement. Roles are hand-checked — not a model guess.

  • How fast a granular packing flows changes the critical behaviour of how rigidity spreads through it.

    In 2D simulations of frictionless sheared disks, force networks at the slowest shear rates showed standard rigidity-percolation exponents (for example nu about 1.21), but above a crossover shear rate the correlation-length exponent rose steadily, while the fractal dimension of the spanning cluster stayed essentially fixed.

    Evidence for the claim as stated.

  • How fast a granular packing flows changes the critical behaviour of how rigidity spreads through it.

    In 2D simulations of frictionless sheared disks, force networks at the slowest shear rates showed standard rigidity-percolation exponents (for example nu about 1.21), but above a crossover shear rate the correlation-length exponent rose steadily, while the fractal dimension of the spanning cluster stayed essentially fixed.

    Scope note — Idealised 2D, frictionless, athermal simulation; percolation defined with a force threshold.

    Limits the claim's scope: a different population, assay, or outcome.

  • Rigidity-percolation results come from frictionless 2D disks, while the intrusion, asteroid and jamming studies involve frictional, cohesive or 3D grains; the authors of the percolation study flag friction and 3D as open questions, so the exponents should not be assumed to apply to real sand.

    Evidence for the claim as stated.

Open questions

Tensions this paper is part of

From concept pages' “where studies disagree.” Disagreement means the same question; scope means different assays, populations, or outcomes.

Discoveries this paper informs or conflicts with

Related papers in this topic

Same topic cluster — not a recommendation engine.