Nonlocal constitutive relation for steady granular flow

Stable small bubble clusters in two-dimensional foams

Key features of the mechanical response of amorphous particulate materials, such as foams, emulsions, and granular media, to applied stress are determined by the frequency and size of particle rearrangements that occur as the system transitions from one mechanically stable state to another. This work describes coordinated experimental and computational studies of bubble rafts, which are quasi-two dimensional systems of bubbles confined to the air-water interface. We focus on small mechanically stable clusters of four, five, six, and seven bubbles with two different sizes with diameter ratio σ_L/σ_S ≃ 1.4. Focusing on small bubble clusters, which can be viewed as subsystems of a larger system, allows us to investigate the full ensemble of clusters that form, measure the respective frequencies with which the clusters occur, and determine the form of the bubble-bubble interactions. We emphasize several important results. First, for clusters with N > 5 bubbles, we find using discrete element simulations that short-range attractive interactions between bubbles give rise to a larger ensemble of distinct mechanically stable clusters compared to that generated by long-range attractive interactions. The additional clusters in systems with short-range attractions possess larger gaps between pairs of neighboring bubbles on the periphery of the clusters. The ensemble of bubble clusters observed in experiments is similar to the ensemble of clusters with long-range attractive interactions. We also compare the frequency with which each cluster occurs in simulations and experiments. We find that the cluster frequencies are extremely sensitive to the protocol used to generate them and only weakly correlated to the energy of the clusters.

General scaling relations for locomotion in granular media

Inspired by dynamic similarity in fluid systems, we have derived a general dimensionless form for locomotion in granular materials, which is validated in experiments and discrete element method (DEM) simulations. The form instructs how to scale size, mass, and driving parameters in order to relate dynamic behaviors of different locomotors in the same granular media. The scaling can be derived by assuming intrusion forces arise from resistive force theory or equivalently by assuming the granular material behaves as a continuum obeying a frictional yield criterion. The scalings are experimentally confirmed using pairs of wheels of various shapes and sizes under many driving conditions in a common sand bed. We discuss why the two models provide such a robust set of scaling laws even though they neglect a number of the complexities of granular rheology. Motivated by potential extraplanetary applications, the dimensionless form also implies a way to predict wheel performance in one ambient gravity based on tests in a different ambient gravity. We confirm this using DEM simulations, which show that scaling relations are satisfied over an array of driving modes even when gravity differs between scaled tests.

Scaling description of non-local rheology

The plastic flow of amorphous materials displays non-local effects, characterized by a cooperativity length scale ξ. We argue that these effects enter in the more general description of surface phenomena near critical points. Using this approach, we obtain a scaling relation between exponents that describe the strain rate profiles in shear driven and pressure driven flow, which we confirm both in numerical models and experimental data. We find empirically that the cooperative length follows closely the characteristic length previously extracted in homogenous bulk flows. This analysis shows that the often used mean field exponents fail to capture quantitatively the non-local effects. Our analysis also explains the unusually large finite size effects previously observed in pressure driven flows.

Well-posed continuum equations for granular flow with compressibility and μ(I)-rheology

Continuum modelling of granular flow has been plagued with the issue of ill-posed dynamic equations for a long time. Equations for incompressible, two-dimensional flow based on the Coulomb friction law are ill-posed regardless of the deformation, whereas the rate-dependent μ(I)-rheology is ill-posed when the non-dimensional inertial number I is too high or too low. Here, incorporating ideas from critical-state soil mechanics, we derive conditions for well-posedness of partial differential equations that combine compressibility with I-dependent rheology. When the I-dependence comes from a specific friction coefficient μ(I), our results show that, with compressibility, the equations are well-posed for all deformation rates provided that μ(I) satisfies certain minimal, physically natural, inequalities.

Nonlocal Elasticity near Jamming in Frictionless Soft Spheres

