Rheology and dynamical heterogeneity in frictionless beads at jamming density

Scale Separation of Shear-Induced Criticality in Glasses

In a sheared steady state, glasses reach a nonequilibrium criticality called yielding criticality. We report that the qualitative nature of this nonequilibrium critical phenomenon depends on the details of the system and that responses and fluctuations are governed by different critical correlation lengths in specific situations. This scale separation of critical lengths arises when the screening of elastic propagation of mechanical signals is not negligible. We also discuss the determinant of the impact of screening effects from the viewpoint of the microscopic dissipation mechanism.

Relaxation times, rheology, and finite size effects for non-Brownian disks in two dimensions

We carry out overdamped simulations in a simple model of jamming-a collection of bidisperse soft core frictionless disks in two dimensions-with the aim to explore the finite size dependence of different quantities, both the relaxation time obtained from the relaxation of the energy and the pressure equivalent of the shear viscosity. The motivation for the paper is the observation [Nishikawa et al., J. Stat. Phys. 182, 37 (2021)0022-471510.1007/s10955-021-02710-8] that there are finite size effects in the relaxation time, τ, that give problems in the determination of the critical divergence, and the claim that this is due to a finite size dependence, τ∼lnN, which makes τ an ill-defined quantity. Beside analyses to determine the relaxation time for the whole system we determine particle relaxation times which allow us to determine both histograms of particle relaxation times and the average particle relaxation times-two quantities that are very useful for the analyses. The starting configurations for the relaxation simulations are of two different kinds-completely random or taken from steady shearing simulations-and we find that the difference between these two cases are bigger than previously noted and that the observed problems in the determination of the critical divergence obtained when starting from random configurations are not present when instead starting the relaxations from shearing configurations. We also argue that the the effect that causes the lnN dependence is not as problematic as asserted. When it comes to the finite size dependence of the pressure equivalent of the shear viscosity we find that our data don't give support for the claimed strong finite size dependence, but also that the finite size dependence is at odds with what one would normally expect for a system with a diverging correlation length, and that this calls for an alternative understanding of the phenomenon of shear-driven jamming.

Instantaneous Normal Modes Reveal Structural Signatures for the Herschel-Bulkley Rheology in Sheared Glasses

The Herschel-Bulkley law, a universal constitutive relation, has been empirically known to be applicable to a vast range of soft materials, including sheared glasses. Although the Herschel-Bulkley law has attracted public attention, its structural origin has remained an open question. In this Letter, by means of atomistic simulation of binary Lennard-Jones glasses, we report that the instantaneous normal modes with negative eigenvalues, or so-called imaginary modes, serve as the structural signatures for the Herschel-Bulkley rheology in sheared glasses.

Translational and rotational velocities in shear-driven jamming of ellipsoidal particles

We study shear-driven jamming of ellipsoidal particles at zero temperature with a focus on the microscopic dynamics. We find that a change from spherical particles to ellipsoids with aspect ratio α=1.02 gives dramatic changes of the microscopic dynamics with much lower translational velocities and a new role for the rotations. Whereas the velocity difference at contacts-and thereby the dissipation-in collections of spheres is dominated by the translational velocities and reduced by the rotations, the same quantity is in collections of ellipsoids instead totally dominated by the rotational velocities. By also examining the effect of different aspect ratios we find that the examined quantities show either a peak or a change in slope at α≈1.2, which thus gives evidence for a crossover between different regions of low and high aspect ratio.

Avalanche Interpretation of the Power-Law Energy Spectrum in Three-Dimensional Dense Granular Flow

Turbulence is ubiquitous in nonequilibrium systems, and it has been noted that even dense granular flows exhibit characteristics that are typical of turbulent flow, such as the power-law energy spectrum. However, studies on the turbulentlike behavior of granular flows are limited to two-dimensional (2D) flow. We demonstrate that the statistics in three-dimensional (3D) flow are qualitatively different from those in 2D flow. We also elucidate that avalanche dynamics can explain this dimensionality dependence. Moreover, we define clusters of collectively moving particles that are equivalent to vortex filaments. The clusters unveil complicated structures in 3D flows that are absent in 2D flows.

Effect of collisional elasticity on the Bagnold rheology of sheared frictionless two-dimensional disks

We carry out constant volume simulations of steady-state, shear-driven flow in a simple model of athermal, bidisperse, soft-core, frictionless disks in two dimensions, using a dissipation law that gives rise to Bagnoldian rheology. Focusing on the small strain rate limit, we map out the rheological behavior as a function of particle packing fraction ϕ and a parameter Q that measures the elasticity of binary particle collisions. We find a Q^{*}(ϕ) that marks the clear crossover from a region characteristic of strongly inelastic collisions, Q<Q^{*}, to a region characteristic of weakly inelastic collisions, Q>Q^{*}, and give evidence that Q^{*}(ϕ) diverges as ϕ→ϕ_{J}, the shear-driven jamming transition. We thus conclude that the jamming transition at any value of Q behaves the same as the strongly inelastic case, provided one is sufficiently close to ϕ_{J}. We further characterize the differing nature of collisions in the strongly inelastic vs weakly inelastic regions, and recast our results into the constitutive equation form commonly used in discussions of hard granular matter.

