Chaos and threshold for irreversibility in sheared suspensions

(Ir)reversibility in dense granular systems driven by oscillating forces

We use computer simulations to study highly dense systems of granular particles that are driven by oscillating forces. We implement different dissipation mechanisms that are used to extract the injected energy. In particular, the action of a simple local Stokes' drag is compared with non-linear and history-dependent frictional forces that act either between particle pairs or between particles and an external container wall. The Stokes' drag leads to particle motion that is periodic with the driving force, even at high densities around close packing where particles undergo frequent collisions. With the introduction of inter-particle frictional forces this "interacting absorbing state" is destroyed and particles start to diffuse around. By reducing the density of the material we go through another transition to a "non-interacting" absorbing state, where particles independently follow the force-induced oscillations without collisions. In the system with particle-wall frictional interactions this transition has signs of a discontinuous phase transition. It is accompanied by a diverging relaxation time, but not by a vanishing order parameter, which rather jumps to zero at the transition.

Anomalous preasymptotic colloid transport by hydrodynamic dispersion in microfluidic capillary flow

The anomalous preasymptotic transport of colloids in a microfluidic capillary flow due to hydrodynamic dispersion is measured by noninvasive nuclear magnetic resonance (NMR). The data indicate a reduced scaling of mean squared displacement with time from the 〈z(t)(2)〉(c) ∼ t(3) behavior for the interaction of a normal diffusion process with a simple shear flow. This nonequilibrium steady-state system is shown to be modeled by a continuous time random walk (CTRW) on a moving fluid. The full propagator of the motion is measured by NMR, providing verification of the assumption of Gaussian jump length distributions in the CTRW model. The connection of the data to microrheology measurements by NMR, in which every particle in a suspension contributes information, is established.

Consequences of anomalous diffusion in disordered systems under cyclic forcing

We study the particle scale response of a 2D frictionless disk system to bulk forcing via cyclic shear with reversal amplitude γ_{r}. We find a subdiffusive γ_{r}-dependent regime, which is consistent with models of anomalous diffusion with scale-invariant cage dynamics, and a crossover to diffusive grain motion at high γ_{r}. Analysis of local displacements of a particle relative to its cage of neighbors reveals a key distinction from thermal systems. Particles are moved by fluctuations of their cage of neighbors rather than rattling in their cage, indicating a distinct cage-breaking mechanism.

Variable-amplitude oscillatory shear response of amorphous materials

Variable-amplitude oscillatory shear tests are emerging as powerful tools to investigate and quantify the nonlinear rheology of amorphous solids, complex fluids, and biological materials. Quite a few recent experimental and atomistic simulation studies demonstrated that at low shear amplitudes, an amorphous solid settles into an amplitude- and initial-conditions-dependent dissipative limit cycle, in which back-and-forth localized particle rearrangements periodically bring the system to the same state. At sufficiently large shear amplitudes, the amorphous system loses memory of the initial conditions, exhibits chaotic particle motions accompanied by diffusive behavior, and settles into a stochastic steady state. The two regimes are separated by a transition amplitude, possibly characterized by some critical-like features. Here we argue that these observations support some of the physical assumptions embodied in the nonequilibrium thermodynamic, internal-variables based, shear-transformation-zone model of amorphous viscoplasticity; most notably that "flow defects" in amorphous solids are characterized by internal states between which they can make transitions, and that structural evolution is driven by dissipation associated with plastic deformation. We present a rather extensive theoretical analysis of the thermodynamic shear-transformation-zone model for a variable-amplitude oscillatory shear protocol, highlighting its success in accounting for various experimental and simulational observations, as well as its limitations. Our results offer a continuum-level theoretical framework for interpreting the variable-amplitude oscillatory shear response of amorphous solids and may promote additional developments.

Flight-crash events in turbulence