We use simulations of frictionless soft sphere packings to identify novel constitutive relations for linear elasticity near the jamming transition. By forcing packings at varying wavelengths, we directly access their transverse and longitudinal compliances. These are found to be wavelength dependent, in violation of conventional (local) linear elasticity. Crossovers in the compliances select characteristic length scales, which signify the appearance of nonlocal effects. Two of these length scales diverge as the pressure vanishes, indicating that critical effects near jamming control the breakdown of local elasticity. We expect these nonlocal constitutive relations to be applicable to a wide range of weakly jammed solids, including emulsions, foams, and granulates.

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.

Microscopic Description of the Granular Fluidity Field in Nonlocal Flow Modeling

A recent granular rheology based on an implicit "granular fluidity" field has been shown to quantitatively predict many nonlocal phenomena. However, the physical nature of the field has not been identified. Here, the granular fluidity is found to be a kinematic variable given by the velocity fluctuation and packing fraction. This is verified with many discrete element simulations, which show that the operational fluidity definition, solutions of the fluidity model, and the proposed microscopic formula all agree. Kinetic theoretical and Eyring-like explanations shed insight into the obtained form.

Why granular media are thermal, and quite normal, after all

Two approaches exist to account for granular dynamics: The athermal one takes grains as elementary, the thermal one considers the total entropy that includes microscopic degrees of freedom such as phonons and electrons. Discrete element method (DEM), granular kinetic theory and athermal statistical mechanics (ASM) belong to the first, granular solid hydrodynamics (GSH) to the second one. A discussion of the conceptual differences between both is given here, leading, among others, to the following insights: 1) While DEM and granular kinetic theory are well justified to take grains as athermal, any entropic consideration is far less likely to succeed. 2) In addition to modeling grains as a gas of dissipative, rigid mass points, it is very helpful take grains as a thermal solid that has been sliced and diced. 3) General principles that appear invalid in granular media are repaired and restored once the true entropy is included. These abnormalities (such as invalidity of the fluctuation-dissipation theorem, granular temperatures failing to equilibrate, and grains at rest unable to explore the phase space) are consequences of the athermal approximation, not properties of granular media.

Structural and cooperative length scales in polymer gels

Understanding the relationship between the material structural details, the geometrical confining constraints, the local dynamical events and the global rheological response is at the core of present investigations on complex fluid properties. In the present article, this problem is addressed on a model yield stress fluid made of highly entangled polymer gels of Carbopol which follows at the macroscopic scale the well-known Herschel-Bulkley rheological law. First, performing local rheology measurements up to high shear rates ([Formula: see text] s^-1)and under confinement, we evidence unambiguously the breakdown of bulk rheology associated with cooperative processes under flow. Moreover, we show that these behaviors are fully captured with a unique cooperativity length [Formula: see text] over the whole range of experimental conditions. Second, we introduce an original optical microscopy method to access structural properties of the entangled polymer gel in the direct space. Performing image correlation spectroscopy of fluorophore-loaded gels, the characteristic size D of carbopol gels microstructure is determined as a function of preparation protocol. Combining both dynamical and structural information shows that the measured cooperative length [Formula: see text] corresponds to 2-5 times the underlying structural size D, thus providing a strong grounding to the "Shear Transformation Zones" modeling approach.

Unified rheology of vibro-fluidized dry granular media: From slow dense flows to fast gas-like regimes

Granular media take on great importance in industry and geophysics, posing a severe challenge to materials science. Their response properties elude known soft rheological models, even when the yield-stress discontinuity is blurred by vibro-fluidization. Here we propose a broad rheological scenario where average stress sums up a frictional contribution, generalizing conventional μ(I)-rheology, and a kinetic collisional term dominating at fast fluidization. Our conjecture fairly describes a wide series of experiments in a vibrofluidized vane setup, whose phenomenology includes velocity weakening, shear thinning, a discontinuous thinning transition, and gaseous shear thickening. The employed setup gives access to dynamic fluctuations, which exhibit a broad range of timescales. In the slow dense regime the frequency of cage-opening increases with stress and enhances, with respect to μ(I)-rheology, the decrease of viscosity. Diffusivity is exponential in the shear stress in both thinning and thickening regimes, with a huge growth near the transition.

