Length scales and self-organization in dense suspension flows

Effects of coordination and stiffness scale separation in disordered elastic networks

Many fibrous materials are modeled as elastic networks featuring a substantial separation between the stiffness scales that characterize different microscopic deformation modes of the network's constituents. This scale separation has been shown to give rise to emergent complexity in these systems' linear and nonlinear mechanical response. Here we study numerically a simple model featuring said stiffness scale separation in two-dimensions and show that its mechanical response is governed by the competition between the characteristic stiffness of collective nonphononic soft modes of the stiff subsystem, and the characteristic stiffness of the soft interactions. We present and rationalize the behavior of the shear modulus of our complex networks across the unjamming transition at which the stiff subsystem alone loses its macroscopic mechanical rigidity. We further establish a relation in the soft-interaction-dominated regime between the shear modulus, the characteristic frequency of nonphononic vibrational modes, and the mesoscopic correlation length that marks the crossover from a disorder-dominated response to local mechanical perturbations in the near field, to a linear, continuumlike response in the far field. The effects of spatial dimension on the observed scaling behavior are discussed, in addition to the interplay between stiffness scales in strain-stiffened networks, which is relevant to understanding the nonlinear mechanics of non-Brownian fibrous biomatter.

Elasticity of self-organized frustrated disordered spring networks

There have been some interesting recent advances in understanding the notion of mechanical disorder in structural glasses and the statistical mechanics of these systems' low-energy excitations. Here we contribute to these advances by studying a minimal model for structural glasses' elasticity in which the degree of mechanical disorder-as characterized by recently introduced dimensionless quantifiers-is readily tunable over a very large range. We comprehensively investigate a number of scaling laws observed for various macro, meso and microscopic elastic properties, and rationalize them using scaling arguments. Interestingly, we demonstrate that the model features the universal quartic glassy vibrational density of states as seen in many atomistic and molecular models of structural glasses formed by cooling a melt. The emergence of this universal glassy spectrum highlights the role of self-organization (toward mechanical equilibrium) in its formation, and elucidates why models featuring structural frustration alone do not feature the same universal glassy spectrum. Finally, we discuss relations to existing work in the context of strain stiffening of elastic networks and of low-energy excitations in structural glasses, in addition to future research directions.

Scaling regimes and fluctuations of observables in computer glasses approaching the unjamming transition

Under decompression, disordered solids undergo an unjamming transition where they become under-coordinated and lose their structural rigidity. The mechanical and vibrational properties of these materials have been an object of theoretical, numerical, and experimental research for decades. In the study of low-coordination solids, understanding the behavior and physical interpretation of observables that diverge near the transition is of particular importance. Several such quantities are length scales (ξ or l) that characterize the size of excitations, the decay of spatial correlations, the response to perturbations, or the effect of physical constraints in the boundary or bulk of the material. Additionally, the spatial and sample-to-sample fluctuations of macroscopic observables such as contact statistics or elastic moduli diverge approaching unjamming. Here, we discuss important connections between all of these quantities and present numerical results that characterize the scaling properties of sample-to-sample contact and shear modulus fluctuations in ensembles of low-coordination disordered sphere packings and spring networks. Overall, we highlight three distinct scaling regimes and two crossovers in the disorder quantifiers χz and χμ as functions of system size N and proximity to unjamming δz. As we discuss, χX relates to the standard deviation σX of the sample-to-sample distribution of the quantity X (e.g., excess coordination δz or shear modulus μ) for an ensemble of systems. Importantly, χμ has been linked to experimentally accessible quantities that pertain to sound attenuation and the density of vibrational states in glasses. We investigate similarities and differences in the behaviors of χz and χμ near the transition and discuss the implications of our findings on current literature, unifying findings in previous studies.

Granular solids transmit stress as two-phase composites