The statistical properties of turbulence differ in an essential way from those of systems in or near thermal equilibrium because of the flux of energy between vastly different scales at which energy is supplied and at which it is dissipated. We elucidate this difference by studying experimentally and numerically the fluctuations of the energy of a small fluid particle moving in a turbulent fluid. We demonstrate how the fundamental property of detailed balance is broken, so that the probabilities of forward and backward transitions are not equal for turbulence. In physical terms, we found that in a large set of flow configurations, fluid elements decelerate faster than accelerate, a feature known all too well from driving in dense traffic. The statistical signature of rare "flight-crash" events, associated with fast particle deceleration, provides a way to quantify irreversibility in a turbulent flow. Namely, we find that the third moment of the power fluctuations along a trajectory, nondimensionalized by the energy flux, displays a remarkable power law as a function of the Reynolds number, both in two and in three spatial dimensions. This establishes a relation between the irreversibility of the system and the range of active scales. We speculate that the breakdown of the detailed balance characterized here is a general feature of other systems very far from equilibrium, displaying a wide range of spatial scales.

Colored-noise Fokker-Planck equation for the shear-induced self-diffusion process of non-Brownian particles

In the literature, it is pointed out that non-Brownian particles tend to show shear-induced diffusive behavior due to hydrodynamic interactions. Several authors indicate a long correlation time of the particle velocities in comparison to Brownian particle velocities modeled by a white noise. This work deals with the derivation of a Fokker-Planck equation both in position and velocity space which describes the process of shear-induced self-diffusion, whereas, so far, this problem has been described by Fokker-Planck equations restricted to position space. The long velocity correlation times actually would necessitate large time-step sizes in the mathematical description of the problem in order to capture the diffusive regime. In fact, time steps of specific lengths pose problems to the derivation of the corresponding Fokker-Planck equation because the whole particle configuration changes during long time-step sizes. On the other hand, small time-step sizes, i.e., in the range of the velocity correlation time, violate the Markov property of the position variable. In this work we regard the problem of shear-induced self-diffusion with respect to the Markov property and reformulate the problem with respect to small time-step sizes. In this derivation, we regard the nondimensionalized Langevin equation and develop a new compact form which allows us to analyze the Langevin equation for all time scales of interest for both Brownian and non-Brownian particles starting from a single equation. This shows that the Fokker-Planck equation in position space should be extended to a colored-noise Fokker-Planck equation in both position and colored-noise velocity space, which we will derive.

Tracking control of colloidal particles through non-homogeneous stationary flows

We consider the problem of controlling the trajectory of a single colloidal particle in a fluid with steady non-homogeneous flow. We use a Langevin equation to describe the dynamics of this particle, where the friction term is assumed to be given by the Faxén's Theorem for the force on a sphere immersed in a stationary flow. We use this description to propose an explicit control force field to be applied on the particle such that it will follow asymptotically any given desired trajectory, starting from an arbitrary initial condition. We show that the dynamics of the controlled particle can be mapped into a set of stochastic harmonic oscillators and that the velocity gradient of the solvent induces an asymmetric coupling between them. We study the particular case of a Brownian particle controlled through a plane Couette flow and show explicitly that the velocity gradient of the solvent renders the dynamics non-stationary and non-reversible in time. We quantify this effect in terms of the correlation functions for the position of the controlled particle, which turn out to exhibit contributions depending exclusively on the non-equilibrium character of the state of the solvent. In order to test the validity of our model, we perform simulations of the controlled particle moving in a simple shear flow, using a hybrid method combining molecular dynamics and multi-particle collision dynamics. We confirm numerically that the proposed guiding force allows for controlling the trajectory of the micro-sized particle by obligating it to follow diverse specific trajectories in fluids with homogeneous shear rates of different strengths. In addition, we find that the non-equilibrium correlation functions in simulations exhibit the same qualitative behavior predicted by the model, thus revealing the presence of the asymmetric non-equilibrium coupling mechanism induced by the velocity gradient.

Connecting the random organization transition and jamming within a unifying model system

While the random organization transition describes the change from reversible to irreversible dynamics in a nonequilibrium system, the athermal jamming transition at zero shear rate occurs when particles can no longer avoid overlaps. Despite the obvious differences between these two transitions, we show that they both occur within the same model packing problem. In this unifying model system the particles are first randomly distributed and then displaced in each step if they overlap. For random displacements we obtain a random organization transition, while jamming occurs in the case of deterministic shifts. We also analyze the critical behavior of random organization. Our results show that random organization and jamming are opposite limits of random sphere packings, and we expect that various equilibrium and nonequilibrium transitions can be formulated as related intermediate packing problems.