Rheology of sediment transported by a laminar flow

Understanding the dynamics of fluid-driven sediment transport remains challenging, as it occurs at the interface between a granular material and a fluid flow. Boyer, Guazzelli, and Pouliquen [Phys. Rev. Lett. 107, 188301 (2011)]PRLTAO0031-900710.1103/PhysRevLett.107.188301 proposed a local rheology unifying dense dry-granular and viscous-suspension flows, but it has been validated only for neutrally buoyant particles in a confined and homogeneous system. Here we generalize the Boyer, Guazzelli, and Pouliquen model to account for the weight of a particle by addition of a pressure P_{0} and test the ability of this model to describe sediment transport in an idealized laboratory river. We subject a bed of settling plastic particles to a laminar-shear flow from above, and use refractive-index-matching to track particles' motion and determine local rheology-from the fluid-granular interface to deep in the granular bed. Data from all experiments collapse onto a single curve of friction μ as a function of the viscous number I_{v} over the range 3×10^{-5}≤I_{v}≤2, validating the local rheology model. For I_{v}<3×10^{-5}, however, data do not collapse. Instead of undergoing a jamming transition with μ→μ_{s} as expected, particles transition to a creeping regime where we observe a continuous decay of the friction coefficient μ≤μ_{s} as I_{v} decreases. The rheology of this creep regime cannot be described by the local model, and more work is needed to determine whether a nonlocal rheology model can be modified to account for our findings.

Intrusion rheology in grains and other flowable materials

The interaction of intruding objects with deformable materials arises in many contexts, including locomotion in fluids and loose media, impact and penetration problems, and geospace applications. Despite the complex constitutive behaviour of granular media, forces on arbitrarily shaped granular intruders are observed to obey surprisingly simple, yet empirical 'resistive force hypotheses'. The physics of this macroscale reduction, and how it might play out in other media, has however remained elusive. Here, we show that all resistive force hypotheses in grains arise from local frictional yielding, revealing a novel invariance within a class of plasticity models. This mechanical foundation, supported by numerical and experimental validations, leads to a general analytical criterion to determine which rheologies can obey resistive force hypotheses. We use it to explain why viscous fluids are observed to perform worse than grains, and to predict a new family of resistive-force-obeying materials: cohesive media such as pastes, gels and muds.

Criticality in Vibrated Frictional Flows at a Finite Strain Rate

We evidence critical fluctuations in the strain rate of granular flows that are weakly vibrated. Strikingly, the critical point arises at finite values of the mean strain rate and vibration strength, far from the yielding critical point at a zero flow rate. We show that the global rheology, as well as the amplitude and correlation time of the fluctuations, are consistent with a mean-field, Landau-like description, where the strain rate and the stress act as conjugated variables. We introduce a general model which captures the observed phenomenology and argue that this type of critical behavior generically arises when self-fluidization competes with friction.

A review on locomotion robophysics: the study of movement at the intersection of robotics, soft matter and dynamical systems

Discovery of fundamental principles which govern and limit effective locomotion (self-propulsion) is of intellectual interest and practical importance. Human technology has created robotic moving systems that excel in movement on and within environments of societal interest: paved roads, open air and water. However, such devices cannot yet robustly and efficiently navigate (as animals do) the enormous diversity of natural environments which might be of future interest for autonomous robots; examples include vertical surfaces like trees and cliffs, heterogeneous ground like desert rubble and brush, turbulent flows found near seashores, and deformable/flowable substrates like sand, mud and soil. In this review we argue for the creation of a physics of moving systems-a 'locomotion robophysics'-which we define as the pursuit of principles of self-generated motion. Robophysics can provide an important intellectual complement to the discipline of robotics, largely the domain of researchers from engineering and computer science. The essential idea is that we must complement the study of complex robots in complex situations with systematic study of simplified robotic devices in controlled laboratory settings and in simplified theoretical models. We must thus use the methods of physics to examine both locomotor successes and failures using parameter space exploration, systematic control, and techniques from dynamical systems. Using examples from our and others' research, we will discuss how such robophysical studies have begun to aid engineers in the creation of devices that have begun to achieve life-like locomotor abilities on and within complex environments, have inspired interesting physics questions in low dimensional dynamical systems, geometric mechanics and soft matter physics, and have been useful to develop models for biological locomotion in complex terrain. The rapidly decreasing cost of constructing robot models with easy access to significant computational power bodes well for scientists and engineers to engage in a discipline which can readily integrate experiment, theory and computation.

