Search

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 phi and a parameter Q that measures the elasticity of binary particle collisions. We find a Q*(phi) 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*(phi) diverges as phi -> phi(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 fJ. 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.

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 sigma. Our analysis determines the critical exponent beta that describes the algebraic divergence of the Bagnold transport coefficients lim((gamma) over dot -> 0) p/(gamma) over dot(2), sigma/(gamma) over dot(2) similar to (phi(J) -phi)(-beta) as the jamming transition phi(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 beta approximate to 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.

We use numerical simulations to investigate the effect that different models of energy dissipation have on the rheology of soft-core frictionless disks, below jamming in two dimensions. We find that it is not necessarily the mass of the particles that determines whether a system has Bagnoldian or Newtonian rheology, but rather the presence or absence of large connected clusters of particles. We demonstrate the key role that tangential dissipation plays in the formation of such clusters and in several models find a transition from Bagnoldian to Newtonian rheology as the packing fraction. is varied. For each model, we show that appropriately scaled rheology curves approach a well defined limit as the mass of the particles decreases and collisions become strongly inelastic.

We numerically study the distributions of global pressure that are found in ensembles of statically jammed and quasistatically sheared systems of bidisperse, frictionless disks at fixed packing fraction phi in two dimensions. We use these distributions to address the question of how pressure increases as phi increases above the jamming point phi(J), p similar to |phi - phi(J) |(y). For statically jammed ensembles, our results are consistent with the exponent y being simply related to the power law of the interparticle soft-core interaction. For sheared systems, however, the value of y is consistent with a nontrivial value, as found previously in rheological simulations., Originally published in dissertation in manuscript form.

In a recent paper [Mandal et al., Phys. Rev. E 88, 022129 (2013)], the nature of spatial correlations of plasticity in hard-sphere glasses was addressed both via computer simulations and in experiments. It was found that the experimentally obtained correlations obey a power law, whereas the correlations from simulations are better fitted by an exponential decay. We here provide direct evidence-via simulations of a hard-sphere glass in two dimensions (2D)-that this discrepancy is a consequence of the finite system size in the 3D simulations. By extending the study to a 2D soft disk model at zero temperature [Durian, Phys. Rev. Lett. 75, 4780 (1995)], the robustness of the power-law decay in sheared amorphous solids is underlined. Deviations from a power law occur when either reducing the packing fraction towards the supercooled regime in the case of hard spheres or changing the dissipation mechanism from contact dissipation to a mean-field-type drag in the case of soft 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.

The jamming of bidisperse soft core disks is considered, using a variety of different protocols to produce the jammed state. In agreement with other works, we find that cooling and compression can lead to a broad range of jamming packing fractions ϕJ, depending on cooling rate and initial configuration; the larger the degree of big particle clustering in the initial configuration, the larger will be the value of ϕJ. In contrast, we find that shearing disrupts particle clustering, leading to a much narrower range of ϕJ as the shear strain rate varies. In the limit of vanishingly small shear strain rate, we find a unique nontrivial value for the jamming density that is independent of the initial system configuration. We conclude that shear driven jamming is a unique and well-defined critical point in the space of shear driven steady states. We clarify the relation between glassy behavior, rigidity, and jamming in such systems and relate our results to recent experiments.

We carry out a finite-size scaling analysis of the jamming transition in frictionless bidisperse soft core disks in two dimensions. We consider two different jamming protocols: (i) quench from random initial positions and (ii) quasistatic shearing. By considering the fraction of jammed states as a function of packing fraction for systems with different numbers of particles, we determine the spatial correlation length critical exponent ν≈1 and show that corrections to scaling are crucial for analyzing the data. We show that earlier numerical results yielding ν<1 are due to the improper neglect of these corrections.