Onset of irreversibility and chaos in amorphous solids under periodic shear

A fundamental problem in the physics of amorphous materials is understanding the transition from reversible to irreversible plastic behavior and its connection to yield. Currently, continuum material modeling relies on phenomenological yield thresholds, however in many cases the transition from elastic to plastic behavior is gradual, which makes it difficult to identify an exact yield criterion. Here we show that under periodic shear, amorphous solids undergo a transition from repetitive, predictable behavior to chaotic, irregular behavior as a function of the strain amplitude. In both the periodic and chaotic regimes, localized particle rearrangements are observed. We associate the point of transition from repetitive to chaotic behavior with the yield strain and suggest that at least for oscillatory shear, yield in amorphous solids is a result of a "transition to chaos."

Particle-scale reversibility in athermal particulate media below jamming

We perform numerical simulations of repulsive, frictionless athermal disks in two and three spatial dimensions undergoing cyclic quasistatic simple shear to investigate particle-scale reversible motion. We identify three classes of steady-state dynamics as a function of packing fraction φ and maximum strain amplitude per cycle γ(max). Point-reversible states, where particles do not collide and exactly retrace their intracycle trajectories, occur at low φ and γ(max). Particles in loop-reversible states undergo numerous collisions and execute complex trajectories but return to their initial positions at the end of each cycle. For sufficiently large φ and γ(max), systems display irreversible dynamics with nonzero self-diffusion. Loop-reversible dynamics enables the reliable preparation of configurations with specified structural and mechanical properties over a broad range of φ.

Contact processes in crowded environments

Periodically sheared colloids at low densities demonstrate a dynamical phase transition from an inactive to active phase as the strain amplitude is increased. The inactive phase consists of no collisions (contacts) between particles in the steady state limit, while in the active phase collisions persist. To investigate this system at higher densities, we construct and study a conserved-particle-number contact process with three-body interactions, which are potentially more likely than two-body interactions at higher densities. For example, consider one active (diffusing) particle colliding with two inactive (nondiffusing) particles such that they become active and consider spontaneous inactivation. In mean field, this system exhibits a continuous dynamical phase transition. Simulations on square lattices also indicate a continuous transition with exponents similar to those measured for the conserved lattice gas (CLG) model. In contrast, the three-body interaction requiring two active particles to activate one inactive particle exhibits a discontinuous transition. Finally, inspired by kinetically constrained models of the glass transition, we investigate the "caging effect" at even higher particle densities to look for a second dynamical phase transition back to an inactive phase. Square lattice simulations suggest a continuous transition with a new set of exponents differing from both the CLG model and what is known as directed percolation, indicating a potentially new universality class for a contact process with a conserved particle number.

Multiple transient memories in sheared suspensions: robustness, structure, and routes to plasticity

Multiple transient memories, originally discovered in charge-density-wave conductors, are a remarkable and initially counterintuitive example of how a system can store information about its driving. In this class of memories, a system can learn multiple driving inputs, nearly all of which are eventually forgotten despite their continual input. If sufficient noise is present, the system regains plasticity so that it can continue to learn new memories indefinitely. Recently, Keim and Nagel [Phys. Rev. Lett. 107, 010603 (2011)] showed how multiple transient memories could be generalized to a generic driven disordered system with noise, giving as an example simulations of a simple model of a sheared non-Brownian suspension. Here, we further explore simulation models of suspensions under cyclic shear, focusing on three main themes: robustness, structure, and overdriving. We show that multiple transient memories are a robust feature independent of many details of the model. The steady-state spatial distribution of the particles is sensitive to the driving algorithm; nonetheless, the memory formation is independent of such a change in particle correlations. Finally, we demonstrate that overdriving provides another means for controlling memory formation and retention.

Shear thickening in non-Brownian suspensions: an excluded volume effect

Shear thickening appears as an increase of the viscosity of a dense suspension with the shear rate, sometimes sudden and violent at high volume fraction. Its origin for noncolloidal suspension with non-negligible inertial effects is still debated. Here we consider a simple shear flow and demonstrate that fluid inertia causes a strong microstructure anisotropy that results in the formation of a shadow region with no relative flux of particles. We show that shear thickening at finite inertia can be explained as an increase of the effective volume fraction when considering the dynamically excluded volume due to these shadow regions.