Eshelby inclusions in granular matter: Theory and simulations

We present a numerical implementation of an active inclusion in a granular material submitted to a biaxial test. We discuss the dependence of the response to this perturbation on two parameters: the intragranular friction coefficient on one hand, and the degree of the loading on the other hand. We compare the numerical results to theoretical predictions taking into account the change of volume of the inclusion as well as the anisotropy of the elastic matrix.

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.

Exponential velocity profile of granular flows down a confined heap

Thick granular flows are essential to many natural and industrial phenomena. Experimentally, it has been well established that the grain velocity profile is linear from the free surface to a certain depth, after which it decreases exponentially in the so-called "creep region". In this paper we obtain an exponential velocity profile based on the force balance of a grain near a wall, where the Janssen effect and the non-locality of interactions between grains are considered. When experimental parameters such as flow angles and friction coefficients are introduced in our model, it is able to reproduce experimental creep velocity profiles previously reported in the literature.

Yielding and flow of cellulose microfibril dispersions in the presence of a charged polymer

The shear flow of microfibrillated cellulose dispersions is still not wholly understood as a consequence of their multi-length-scale heterogeneity. We added carboxymethyl cellulose, a charged polymer, that makes cellulose microfibril dispersions more homogeneous at the submicron and macro scales. We then compared the yielding and flow behavior of these dispersions to that of typical thixotropic yield-stress fluids. Despite the apparent homogeneity of the dispersions, their flow velocity profiles in cone-plate geometry, as measured by rheo-MRI velocimetry, differ strongly from those observed for typical thixotropic model systems: the viscosity across the gap is not uniform, despite a flat stress field across the gap. We describe these velocity profiles with a nonlocal model, and attribute the non-locality to persistent micron-scale structural heterogeneity.

Steady flow dynamics during granular impact

We study experimentally and computationally the dynamics of granular flow during impacts where intruders strike a collection of disks from above. In the regime where granular force dynamics are much more rapid than the intruder motion, we find that the particle flow near the intruder is proportional to the instantaneous intruder speed; it is essentially constant when normalized by that speed. The granular flow is nearly divergence free and remains in balance with the intruder, despite the latter's rapid deceleration. Simulations indicate that this observation is insensitive to grain properties, which can be explained by the separation of time scales between intergrain force dynamics and intruder dynamics. Assuming there is a comparable separation of time scales, we expect that our results are applicable to a broad class of dynamic or transient granular flows. Our results suggest that descriptions of static-in-time granular flows might be extended or modified to describe these dynamic flows. Additionally, we find that accurate grain-grain interactions are not necessary to correctly capture the granular flow in this regime.

Unsteady granular flows down an inclined plane

The continuum description of granular flows is still a challenge despite their importance in many geophysical and industrial applications. We extend previous works, which have explored steady flow properties, by focusing on unsteady flows accelerating or decelerating down an inclined plane in the simple shear configuration. We solve the flow kinematics analytically, including predictions of evolving velocity and stress profiles and the duration of the transient stage. The solution shows why and how granular materials reach steady flow on slopes steeper than the angle of repose and how they decelerate on shallower slopes. The model might facilitate development of natural hazard assessment and may be modified in the future to explore unsteady granular flows in different configurations.

Anomalous energy cascades in dense granular materials yielding under simple shear deformations