A basic problem in the science of realistic granular matter is the plethora of heuristic models of the stress field in the absence of a first-principles theory. Such a theory is formulated here, based on the idea that static granular assemblies can be regarded as two-phase composites. A thought experiment is described, demonstrating that the state of such materials can be varied continuously from marginal stability, via a two-phase granular assembly, then porous structure, and finally be made perfectly elastic. For completeness, I review briefly the condition for marginal stability in infinitely large assemblies. The general solution for the stress equations in d=2 is reviewed in detail and shown to be consistent with the two-phase idea. A method for identifying the phases of finite regions in larger systems is constructed, providing a stability parameter that quantifies the "proximity" to the marginally stable state. The difficulty involved in deriving stress fields in such composites is a unique constraint on the boundary between phases, and, to highlight it, a simple case of a stack of plates of alternating phase is solved explicitly. An effective medium approximation, which satisfies this constraint, is then developed and analyzed in detail. This approach forms a basis for the extension of the stress theory to general granular solids that are not marginally stable or at the yield threshold.

Strain-Controlled Critical Slowing Down in the Rheology of Disordered Networks

Networks and dense suspensions frequently reside near a boundary between soft (or fluidlike) and rigid (or solidlike) regimes. Transitions between these regimes can be driven by changes in structure, density, or applied stress or strain. In general, near the onset or loss of rigidity in these systems, dissipation-limiting heterogeneous nonaffine rearrangements dominate the macroscopic viscoelastic response, giving rise to diverging relaxation times and power-law rheology. Here, we describe a simple quantitative relationship between nonaffinity and the excess viscosity. We test this nonaffinity-viscosity relationship computationally and demonstrate its rheological consequences in simulations of strained filament networks and dense suspensions. We also predict critical signatures in the rheology of semiflexible and stiff biopolymer networks near the strain stiffening transition.

Anomalous linear elasticity of disordered networks

Continuum elasticity is a powerful tool applicable in a broad range of physical systems and phenomena. Yet, understanding how and on what scales material disorder may lead to the breakdown of continuum elasticity is not fully understood. We show, based on recent theoretical developments and extensive numerical computations, that disordered elastic networks near a critical rigidity transition, such as strain-stiffened fibrous biopolymer networks that are abundant in living systems, reveal an anomalous long-range linear elastic response below a correlation length. This emergent anomalous elasticity, which is non-affine in nature, is shown to feature a qualitatively different multipole expansion structure compared to ordinary continuum elasticity, and a slower spatial decay of perturbations. The potential degree of universality of these results, their implications (e.g. for cell-cell communication through biological extracellular matrices) and open questions are briefly discussed.

Stiffening of under-constrained spring networks under isotropic strain

Disordered spring networks are a useful paradigm to examine macroscopic mechanical properties of amorphous materials. Here, we study the elastic behavior of under-constrained spring networks, i.e. networks with more degrees of freedom than springs. While such networks are usually floppy, they can be rigidified by applying external strain. Recently, an analytical formalism has been developed to predict the scaling behavior of the elastic network properties close to this rigidity transition. Here we numerically show that these predictions apply to many different classes of spring networks, including phantom triangular, Delaunay, Voronoi, and honeycomb networks. The analytical predictions further imply that the shear modulus G scales linearly with isotropic stress T close to the rigidity transition. However, this seems to be at odds with recent numerical studies suggesting an exponent between G and T that is smaller than one for some network classes. Using increased numerical precision and shear stabilization, we demonstrate here that close to the transition a linear scaling, G ∼ T, holds independent of the network class. Finally, we show that our results are not or only weakly affected by finite-size effects, depending on the network class.

Nonlocal Effects Reflect the Jamming Criticality in Frictionless Granular Flows Down Inclines

The jamming transition is accompanied by a rich phenomenology such as hysteresis or nonlocal effects that is still not well understood. Here, we experimentally investigate a model frictionless granular layer flowing down an inclined plane as a way to disentangle generic collective effects from those arising from frictional interactions. We find that thin frictionless granular layers are devoid of hysteresis of the avalanche angle, yet the layer stability increases as it gets thinner. Steady rheological laws obtained for different layer thicknesses can be collapsed into a unique master curve, supporting the idea that nonlocal effects are the consequence of the usual finite-size effects associated with the presence of a critical point. This collapse indicates that the so-called isostatic length l^{*}, the scale on which pinning a boundary freezes all remaining floppy modes, governs the effect of boundaries on flow and rules out other propositions made in the past.