Suppression and emergence of granular segregation under cyclic shear

While convective flows are implicated in many granular segregation processes, the associated particle-scale rearrangements are not well understood. A three-dimensional bidisperse mixture segregates under steady shear, but the cyclically driven system either remains mixed or segregates slowly. Individual grain motion shows no signs of particle-scale segregation dynamics that precede bulk segregation. Instead, we find that the transition from nonsegregating to segregating flow is accompanied by significantly less reversible particle trajectories and the emergence of a convective flow field.

Oscillatory athermal quasistatic deformation of a model glass

We report computer simulations of oscillatory athermal quasistatic shear deformation of dense amorphous samples of a three-dimensional model glass former. A dynamical transition is observed as the amplitude of the deformation is varied: For large values of the amplitude the system exhibits diffusive behavior and loss of memory of the initial conditions, whereas localization is observed for small amplitudes. Our results suggest that the same kind of transition found in driven colloidal systems is present in the case of amorphous solids (e.g., metallic glasses). The onset of the transition is shown to be related to the onset of energy dissipation. Shear banding is observed for large system sizes, without, however, affecting qualitative aspects of the transition.

Translational and rotational dynamics of colloidal particles in suspension: effect of shear

We report a generalization of a nonequilibrium thermodynamic theory for the mesoscopic dynamics of radially symmetric interacting particles to anisotropic pairwise interactions and attain the one- and two-particle Fokker-Planck kinetics equations at a low-density limit that provides the translational-rotational coupling of their motion due to hydrodynamic interactions, from which we derived the balance equations of linear, angular momentum, and energy dissipation due to particle interactions and energy interchange with heat bath. In this low-density approximation, an already-known virial expression for the long-time translational collective diffusion coefficient of an orientational isotropic suspension in terms of the fluid equilibrium microstructure is recovered. An external shear flow induces, in the diffusive regime, vorticity effects into the rotational diffusion property of the colloidal particles. They manifest in the appearance of the particle's rotational viscosity due to vortex flow. The Smoluchowski equation that governs the dynamical relaxation of colloid microstructure due to particle's Brownian motion under stationary flow is provided.

Irreversibility and chaos: role of lubrication interactions in sheared suspensions

We investigate non-Brownian particles suspended in a periodic shear-flow using simulations. Following Metzger and Butler [Phys. Rev. E 82, 051406 (2010)], we show that the chaotic dynamics arising from lubrication interactions are too weak to generate an observable particle dispersion. The irreversibility observed in periodic flow is dominated by contact interactions. Nonetheless, we show that lubrication interactions must be included in the calculation to obtain results that agree with experiments.

Complex oscillatory yielding of model hard-sphere glasses

The yielding behavior of hard sphere glasses under large-amplitude oscillatory shear has been studied by probing the interplay of Brownian motion and shear-induced diffusion at varying oscillation frequencies. Stress, structure and dynamics are followed by experimental rheology and Browian dynamics simulations. Brownian-motion-assisted cage escape dominates at low frequencies while escape through shear-induced collisions at high ones, both related with a yielding peak in G''. At intermediate frequencies a novel, for hard sphere glasses, double peak in G'' is revealed reflecting both mechanisms. At high frequencies and strain amplitudes a persistent structural anisotropy causes a stress drop within the cycle after strain reversal, while higher stress harmonics are minimized at certain strain amplitudes indicating an apparent harmonic response.

Melting a granular glass by cooling

Driven granular systems readily form glassy phases at high particle volume fractions and low driving amplitudes. We use computer simulations of a driven granular glass to evidence a reentrance melting transition into a fluid state, which, contrary to intuition, occurs by reducing the amplitude of the driving. This transition is accompanied by anomalous particle dynamics and superdiffusive behavior on intermediate time scales. We highlight the special role played by frictional interactions, which help particles to escape their glassy cages. Such an effect is in striking contrast to what friction is expected to do: reduce particle mobility by making them stick.

Impact of chaos and Brownian diffusion on irreversibility in Stokes flows