By using molecular dynamics (MD) simulations of dense granular particles in two dimensions, we study turbulent-like structures of their non-affine velocities under simple shear deformations. We find that the spectrum of non-affine velocities, introduced as an analog of the energy spectrum for turbulent flows, exhibits the power-law decay if the system is yielding in a quasi-static regime, where large-scale collective motions and inelastic interactions of granular particles are crucial for the anomalous cascade of kinetic energy. Based on hydrodynamic equations of dense granular materials, which include both kinetic and contact contributions in constitutive relations, we derive a theoretical expression for the spectrum, where a good agreement between the result of MD simulations and theoretical prediction is established over a wide range of length scales.

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.

Non-local rheology in dense granular flows: Revisiting the concept of fluidity

The aim of this article is to discuss the concepts of non-local rheology and fluidity, recently introduced to describe dense granular flows. We review and compare various approaches based on different constitutive relations and choices for the fluidity parameter, focusing on the kinetic elasto-plastic model introduced by Bocquet et al. (Phys. Rev. Lett 103, 036001 (2009)) for soft matter, and adapted for granular matter by Kamrin et al. (Phys. Rev. Lett. 108, 178301 (2012)), and the gradient expansion of the local rheology μ(I) that we have proposed (Phys. Rev. Lett. 111, 238301 (2013)). We emphasise that, to discriminate between these approaches, one has to go beyond the predictions derived from linearisation around a uniform stress profile, such as that obtained in a simple shear cell. We argue that future tests can be based on the nature of the chosen fluidity parameter, and the related boundary conditions, as well as the hypothesis made to derive the models and the dynamical mechanisms underlying their dynamics.

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.

Effective Wall Friction in Wall-Bounded 3D Dense Granular Flows

We report numerical simulations on granular shear flows confined between two flat but frictional sidewalls. Novel regimes differing by their strain localization features are observed. They originate from the competition between dissipation at the sidewalls and dissipation in the bulk of the flow. The effective friction at sidewalls is characterized (effective friction coefficient and orientation of the friction force) for each regime, and its interdependence with slip and force fluctuations is pointed out. We propose a simple scaling law linking the slip velocity to the granular temperature in the main flow direction which leads naturally to another scaling law for the effective friction.

Anisotropy of weakly vibrated granular flows

We experimentally probe the anisotropy of weakly vibrated flowing granular media. Depending on the driving parameters-flow rate and vibration strength-this anisotropy varies significantly. We show how the anisotropy collapses when plotted as a function of the driving stresses, uncovering a direct link between stresses and anisotropy. Moreover, our data suggest that for small anisotropies, the shear stresses vanish. Anisotropy of the fabric of granular media thus plays a crucial role in determining the rheology of granular flows.

Measurement of Stress Redistribution in Flowing Emulsions

We study how local rearrangements alter droplet stresses within flowing dense quasi-two-dimensional emulsions at area fractions ϕ≥0.88. Using microscopy, we measure droplet positions while simultaneously using their deformed shape to measure droplet stresses. We find that rearrangements alter nearby stresses in a quadrupolar pattern: stresses on neighboring droplets tend to either decrease or increase depending on location. The stress redistribution is more anisotropic with increasing ϕ. The spatial character of the stress redistribution influences where subsequent rearrangements occur. Our results provide direct quantitative support for rheological theories of dense amorphous materials that connect local rearrangements to changes in nearby stress.

Hydrodynamic instabilities in shear flows of dry cohesive granular particles

We extend the dynamic van der Waals model introduced by A. Onuki [Phys. Rev. Lett., 2005, 94, 054501] to the description of cohesive granular flows under a plane shear to study their hydrodynamic instabilities. By numerically solving the dynamic van der Waals model, we observed various heterogeneous structures of density fields in steady states, where the viscous heating is balanced with the energy dissipation caused by inelastic collisions. Based on the linear stability analysis, we found that the spatial structures are determined by the mean volume fraction, the applied shear rate, and the inelasticity, where the instability is triggered if the system is thermodynamically unstable, i.e. the pressure, p, and the volume fraction, ϕ, satisfy ∂p/∂ϕ < 0.

Shear-Rate-Independent Diffusion in Granular Flows

