Transition from the viscous to inertial regime in dense suspensions

Nonmonotonic Constitutive Curves and Shear Banding in Dry and Wet Granular Flows

We use particle simulations to map comprehensively the shear rheology of dry and wet granular matter comprising particles of finite stiffness, in both fixed pressure and fixed volume protocols. At fixed pressure we find nonmonotonic constitutive curves that are shear thinning, whereas at fixed volume we find nonmonotonic constitutive curves that are shear thickening. We show that the presence of one nonmonotonicity does not imply the other. Instead, there exists a signature in the volume fraction measured under fixed pressure that, when present, ensures nonmonotonic constitutive curves at fixed volume. In the context of dry granular flow we show that gradient and vorticity bands arise under fixed pressure and volume, respectively, as implied by the constitutive curves. For wet systems our results are consistent with a recent experimental observation of shear thinning at fixed pressure. We furthermore predict discontinuous shear thickening in the absence of critical load friction.

Compaction Evolution and Mechanisms of Granular Materials Due to Gyratory Shearing

Granular systems, no matter whether they are dry or saturated, are commonly encountered in both natural scenarios and engineering applications. In this work, we tackle the compaction problem of both dry and saturated granular systems under gyratory shearing compaction, where particles are subjected to constant pressure and continuous shear rate, which is quite different from the traditional cyclic shearing compaction. Such phenomena are crucial to the compaction of asphalt mixtures or soils in civil engineering and can be extended to other areas, such as powder processing and pharmaceutical engineering. In this study, we investigated the behavior of both dry and fully saturated mono-dispersed granular materials under gyratory shearing compaction using the discrete element method (DEM) and found that the gyratory speed or interstitial fluid viscosity has almost no impact on the compaction behavior, while the pressure and the particle size play more important roles. Additionally, it is the inertial time scale which dictates the compaction behavior under gyratory shearing in most cases; meanwhile, the viscous time scale can also have influence in some conditions. This work determines the similarity and unity between the granular gyratory compaction and the rheology of granular systems, which has direct relevance to various natural and engineering systems.

Rheology of Suspensions of Non-Brownian Soft Spheres across the Jamming and Viscous-to-Inertial Transitions

The rheology of suspensions of non-Brownian soft spheres is studied across jamming but also across the viscous and inertial regimes using a custom pressure- and volume-imposed rheometer. The study shows that the granular rheology found for suspensions of hard spheres can be extended to a soft granular rheology (SGranR) by renormalizing the critical volume fraction and friction coefficient to pressure-dependent values and using the addition of the viscous and inertial stress scales. This SGranR encompasses rheological behaviors on both sides of the jamming transition, resulting in an approximate collapse of the rheological data into two branches when scaled with the distance to jamming, as observed for soft colloids. This research suggests that suspensions of soft particles across flow regimes can be described by a unified SGranR framework around the jamming transition.

Order-disorder transition during shear thickening in bidisperse dense suspensions

Shear thickening of multimodal suspensions has proven difficult to understand because the rheology depends largely on the microscopic details of stress-induced frictional contacts at different particle size distributions (PSDs). Our discrete particle simulations below a critical volume fraction ϕc over a broad range of shear rates and PSDs elucidate the basic mechanism of order-disorder transition. Around the theoretical optimal PSD (relative content of small particles ζ1= 0.26), particles order into a layered structure in the Newtonian regime. At the onset of shear thickening, this layered structure transforms to a disordered one, accompanied by an abrupt viscosity jump. Minor increase in large-large particle contacts after the order-disorder transition causes apparent increase in radial force along the compressional axis. Bidisperse suspensions with less regular but stable layered structure at ζ1= 0.50 show good fluidity in the shear thickening regime. This work shows that in inertial flows where particle collisions dominate, order-disorder transition could play an essential role in shear thickening for bidisperse suspensions.

Collapse dynamics and deposition morphology of low-viscocohesive granular columns on a rough horizontal surface

