A predictive, size-dependent continuum model for dense granular flows

Soft matter physics of the ground beneath our feet

The soft part of the Earth's surface - the ground beneath our feet - constitutes the basis for life and natural resources, yet a general physical understanding of the ground is still lacking. In this critical time of climate change, cross-pollination of scientific approaches is urgently needed to better understand the behavior of our planet's surface. The major topics in current research in this area cross different disciplines, spanning geosciences, and various aspects of engineering, material sciences, physics, chemistry, and biology. Among these, soft matter physics has emerged as a fundamental nexus connecting and underpinning many research questions. This perspective article is a multi-voice effort to bring together different views and approaches, questions and insights, from researchers that work in this emerging area, the soft matter physics of the ground beneath our feet. In particular, we identify four major challenges concerned with the dynamics in and of the ground: (I) modeling from the grain scale, (II) near-criticality, (III) bridging scales, and (IV) life. For each challenge, we present a selection of topics by individual authors, providing specific context, recent advances, and open questions. Through this, we seek to provide an overview of the opportunities for the broad Soft Matter community to contribute to the fundamental understanding of the physics of the ground, strive towards a common language, and encourage new collaborations across the broad spectrum of scientists interested in the matter of the Earth's surface.

Shear zones in granular mixtures of hard and soft particles with high and low friction

Granular materials show inhomogeneous flows characterized by strain localization. When strain is localized in a sheared granular material, rigid regions of a nearly undeformed state are separated by shear bands, where the material yields and flows. The characteristics of the shear bands are determined by the geometry of the system, the micromechanical material properties, and the kinematics at the particle level. For a split-bottom shear cell, recent experimental work has shown that mixtures of hard, frictional and soft, nearly frictionless particles exhibit wider shear zones than samples with only one of the two components. To explain this finding, we investigate the shear zone properties and the stress response of granular mixtures using discrete element simulations. We show that both interparticle friction and elastic modulus determine the shear-band properties and packing density of granular mixtures of various mixing ratios, but their stress response depends strongly on the interparticle friction. Our study provides a fundamental understanding of the micromechanics of shear band formation in granular mixtures.

Dynamic imaging of force chains in 3D granular media

Granular media constitute the most abundant form of solid matter on Earth and beyond. When external forces are applied to a granular medium, the forces are transmitted through it via chains of contacts among grains-force chains. Understanding the spatial structure and temporal evolution of force chains constitutes a fundamental goal of granular mechanics. Here, we introduce an experimental technique, interference optical projection tomography, to study force chains in three-dimensional (3D) granular packs under triaxial shear loads and illustrate the technique with random assemblies of spheres and icosahedra. We find that, in response to an increasing vertical load, the pack of spheres forms intensifying vertical force chains, while the pack of icosahedra forms more interconnected force-chain networks. This provides microscopic insights into why particles with more angularity are more resistant to shear failure-the interconnected force-chain network is stronger (that is, more resilient to topological collapse) than the isolated force chains in round particles. The longer force chains with less branching in the pack of round particles are more likely to buckle, which leads to the macroscopic failure of the pack. This work paves the way for understanding the grain-scale underpinning of localized failure of 3D granular media, such as shear localization in landslides and stick-slip frictional motion in tectonic and induced earthquakes.

The combined effect of cohesion and finite size on the collapse of wet granular columns

The collapse of low-saturation liquid-containing granular materials is prevalent in nature and industrial processes, and understanding the associated transient dynamics is extremely important for exploring such complex flow processes. In this paper, the collapse of a finite-size wet granular column is systematically studied and the determinants affecting its dynamics are analyzed based on the discrete element model for wet particles and the corresponding small-scale experiments. With the aid of parametric analysis, the dimensionless cohesion parameter containing the system size and grain-scale bond number is proposed, and its relevance in characterizing column stability and collapse dynamics of wet granular materials is further confirmed. For the collapse of wet granular columns with a fixed aspect ratio, the initial height contained in the cohesion parameter is verified to be a manifestation of the finite size effect, which is present in a wet granular collapse and is coupled with the cohesive effect. Such a coupling effect is taken into account in our proposed scaling laws that can be applied to uniformly describe the deposit morphology of wet granular columns after collapse.

Microscopic origin of granular fluidity: An experimental investigation