Dynamic length scales in athermal, shear-driven jamming of frictionless disks in two dimensions

We carry out numerical simulations of athermally sheared, bidisperse, frictionless disks in two dimensions. From an appropriately defined velocity correlation function, we determine that there are two diverging length scales, ξ and ℓ, as the jamming transition is approached. We analyze our results using a critical scaling ansatz for the correlation function and argue that the more divergent length ℓ is a consequence of a dangerous irrelevant scaling variable and that it is ξ, which is the correlation length that determines the divergence of the system viscosity as jamming is approached from below in the liquid phase. We find that ξ∼(ϕ_{J}-ϕ)^{-ν} diverges with the critical exponent ν=1. We provide evidence that ξ measures the length scale of fluctuations in the rotation of the particle velocity field, while ℓ measures the length scale of fluctuations in the divergence of the velocity field.

Universal Relaxation Dynamics of Sphere Packings below Jamming

We show that non-Brownian suspensions of repulsive spheres below jamming display a slow relaxational dynamics with a characteristic timescale that diverges at jamming. This slow timescale is fully encoded in the structure of the unjammed packing and can be readily measured via the vibrational density of states. We show that the corresponding dynamic critical exponent is the same for randomly generated and sheared packings. Our results show that a wide variety of physical situations, from suspension rheology to algorithmic studies of the jamming transition are controlled by a unique diverging timescale, with a universal critical exponent.

Self-stresses control stiffness and stability in overconstrained disordered networks

We investigate the interplay between prestress and mechanical properties in random elastic networks. To do this in a controlled fashion, we introduce an algorithm for creating random free-standing frames that support exactly one state of self-stress. By multiplying all the bond tensions in this state of self-stress by the same number-which with the appropriate normalization corresponds to the physical prestress inside the frame-we systematically evaluate the linear mechanical response of the frame as a function of prestress. After proving that the mechanical moduli of affinely deforming frames are rigorously independent of prestress, we turn to nonaffinely deforming frames. In such frames, prestress has a profound effect on linear response: not only can it change the values of the linear modulus-an effect we demonstrate to be related to a suppressive effect of prestress on nonaffinity-but prestresses also generically trigger a bistable mechanical response. Thus, prestress can be leveraged to both augment the mechanical response of network architectures on the fly, and to actuate finite deformations. These control modalities may be of use in the design of both novel responsive materials and soft actuators.

A minimal-length approach unifies rigidity in underconstrained materials

We present an approach to understand geometric-incompatibility-induced rigidity in underconstrained materials, including subisostatic 2D spring networks and 2D and 3D vertex models for dense biological tissues. We show that in all these models a geometric criterion, represented by a minimal length [Formula: see text], determines the onset of prestresses and rigidity. This allows us to predict not only the correct scalings for the elastic material properties, but also the precise magnitudes for bulk modulus and shear modulus discontinuities at the rigidity transition as well as the magnitude of the Poynting effect. We also predict from first principles that the ratio of the excess shear modulus to the shear stress should be inversely proportional to the critical strain with a prefactor of 3. We propose that this factor of 3 is a general hallmark of geometrically induced rigidity in underconstrained materials and could be used to distinguish this effect from nonlinear mechanics of single components in experiments. Finally, our results may lay important foundations for ways to estimate [Formula: see text] from measurements of local geometric structure and thus help develop methods to characterize large-scale mechanical properties from imaging data.

Quasilocalized states of self stress in packing-derived networks

States of self stress (SSS) are assignments of forces on the edges of a network that satisfy mechanical equilibrium in the absence of external forces. In this work we show that a particular class of quasilocalized SSS in packing-derived networks, first introduced by D.M. Sussman, C.P. Goodrich, A.J. Liu (Soft Matter 12, 3982 (2016)), are characterized by a decay length that diverges as [Formula: see text] , where [Formula: see text] is the mean connectivity of the network, and [Formula: see text] is the Maxwell threshold in two dimensions, at odds with previous claims. Our results verify the previously proposed analogy between quasilocalized SSS and the mechanical response to a local dipolar force in random networks of relaxed Hookean springs. We show that the normalization factor that distinguishes between quasilocalized SSS and the response to a local dipole constitutes a measure of the mechanical coupling of the forced spring to the elastic network in which it is embedded. We further demonstrate that the lengthscale that characterizes quasilocalized SSS does not depend on its associated degree of mechanical coupling, but instead only on the network connectivity.