Using the three-dimensional discrete element method, we numerically investigate the collapse dynamics and deposition morphology of low-viscocohesive granular columns on a rough-horizontal plane by systematically varying a broad range of values of the initial column aspect ratio, cohesive stress, and liquid viscosity. The results show that the kinetic energy, half runout time, and runout distance increase with increasing the initial column aspect ratio but decrease with increasing the cohesive and viscous effects of the binding liquid, while the toe angle and deposit height decrease with increasing the aspect ratio and increase with increasing cohesive stress and liquid viscosity. Remarkably, by defining a dimensionless scaling number that incorporates the Bond number and initial column aspect ratio, this allows us to nicely describe the kinetic energy, half runout time, deposition height, runout distance, and toe angle. These unified descriptions may provide insights into the physical properties of the collapse dynamics and deposition morphology of low-viscocohesive granular columns, leading to good explanations of the complex properties of natural disaster events.

Slow and fast particles in shear-driven jamming: Critical behavior

We do extensive simulations of a simple model of shear-driven jamming in two dimensions to determine and analyze the velocity distribution at different densities ϕ around the jamming density ϕ_{J} and at different low shear strain rates, γ[over ̇]. We then find that the velocity distribution is made up of two parts which are related to two different physical processes which we call the slow process and the fast process as they are dominated by the slower and the faster particles, respectively. Earlier scaling analyses have shown that the shear viscosity η, which diverges as the jamming density is approached from below, consists of two different terms, and we present strong evidence that these terms are related to the two different processes: the leading divergence is due to the fast process, whereas the correction-to-scaling term is due to the slow process. The analysis of the slow process is possible thanks to the observation that the velocity distribution for different γ[over ̇] and ϕ at and around the shear-driven jamming transition has a peak at low velocities and that the distribution has a constant shape up to and slightly above this peak. We then find that it is possible to express the contribution to the shear viscosity due to the slow process in terms of height and position of the peak in the velocity distribution and find that this contribution matches the correction-to-scaling term, determined through a standard critical scaling analysis. A further observation is that the collective particle motion is dominated by the slow process. In contrast to the usual picture in critical phenomena with a direct link between the diverging correlation length and a diverging order parameter, we find that correlations and shear viscosity decouple since they are controlled by different sets of particles and that shear-driven jamming is thus an unusual kind of critical phenomenon.

Non-Newtonian fluid gating membranes with acoustically responsive and self-protective gas transport control

Control of gas transport through porous media is desired in multifarious processes such as chemical reactions, interface absorption, and medical treatment. Liquid gating technology, based on dynamically adaptive interfaces, has been developed in recent years and has shown excellent control capability in gas manipulation-the reversible opening and closing of a liquid gate for gas transport as the applied pressure changes. Here, we report a new strategy to achieve self-protective gas transport control by regulating the dynamic porous interface in a non-Newtonian fluid gating membrane based on the shear thickening fluid. The gas transport process can be suspended and restored via modulation of the acoustic field, owing to the transition of particle-to-particle interactions in a confined geometry. Our experimental and theoretical results support the stability and tunability of the gas transport control. In addition, relying on the shear thickening behaviour of the gating fluid, the transient response can be achieved to resist high-impact pressure. This strategy could be utilized to design integrated smart materials used in complex and extreme environments such as hazardous and explosive gas transportation.

Influence of fine particle content in debris flows on alluvial fan morphology

Alluvial fans are large-scale depositional structures commonly found at the base of mountain ranges. They are relatively soil-rich compared to the rocky terrains, or catchment areas, from which their material originates. When frequented by debris flows (massive, muddy, rocky flows) they contribute significantly to local hazards as they carry focused, collisional, fast-moving materials across alluvial fans, unpredictable in size, speed, and direction. We research how fine particle content in debris flows correlates with directional changes, i.e., debris flow avulsions. Toward this, we analyzed field data from two neighboring alluvial fans in the White Mountains (California, USA) that exhibit dramatically different topographies despite their proximity and associated similar long-term climates. Informed by these measurements, we performed long-term and incremental alluvial fan experiments built by debris flows with systematically-varied fine particle content. We found that (1) decreasing fine particle content increases the variability of fan slopes and associated channelization dynamics, and (2) for all mixtures longer-term continuous alluvial fan experiments form more complex surface channelizations than repeated flows for the same total time, indicating the importance of both particle sizes and timescales on alluvial fan surface morphology.

Viscous to Inertial Transition in Dense Granular Suspension

