Effect of inertia on sheared disordered solids: critical scaling of avalanches in two and three dimensions

Effects of Inertia on the Steady-Shear Rheology of Disordered Solids

We study the finite-shear-rate rheology of disordered solids by means of molecular dynamics simulations in two dimensions. By systematically varying the damping strength ζ in the low-temperature limit, we identify two well-defined flow regimes, separated by a thin (temperature-dependent) crossover region. In the overdamped regime, the athermal rheology is governed by the competition between elastic forces and viscous forces, whose ratio gives the Weissenberg number Wi∝ζγ[over ˙]; the macroscopic stress Σ follows the frequently encountered Herschel-Bulkley law Σ=Σ_{0}+ksqrt[Wi], with yield stress Σ_{0}>0. In the underdamped (inertial) regime, dramatic changes in the rheology are observed for low damping: the flow curve becomes nonmonotonic. This change is not caused by longer-lived correlations in the particle dynamics at lower damping; instead, for weak dissipation, the sample heats up considerably due to, and in proportion to, the driving. By thermostating more or less underdamped systems, we are able to link quantitatively the rheology to the kinetic temperature and the shear rate, rescaled with Einstein's vibration frequency.

Role of inertia in the rheology of amorphous systems: A finite-element-based elastoplastic model

A simple finite-element analysis with varying damping strength is used to model the athermal shear rheology of densely packed glassy systems at a continuum level. We focus on the influence of dissipation on bulk rheological properties. Our numerical studies, done over a wide range of damping coefficients, identify two well-separated rheological regimes along with a crossover region controlled by a critical damping. In the overdamped limit, inertial effects are negligible and the rheological response is well described by the commonly observed Herschel-Bulkley equation. In stark contrast, inertial vibrations in the underdamped regime prompt a significant drop in the mean-stress level, leading to a nonmonotonic constitutive relation. The observed negative slope in the flow curve, which is a signature of mechanical instability and thus permanent shear banding, arises from the sole influence of inertia, in qualitative agreement with the recent molecular dynamics study of Nicolas et al., Phys. Rev. Lett. 116, 058303 (2016).

Cooperativity flows and shear-bandings: a statistical field theory approach

Cooperativity effects have been proposed to explain the non-local rheology in the dynamics of soft jammed systems. Based on the analysis of the free-energy model proposed by L. Bocquet, A. Colin and A. Ajdari, Phys. Rev. Lett., 2009, 103, 036001, we show that cooperativity effects resulting from the non-local nature of the fluidity (inverse viscosity) are intimately related to the emergence of shear-banding configurations. This connection materializes through the onset of inhomogeneous compact solutions (compactons), wherein the fluidity is confined to finite-support subregions of the flow and strictly zero elsewhere. The compacton coexistence with regions of zero fluidity ("non-flowing vacuum") is shown to be stabilized by the presence of mechanical noise, which ultimately shapes up the equilibrium distribution of the fluidity field, the latter acting as an order parameter for the flow-noflow transitions occurring in the material.

The effect of a reversible shear transformation on plastic deformation of an amorphous solid

Langevin dynamics simulations are performed to investigate the plastic response of a model glass to a local shear transformation in a quiescent system. The deformation of the material is induced by a spherical inclusion that is gradually strained into an ellipsoid of the same volume and then reverted back into the sphere. We show that the number of cage-breaking events increases with increasing strain amplitude of the shear transformation. The results of numerical simulations indicate that the density of cage jumps is larger in the cases of weak damping or slow shear transformation. Remarkably, we also found that, for a given strain amplitude, the peak value of the density profiles is a function of the ratio of the damping coefficient and the time scale of the shear transformation.

Criticality in the Approach to Failure in Amorphous Solids