Granular fluidity has been central to the development of nonlocal constitutive equations, which are necessary for characterizing nonlocal effects observed in experimental granular flow data. However, validation of these equations has been largely computational due to challenges in laboratory experiments. Specifically, the origin of the fluidity on a microscopic, single-particle level is still unproven. In this work, we present an experimental validation of a microscopic definition of granular fluidity, and show the importance of basal boundary conditions to the validity of the theory.

Force on a sphere suspended in flowing granulate

We investigate the force of flowing granular material on an obstacle. A sphere suspended in a discharging silo experiences both the weight of the overlaying layers and drag of the surrounding moving grains. In experiments with frictional hard glass beads, the force on the obstacle was practically flow-rate independent. In contrast, flow of nearly frictionless soft hydrogel spheres added drag to the gravitational force. The dependence of the total force on the obstacle diameter is qualitatively different for the two types of material: It grows quadratically with the obstacle diameter in the soft, low-friction material, while it grows much weaker, nearly linearly with the obstacle diameter, in the bed of glass spheres. In addition to the drag, the obstacle embedded in flowing low-friction soft particles experiences a total force from the top as if immersed in a hydrostatic pressure profile, but a much lower counterforce acting from below. In contrast, when embedded in frictional, hard particles, a strong pressure gradient forms near the upper obstacle surface.

Controlling rheology via boundary conditions in dense granular flows

Boundary shape, particularly roughness, strongly controls the amount of wall slip in dense granular flows. In this paper, we aim to quantify and understand which aspects of a dense granular flow are controlled by the boundary conditions, and to incorporate these observations into a cooperative nonlocal model characterizing slow granular flows. To examine the influence of boundary properties, we perform experiments on a quasi-2D annular shear cell with a rotating inner wall and a fixed outer wall; the latter is selected among 6 walls with various roughnesses, local concavity, and compliance. We find that we can successfully capture the full flow profile using a single set of empirically determined model parameters, with only the wall slip velocity set by direct observation. Through the use of photoelastic particles, we observe how the internal stresses fluctuate more for rougher boundaries, corresponding to a lower wall slip, and connect this observation to the propagation of nonlocal effects originating from the wall. Our measurements indicate a universal relationship between dimensionless fluidity and velocity.

Frictional Weakening of Vibrated Granular Flows

We computationally study the frictional properties of sheared granular media subjected to harmonic vibration applied at the boundary. Such vibrations are thought to play an important role in weakening flows, yet the independent effects of amplitude, frequency, and pressure on the process have remained unclear. Based on a dimensional analysis and DEM simulations, we show that, in addition to a previously proposed criterion for peak acceleration that leads to breaking of contacts, weakening requires the absolute amplitude squared of the displacement to be sufficiently large relative to the confining pressure. The analysis provides a basis for predicting flows subjected to arbitrary external vibration and demonstrates that a previously unrecognized second process that is dependent on dissipation contributes to shear weakening under vibrations.

Quantitative Understanding of the Onset of Dense Granular Flows

The question of when and how dense granular materials start to flow under stress, despite many industrial and geophysical applications, remains largely unresolved. We develop and test a simple equation for the onset of quasistatic flows of granular materials which is based on the frictional aging of the granular packing. The result is a nonmonotonic stress-strain relation which-akin to classical friction-is independent of the shear rate. This relation suffices to understand the quasistatic deformations of aging granular media, and its solid-to-liquid transition. Our results also elucidate the (flow) history dependence of the mechanical properties, and the sensitivity to the initial preparation of granular media.

The effect of grain shape and material on the nonlocal rheology of dense granular flows

Nonlocal rheologies allow for the modeling of granular flows from the creeping to intermediate flow regimes, using a small number of parameters. In this paper, we report on experiments testing how particle properties affect the model parameters used in the Kamrin & Koval cooperative nonlocal model, using particles of three different shapes (circles, ellipses, and pentagons) and three different materials, including one which allows for the measurement of stresses via photoelasticity. Our experiments are performed on a quasi-2D annular shear cell with a rotating inner wall and a fixed outer wall. Each type of particle is found to exhibit flows which are well-fit by nonlocal rheology, with each particle having a distinct triad of the local, nonlocal, and frictional parameters. While the local parameter b is always approximately unity, the nonlocal parameter A depends sensitively on both the particle shape and material. The critical stress ratio μ_s, above which Coulomb failure occurs, varies for particles with the same material but different shape, indicating that geometric friction can dominate over material friction.

Turbulent-like velocity fluctuations in two-dimensional granular materials subject to cyclic shear