Granular suspensions present a transition from a Newtonian rheology in the Stokes limit to a Bagnoldian rheology when inertia is increased. A custom rheometer that can be run in a pressure- or a volume-imposed mode is used to examine this transition in the dense regime close to jamming. By varying systematically the interstitial fluid, shear rate, and packing fraction in volume-imposed measurements, we show that the transition takes place at a Stokes number of 10 independent of the packing fraction. Using pressure-imposed rheometry, we investigate whether the inertial and viscous regimes can be unified as a function of a single dimensionless number based on stress additivity.

Mass Changes the Diffusion Coefficient of Particles with Ligand-Receptor Contacts in the Overdamped Limit

Inertia does not generally affect the long-time diffusion of passive overdamped particles in fluids. Yet a model starting from the Langevin equation predicts a surprising property of particles coated with ligands that bind reversibly to surface receptors: heavy particles diffuse more slowly than light ones of the same size. We show this by simulation and by deriving an analytic formula for the mass-dependent diffusion coefficient in the overdamped limit. We estimate the magnitude of this effect for a range of biophysical ligand-receptor systems, and find it is potentially observable for tailored micronscale DNA-coated colloids.

Granular porous landslide tsunami modelling - the 2014 Lake Askja flank collapse

Subaerial landslides and volcano flank collapses can generate tsunamis with devastating consequences. The lack of comprehensive models incorporating both the landslide and the wave mechanics represents a gap in providing consistent predictions of real events. Here, we present a novel three-dimensional granular landslide and tsunami model and apply it to the 2014 Lake Askja landslide tsunami. For the first time, we consistently simulate small-scale laboratory experiments as well as full scale catastrophic events with the same model. The model captures the complete event chain from the landslide dynamics to the wave generation and inundation. Unique and complete field data, along with the limited geographic extent of Lake Askja enabled a rigorous validation. The model gives deep insights into the physical landslide processes and improves our understanding and prediction capabilities of frequent and catastrophic landslide tsunamis.

Tsunamis are devastating events. They are especially difficult to predict, when generated by landslides. In this paper, the authors overcome this issue by modelling the landslide and the tsunami in a unified framework in unprecedented detail.

Directional shear jamming of frictionless ellipses

In this work we study shear reversals of dense non-Brownian suspensions composed of cohesionless elliptical particles. By numerical simulations, we show that a new fragility appears for frictionless ellipses in the flowing states, where particles can flow indefinitely in one direction at applied shear stresses but shear jam in the other direction upon shear stress reversal. This new fragility, absent in the isotropic particle case, is linked to the directional order of the elongated particles at steady shear and its reorientation at shear stress reversal, which forces the suspensions to pass through a more disordered state with an increased number of contacts in which it might get arrested.

Granular-slurry rheology and asphalt compaction

Hot mixed asphalt (HMA) is a mixture of particles (coarse and fine aggregates) and interstitial fluid (asphalt binder) designed to compact and harden for long-lasting roads. In this study, we implement a two-scale approach to capture the compaction behaviour of hot asphalt mixtures using both a granular-slurry rheology (GSR) at a smaller scale and a discrete element method (DEM) simulation at the scale of a compactor. We show that this modelling effort captures the compaction of HMA with different binder viscosities modified by adding graphene nano-platelets (GNP). This research has the capacity to shed light on how the properties of mixture components can influence compaction efficiency and effectiveness.

Scaling behavior of the tensile strength of viscocohesive granular aggregates

We numerically analyze the tensile strength of a single wet agglomerate modeled as a viscocohesive aggregate impacting a flat surface by using the discrete-element simulations. The viscocohesive agglomerate composed of primary spherical particles with the inclusion of the interstitial liquid in the form of the capillary bridges characterized by the cohesive and viscous forces between particles is extracted from a cuboidal sample of granular materials by applying a spherical probe. The tensile strength is measured from the impact test of a wet agglomerate by systematically varying different values of the surface tension of the interstitial liquid, the liquid viscosity, and the impact speed. We show that the tensile stress increases immediately when the collision occurs between the agglomerate and the flat surface. The peak of the tensile stress obtained after the collision, then decreases smoothly with increasing the particle movement. The maximum tensile stress is defined to be the tensile strength of such agglomerate. It is remarkable that the normalized tensile strength of such agglomerate can be well described as a function of a dimensionless impact number that incorporates the capillary number and Stokes number (calculated from the surface tension and the viscosity of the liquid and the impact rate of the agglomerate), thus providing the confirmation for the unified representation of the liquid properties and the impact rate of wet granular media.

