Quasistatic to inertial transition in granular materials and the role of fluctuations

Rheology of nonconvex granular flows based on particle rotational characteristics

Particle shape has a profound impact on the flow behaviors of granular materials, yet effectively incorporating the role of particle shape into granular rheology remains challenging. In this study, we employ three representative types of nonconvex particles generated through the multisphere approach and identify a consistent one-to-one relationship between the rescaled friction coefficient and the inertial number I across both inertial and quasistatic flow regimes. However, variations in particle shape cause notable deviations in rheological data compared to their spherical counterparts. Based on the observed dependence of rheological data on I for various nonconvex particles and their convergence at high I, we propose an inertial number I_{s} to effectively capture the impact of particle shape on flow states. The model parameters defining I_{s} are shown to be nearly independent of flow states and configurations, with physical interpretations related to particle rotational characteristics during shear. For practical application, we propose an empirical formula to capture the dependence of model parameters on particle geometrical shapes. The robustness of the proposed model is validated by predicting flow in an inclined flow configuration and applying it to additional nonconvex particles with more irregular and asymmetric features. This establishes a crucial foundation for extending the application of this generalized rheological model to other complex granular flows.

Granular effects in sea ice rheology in the marginal ice zone

Sea ice in the marginal ice zone (MIZ) consists of relatively small floes with a wide size span. In response to oceanic and atmospheric forcing, it behaves as an approximately two-dimensional, highly polydisperse granular material. The established viscous-plastic rheologies used in continuum sea ice models are not suitable for the MIZ; the collisional rheology, in which sea ice is treated as a granular gas, captures only one aspect of the granular behaviour, typical for a narrow range of conditions when dynamics is dominated by binary floe collisions. This paper reviews rheology models and concepts from research on granular materials relevant for MIZ dynamics (average stress as a result of 'microscopic' interactions of grains; [Formula: see text] and collisional rheologies). Idealized discrete-element simulations are used to illustrate granular effects and strong influence of the floe size distribution on strain-stress relationships in sheared sea ice, demonstrating the need for an MIZ rheology model capturing the whole range of 'regimes', from quasi-static/dense flow in the inner MIZ to the inertial flow in the outer MIZ. This article is part of the theme issue 'Theory, modelling and observations of marginal ice zone dynamics: multidisciplinary perspectives and outlooks'.

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.

Three-dimensional and real-scale modeling of flow regimes in dense snow avalanches

Snow avalanches cause fatalities and economic loss worldwide and are one of the most dangerous gravitational hazards in mountainous regions. Various flow behaviors have been reported in snow avalanches, making them challenging to be thoroughly understood and mitigated. Existing popular numerical approaches for modeling snow avalanches predominantly adopt depth-averaged models, which are computationally efficient but fail to capture important features along the flow depth direction such as densification and granulation. This study applies a three-dimensional (3D) material point method (MPM) to explore snow avalanches in different regimes on a complex real terrain. Flow features of the snow avalanches from release to deposition are comprehensively characterized for identification of the different regimes. In particular, brittle and ductile fractures are identified in the different modeled avalanches shortly after their release. During the flow, the analysis of local snow density variation reveals that snow granulation requires an appropriate combination of snow fracture and compaction. In contrast, cohesionless granular flows and plug flows are mainly governed by expansion and compaction hardening, respectively. Distinct textures of avalanche deposits are characterized, including a smooth surface, rough surfaces with snow granules, as well as a surface showing compacting shear planes often reported in wet snow avalanche deposits. Finally, the MPM modeling is verified with a real snow avalanche that occurred at Vallée de la Sionne, Switzerland. The MPM framework has been proven as a promising numerical tool for exploring complex behavior of a wide range of snow avalanches in different regimes to better understand avalanche dynamics. In the future, this framework can be extended to study other types of gravitational mass movements such as rock/glacier avalanches and debris flows with implementation of modified constitutive laws.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10346-021-01692-8.

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.

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.

Two mechanisms of momentum transfer in granular flows

This Rapid Communication highlights the physical processes at the origin of the constitutive law of dense granular flows. In simulated plane shear flows, we present a micro-mechanical expression for the phenomenological friction law μ(I). The expression highlights two distinct pathways for momentum transport-through either balanced contact forces or grain micro-acceleration. We show that these two rate-dependent processes control and explain the friction law. This understanding may help advance rheological models for granular materials and other soft materials such as emulsions and suspensions.

Rheological transition in simple shear of moderately dense assemblies of dry cohesive granules

The rheology of homogeneous cohesive granular assemblies under shear at moderate volume fractions is investigated using the discrete element method for both frictionless and frictional granules. A transition in rheology from inertial to quasistatic scaling is observed at volume fractions below the jamming point of noncohesive systems, which is a function of the granular temperature, energy dissipation, and cohesive potential. The transition is found to be the result of growing clusters, which eventually percolate the domain, and change the mode of momentum transport in the system. Differences in the behavior of the shear stress normalized by the pressure are observed when frictionless and frictional cases are compared. These differences are explained through contact anisotropy after percolation occurs. Both frictionless and frictional systems are found to be vulnerable to instabilities after full system percolation has occurred, where the former becomes thermodynamically unstable and the latter may form shear bands. Finally, implications for constitutive modeling are discussed.

Scaling laws for the mechanics of loose and cohesive granular materials based on Baxter's sticky hard spheres

We have conducted discrete element simulations (pfc3d) of very loose, cohesive, granular assemblies with initial configurations which are drawn from Baxter's sticky hard sphere (SHS) ensemble. The SHS model is employed as a promising auxiliary means to independently control the coordination number z_{c} of cohesive contacts and particle volume fraction ϕ of the initial states. We focus on discerning the role of z_{c} and ϕ for the elastic modulus, failure strength, and the plastic consolidation line under quasistatic, uniaxial compression. We find scaling behavior of the modulus and the strength, which both scale with the cohesive contact density ν_{c}=z_{c}ϕ of the initial state according to a power law. In contrast, the behavior of the plastic consolidation curve is shown to be independent of the initial conditions. Our results show the primary control of the initial contact density on the mechanics of cohesive granular materials for small deformations, which can be conveniently, but not exclusively explored within the SHS-based assembling procedure.

Force fluctuations on a wall in interaction with a granular lid-driven cavity flow

The force fluctuations experienced by a boundary wall subjected to a lid-driven cavity flow are investigated by means of numerical simulations based on the discrete-element method. The time-averaged dynamics inside the cavity volume and the resulting steady force on the wall are governed by the boundary macroscopic inertial number, the latter being derived from the shearing velocity and the confinement pressure imposed at the top. The force fluctuations are quantified through measuring both the autocorrelation of force time series and the distributions of grain-wall forces, at distinct spatial scales from particle scale to wall scale. A key result is that the grain-wall force distributions are entirely driven by the boundary macroscopic inertial number, whatever the spatial scale considered. In particular, when the wall scale is considered, the distributions are found to evolve from nearly exponential to nearly Gaussian distributions by decreasing the macroscopic inertial number. The transition from quasistatic to dense inertial flow is well identified through remarkable changes in the shapes of the distributions of grain-wall forces, accompanied by a loss of system memory in terms of the mesoscale force transmitted toward the wall.

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.

Phase diagram for inertial granular flows

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