We perform a systematic experimental study to investigate the velocity fluctuations in the two-dimensional granular matter of low and high friction coefficients subjected to cyclic shear of a range of shear amplitudes, whose velocity fields are strikingly turbulent-like with vortices of different scales. The scaling behaviors of both the transverse velocity power spectra E_T(k) ∝ k^-α_T and, more severely, the longitudinal velocity power spectra E_L(k) ∝ k^-α_L are affected by the prominent peak centered around k ≈ 2π of the inter-particle distance due to the static structure factor of the hard-particle nature in contrast to the real turbulence. To reduce the strong peak effect to the actual values of α_ν (the subscript 'ν' refers to either T or L), we subsequently analyze the second-order velocity structure functions of S(2)ν(r) in real space, which show the power-law scalings of S(2)ν(r) ∝ r^β_ν for both modes. From the values of β_ν, we deduce the corresponding α_ν from the scaling relations of α_ν = β_ν + 2. The deduced values of α_ν increase continuously with the shear amplitude γ_m, showing no signature of yielding transition, and are slightly larger than α_ν = 2.0 at the limit of γ_m → 0, which corresponds to the elastic limit of the system, for all γ_m. The inter-particle friction coefficients show no significant effect on the turbulent-like velocity fluctuations. Our findings suggest that the turbulent-like collective particle motions are governed by both the elasticity and plasticity in cyclically sheared granular materials.

The Development and Application of a TFM for Dense Particle Flow and Mixing in Rotating Drums

The two-fluid model (TFM) coupled with the kinetic theory of granular flow (KTGF) has gradually been used for modeling dense granular flows and mixing in rotating drums in recent years. In the present paper, a review is made from the perspective of model development and model application. It is found that several frictional viscosity models were proposed to consider the enduring contact of dense particles for the specific rotating studied, but there is still a lack of a universal model. The model is validated by various experiment results and the applicability is indicated. The model is used for investigating dynamic particle flow, and the effects of the parameters on granular flow behavior and flight design. Although the model theoretically has the advantage of saving computing resources, and is suitable for industrial-scale modeling, it is found that the model is used for the research of laboratory-scale rotating drums (diameter less than 0.5 m) and has not been used for industrial rotating drum analysis. Moreover, recommendations for future work are provided.

Gravity Governs Shear Localization in Confined Dense Granular Flows

The prediction of flow profiles of slowly sheared granular materials is a major geophysical and industrial challenge. Understanding the role of gravity is particularly important for future planetary exploration in varying gravitational environments. Using the principle of minimization of energy dissipation, and combining experiments and variational analysis, we disentangle the contributions of the gravitational acceleration, confining pressure, and layer thickness on shear strain localization induced by moving fault boundaries at the bottom of a granular layer. The flow profile is independent of the gravity for geometries with a free top surface. However, under a confining pressure or if the sheared layer withstands the weight of the upper layers, increasing gravity promotes the transition from closed shear zones buried in the bulk to open ones that intersect the top surface. We show that the center position and width of the shear zone and the axial angular velocity at the top surface follow universal scaling laws when properly scaled by the gravity, applied pressure, and layer thickness. Our finding that the flow profiles lie on a universal master curve opens the possibility to predict the quasistatic shear flow of granular materials in varying gravitational environments.

Rotational diffusion and rotational correlations in frictional amorphous disk packings under shear

We show here that rotations of round particles in amorphous disk packing reveal various nontrivial microscopic features when the packing is close to rigidification. We analyze experimental measurements on disk packing subjected to simple shear deformation with various inter-particle friction coefficients and across a range of volume fractions where the system is known to stiffen. The analysis of measurements indicates that shear induces diffusive microrotation, that can be both enhanced and suppressed depending upon the volume fraction as well as the inter-particle friction. Rotations also display persistent anticorrelated motion. Spatial correlations in microrotation are observed to be directly correlated with system pressure. These observations point towards the broader mechanical relevance of collective dynamics in the rotational degree of freedom of particles.

Numerical investigation on effect of particle aspect ratio on the dynamical behaviour of ellipsoidal particle flow