Time scales and rheology of visco-cohesive granular flows

In the presence of viscous and cohesive interactions between particles, a granular flow is governed by several characteristic time and stress scales that determine its rheological properties (shear stress, packing fraction, effective viscosities). In this paper, we revisit and extend the scaling arguments used previously for dry cohesionless granular flows and suspensions. We show that the rheology can be in principle described by a single dimensionless control parameter that includes all characteristic times. We also briefly present simulation results for 2D sheared suspensions and 3D wet granular flows where the effective friction coefficient and packing fraction are consistently described as functions of this unique control parameter.

Transition from steady shear to oscillatory shear rheology of dense suspensions

Recent studies have highlighted that oscillatory and time-dependent shear flows might help increase the flowability of dense suspensions. While most focus has been on cross-flows we here study a simple two-dimensional suspensions where we apply simultaneously oscillatory and stationary shear along the same direction. We first show that the dissipative viscosities in this set-up significantly decrease with an increasing shear-rate magnitude of the oscillations and given that the oscillatory strain is small, in a similar fashion as found previously for cross-flow oscillations. As for cross-flow oscillations, the decrease can be attributed to the large decrease in the number of contacts and an altered microstructure as one transitions from a steady shear to an oscillatory shear dominated rheology. As subresults we find both an extension to the μ(J) rheology, a constitutive relationship between the shear stresses and the shear rate, valid for oscillatory shear flows and that shear-jamming of frictional particles at oscillatory shear dominated flows occurs at higher packing fractions compared to steady shear dominated flows.

The role of inter-particle friction on rheology and texture of wet granular flows

In order to get insight into the rheology and texture of rough unsaturated granular flows, we study the effects of the inter-particle friction coefficient on the macroscopic attributes and the texture variables of steady-state shearing flow of wet granular materials by relying on three-dimensional (3D) particle dynamics simulations. The macroscopic attributes are characterized by the macroscopic friction coefficient, macroscopic cohesion, and the packing fraction. The microstructural variables are characterized by the fabric and force anisotropies, the coordination number, and the stress transmission ratio. We show that the macroscopic observables behave as a function of the inertial number as a dry case for different values of the inter-particle friction coefficient. In particular, the macroscopic friction coefficient increases and the packing fraction decreases rapidly for small values of the friction coefficient, then they almost reach plateaus for higher values of the friction coefficient. Interestingly, all the macroscopic observables nicely behave as a function of the small values of the friction coefficient. Similarly, we also observe these characteristics for the fabric and force anisotropies and the coordination number as well as the stress transmission ratio which reflects the intermediate relationship between the microstructure and the mechanical behavior of such flows.

Shear-Induced Heteroaggregation of Oppositely Charged Colloidal Particles

This paper investigates numerically the shear-induced aggregation of mixed populations of colloidal particles leading to the formation of clusters. Suspensions with different amounts of positively and negatively charged colloidal particles are simulated. To resolve the aggregation kinetics and structural properties of the formed clusters, we resort to a mixed deterministic-stochastic simulation method. The method is built on a combination of a Monte Carlo algorithm to sample a statistically expected sequence of encounter events between the suspended particles and a discrete element method built in the framework of Stokesian dynamics to simulate the encounters in a fully predictive manner. Results reveal a strong influence of the composition of the population on both the aggregation kinetics and the aggregate structure. In particular, we observe a size-stabilization phenomenon taking place in the suspension when the relative concentration of the majority particles lies in the range 80-85%; i.e., starting from primary particles, after a short growth period, we observed a cessation of aggregation. Inspection of the aggregate morphology shows that the formed clusters are composed of few minority particles placed in the inner region, while the aggregate surface is covered by majority particles, acting to provide a shielding effect against further growth.

Soft granular particles sheared at a controlled volume: rate-dependent dynamics and the solid-fluid transition