We computationally study the behavior of the diffusion coefficient D in granular flows of monodisperse and bidisperse particles spanning regions of relatively high and low shear rate in open and closed laterally confined heaps. Measurements of D at various flow rates, streamwise positions, and depths collapse onto a single curve when plotted as a function of γd2, where d is the local mean particle diameter and γ is the local shear rate. When γ is large, D is proportional to γd2, as in previous studies. However, for γd2 below a critical value, D is independent of γd2. The acceleration due to gravity g and particle stiffness (or, equivalently, the binary collision time t(c)) together determine the transition in D between regimes. This suggests that while shear rate and particle size determine diffusion at relatively high shear rates in surface-driven flows, diffusion at low shear rates is an elastic phenomenon with time and length scales dependent on gravity (sqrt d/g) and particle stiffness (t(c)sqrt(dg), respectively.

Mechanical fluctuations suppress the threshold of soft-glassy solids: The secular drift scenario

We propose a dynamical mechanism leading to the fluidization by external mechanical fluctuations of soft-glassy amorphous material driven below the yield stress. The model is based on the combination of memory effect and nonlinearity, leading to an accumulation of tiny effects over a long term. We test this scenario on a granular packing driven mechanically below the Coulomb threshold. We provide evidence for an effective viscous response directly related to small stress modulations in agreement with the theoretical prediction of a generic secular drift. We propose to extend this result more generally to a large class of glassy systems.

Steady shearing flows of deformable, inelastic spheres

We extend models for granular flows based on the kinetic theory beyond the critical volume fraction at which a rate-independent contribution to the stresses develops. This involves the incorporation of a measure of the duration of the particle interaction before and after this volume fraction. At volume fractions less than the critical, the stress components contain contributions from momentum exchanged in collisions that are influenced by the particle elasticity. At volume fractions greater than the critical, the stress components contain both static contributions from particle elasticity and dynamic contributions from the momentum transfer associated with the release of elastic energy by the breaking of force chains. A simple expression for the duration of a collision before and after the critical volume fraction permits a smooth transition between the two regimes and predictions for the components of the stress in steady, homogeneous shearing that are in good agreement with the results of numerical simulations. Application of the theory to steady, inhomogeneous flows reproduces the features of such flows seen in numerical simulations and physical experiments.

Experimental observation of local rearrangements in dense quasi-two-dimensional emulsion flow

We experimentally study rearranging regions in slow athermal flow by observing the flow of a concentrated oil-in-water emulsion in a thin chamber with a constricting hopper shape. The gap of the chamber is smaller than the droplet diameters, so that the droplets are compressed into quasi-two-dimensional pancakes. We focus on localized rearrangements known as "T1 events" where four droplets exchange neighbors. Flowing droplets are deformed due to forces from neighboring droplets, and these deformations are decreased by nearby T1 events, with a spatial dependence related to the local structure. We see a tendency of the T1 events to occur in small clusters.

Nonlocal effects in sand flows on an inclined plane

The flow of sand on a rough inclined plane is investigated experimentally. We directly show that a jammed layer of grains spontaneously forms below the avalanche. Its properties and its relation with the rheology of the flowing layer of grains are presented and discussed. In a second part, we study the dynamics of erosion and deposition solitary waves in the domain where they are transversally stable. We characterize their shapes and velocity profiles. We relate their translational velocity to the stopping height and to the mass trapped in the avalanche. Finally, we use the velocity profile to get insight into the rheology very close to the jamming limit.

Applying GSH to a wide range of experiments in granular media

Granular solid hydrodynamics (GSH) is a continuum-mechanical theory for granular media, whose wide range of applicability is shown in this paper. Simple, frequently analytic solutions are related to classic observations at different shear rates, including: i) static stress distribution, clogging; ii) elasto-plastic motion: loading and unloading, approach to the critical state, angle of stability and repose; iii) rapid dense flow: the μ-rheology, Bagnold scaling and the stress minimum; iv) elastic waves, compaction, wide and narrow shear band. Less conventional experiments have also been considered: shear jamming, creep flow, visco-elastic behavior and non-local fluidization. With all these phenomena ordered, related, explained and accounted for, though frequently qualitatively, we believe that GSH may be taken as a unifying framework, providing the appropriate macroscopic vocabulary and mindset that help one coming to terms with the breadth of granular physics.