Flow of ellipsoidal particles in a modal shear cell was investigated at the microdynamic level based on discrete element method simulations. In a stress-controlled double-shear condition, the flow was studied by varying the aspect ratio of ellipsoidal particles and comparing with the flow of spherical particle assembly in terms of some key properties, including particle alignment, linear velocity, angular velocity, porosity, contact force and contact energy. It was found that particle elongation impacts the rotational displacement around the axis perpendicular to the shear direction, which causes that the ellipsoidal particles with higher elongation are more aligned with the direction of the shear velocity, with more uniform force network. This then affects other particle properties. The fluctuation of linear velocity and the angular velocity decreases with an increase in particle aspect ratio, although the particle elongation does not significantly affect the flow velocity gradient. There is a reduction in both normal and tangential forces per contact with an increase of particle elongation. Due to the variation of the particle alignment with elongation, the standard deviation of the contact energies increases and then reduces when an increase in particle aspect ratio occurs, and on contrary, the porosity has an opposite variation trend.

Interplay between hysteresis and nonlocality during onset and arrest of flow in granular materials

The jamming transition in granular materials is well-known for exhibiting hysteresis, wherein the level of shear stress required to trigger flow is larger than that below which flow stops. Although such behavior is typically modeled as a simple non-monotonic flow rule, the rheology of granular materials is also nonlocal due to cooperativity at the grain scale, leading for instance to increased strengthening of the flow threshold as system size is reduced. We investigate how these two effects - hysteresis and nonlocality - couple with each other by incorporating non-monotonicity of the flow rule into the nonlocal granular fluidity (NGF) model, a nonlocal constitutive model for granular flows. By artificially tuning the strength of nonlocal diffusion, we demonstrate that both ingredients are key to explaining certain features of the hysteretic transition between flow and arrest. Finally, we assess the ability of the NGF model to quantitatively predict material behavior both around the transition and in the flowing regime, through stress-driven discrete element method (DEM) simulations of flow onset and arrest in various geometries. Along the way, we develop a new methodology to compare deterministic model predictions with the stochastic behavior exhibited by the DEM simulations around the jamming transition.

Efficacy of simple continuum models for diverse granular intrusions

Granular intrusion is commonly observed in natural and human-made settings. Unlike typical solids and fluids, granular media can simultaneously display fluid-like and solid-like characteristics in a variety of intrusion scenarios. This multi-phase behavior increases the difficulty of accurately modeling these and other yielding (or flowable) materials. Micro-scale modeling methods, such as DEM (Discrete Element Method), capture this behavior by modeling the media at the grain scale, but there is often interest in the macro-scale characterizations of such systems. We examine the efficacy of a macro-scale continuum approach in modeling and understanding the physics of various macroscopic phenomena in a variety of granular intrusion cases using two basic frictional yielding constitutive models. We compare predicted granular force response and material flow to experimental data in four quasi-2D intrusion cases: (1) depth-dependent force response in horizontal submerged-intruder motion; (2) separation-dependent drag variation in parallel-plate vertical-intrusion; (3) initial-density-dependent drag fluctuations in free surface plowing, and (4) flow zone development during vertical plate intrusions in under-compacted granular media. Our continuum modeling approach captures the flow process and drag forces while providing key meso- and macro-scopic insights. The modeling results are then compared to experimental data. Our study highlights how continuum modeling approaches provide an alternative for efficient modeling as well as a conceptual understanding of various granular intrusion phenomena.

The effect of boundary roughness on dense granular flows

In the field of granular rheology, an important open question is to understand the influence of boundary conditions on granular flows. We perform experiments in a quasi-2D annular shear cell subject to 6 different boundaries with controlled roughness/compliance. We characterize the granular slip at the boundaries to investigate which aspects of a dense granular flow can be controlled by the choice of boundary condition. Photoelastic techniques are implemented to measure the stress fields P(r) and τ(r) throughout the material. A full inverse-analysis of the fringes within each disk provides the vector force at each contact. This allows us to measure the continuum stress field by coarse-graining internal forces. We have observed that boundary roughness and compliance strongly controls the flow profile v(r) and shear rate profile γ˙(r). We also observed that boundary roughness and compliance play a significant role in the pressure profile P(r) and shear stress profile τ(r).

Coupling non-local rheology and volume of fluid (VOF) method: a finite volume method (FVM) implementation

Additional to a behavior switching between solid-like and liquid-like, dense granular flows also present propagating grain size-dependent effects also called non-local effects. Such behaviors cannot be efficiently modeled by standard rheologies such as µ(I)-rheology but have to be dealt with advanced non-local models. Unfortunately, these models are still new and cannot be used easily nor be used for various configurations. We propose in this work a FVM implementation of the recently popular NGF model coupled with the VOF method in order to both make non-local modeling more accessible to everyone and suitable not only for single-phase flows but also for two-phase flows. The proposed implementation has the advantage to be extremely straightforward and to only require a supplementary stabilization loop compared to the theoretical equations. We then applied our new framework to both single and two-phase flows for validation.