We study the responses of fluid-immersed soft hydrogel spheres that are sheared under controlled volume fractions. Slippery, deformable particles along with the density-matched interstitial fluid are sandwiched between two opposing rough cones, allowing studies for a wide range of volume fraction φ both above and below the jamming of granular suspension. We utilize sudden cessations of shearing, accompanied by refraction-matched internal imaging, to supplement the conventional flow-curve measurements. At sufficiently high volume fractions, the settling of particles after the cessations exhibits a continuous yet distinct transition over the change of the shear rate. Such changes back out the qualitative difference in the state of flowing prior to the cessations: the quasi-static yielding of a tightly packed network, as opposed to the rapid sliding of particles mediated by the interstitial fluid whose dynamics depends on the driving rate. In addition, we determine the solid-fluid transition using two independent methods: the extrapolation of stress residues and the estimated yield stress from high values of φ, and the settling of particles upon shear cessations as φ goes across the transition. We also verify the power law on values of characteristic stress with respect to the distance from jamming φ - φc, with an exponent close to 2. These results demonstrate a multitude of relaxation timescales behind the dynamics of soft particles, and raise questions on how we extend the existing paradigms of the flow of a densely packed system when the softness is actively involved.

Mobility in immersed granular materials upon cyclic loading

We study the mobility of objects embedded in an immersed granular packing and subjected to cyclic loadings. With this aim, we conducted uplift experiments whereby a horizontal plate is embedded in the packing and subjected to a vertical cyclic force oscillating between zero and a maximum amplitude. Tests performed at different cyclic force frequencies and amplitudes evidence the development of three mobility regimes whereby the plate stays virtually immobile, moves up steadily, or slowly creeps upwards. Results show that steady plate uplift can occur at lower force magnitudes when the frequency is increased. We propose an interpretation of this frequency-weakening behavior based on force relaxation experiments and on the analysis of the mobility response of theoretical viscoelastoplastic mechanical analog. These results and analysis point out inherent differences in mobility response between steady and cyclic loadings in immersed granular materials.

Shear jamming, discontinuous shear thickening, and fragile states in dry granular materials under oscillatory shear