Vortices Enhance Diffusion in Dense Granular Flows

This Letter introduces unexpected diffusion properties in dense granular flows and shows that they result from the development of partially jammed clusters of grains, or granular vortices. Transverse diffusion coefficients D and average vortex sizes ℓ are systematically measured in simulated plane shear flows at differing inertial numbers I revealing (i) a strong deviation from the expected scaling D∝d^{2}γ[over ˙] involving the grain size d and shear rate γ[over ˙] and (ii) an increase in average vortex size ℓ at low I, following ℓ∝dI^{-1/2} but limited by the system size. A general scaling D∝ℓdγ[over ˙] is introduced that captures all the measurements and highlights the key role of vortex size. This leads to establishing a scaling for the diffusivity in dense granular flow as D∝d^{2}sqrt[γ[over ˙]/t_{i}] involving the geometric average of shear time 1/γ[over ˙] and inertial time t_{i} as the relevant time scale. Analysis of grain trajectories is further evidence that this diffusion process arises from a vortex-driven random walk.

Comment on "Spatial structure of states of self stress in jammed systems" by D. M. Sussman, C. P. Goodrich, and A. J. Liu, Soft Matter, 2016, 12, 3982

Sussman, Goodrich and Liu recently introduced a novel definition of states of self stress in packing-derived networks, and reported that the lengthscale that characterizes these states depends on the network connectivity z and spatial dimension đ as (z - 2đ)^-0.8 in two dimensions, and as (z - 2đ)^-0.6 in three dimensions. Here we derive an explicit expression for these particular states of self stress, and show that they are equivalent to the force response to a local dipolar force in random networks of relaxed Hookean springs, previously shown to be characterized by the lengthscale l_c ∼ (z - 2đ)^-1/2. We conclude that the systems studied by Sussman et al. are insufficient in size to observe the correct scaling with connectivity of the characteristic lengthscale of states of self stress.

Effect of friction on dense suspension flows of hard particles

We use numerical simulations to study the effect of particle friction on suspension flows of non-Brownian hard particles. By systematically varying the microscopic friction coefficient μ_{p} and the viscous number J, we build a phase diagram that identifies three regimes of flow: frictionless, frictional sliding, and rolling. Using energy balance in flow, we predict relations between kinetic observables, confirmed by numerical simulations. For realistic friction coefficients and small viscous numbers (below J∼10^{-3}), we show that the dominating dissipative mechanism is sliding of frictional contacts, and we characterize asymptotic behaviors as jamming is approached. Outside this regime, our observations support the idea that flow belongs to the universality class of frictionless particles. We discuss recent experiments in the context of our phase diagram.

Clogging of microfluidic systems

The transport of suspensions of microparticles in confined environments is associated with complex phenomena at the interface of fluid mechanics and soft matter. Indeed, the deposition and assembly of particles under flow involve hydrodynamic, steric and colloidal forces, and can lead to the clogging of microchannels. The formation of clogs dramatically alters the performance of both natural and engineered systems, effectively limiting the use of microfluidic technology. While the fouling of porous filters has been studied at the macroscopic level, it is only recently that the formation of clogs has been considered at the pore-scale, using microfluidic devices. In this review, we present the clogging mechanisms recently reported for suspension flows of colloidal particles and for biofluids in microfluidic channels, including sieving, bridging and aggregation of particles. We discuss the technological implications of the clogging of microchannels and the schemes that leverage the formation of clogs. We finally consider some of the outstanding challenges involving clogging in human health, which could be tackled with microfluidic methods.

Effect of particle collisions in dense suspension flows