Signature of Transition between Granular Solid and Fluid: Rate-Dependent Stick Slips in Steady Shearing

Despite extensive studies on either smooth granular-fluid flow or the solidlike deformation at the slow limit, the change between these two extremes remains largely unexplored. By systematically investigating the fluctuations of tightly packed grains under steady shearing, we identify a transition zone with prominent stick-slip avalanches. We establish a state diagram, and propose a new dimensionless shear rate based on the speed dependence of interparticle friction and particle size. With fluid-immersed particles confined in a fixed volume and forced to "flow" at viscous numbers J decades below reported values, we answer how a granular system can transition to the regime sustained by solid-to-solid friction that goes beyond existing paradigms based on suspension rheology.

Evolution of shear zones in granular packings under pressure

Stress transmission in realistic granular media often occurs under external load and in the presence of boundary slip. We investigate shear localization in a split-bottom Couette cell with smooth walls subject to a confining pressure experimentally and by means of numerical simulations. We demonstrate how the characteristics of the shear zone, such as its center position and width, evolve as the confining pressure and wall slip modify the local effective friction coefficient of the material. For increasing applied pressure, the shear zone evolves toward the center of the cylinder and grows wider and the angular velocity reduces compared to the driving rate of the bottom disk. Moreover, the presence of slip promotes the transition from open shear zones at the top surface to closed shear zones inside the bulk. We also systematically vary the ratio of the effective friction near the bottom plate and in the bulk in simulations and observe the resulting impact on the surface flow profile. Besides the boundary conditions and external load, material properties such as grain size are also known to influence the effective friction coefficient. However, our numerical results reveal that the center position and width of the shear zone are insignificantly affected by the choice of the grain size as far as it remains small compared to the radius of the rotating bottom disk.

Heaps of sand in flows within a split-bottom Couette cell

In this paper, we study the flow of angular grains in a split-bottom Couette cell. Grains departing from a spherical shape result in collective flow fields that form a heap on the free surface. Here we extend on previous observations in split-bottom cells, exploring a wider range of flows within the inertial regime and finding a richness collection of behaviours. Surface height profiles and velocity profiles are accurately measured with digital image analysis. These measurements allow the characterization of the flow regimes within the cell and the heap morphology. We show that the known flow regimes in split-bottom geometries, like the universal and wall-collapsed regimes, can also be observed in moderately high inertial flows, extending the range for studying universal shear banding. The heap morphology is amplified by the flow inertia, with a partial collapse when the cell comes to a halt. Moreover, at high angular velocities, flows under low confinement will spread radially outwards, while flows under high confinement will develop localized particle ejections. Our results complement the observation of free-surface deformations of flows of nonspherical grains. These observations suggest a need for considering deformable free surface boundary conditions in the simulation of angular grains during shear, with repercussions in the characterization and prediction of natural mass flows.

Analytical solutions for dense, inclined, granular flow over a rigid, bumpy base

We first phrase a boundary-value problem for a dense, steady, fully-developed, gravitational flow of identical inelastic spheres over in inclined bumpy base in the absence of sidewalls. We then obtain approximate analytical solutions for the profiles of the solid volume fraction, the strength of the velocity fluctuations, and the mean velocity of the flow. We compare these with those obtained in numerical solutions of the exact equations.

Dilation as a precursor in a continuous granular fault

We analyze the dilation of the system in a cylindrical granular fault consisting of one single layer of disks submitted to both normal pressure and continuous and slow shear, which results in intermittent and sudden energy release events that reproduce the main laws of seismicity. The dilation of the system can be separated into two parts: a smooth increase of dilation, plus sudden changes both contracting and dilating the medium, which are correlated to abrupt jumps -both positive and negative- in the measured resisting torque. We explain the four possible (and existing) general scenarios combining those two variables: dilation jumps and torque jumps, thanks to the assumption of an optimal local angle in the direction of force chains, and each reorganization of the structure as a replacement of the force chain holding most of the applied stress. The average rate of increase of global dilation varies monotonically with the size of the energy release event, making dilation a plausible candidate to predict catastrophic events in such earthquake-like systems.

Evidence of a non-local ø(I) response

Granular dilatancy has been previously characterised through a simple linear relationship between the packing fraction and dimensionless shear rate. However, this relationship was developed for granular flows in a simple shear cell geometry. Here we examine inertial volume changes in a shear cell with gravity, a vertical chute, and a pseudo-2D hopper. In so doing, we show that the packing fraction displays both a local and non-local response, analogous to what is typically observed for the stress ratio µ.