We numerically study the linear response of two-dimensional frictional granular materials under oscillatory shear. The storage modulus G^{'} and the loss modulus G^{''} in the zero strain rate limit depend on the initial strain amplitude of the oscillatory shear before measurement. The shear jammed state (satisfying G^{'}>0) can be observed at an amplitude greater than a critical initial strain amplitude. The fragile state is defined by the emergence of liquid-like and solid-like states depending on the form of the initial shear. In this state, the observed G^{'} after the reduction of the strain amplitude depends on the phase of the external shear strain. The loss modulus G^{''} exhibits a discontinuous jump corresponding to discontinuous shear thickening in the fragile state.

Stress Relaxation above and below the Jamming Transition

We numerically investigate stress relaxation in soft athermal disks to reveal critical slowing down when the system approaches the jamming point. The exponents describing the divergence of the relaxation time differ dramatically depending on whether the transition is approached from the jammed or unjammed phase. This contrasts sharply with conventional dynamic critical scaling scenarios, where a single exponent characterizes both sides. We explain this surprising difference in terms of the vibrational density of states, which is a key ingredient of linear viscoelastic theory. The vibrational density of states exhibits an extra slow mode that emerges below jamming, which we utilize to demonstrate the anomalous exponent below jamming.

Additive rheology of complex granular flows

Granular flows are omnipresent in nature and industrial processes, but their rheological properties such as apparent friction and packing fraction are still elusive when inertial, cohesive and viscous interactions occur between particles in addition to frictional and elastic forces. Here we report on extensive particle dynamics simulations of such complex flows for a model granular system composed of perfectly rigid particles. We show that, when the apparent friction and packing fraction are normalized by their cohesion-dependent quasistatic values, they are governed by a single dimensionless number that, by virtue of stress additivity, accounts for all interactions. We also find that this dimensionless parameter, as a generalized inertial number, describes the texture variables such as the bond network connectivity and anisotropy. Encompassing various stress sources, this unified framework considerably simplifies and extends the modeling scope for granular dynamics, with potential applications to powder technology and natural flows.

Granular materials are abundant in nature, but we haven’t fully understood their rheological properties as complex interactions between particles are involved. Here, Vo et al. show that granular flows can be described by a generalized dimensionless number based on stress additivity.

Discrete Element Simulation and Validation of a Mixing Process of Granular Materials

The mixing processes of granular materials have gained wide interest among various fields of science and engineering. In this study, our focus is a mixing process for offshore mining. We conducted numerical simulations using the discrete element method (DEM) in comparison with experimental works on mixing color sand. Careful calibration of initial packing densities has been performed for the simulations. For validation, the steady-state torques on the mixer head, the maximal increase of surface height after mixing, and the surface mixing patterns have been compared. The effect of particle size on the simulation results has been clarified. With the particle size approaching the actual particle size, consistent torques and mixing patterns indicate the capability of the DEM code for studying the particular mixing process, while the results for the maximal increase of surface height should be interpreted with more caution.

Shear-driven flow of athermal, frictionless, spherocylinder suspensions in two dimensions: Particle rotations and orientational ordering

We use numerical simulations to study the flow of a bidisperse mixture of athermal, frictionless, soft-core two-dimensional spherocylinders driven by a uniform steady-state simple shear applied at a fixed volume and a fixed finite strain rate γ[over ̇]. Energy dissipation is via a viscous drag with respect to a uniformly sheared host fluid, giving a simple model for flow in a non-Brownian suspension with Newtonian rheology. Considering a range of packing fractions ϕ and particle asphericities α at small γ[over ̇], we study the angular rotation θ[over ̇]_{i} and the nematic orientational ordering S_{2} of the particles induced by the shear flow, finding a nonmonotonic behavior as the packing ϕ is varied. We interpret this nonmonotonic behavior as a crossover from dilute systems at small ϕ, where single-particle-like behavior occurs, to dense systems at large ϕ, where the geometry of the dense packing dominates and a random Poisson-like process for particle rotations results. We also argue that the finite nematic ordering S_{2} is a consequence of the shearing serving as an ordering field, rather than a result of long-range cooperative behavior among the particles. We arrive at these conclusions by consideration of (i) the distribution of waiting times for a particle to rotate by π, (ii) the behavior of the system under pure, as compared to simple, shearing, (iii) the relaxation of the nematic order parameter S_{2} when perturbed away from the steady state, and (iv) by construction, a numerical mean-field model for the rotational motion of a particle. Our results also help to explain the singular behavior observed when taking the α→0 limit approaching circular disks.

Evolution of wet agglomerates inside inertial shear flow of dry granular materials

We use particle dynamics simulations to investigate the evolution of a wet agglomerate inside homogeneous shear flows of dry particles. The agglomerate is modeled by introducing approximate analytical expressions of capillary and viscous forces between particles in addition to frictional contacts. During shear flow, the agglomerate may elongate, break, or be eroded by loss of its capillary bonds and primary particles. By systematically varying the shear rate and surface tension of the binding liquid, we characterize the rates of these dispersion modes. All the rates increase with increasing inertial number of the flow and decreasing cohesion index of the agglomerate. We show that the data points for each mode collapse on a master curve for a dimensionless scaling parameter that combines the inertial number and the cohesion index. The erosion rate vanishes below a cutoff value of the scaling parameter. This leads to a power-law borderline between the vanishing erosion states and erosion states in the phase space defined by the inertial number and the cohesion index.

Effective drag of a rod in fluid-saturated granular beds

We measure the drag encountered by a vertically oriented rod moving across a sedimented granular bed immersed in a fluid under steady-state conditions. At low rod speeds, the presence of the fluid leads to a lower drag because of buoyancy, whereas a significantly higher drag is observed with increasing speeds. The drag as a function of the depth is observed to decrease from being quadratic at low speeds to appearing more linear at higher speeds. By scaling the drag with the average weight of the grains acting on the rod, we obtain the effective friction μ_{e} encountered over six orders of magnitude of speeds. While a constant μ_{e} is found when the grain size, rod depth, and fluid viscosity are varied at low speeds, a systematic increase is observed as the speed is increased. We analyze μ_{e} in terms of the inertial number I and viscous number J to understand the relative importance of inertia and viscous forces, respectively. For sufficiently high fluid viscosities, we find that the effect of varying the speed, depth, and viscosity can be described by the empirical function μ_{e}=μ_{o}+kJ^{n}, where μ_{o} is the effective friction measured in the quasistatic limit, and k and n are material constants. The drag is then analyzed in terms of the effective viscosity η_{e} and found to decrease systematically as a function of J. We further show that η_{e} as a function of J is directly proportional to the fluid viscosity and the μ_{e} encountered by the rod.

Shear-driven flow of athermal, frictionless, spherocylinder suspensions in two dimensions: Stress, jamming, and contacts

We use numerical simulations to study the flow of a bidisperse mixture of athermal, frictionless, soft-core two-dimensional spherocylinders driven by a uniform steady-state shear strain applied at a fixed finite rate. Energy dissipation occurs via a viscous drag with respect to a uniformly sheared host fluid, giving a simple model for flow in a non-Brownian suspension and resulting in a Newtonian rheology. We study the resulting pressure p and deviatoric shear stress σ of the interacting spherocylinders as a function of packing fraction ϕ, strain rate γ[over ̇], and a parameter α that measures the asphericity of the particles; α is varied to consider the range from nearly circular disks to elongated rods. We consider the direction of anisotropy of the stress tensor, the macroscopic friction μ=σ/p, and the divergence of the transport coefficient η_{p}=p/γ[over ̇] as ϕ is increased to the jamming transition ϕ_{J}. From a phenomenological analysis of Herschel-Bulkley rheology above jamming, we estimate ϕ_{J} as a function of asphericity α and show that the variation of ϕ_{J} with α is the main cause for differences in rheology as α is varied; when plotted as ϕ/ϕ_{J}, rheological curves for different α qualitatively agree. However, a detailed scaling analysis of the divergence of η_{p} for our most elongated particles suggests that the jamming transition of spherocylinders may be in a different universality class than that of circular disks. We also compute the number of contacts per particle Z in the system and show that the value at jamming Z_{J} is a nonmonotonic function of α that is always smaller than the isostatic value. We measure the probability distribution of contacts per unit surface length P(ϑ) at polar angle ϑ with respect to the spherocylinder spine and find that as α→0 this distribution seems to diverge at ϑ=π/2, giving a finite limiting probability for contacts on the vanishingly small flat sides of the spherocylinder. Finally, we consider the variation of the average contact force as a function of location on the particle surface.

Local Rheology Relation with Variable Yield Stress Ratio across Dry, Wet, Dense, and Dilute Granular Flows

Dry, wet, dense, and dilute granular flows have been previously considered fundamentally different and thus described by distinct, and in many cases incompatible, rheologies. We carry out extensive simulations of granular flows, including wet and dry conditions, various geometries and driving mechanisms (boundary driven, fluid driven, and gravity driven), many of which are not captured by standard rheology models. For all simulated conditions, except for fluid-driven and gravity-driven flows close to the flow threshold, we find that the Mohr-Coulomb friction coefficient μ scales with the square root of the local Péclet number Pe provided that the particle diameter exceeds the particle mean free path. With decreasing Pe and granular temperature gradient M, this general scaling breaks down, leading to a yield condition with a variable yield stress ratio characterized by M.

Dynamic Vorticity Banding in Discontinuously Shear Thickening Suspensions

It has recently been argued that steady-state vorticity bands cannot arise in shear thickening suspensions because the normal stress imbalance across the interface between the bands will set up particle migrations. In this Letter, we develop a simple continuum model that couples shear thickening to particle migration. We show by linear stability analysis that homogeneous flow is unstable towards vorticity banding, as expected, in the regime of negative constitutive slope. In full nonlinear computations, we show, however, that the resulting vorticity bands are unsteady, with spatiotemporal patterns governed by stress-concentration coupling. We furthermore show that these dynamical bands also arise in direct particle simulations, in good agreement with the continuum model.

Discontinuous shear thickening in Brownian suspensions

Discontinuous shear thickening in dense suspensions naturally emerges from the activation of frictional forces by shear flow in non-Brownian systems close to jamming. Yet, this physical picture is incomplete as most experiments study soft colloidal particles subject to thermal fluctuations. To characterize discontinuous shear thickening in colloidal suspensions, we use computer simulations to provide a complete description of the competition between athermal jamming, frictional forces, thermal motion, particle softness, and shear flow. We intentionally neglect hydrodynamics, electrostatics, lubrication, and inertia, but can nevertheless achieve quantitative agreement with experimental findings. In particular, shear thickening corresponds to a crossover between frictionless and frictional jamming regimes which is controlled by thermal fluctuations and particle softness and occurs at a softness dependent Péclet number. We also explore the consequences of our findings for constant pressure experiments, and critically discuss the reported emergence of "S-shaped" flow curves. Our work provides the minimal ingredients to quantitatively interpret a large body of experimental work on discontinuous shear thickening in colloidal suspensions.

Viscoinertial regime of immersed granular flows

By means of extensive coupled molecular dynamics-lattice Boltzmann simulations, accounting for grain dynamics and subparticle resolution of the fluid phase, we analyze steady inertial granular flows sheared by a viscous fluid. We show that, for a broad range of system parameters (shear rate, confining stress, fluid viscosity, and relative fluid-grain density), the frictional strength and packing fraction can be described by a modified inertial number incorporating the fluid effect. In a dual viscous description, the effective viscosity diverges as the inverse square of the difference between the packing fraction and its jamming value, as observed in experiments. We also find that the fabric and force anisotropies extracted from the contact network are well described by the modified inertial number, thus providing clear evidence for the role of these key structural parameters in dense suspensions.

Rheology of dense granular flows for elongated particles

We study the rheology of dense granular flows for frictionless spherocylinders by means of 3D numerical simulations. As in the case of spherical particles, the effective friction μ is an increasing function of the inertial number I, and we systematically investigate the dependence of μ on the particle aspect ratio Q, as well as that of the normal stress differences, the volume fraction, and the coordination number. We show in particular that the quasistatic friction coefficient is nonmonotonic with Q: from the spherical case Q=1, it first sharply increases, reaches a maximum around Q≃1.05, and then gently decreases until Q=3, passing its initial value at Q≃2. We provide a microscopic interpretation for this unexpected behavior through the analysis of the distribution of dissipative contacts around the particles: as compared to spheres, slightly elongated grains enhance contacts in their central cylindrical band, whereas at larger aspect ratios particles tend to align and dissipate by preferential contacts at their hemispherical caps.

Discontinuous change of shear modulus for frictional jammed granular materials

The shear modulus of jammed frictional granular materials with harmonic repulsive interaction under an oscillatory shear is numerically investigated. It is confirmed that the storage modulus, the real part of the shear modulus, for frictional grains with sufficiently small strain amplitude γ_{0} discontinuously emerges at the jamming transition point. The storage modulus for small γ_{0} differs from that of frictionless grains even in the zero friction limit, whereas they are almost identical with each other for sufficiently large γ_{0}, where the transition becomes continuous. The stress-strain curve exhibits a hysteresis loop even for a small strain, which connects a linear region for sufficiently small strain to another linear region for larger strain. We propose a scaling law to interpolate between the states of small and large γ_{0}.

Drag law for an intruder in granular sediments

We investigate the drag experienced by a spherical intruder moving through a medium consisting of granular hydrogels immersed in water as a function of its depth, size, and speed. The medium is observed to display a yield stress with a finite force required to move the intruder in the quasistatic regime at low speeds before rapidly increasing at high speeds. In order to understand the relevant time scales that determine drag, we estimate the inertial number I given by the ratio of the time scales required to rearrange grains due to the overburden pressure and imposed shear and the viscous number J given by the ratio of the time scales required to sediment grains in the interstitial fluid and imposed shear. We find that the effective friction μ_{e} encountered by the intruder can be parametrized by I=sqrt[ρ_{g}/P_{p}]v_{i}, where ρ_{g} is the density of the granular hydrogels, v_{i} is the intruder speed, and P_{p} is the overburden pressure due to the weight of the medium, over a wide range of I where the Stokes number St=I^{2}/J≫1. We then show that μ_{e} can be described by the function μ_{e}(I)=μ_{0}+αI^{β}, where μ_{0}, α, and β are constants that depend on the medium. This formula can be used to predict the drag experienced by an intruder of a different size at a different depth in the same medium as a function of its speed.

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.

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.