We study a reversal process in Stokes flows in the presence of weak diffusion in order to clarify the distinct effects that chaotic flows have on the loss of reversibility relative to nonchaotic flows. In all linear flows, including a representation of the baker's map, we show that the decay of reversibility presents universal properties. In nonlinear chaotic and nonchaotic flows, we show that this universality breaks down due to the distribution of strain rates. In the limit of infinitesimal diffusivity, we predict qualitatively distinct behavior in the chaotic case.

Hysteresis and return-point memory in colloidal artificial spin ice systems

Using computer simulations, we investigate hysteresis loops and return-point memory for artificial square and kagome spin ice systems by cycling an applied bias force and comparing microscopic effective spin configurations throughout the hysteresis cycle. Return-point memory loss is caused by motion of individual defects in kagome ice or of grain boundaries in square ice. In successive cycles, return-point memory is recovered rapidly in kagome ice. Memory is recovered more gradually in square ice due to the extended nature of the grain boundaries. Increasing the amount of quenched disorder increases the defect density but also enhances the return-point memory since the defects become trapped more easily.

Continuous separation of microparticles in a microfluidic channel via the elasto-inertial effect of non-Newtonian fluid

Pure separation and sorting of microparticles from complex fluids are essential for biochemical analyses and clinical diagnostics. However, conventional techniques require highly complex and expensive labeling processes for high purity separation. In this study, we present a simple and label-free method for separating microparticles with high purity using the elasto-inertial characteristic of a non-Newtonian fluid in microchannel flow. At the inlet, particle-containing sample flow was pushed toward the side walls by introducing sheath fluid from the center inlet. Particles of 1 μm and 5 μm in diameter, which were suspended in viscoelastic fluid, were successfully separated in the outlet channels: larger particles were notably focused on the centerline of the channel at the outlet, while smaller particles continued flowing along the side walls with minimal lateral migration towards the centerline. The same technique was further applied to separate platelets from diluted whole blood. Through cytometric analysis, we obtained a purity of collected platelets of close to 99.9%. Conclusively, our microparticle separation technique using elasto-inertial forces in non-Newtonian fluid is an effective method for separating and collecting microparticles on the basis of size differences with high purity.

Onset of irreversibility in cyclic shear of granular packings

We investigate the onset of irreversibility in a dense granular medium subjected to cyclic shear in a split-bottom geometry. To probe the micro- and mesoscale, we image bead trajectories in three dimensions throughout a series of shear strain oscillations. Although beads lose and regain contact with neighbors during a cycle, the global topology of the contact network exhibits reversible properties for low oscillation amplitudes. With increasing reversal amplitude, a transition to an irreversible diffusive regime occurs.

Assembly of vorticity-aligned hard-sphere colloidal strings in a simple shear flow

Colloidal suspensions self-assemble into equilibrium structures ranging from face- and body-centered cubic crystals to binary ionic crystals, and even kagome lattices. When driven out of equilibrium by hydrodynamic interactions, even more diverse structures can be accessed. However, mechanisms underlying out-of-equilibrium assembly are much less understood, though such processes are clearly relevant in many natural and industrial systems. Even in the simple case of hard-sphere colloidal particles under shear, there are conflicting predictions about whether particles link up into string-like structures along the shear flow direction. Here, using confocal microscopy, we measure the shear-induced suspension structure. Surprisingly, rather than flow-aligned strings, we observe log-rolling strings of particles normal to the plane of shear. By employing Stokesian dynamics simulations, we address the mechanism leading to this out-of-equilibrium structure and show that it emerges from a delicate balance between hydrodynamic and interparticle interactions. These results demonstrate a method for assembling large-scale particle structures using shear flows.

Invited Article: Refractive index matched scanning of dense granular materials

We review an experimental method that allows to probe the time-dependent structure of fully three-dimensional densely packed granular materials and suspensions by means of particle recognition. The method relies on submersing a granular medium in a refractive index matched fluid. This makes the resulting suspension transparent. The granular medium is then visualized by exciting, layer by layer, the fluorescent dye in the fluid phase. We collect references and unreported experimental know-how to provide a solid background for future development of the technique, both for new and experienced users.

Transverse alignment of fibers in a periodically sheared suspension: an absorbing phase transition with a slowly varying control parameter