Microscopic Origin of Nonlocal Rheology in Dense Granular Materials

We study the microscopic origin of nonlocality in dense granular media. Discrete element simulations reveal that macroscopic shear results from a balance between microscopic elementary rearrangements occurring in opposite directions. The effective macroscopic fluidity of the material is controlled by these velocity fluctuations, which are responsible for nonlocal effects in quasistatic regions. We define a new micromechanically based unified constitutive law describing both quasistatic and inertial regimes, valid for different system configurations.

Nonlocal continuum modeling of dense granular flow in a split-bottom cell with a vane-shaped intruder

Shear flow in one spatial region of a dense granular material-induced, for example, through the motion of a boundary-fluidizes the entire granular material. One consequence is that the yield condition vanishes throughout the granular material-even in regions that are very far from the "primary," boundary-driven shear flow. This phenomenon may be characterized through the mechanics of intruders embedded in the granular medium. When there is no primary flow, a critical load must be exceeded to move the intruder; however, in the presence of a primary flow, intruder motion occurs even when an arbitrarily small external load is applied to an intruder embedded in a region far from the sheared zone. In this paper, we apply the nonlocal granular fluidity (NGF) model-a continuum model that involves higher-order flow gradients-to simulate the specific case of dense flow in a split-bottom cell with a vane-shape intruder. Our simulations quantitatively capture the key features of the experimentally observed phenomena: (1) the vanishing of the yield condition, (2) an exponential-type relationship between the applied torque and the rotation rate, (3) the effect of the distance between the intruder and the primary flow zone, and (4) the direction-dependence of the torque/rotation-rate relation, in which the observed relation changes depending on whether the intruder is forced to rotate along with or counter to the primary flow. Importantly, this represents the first fully three-dimensional validation test for a nonlocal model for dense granular flow in general and for the NGF model in particular.

Integration through transients approach to the μ(I) rheology

This work generalizes the granular integration through transients formalism introduced by Kranz et al. [Phys. Rev. Lett. 121, 148002 (2018)10.1103/PhysRevLett.121.148002] to the determination of the pressure. We focus on the Bagnold regime and provide theoretical support to the empirical μ(I) rheology laws that have been successfully applied in many granular flow problems. In particular, we confirm that the interparticle friction is irrelevant in the regime where the μ(I) laws apply.

Power-Law Scaling in Granular Rheology across Flow Geometries

Based on discrete element method simulations, we propose a new form of the constitutive equation for granular flows independent of packing fraction. Rescaling the stress ratio μ by a power of dimensionless temperature Θ makes the data from a wide set of flow geometries collapse to a master curve depending only on the inertial number I. The basic power-law structure appears robust to varying particle properties (e.g., surface friction) in both 2D and 3D systems. We show how this rheology fits and extends frameworks such as kinetic theory and the nonlocal granular fluidity model.

Characterization of microcrystalline cellulose spheres and prediction of hopper flow based on a μ(I)-rheology model

The objective of this study was to characterize the rheology of a pharmaceutical material in the context of the µ(I)-rheology model and to use this model to predict powder flow in a manufacturing operation that is relevant to pharmaceutical manufacturing. The rheology of microcrystalline cellulose spheres was therefore characterized in terms of the μ(I)-rheology model using a modified Malvern Kinexus rheometer. As an example of an important problem in pharmaceutical manufacturing, the flow of these particles from a hopper was studied experimentally and numerically using a continuum Navier-Stokes solver based on the Volume-Of-Fluid (VOF) interface-capturing numerical method. The work shows that the rheology of this typical pharmaceutical material can be measured using a modified annular shear rheometer and that the results can be interpreted in terms of the μ(I)-rheology model. It is demonstrated that both the simulation results and the experimental data show a constant hopper discharge rate. It is noted that the model can suffer from ill-posedness and it is shown how an increasingly fine grid resolution can result in predictions that are not entirely physically realistic. This shortcoming of the numerical framework implies that caution is required when making a one-to-one comparison with experimental data.

Continuously Sheared Granular Matter Reproduces in Detail Seismicity Laws