We study nonlocal effects associated with particle collisions in dense suspension flows, in the context of the Affine Solvent Model, known to capture various aspects of the jamming transition. We show that an individual collision changes significantly the velocity field on a characteristic volume Ω_{c}∼1/δz that diverges as jamming is approached, where δz is the deficit in coordination number required to jam the system. Such an event also affects the contact forces between particles on that same volume Ω_{c}, but this change is modest in relative terms, of order f_{coll}∼f[over ¯]^{0.8}, where f[over ¯] is the typical contact force scale. We then show that the requirement that coordination is stationary (such that a collision has a finite probability to open one contact elsewhere in the system) yields the scaling of the viscosity (or equivalently the viscous number) with coordination deficit δz. The same scaling result was derived [E. DeGiuli, G. Düring, E. Lerner, and M. Wyart, Phys. Rev. E 91, 062206 (2015)PLEEE81539-375510.1103/PhysRevE.91.062206] via different arguments making an additional assumption. The present approach gives a mechanistic justification as to why the correct finite size scaling volume behaves as 1/δz and can be used to recover a marginality condition known to characterize the distributions of contact forces and gaps in jammed packings.

Phase diagram for inertial granular flows

Flows of hard granular materials depend strongly on the interparticle friction coefficient μ_{p} and on the inertial number I, which characterizes proximity to the jamming transition where flow stops. Guided by numerical simulations, we derive the phase diagram of dense inertial flow of spherical particles, finding three regimes for 10^{-4}≲I≲10^{-1}: frictionless, frictional sliding, and rolling. These are distinguished by the dominant means of energy dissipation, changing from collisional to sliding friction, and back to collisional, as μ_{p} increases from zero at constant I. The three regimes differ in their kinetics and rheology; in particular, the velocity fluctuations and the stress ratio both display nonmonotonic behavior with μ_{p}, corresponding to transitions between the three regimes of flow. We rationalize the phase boundaries between these regimes, show that energy balance yields scaling relations between microscopic properties in each of them, and derive the strain scale at which particles lose memory of their velocity. For the frictional sliding regime most relevant experimentally, we find for I≥10^{-2.5} that the growth of the macroscopic friction μ(I) with I is induced by an increase of collisional dissipation. This implies in that range that μ(I)-μ(0)∼I^{1-2b}, where b≈0.2 is an exponent that characterizes both the dimensionless velocity fluctuations L∼I^{-b} and the density of sliding contacts χ∼I^{b}.

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.

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.

Unified theory of inertial granular flows and non-Brownian suspensions

Rheological properties of dense flows of hard particles are singular as one approaches the jamming threshold where flow ceases both for aerial granular flows dominated by inertia and for over-damped suspensions. Concomitantly, the length scale characterizing velocity correlations appears to diverge at jamming. Here we introduce a theoretical framework that proposes a tentative, but potentially complete, scaling description of stationary flows. Our analysis, which focuses on frictionless particles, applies both to suspensions and inertial flows of hard particles. We compare our predictions with the empirical literature, as well as with novel numerical data. Overall, we find a very good agreement between theory and observations, except for frictional inertial flows whose scaling properties clearly differ from frictionless systems. For overdamped flows, more observations are needed to decide if friction is a relevant perturbation. Our analysis makes several new predictions on microscopic dynamical quantities that should be accessible experimentally.

Diverging viscosity and soft granular rheology in non-Brownian suspensions

We use large scale computer simulations and finite-size scaling analysis to study the shear rheology of dense three-dimensional suspensions of frictionless non-Brownian particles in the vicinity of the jamming transition. We perform simulations of soft repulsive particles at constant shear rate, constant pressure, and finite system size and carefully study the asymptotic limits of large system sizes and infinitely hard particle repulsion. We first focus on the asymptotic behavior of the shear viscosity in the hard particle limit. By measuring the viscosity increase over about 5 orders of magnitude, we are able to confirm its asymptotic power law divergence close to the jamming transition. However, a precise determination of the critical density and critical exponent is difficult due to the "multiscaling" behavior of the viscosity. Additionally, finite-size scaling analysis suggests that this divergence is accompanied by a growing correlation length scale, which also diverges algebraically. Finally, we study the effect of particle softness and propose a natural extension of the standard granular rheology, which we test against our simulation data. Close to the jamming transition, this "soft granular rheology" offers a detailed description of the nonlinear rheology of soft particles, which differs from earlier empirical scaling forms.

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.