Shearing solutions of fibers or polymers tends to align fiber or polymers in the flow direction. Here, non-Brownian rods subjected to oscillatory shear align perpendicular to the flow while the system undergoes a nonequilibrium absorbing phase transition. The slow alignment of the fibers can drive the system through the critical point and thus promote the transition to an absorbing state. This picture is confirmed by a universal scaling relation that collapses the data with critical exponents that are consistent with conserved directed percolation.

Generic transient memory formation in disordered systems with noise

Out-of-equilibrium disordered systems may form memories of external driving in a remarkable fashion. The system "remembers" multiple values from a series of training inputs yet "forgets" nearly all of them at long times despite the inputs being continually repeated. Here, learning and forgetting are inseparable aspects of a single process. The memory loss may be prevented by the addition of noise. We identify a class of systems with this behavior, giving as an example a model of non-Brownian suspensions under cyclic shear.

Irreversibility and chaos: role of long-range hydrodynamic interactions in sheared suspensions

Non-Brownian particles suspended in an oscillatory shear flow are studied numerically. In these systems it is often assumed that chaos (due to the long-range nature of the hydrodynamic interaction between particles) plus noise (contact or roughness) lead to irreversible behavior. However, we demonstrate that the long-range hydrodynamic interactions are not a source, nor even a magnifier, of irreversibility when coupled with nonhydrodynamic interactions. Additionally, analysis reveals that the apparent anisotropy of the particle diffusion is due to coupling of the shear flow and transverse diffusion.

Hydrodynamic irreversibility in particle suspensions with nonuniform strain

A dynamical phase transition from reversible to irreversible behavior occurs when particle suspensions are subjected to uniform oscillatory shear, even in the Stokes flow limit. We consider a more general situation with nonuniform strain (e.g., oscillatory channel flow), which is observed to exhibit markedly different dynamics. The onset of irreversibility is delayed, and occurs simultaneously across the entire channel. This behavior is only partially explained by self-organization and shear-induced migration. The onset is accompanied by long-range correlated particle motion even at the channel center, where the strain is negligible; this motion prevents the system from evolving into a reversible state.

Shear melting of a colloidal glass

We use confocal microscopy to explore shear melting of colloidal glasses, which occurs at strains of approximately 0.08, coinciding with a strongly non-Gaussian step size distribution. For larger strains, the particle mean square displacement increases linearly with strain and the step size distribution becomes Gaussian. The effective diffusion coefficient varies approximately linearly with shear rate, consistent with a modified Stokes-Einstein relationship in which thermal energy is replaced by shear energy and the length scale is set by the size of cooperatively moving regions consisting of approximately 3 particles.

Nonequilibrium fluctuation-dissipation relations of interacting Brownian particles driven by shear

We present a detailed analysis of the fluctuation-dissipation theorem (FDT) close to the glass transition in colloidal suspensions under steady shear using mode coupling approximations. Starting point is the many-particle Smoluchowski equation. Under shear, detailed balance is broken and the response functions in the stationary state are smaller at long times than estimated from the equilibrium FDT. An asymptotically constant relation connects response and fluctuations during the shear driven decay, restoring the form of the FDT with, however, a ratio different from the equilibrium one. At short times, the equilibrium FDT holds. We follow two independent approaches whose results are in qualitative agreement. To discuss the derived fluctuation-dissipation ratios, we show an exact reformulation of the susceptibility which contains not the full Smoluchowski operator as in equilibrium, but only its well defined Hermitian part. This Hermitian part can be interpreted as governing the dynamics in the frame comoving with the probability current. We present a simple toy model which illustrates the FDT violation in the sheared colloidal system.

Microscopic structural relaxation in a sheared supercooled colloidal liquid

The rheology of dense amorphous materials under large shear strain is not fully understood, partly due to the difficulty of directly viewing the microscopic details of such materials. We use a colloidal suspension to simulate amorphous materials and study the shear-induced structural relaxation with fast confocal microscopy. We quantify the plastic rearrangements of the particles in several ways. Each of these measures of plasticity reveals spatially heterogeneous dynamics, with localized regions where many particles are strongly rearranging by these measures. We examine the shapes of these regions and find them to be essentially isotropic, with no alignment in any particular direction. Furthermore, individual particles are equally likely to move in any direction other than the overall bias imposed by the strain.

