Eddy viscosity in dense granular flows

Inelastic rotations and pseudoturbulent plastic avalanches in crystals

Plastic deformations in crystals produce microstructures with randomly oriented patches of unstressed lattice forming complex textures. We use a mesoscopic Landau-type tensorial model of crystal plasticity to show that in such textures rotations can originate from crystallographically exact microslips which self organize in the form of laminates of a pseudotwin type. The formation of such laminates can be viewed as an effective internal "wrinkling" of the crystal lattice. While such "wrinkling" disguises itself as an elastically neutral rotation, behind it is inherently dissipative, dislocation-mediated process. Our numerical experiments reveal pseudoturbulent effective rotations with power-law distributed spatial correlations which suggests that the process of dislocational self-organization is inherently unstable and points toward the necessity of a probabilistic description of crystal plasticity.

Classifying grains using behaviour-informed machine learning

Sorting granular materials such as ores, coffee beans, cereals, gravels and pills is essential for applications in mineral processing, agriculture and waste recycling. Existing sorting methods are based on the detection of contrast in grain properties including size, colour, density and chemical composition. However, many grain properties cannot be directly detected in-situ, which significantly impairs sorting efficacy. We show here that a simple neural network can infer contrast in a wide range of grain properties by detecting patterns in their observable kinematics. These properties include grain size, density, stiffness, friction, dissipation and adhesion. This method of classification based on behaviour can significantly widen the range of granular materials that can be sorted. It can similarly be applied to enhance the sorting of other particulate materials including cells and droplets in microfluidic devices.

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.

Self-diffusion scalings in dense granular flows

We report on measurements of self-diffusion coefficients in discrete numerical simulations of steady, homogeneous, collisional shearing flows of nearly identical, frictional, inelastic spheres. We focus on a range of relatively high solid volume fractions that are important in those terrestrial gravitational shearing flows that are dominated by collisional interactions. Diffusion over this range of solid fraction has not been well characterized in previous studies. We first compare the measured values with an empirical scaling based on shear rate previously proposed in the literature, and highlight the presence of anisotropy and the solid fraction dependence. We then compare the numerical measurements with those predicted by the kinetic theory for shearing flows of inelastic spheres and offer an explanation for why the measured and predicted values differ.

Shear-induced diffusion: the role of granular clusters

This paper is concerned with the physical mechanisms controlling shear-induced diffusion in dense granular flows. The starting point is that of the granular random walk occurring in diluted granular flows, which underpins Bagnold’s scaling relating the coefficient of self-diffusion to the grain size and shear rate. By means of DEM simulations of plane shear flows, we measure some deviations from this scaling in dense granular flows with and without contact adhesion. We propose to relate these deviations to the development of correlated motion of grains in these flows, which impacts the magnitude of grain velocity fluctuations and their time persistence.

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.

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.

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.

Avalanche Interpretation of the Power-Law Energy Spectrum in Three-Dimensional Dense Granular Flow

Turbulence is ubiquitous in nonequilibrium systems, and it has been noted that even dense granular flows exhibit characteristics that are typical of turbulent flow, such as the power-law energy spectrum. However, studies on the turbulentlike behavior of granular flows are limited to two-dimensional (2D) flow. We demonstrate that the statistics in three-dimensional (3D) flow are qualitatively different from those in 2D flow. We also elucidate that avalanche dynamics can explain this dimensionality dependence. Moreover, we define clusters of collectively moving particles that are equivalent to vortex filaments. The clusters unveil complicated structures in 3D flows that are absent in 2D flows.

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.

S-shaped velocity deformation induced by ionic hydration in aqueous salt solution flow

Ionic hydration shells are the most noticeable microscopic feature in an aqueous salt solution, and have attracted attention due to their possible contribution to its flow behavior. In this paper, we find by molecular dynamic simulations that an S-shaped velocity profile is induced by the ionic hydration shells in the nano channel flow. Our theoretical analysis implies a linear relationship between the energy density inside the first hydration shell of the ions and the deformation strength of the velocity profiles of aqueous salt solutions, where the deformation strength is quantified by the curvature length defined by the linear deviation extended from the velocity profile. Our simulation results confirm that such a linear relationship holds for chloride salt solutions with monatomic cations, e.g., K, Na, Ca, Mg, Al and the Na/Ca models by varying the valence number of Na and Ca in the salt solutions. Furthermore, the influence of the flow velocity and the channel width upon the velocity deformation strength are also investigated. Our results indicate that the calculated curvature length provides a numerical evaluation for nano flow behavior and would be helpful in nanofluidic device design.