Failure of amorphous solids is fundamental to various phenomena, including landslides and earthquakes. Recent experiments indicate that highly plastic regions form elongated structures that are especially apparent near the maximal shear stress Σmax where failure occurs. This observation suggested that Σmax acts as a critical point where the length scale of those structures diverges, possibly causing macroscopic transient shear bands. Here, we argue instead that the entire solid phase (Σ<Σmax) is critical, that plasticity always involves system-spanning events, and that their magnitude diverges at Σmax independently of the presence of shear bands. We relate the statistics and fractal properties of these rearrangements to an exponent θ that captures the stability of the material, which is observed to vary continuously with stress, and we confirm our predictions in elastoplastic models.

Probing relevant ingredients in mean-field approaches for the athermal rheology of yield stress materials

Although the notion of mechanical noise is expected to play a key role in the non-linear rheology of athermally sheared amorphous systems, its characterization has so far remained elusive. Here, we show using molecular dynamic simulations that in spite of the presence of strong spatio-temporal correlations in the system, the local stress exhibits normal diffusion under the effect of the mechanical noise in the finite driving regime. The diffusion constant appears to be proportional to the mean plastic activity. Our data suggests that the corresponding proportionality constant is density independent, and can be directly related to the specific form of the rheological flow curve, pointing the way to a generic way of modeling mechanical noise in mean-field equations.

Shear-transformation-zone theory of yielding in athermal amorphous materials

Yielding transitions in athermal amorphous materials undergoing steady-state shear flow resemble critical phenomena. Historically, they have been described by the Herschel-Bulkley rheological formula, which implies singular behaviors at yield points. In this paper, I examine this class of phenomena using an elementary version of the thermodynamic shear-transformation-zone (STZ) theory, focusing on the role of the effective disorder temperature, and paying special attention to scaling and dimensional arguments. I find a wide variety of Herschel-Bulkley-like rheologies but, for fundamental reasons not specific to the STZ theory, conclude that the yielding transition is not truly critical. In particular, for realistic many-body models with short-range interactions, there is a correlation length that grows rapidly but ultimately saturates near the yield point.

Density scaling and quasiuniversality of flow-event statistics for athermal plastic flows

Athermal steady-state plastic flows were simulated for the Kob-Andersen binary Lennard-Jones system and its repulsive version in which the sign of the attractive terms is changed to a plus. Properties evaluated include the distributions of energy drops, stress drops, and strain intervals between the flow events. We show that simulations at a single density in conjunction with an equilibrium-liquid simulation at the same density allow one to predict the plastic flow-event statistics at other densities. This is done by applying the recently established "hidden scale invariance" of simple liquids to the glass phase. The resulting scaling of flow-event properties reveals quasiuniversality, i.e., that the probability distributions of energy drops, stress drops, and strain intervals in properly reduced units are virtually independent of the microscopic pair potentials.

Scaling description of the yielding transition in soft amorphous solids at zero temperature

Yield stress materials flow if a sufficiently large shear stress is applied. Although such materials are ubiquitous and relevant for industry, there is no accepted microscopic description of how they yield, even in the simplest situations in which temperature is negligible and in which flow inhomogeneities such as shear bands or fractures are absent. Here we propose a scaling description of the yielding transition in amorphous solids made of soft particles at zero temperature. Our description makes a connection between the Herschel-Bulkley exponent characterizing the singularity of the flow curve near the yield stress Σc, the extension and duration of the avalanches of plasticity observed at threshold, and the density P(x) of soft spots, or shear transformation zones, as a function of the stress increment x beyond which they yield. We argue that the critical exponents of the yielding transition may be expressed in terms of three independent exponents, θ, df, and z, characterizing, respectively, the density of soft spots, the fractal dimension of the avalanches, and their duration. Our description shares some similarity with the depinning transition that occurs when an elastic manifold is driven through a random potential, but also presents some striking differences. We test our arguments in an elasto-plastic model, an automaton model similar to those used in depinning, but with a different interaction kernel, and find satisfying agreement with our predictions in both two and three dimensions.