We introduce a shear experiment that quantitatively reproduces the main laws of seismicity. By continuously and slowly shearing a compressed monolayer of disks in a ringlike geometry, our system delivers events of frictional failures with energies following a Gutenberg-Richter law. Moreover, foreshocks and aftershocks are described by Omori laws and interevent times also follow exactly the same distribution as real earthquakes, showing the existence of memory of past events. Other features of real earthquakes qualitatively reproduced in our system are both the existence of a quiescence preceding some main shocks, as well as magnitude correlations linked to large quakes. The key ingredient of the dynamics is the nature of the force network, governing the distribution of frictional thresholds.

Nonlocal Effects in Inhomogeneous Flows of Soft Athermal Disks

We numerically investigate nonlocal effects on inhomogeneous flows of soft athermal disks close to but below their jamming transition. We employ molecular dynamics to simulate Kolmogorov flows, in which a sinusoidal flow profile with fixed wave number is externally imposed, resulting in a spatially inhomogeneous shear rate. We find that the resulting rheology is strongly wave-number-dependent, and that particle migration, while present, is not sufficient to describe the resulting stress profiles within a conventional local model. We show that, instead, stress profiles can be captured with nonlocal constitutive relations that account for gradients to fourth order. Unlike nonlocal flow in yield stress fluids, we find no evidence of a diverging length scale.

Modeling Segregation in Granular Flows

Accurate continuum models of flow and segregation of dense granular flows are now possible. This is the result of extensive comparisons, over the last several years, of computer simulations of increasing accuracy and scale, experiments, and continuum models, in a variety of flows and for a variety of mixtures. Computer simulations-discrete element methods (DEM)-yield remarkably detailed views of granular flow and segregation. Conti-nuum models, however, offer the best possibility for parametric studies of outcomes in what could be a prohibitively large space resulting from the competition between three distinct driving mechanisms: advection, diffusion, and segregation. We present a continuum transport equation-based framework, informed by phenomenological constitutive equations, that accurately predicts segregation in many settings, both industrial and natural. Three-way comparisons among experiments, DEM, and theory are offered wherever possible to validate the approach. In addition to the flows and mixtures described here, many straightforward extensions of the framework appear possible.

Surface flow profiles for dry and wet granular materials by Particle Tracking Velocimetry; the effect of wall roughness

Two-dimensional Particle Tracking Velocimetry (PTV) is a promising technique to study the behaviour of granular flows. The aim is to experimentally determine the free surface width and position of the shear band from the velocity profile to validate simulations in a split-bottom shear cell geometry. The position and velocities of scattered tracer particles are tracked as they move with the bulk flow by analyzing images. We then use a new technique to extract the continuum velocity field, applying coarse-graining with the postprocessing toolbox MercuryCG on the discrete experimental PTV data. For intermediate filling heights, the dependence of the shear (or angular) velocity on the radial coordinate at the free surface is well fitted by an error function. From the error function, we get the width and the centre position of the shear band. We investigate the dependence of these shear band properties on filling height and rotation frequencies of the shear cell for dry glass beads for rough and smooth wall surfaces. For rough surfaces, the data agrees with the existing experimental results and theoretical scaling predictions. For smooth surfaces, particle-wall slippage is significant and the data deviates from the predictions. We further study the effect of cohesion on the shear band properties by using small amount of silicon oil and glycerol as interstitial liquids with the glass beads. While silicon oil does not lead to big changes, glycerol changes the shear band properties considerably. The shear band gets wider and is situated further inward with increasing liquid saturation, due to the correspondingly increasing trend of particles to stick together.

The dynamics of granular flow from a silo with two symmetric openings

The dynamics of granular flow in a rectangular silo with two symmetrically placed exit openings is investigated using particle image velocimetry (PIV), flow rate measurements and discrete element modelling (DEM). The flow of mustard seeds in a Perspex silo is recorded using a high-speed camera and the resulting image frames are analysed using PIV to obtain velocity, velocity divergence and shear rate plots. A change in flow structure is observed as the distance L between the two openings is varied. The mass flow rate is shown to be at a maximum at zero opening separation, decreasing as L is increased; it then reaches a minimum before rising to an equilibrium rate close to two times that of an isolated (non-interacting) opening. The flow rate experiment is repeated using amaranth and screened sand and similar behaviour is observed. Although this result is in contrast with some recent DEM and physical experiments in silo systems, this effect has been reported in an analogous system: the evacuation of pedestrians from a room through two doors. Our experimental results are replicated using DEM and we show that inter-particle friction controls the flow rate behaviour and explains the discrepancies in the literature results.

Energy Fluctuations in Slowly Sheared Granular Materials