Vortices Enhance Diffusion in Dense Granular Flows

This Letter introduces unexpected diffusion properties in dense granular flows and shows that they result from the development of partially jammed clusters of grains, or granular vortices. Transverse diffusion coefficients D and average vortex sizes ℓ are systematically measured in simulated plane shear flows at differing inertial numbers I revealing (i) a strong deviation from the expected scaling D∝d^{2}γ[over ˙] involving the grain size d and shear rate γ[over ˙] and (ii) an increase in average vortex size ℓ at low I, following ℓ∝dI^{-1/2} but limited by the system size. A general scaling D∝ℓdγ[over ˙] is introduced that captures all the measurements and highlights the key role of vortex size. This leads to establishing a scaling for the diffusivity in dense granular flow as D∝d^{2}sqrt[γ[over ˙]/t_{i}] involving the geometric average of shear time 1/γ[over ˙] and inertial time t_{i} as the relevant time scale. Analysis of grain trajectories is further evidence that this diffusion process arises from a vortex-driven random walk.

Dynamic X-ray radiography reveals particle size and shape orientation fields during granular flow

When granular materials flow, the constituent particles segregate by size and align by shape. The impacts of these changes in fabric on the flow itself are not well understood, and thus novel non-invasive means are needed to observe the interior of the material. Here, we propose a new experimental technique using dynamic X-ray radiography to make such measurements possible. The technique is based on Fourier transformation to extract spatiotemporal fields of internal particle size and shape orientation distributions during flow, in addition to complementary measurements of velocity fields through image correlation. We show X-ray radiography captures the bulk flow properties, in contrast to optical methods which typically measure flow within boundary layers, as these are adjacent to any walls. Our results reveal the rich dynamic alignment of particles with respect to streamlines in the bulk during silo discharge, the understanding of which is critical to preventing destructive instabilities and undesirable clogging. The ideas developed in this paper are directly applicable to many other open questions in granular and soft matter systems, such as the evolution of size and shape distributions in foams and biological materials.

Planar granular shear flow under external vibration

We present results from a planar shear experiment in which a two-dimensional horizontal granular assembly of pentagonal particles sheared between two parallel walls is subjected to external vibration. Particle tracking and photoelastic measurements are used to quantify both grain scale motion and interparticle stresses with and without imposed vibrations. We characterize the particle motion in planar shear and find that flow of these strongly interlocking particles consists of transient vortex motion with a mean flow given by the sum of exponential profiles imposed by the shearing walls. Vibration is applied either through the shearing surface or as bulk vertical vibration of the entire shearing region with dimensionless accelerations Γ=A(2πf)^{2}/g≈0-2. In both cases, increasing amplitude of vibration A at fixed frequency f leads to failure of the force network, reduction in mean stress, and a corresponding reduction in imposed strain. Vibration of the shearing surface is shown to induce the preferential slipping of large-angle force chains. These effects are insensitive to changes in frequency in the range studied (f=30-120 Hz), as sufficiently large displacements are required to relieve the geometrical frustration of the jammed states.

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.

Experimental Validation of a Nonextensive Scaling Law in Confined Granular Media

In this Letter, we address the relationship between the statistical fluctuations of grain displacements for a full quasistatic plane shear experiment, and the corresponding anomalous diffusion exponent α. We experimentally validate a particular case of the Tsallis-Bukman scaling law, α=2/(3-q), where q is obtained by fitting the probability density function (PDF) of the displacement fluctuations with a q-Gaussian distribution, and the diffusion exponent is measured independently during the experiment. Applying an original technique, we are able to evince a transition from an anomalous diffusion regime to a Brownian behavior as a function of the length of the strain window used to calculate the displacements of the grains. The outstanding conformity of fitting curves to a massive amount of experimental data shows a clear broadening of the fluctuation PDFs as the length of the strain window decreases, and an increment in the value of the diffusion exponent-anomalous diffusion. Regardless of the size of the strain window considered in the measurements, we show that the Tsallis-Bukman scaling law remains valid, which is the first experimental verification of this relationship for a classical system at different diffusion regimes. We also note that the spatial correlations show marked similarities to the turbulence in fluids, a promising indication that this type of analysis can be used to explore the origins of the macroscopic friction in confined granular materials.