Self-organized criticality in sheared suspensions

Recent studies reveal that suspensions of neutrally buoyant non-brownian particles driven by slow periodic shear can undergo a dynamical phase transition between a fluctuating irreversible steady state and an absorbing reversible state. Using a computer model, we show that such systems exhibit self-organized criticality when a finite particle sedimentation velocity v(s) is introduced. Under periodic shear, these systems evolve, without external intervention, towards the shear-dependent critical concentration phi(c) as v(s) is reduced. This state is characterized by power-law distributions in the lifetime and size of fluctuating clusters. Experiments exhibit similar behavior and, as v(s) is reduced, yield steady-state values of phi that tend towards the phi(c) corresponding to the applied shear.

Random organization and plastic depinning

We provide evidence that the general phenomenon of plastic depinning can be described as an absorbing phase transition, and shows the same features as the random organization which was recently studied in periodically driven particle systems [L. Corte, Nature Phys. 4, 420 (2008)]. In the plastic flow system, the pinned regime corresponds to the absorbing state and the moving state corresponds to the fluctuating state. When an external force is suddenly applied, the system eventually organizes into one of these two states with a time scale that diverges as a power law at a nonequilibrium transition. We propose a simple experiment to test for this transition in systems with random disorder.

Universality class of the reversible-irreversible transition in sheared suspensions

Collections of non-Brownian particles suspended in a viscous fluid and subjected to oscillatory shear at very low Reynolds number have recently been shown to exhibit a remarkable dynamical phase transition separating reversible from irreversible behavior as the strain amplitude or volume fraction are increased. We present a simple model for this phenomenon, based on which we argue that this transition lies in the universality class of the conserved directed percolation models. This leads to predictions for the scaling behavior of a large number of experimental observables. Non-Brownian suspensions under oscillatory shear may thus constitute the first experimental realization of an inactive-active phase transition which is not in the universality class of conventional directed percolation.

Growing length and time scales in a suspension of athermal particles

We simulate a relaxation process of non-Brownian particles in a sheared viscous medium; the system is subject to the small shear strain and then undergoes relaxation. We estimate the exponents with which the relaxation time and the correlation length diverge as the density approaches the jamming density from below. In particular, the dynamic critical exponent is estimated as 4.6(2). It is also found that shear stress undergoes power-law decay at the jamming density, which is reminiscent of critical slowing down.

Particle pressure in a sheared suspension: a bridge from osmosis to granular dilatancy

The normal stress exerted by particles in a sheared suspension is measured by analogy with a method used to measure osmotic pressure in solutions. Particles in a liquid are confined by a fine screen to a gap between two vertical concentric cylinders, the inner of which rotates. Pressure in the liquid is sensed either by a manometer or by a pressure transducer across the screen. The particles are large enough so that Brownian motion and equilibrium osmotic pressure are vanishingly small. The measured pressure yields the shear-induced particle pressure Pi, the nonequilibrium continuation of equilibrium osmotic pressure. For volume fractions 0.3< or =varphi< or =0.5, Pi is strongly dependent on varphi, and linear in shear rate. Comparisons of the measured particle pressure with modeling and simulation show good agreement.

Stochastic transport of particles in straining flows

Important features associated with the segregration of particles in turbulent flow are investigated by considering the statistical distribution (phase-space number density) of particles subject to the combined effects of straining flow and stochastic forcing. A Fokker-Planck model is used to obtain results for the phase-space distributions of particles that are entrained into straining flow fields. The analysis shows that, in marked contrast to the zero strain case, nonsingular steady-state distributions are generated, and also confirms that the diffusional effect resulting from stochastic forcing is sufficient to offset the otherwise singular distributions that would result from the indefinite accumulation of particles along stagnation lines. The influence of particle inertia (Stokes number) on the form of the resulting distributions is considered and several significant results are observed. The influence of strain rate on the attenuation of particle kinetic stresses is quantified and explained. The development of large third-order velocity moments is observed for Stokes numbers above a critical value. The mechanism underlying this phenomenon is seen to be a generic feature of particle transport in flows where vortex structures induce local counterflows of particles. The system therefore provides an ideal test for closure models for third-order moments of particle velocities, and here the standard Chapman-Enskog approximation is assessed.