Critical scaling near the yielding transition in granular media

We show that the yielding transition in granular media displays second-order critical-point scaling behavior. We carry out discrete element simulations in the low-inertial-number limit for frictionless, purely repulsive spherical grains undergoing simple shear at fixed nondimensional shear stress Σ in two and three spatial dimensions. To find a mechanically stable (MS) packing that can support the applied Σ, isotropically prepared states with size L must undergo a total strain γ_{ms}(Σ,L). The number density of MS packings (∝γ_{ms}^{-1}) vanishes for Σ>Σ_{c}≈0.11 according to a critical scaling form with a length scale ξ∝|Σ-Σ_{c}|^{-ν}, where ν≈1.7-1.8. Above the yield stress (Σ>Σ_{c}), no MS packings that can support Σ exist in the large-system limit L/ξ≫1. MS packings generated via shear possess anisotropic force and contact networks, suggesting that Σ_{c} is associated with an upper limit in the degree to which these networks can be deformed away from those for isotropic packings.