Nonlocal modeling of granular flows down inclines

Flows of granular media down a rough inclined plane demonstrate a number of nonlocal phenomena. We apply the recently proposed nonlocal granular fluidity model to this geometry and find that the model captures many of these effects. Utilizing the model's dynamical form, we obtain a formula for the critical stopping height of a layer of grains on an inclined surface. Using an existing parameter calibration for glass beads, the theoretical result compares quantitatively to existing experimental data for glass beads. This provides a stringent test of the model, whose previous validations focused on driven steady-flow problems. For layers thicker than the stopping height, the theoretical flow profiles display a thickness-dependent shape whose features are in agreement with previous discrete particle simulations. We also address the issue of the Froude number of the flows, which has been shown experimentally to collapse as a function of the ratio of layer thickness to stopping height. While the collapse is not obvious, two explanations emerge leading to a revisiting of the history of inertial rheology, which the nonlocal model references for its homogeneous flow response.

Fragmentation processes in two-phase materials

We investigate the fragmentation process of solid materials with crystalline and amorphous phases using the the discrete element method. Damage initiates inside spherical samples above the contact zone in a region where the circumferential stress field is tensile. Cracks initiated in this region grow to form meridional planes. If the collision energy exceeds a critical value which depends on the material's internal structure, cracks reach the sample surface resulting in fragmentation. We show that this primary fragmentation mechanism is very robust with respect to the internal structure of the material. For all configurations, a sharp transition from the damage to the fragmentation regime is observed, with smaller critical collision energies for crystalline samples. The mass distribution of the fragments follows a power law for small fragments with an exponent that is characteristic for the branching merging process of unstable cracks. Moreover this exponent depends only on the dimensionally of the system and not on the microstructure.

Continuum modeling of secondary rheology in dense granular materials

Recent dense granular flow experiments have shown that shear deformation in one region of a granular medium fluidizes its entirety, including regions far from the sheared zone, effectively erasing the yield condition everywhere. This enables slow creep deformation to occur when an external force is applied to a probe in the nominally static regions of the material. The apparent change in rheology induced by far-away motion is termed the "secondary rheology," and a theoretical rationalization of this phenomenon is needed. Recently, a new nonlocal granular rheology was successfully used to predict steady granular flow fields, including grain-size-dependent shear-band widths in a wide variety of flow configurations. We show that the nonlocal fluidity model is also capable of capturing secondary rheology. Specifically, we explore creep of a circular intruder in a two-dimensional annular Couette cell and show that the model captures all salient features observed in experiments, including both the rate-independent nature of creep for sufficiently slow driving rates and the faster-than-linear increase in the creep speed with the force applied to the intruder.

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.

Poiseuille flow of soft glasses in narrow channels: from quiescence to steady state

Using numerical simulations, the onset of Poiseuille flow in a confined soft glass is investigated. Starting from the quiescent state, steady flow sets in at a time scale which increases with a decrease in applied forcing. At this onset time scale, a rapid transition occurs via the simultaneous fluidization of regions having different local stresses. In the absence of steady flow at long times, creep is observed even in regions where the local stress is larger than the bulk yielding threshold. Finally, we show that the time scale to attain steady flow depends strongly on the history of the initial state.

Boundary conditions for soft glassy flows: slippage and surface fluidization

We explore the question of surface boundary conditions for the flow of a dense emulsion. We make use of microlithographic tools to create surfaces with well controlled roughness patterns and measure using dynamic confocal microscopy both the slip velocity and the shear rate close to the wall, which we relate to the notion of surface fluidization. Both slippage and wall fluidization depend non-monotonously on the roughness. We interpret this behavior within a simple model in terms of the building of a stratified layer and the activation of plastic events by the surface roughness.