Here we show the first experimental measurement of the particle-scale energy fluctuations ΔE in a slowly sheared layer of photoelastic disks. Starting from an isotropically jammed state, applying shear causes the shear-induced stochastic strengthening and weakening of particle-scale energies, whose statistics and dynamics govern the evolution of the macroscopic stress-strain curve. We find that the ΔE behave as a temperaturelike noise field, showing a novel, Boltzmann-type, double-exponential distribution at any given shear strain γ. Following the framework of the soft glassy rheology theory, we extract an effective temperature χ from the statistics of the energy fluctuations to interpret the slow startup shear (shear starts from an isotropically jammed state) of granular materials as an "aging" process: Starting below one, χ gradually approaches one as γ increases, similar to those of spin glasses, thermal glasses, and bulk metallic glasses.

Geological implication of grain-size segregation in dense granular matter

To the current common belief, grain size segregation in granular matter requires sufficient porosity. Therefore, grain size segregation found in a natural fault gouge could imply elevated fluid pressure and the reduced normal stress on fault, possibly caused by the frictional heat during an earthquake. To clarify whether fluidization is essential to grain size segregation, we conduct numerical simulation on a simple model of fault gouge in a plane shear geometry under constant volume condition: the volume fraction is fixed at 0.6, at which the granular system possesses yield stress. We observe apparent grain size segregation at this volume fraction, meaning that grain size segregation alone does not imply fluidization of granular matter. We also show that segregation is driven by the nonlinear velocity profile, and that the gravity is not essential to segregation. The physical condition tested here may be relevant to earthquake faults: the normal stress of 1 MPa, the sliding velocity of 1 m s^-1, and the duration of 0.1 s.This article is part of the theme issue 'Statistical physics of fracture and earthquakes'.

Strain localization in dry sheared granular materials: A compactivity-based approach

Shear banding is widely observed in natural fault zones as well as in laboratory experiments on granular materials. Understanding the dynamics of strain localization under different loading conditions is essential for quantifying strength evolution of fault gouge and energy partitioning during earthquakes and characterizing rheological transitions and fault zone structure changes. To that end, we develop a physics-based continuum model for strain localization in sheared granular materials. The grain-scale dynamics is described by the shear transformation zone (STZ) theory, a nonequilibrium statistical thermodynamic framework for viscoplastic deformation in amorphous materials. Using a finite strain computational framework, we investigate the initiation and growth of complex shear bands under a variety of loading conditions and identify implications for strength evolution and the ductile to brittle transition. Our numerical results show similar localization patterns to field and laboratory observations and suggest that shear zones show more ductile response at higher confining pressures, lower dilatancy, and loose initial conditions. Lower pressures, higher dilatancy, and dense initial conditions favor a brittle response and larger strength drops. These findings shed light on a range of mechanisms for strength evolution in dry sheared granular materials and provide a critical input to physics-based multiscale models of fault zone instabilities.

Plastic flow and localization in an amorphous material: Experimental interpretation of the fluidity

We present a thorough study of the plastic response of a granular material progressively loaded. We study experimentally the evolution of the plastic field from a homogeneous one to a heterogeneous one and its fluctuations in terms of incremental strain. We show that the plastic field can be decomposed in two components evolving on two decoupled strain increment scales. We argue that the slowly varying part of the field can be identified with the so-called fluidity field introduced recently to interpret the rheological behavior of amorphous materials. This fluidity field progressively concentrates along a macroscopic direction corresponding to the Mohr-Coulomb angle.

Hydrodynamics of shape-driven rigidity transitions in motile tissues

In biological tissues, it is now well-understood that mechanical cues are a powerful mechanism for pattern regulation. While much work has focused on interactions between cells and external substrates, recent experiments suggest that cell polarization and motility might be governed by the internal shear stiffness of nearby tissue, deemed "plithotaxis". Meanwhile, other work has demonstrated that there is a direct relationship between cell shapes and tissue shear modulus in confluent tissues. Joining these two ideas, we develop a hydrodynamic model that couples cell shape, and therefore tissue stiffness, to cell motility and polarization. Using linear stability analysis and numerical simulations, we find that tissue behavior can be tuned between largely homogeneous states and patterned states such as asters, controlled by a composite "morphotaxis" parameter that encapsulates the nature of the coupling between shape and polarization. The control parameter is in principle experimentally accessible, and depends both on whether a cell tends to move in the direction of lower or higher shear modulus, and whether sinks or sources of polarization tend to fluidize the system.