Enstrophy cascades in two-dimensional dense granular flows

Employing two-dimensional molecular dynamics simulations of dense granular materials under simple shear deformations, we investigate vortex structures of particle rearrangements. Introducing vorticity fields as a measure of local spinning motions of the particles, we observe their heterogeneous distributions, where statistics of vorticity fields exhibit the highly non-Gaussian behavior and typical domain sizes of vorticity fields significantly increase if the system is yielding under quasistatic deformations. In such dense granular flows, a power-law decay of vorticity spectra can be observed at mesoscopic scale, implying anomalous local structures of kinetic energy dissipation. We explain the power-law decay, or enstrophy cascades in dense granular materials, by a dimensional analysis, where the dependence of vorticity spectra not only on the wave number, but also on the shear rate, is well explained. From our dimensional analyses, the scaling of granular temperature and rotational kinetic energy is also predicted.

Critical scaling of Bagnold rheology at the jamming transition of frictionless two-dimensional disks

We carry out constant volume simulations of steady-state shear-driven rheology in a simple model of bidisperse soft-core frictionless disks in two dimensions, using a dissipation law that gives rise to Bagnoldian rheology. We discuss in detail the critical scaling ansatz for the shear-driven jamming transition and carry out a detailed scaling analysis of our resulting data for pressure p and shear stress σ. Our analysis determines the critical exponent β that describes the algebraic divergence of the Bagnold transport coefficients lim_{γ[over ̇]→0}p/γ[over ̇]^{2},σ/γ[over ̇]^{2}∼(ϕ_{J}-ϕ)^{-β} as the jamming transition ϕ_{J} is approached from below. For the low strain rates considered in this work, we show that it is still necessary to consider the leading correction-to-scaling term in order to achieve a self-consistent analysis of our data, in which the critical parameters become independent of the size of the window of data used in the analysis. We compare our resulting value β≈5.0±0.4 against previous numerical results and competing theoretical models. Our results confirm that the shear-driven jamming transition in Bagnoldian systems is well described by a critical scaling theory and we relate this scaling theory to the phenomenological constituent laws for dilatancy and friction.

Universal rescaling of flow curves for yield-stress fluids close to jamming

The experimental flow curves of four different yield-stress fluids with different interparticle interactions are studied near the jamming concentration. By appropriate scaling with the distance to jamming all rheology data can be collapsed onto master curves below and above jamming that meet in the shear-thinning regime and satisfy the Herschel-Bulkley and Cross equations, respectively. In spite of differing interactions in the different systems, master curves characterized by universal scaling exponents are found for the four systems. A two-state microscopic theory of heterogeneous dynamics is presented to rationalize the observed transition from Herschel-Bulkley to Cross behavior and to connect the rheological exponents to microscopic exponents for the divergence of the length and time scales of the heterogeneous dynamics. The experimental data and the microscopic theory are compared with much of the available literature data for yield-stress systems.

Relaxation times and rheology in dense athermal suspensions

We study the jamming transition in a model of elastic particles under shear at zero temperature. The key quantity is the relaxation time τ which is obtained by stopping the shearing and letting energy and pressure decay to zero. At many different densities and initial shear rates we do several such relaxations to determine the average τ. We establish that τ diverges with the same exponent as the viscosity and determine another exponent from the relation between τ and the coordination number. Though most of the simulations are done for the model with dissipation due to the motion of particles relative to an affinely shearing substrate, we also examine a model, where the dissipation is instead due to velocity differences of disks in contact, and confirm that the above-mentioned exponent is the same for these two models. We also consider finite size effects on both τ and the coordination number.

Universality of jamming criticality in overdamped shear-driven frictionless disks

We investigate the criticality of the jamming transition for overdamped shear-driven frictionless disks in two dimensions for two different models of energy dissipation: (i) Durian's bubble model with dissipation proportional to the velocity difference of particles in contact, and (ii) Durian's "mean-field" approximation to (i), with dissipation due to the velocity difference between the particle and the average uniform shear flow velocity. By considering the finite-size behavior of pressure, the pressure analog of viscosity, and the macroscopic friction σ/p, we argue that these two models share the same critical behavior.

Collective motion of macroscopic spheres floating on capillary ripples: dynamic heterogeneity and dynamic criticality

When a densely packed monolayer of macroscopic spheres floats on chaotic capillary Faraday waves, a coexistence of large scale convective motion and caging dynamics typical for glassy systems is observed. We subtract the convective mean flow using a coarse graining (homogenization) method and reveal subdiffusion for the caging time scales followed by a diffusive regime at later times. We apply the methods developed to study dynamic heterogeneity and show that the typical time and length scales of the fluctuations due to rearrangements of observed particle groups significantly increase when the system approaches its largest experimentally accessible packing concentration. To connect the system to the dynamic criticality literature, we fit power laws to our results. The resultant critical exponents are consistent with those found in densely packed suspensions of colloids.