Velocity correlations in dense granular flows observed with internal imaging

Multilayered noise model for transport in complex environments

Transport in complex fluidic environments often exhibits transient subdiffusive dynamics accompanied by non-Gaussian probability density profiles featuring a nonmonotonic non-Gaussian parameter. Such properties cannot be adequately explained by the original theory of Brownian motion. Based on an extension of kinetic theory, this study introduces a chain of hierarchically coupled random walks approach that effectively captures all these intriguing characteristics. If the environment consists of a series of independent white noise sources, then the problem can be expressed as a system of hierarchically coupled Ornstein-Uhlenbech equations. Due to the linearity of the system, the most essential transport properties have a closed analytical form.

Subdiffusion model for granular discharge in a submerged silo

Silo discharge has been extensively studied for decades although questions remain regarding the nature of the velocity field, particularly for submerged systems. In this work, fluid-driven granular drainage was performed in a quasi-two-dimensional silo with grains submerged in fluid. While the observed Gaussian velocity profiles were generally consistent with current diffusion models, the diffusion length was found to significantly decrease with height in contrast to the increases previously seen in dry silos. We propose a phenomenological anomalous diffusion model for the spreading of the flow upwards in the cell, with the fluid-driven flows we study here falling in the category of subdiffusive behavior. As the viscous characteristics of the system were amplified, the diffusion length increased and the shape of the flowing zone in the silo changed, deviating further from the parabolic form predicted by traditional normal diffusion models, in effect becoming more subdiffusive as quantified by a decreasing diffusion exponent.

Strain dependent vorticity in sheared granular media

Displacement fields in sheared particle packings often display vortex-like structures that reveal essential features about the mechanical state of the collection of particles. There are several metrics to quantify these flow field features, yet extracting such quantitative metrics from flow field or particle tracking data involves making numerous choices on the time and length scales over which to average. Here we employ a much used experimental data set on sheared disk packings to explore how such arbitrary data mining choices affect the obtained results. We focus on calculating the strain dependent vorticity, as this metric is a differential method hence potentially sensitive to the way it is computed. We find that the total surface area with an absolute vorticity above a certain threshold approaches a plateau value as shear progresses. This plateau value exhibits a non-monotonic dependence on packing fraction. We also show which range of choices yields results that can support an analysis method independent, physical interpretation of the flow field data.

Shear-induced diffusion in dense granular fluids

Granular materials are comprised of solid, athermal grains. Whilst immune to thermal motion, these grains move and diffuse when they undergo shear deformation. Here we introduce this process of shear-induced diffusion with a focus on dense flows. The goal is to present the established scaling laws for continuum diffusivity and to relate them to the micro-mechanisms of a granular random walk. We then suggest how this knowledge may help advance our understanding of granular rheology and diffusion in other soft-materials.

Poiseuille Flow of Soft Polycrystals in 2D Rough Channels

Polycrystals are partially ordered solids where crystalline order extends over mesoscopic length scales, namely, the grain size. We study the Poisuielle flow of such materials in a rough channel. In general, similar to yield stress fluids, three distinct dynamical states, namely, flowing, stick-slip, and jammed can be observed, with a yield threshold dependent on channel width. Importantly, the interplay between the finite channel width, and the intrinsic ordering scale (the grain size) leads to a new type of spatiotemporal heterogeneity. In wide channels, although the average flow profile remains pluglike, at the underlying granular level, there is vigorous grain remodeling activity resulting from the velocity heterogeneity among the grains. As the channel width approaches typical grain size, the flowing polycrystalline state breaks up into a spatially heterogeneous mixture of flowing liquid like patches and chunks of nearly static grains. Despite these static grains, the average velocity still shows a parabolic profile, dominated by the moving liquidlike patches. However, the solid-liquid front moves at nearly constant speed in the opposite direction of the external drive.

Reticulate collisional structure in boundary-driven granular gases

We report a peculiar head-on collision network between two vibrating boundaries in experiments performed during a parabolic flight and in a laboratory using horizontal vibration. This structure is a new ordering, which is due to an orientation correlation between the relative position and velocity of any particle pair. It weakens the collision frequency and produces a long-range boundary effect. Moreover, we find the molecular chaos assumption is violated in a larger portion of the phase space. Using an anisotropic distribution model, we modify angular integration results and compare them to the results of the kinetic theory.

Heterogeneities Shape Passive Intracellular Transport

A living cell's interior is one of the most complex and intrinsically dynamic systems, providing an elaborate interplay between cytosolic crowding and ATP-driven motion that controls cellular functionality. Here, we investigated two distinct fundamental features of the merely passive, non-biomotor-shuttled material transport within the cytoplasm of Dictyostelium discoideum cells: the anomalous non-linear scaling of the mean-squared displacement of a 150-nm-diameter particle and non-Gaussian distribution of increments. Relying on single-particle tracking data of 320,000 data points, we performed a systematic analysis of four possible origins for non-Gaussian transport: 1) sample-based variability, 2) rarely occurring strong motion events, 3) ergodicity breaking/aging, and 4) spatiotemporal heterogeneities of the intracellular medium. After excluding the first three reasons, we investigated the remaining hypothesis of a heterogeneous cytoplasm as cause for non-Gaussian transport. A, to our knowledge, novel fit model with randomly distributed diffusivities implementing medium heterogeneities suits the experimental data. Strikingly, the non-Gaussian feature is independent of the cytoskeleton condition and lag time. This reveals that efficiency and consistency of passive intracellular transport and the related anomalous scaling of the mean-squared displacement are regulated by cytoskeleton components, whereas cytoplasmic heterogeneities are responsible for the generic, non-Gaussian distribution of increments.

Note: Eliminating stripe artifacts in light-sheet fluorescence imaging

We report two techniques to mitigate stripe artifacts in light-sheet fluorescence imaging. The first uses an image processing algorithm called the multidirectional stripe remover method to filter stripes from an existing image. The second uses an elliptical holographic diffuser with strong scattering anisotropy to prevent stripe formation during image acquisition. These techniques facilitate accurate interpretation of image data, especially in denser samples. They are also facile and cost-effective.

Analyzing X-ray tomographies of granular packings

Starting from three-dimensional volume data of a granular packing, as, e.g., obtained by X-ray Computed Tomography, we discuss methods to first detect the individual particles in the sample and then analyze their properties. This analysis includes the pair correlation function, the volume and shape of the Voronoi cells, and the number and type of contacts formed between individual particles. We mainly focus on packings of monodisperse spheres, but we will also comment on other monoschematic particles such as ellipsoids and tetrahedra. This paper is accompanied by a package of free software containing all programs (including source code) and an example three-dimensional dataset which allows the reader to reproduce and modify all examples given.

Granular rotor as a probe for a nonequilibrium bath

This study numerically and analytically investigates the dynamics of a rotor under viscous or dry friction as a nonequilibrium probe of a granular gas. In order to demonstrate the role of the rotor as a probe for a nonequilibrium bath, the molecular dynamics (MD) simulation of the rotor is performed under viscous or dry friction surrounded by a steady granular gas under gravity. A one-to-one map between the velocity distribution function (VDF) of the granular gas and the angular distribution function for the rotor is theoretically derived. The MD simulation demonstrates that the one-to-one map accurately infers the local VDF of the granular gas from the angular VDF of the rotor, and vice versa.

Nonequilibrium Brownian motion beyond the effective temperature

The condition of thermal equilibrium simplifies the theoretical treatment of fluctuations as found in the celebrated Einstein’s relation between mobility and diffusivity for Brownian motion. Several recent theories relax the hypothesis of thermal equilibrium resulting in at least two main scenarios. With well separated timescales, as in aging glassy systems, equilibrium Fluctuation-Dissipation Theorem applies at each scale with its own “effective” temperature. With mixed timescales, as for example in active or granular fluids or in turbulence, temperature is no more well-defined, the dynamical nature of fluctuations fully emerges and a Generalized Fluctuation-Dissipation Theorem (GFDT) applies. Here, we study experimentally the mixed timescale regime by studying fluctuations and linear response in the Brownian motion of a rotating intruder immersed in a vibro-fluidized granular medium. Increasing the packing fraction, the system is moved from a dilute single-timescale regime toward a denser multiple-timescale stage. Einstein’s relation holds in the former and is violated in the latter. The violation cannot be explained in terms of effective temperatures, while the GFDT is able to impute it to the emergence of a strong coupling between the intruder and the surrounding fluid. Direct experimental measurements confirm the development of spatial correlations in the system when the density is increased.

Still water: dead zones and collimated ejecta from the impact of granular jets

When a dense granular jet hits a target, it forms a large dead zone and ejects a highly collimated conical sheet with a well-defined opening angle. Using experiments, simulations, and continuum modeling, we find that this opening angle is insensitive to the precise target shape and the dissipation mechanisms in the flow. We show that this surprising insensitivity arises because dense granular jet impact, though highly dissipative, is nonetheless controlled by the limit of perfect fluid flow.

Eddy viscosity in dense granular flows

We present a seminal set of experiments on dense granular flows in the stadium shear geometry. The advantage of this geometry is that it produces steady shear flow over large deformations, in which the shear stress is constant. The striking result is that the velocity profiles exhibit an S shape, and are not linear as local constitutive laws would predict. We propose a model that suggests this is a result of wall perturbations which span through the system due to the nonlocal behavior of the material. The model is analogous to that of eddy viscosity in turbulent boundary layers, in which the distance to the wall is introduced to predict velocity profiles. Our findings appear pivotal in a number of experimental and practical situations involving dense granular flows next to a boundary. They could further be adapted to other similar materials such as dense suspensions, foams, or emulsions.

Invited Article: Refractive index matched scanning of dense granular materials

We review an experimental method that allows to probe the time-dependent structure of fully three-dimensional densely packed granular materials and suspensions by means of particle recognition. The method relies on submersing a granular medium in a refractive index matched fluid. This makes the resulting suspension transparent. The granular medium is then visualized by exciting, layer by layer, the fluorescent dye in the fluid phase. We collect references and unreported experimental know-how to provide a solid background for future development of the technique, both for new and experienced users.

Fluctuations of grains inside a discharging two-dimensional silo

We present experimental data corresponding to a two-dimensional dense granular flow, namely, the gravity-driven discharge of grains from a small opening in a silo. We study the local velocity field at the scale of single grains at different places with the help of particle-tracking techniques. From these data, the velocity profiles can be obtained and the validity of some long-standing approaches can be assessed. Moreover, the fluctuations of the velocities are taken into consideration to characterize the features of the advective motion (due to the gravity force) and the diffusive motion, which shows nontrivial behavior.

Packings of deformable spheres

We present an experimental study of disordered packings of deformable spheres. Fluorescent hydrogel spheres immersed in water together with a tomography technique enabled the imaging of the three-dimensional arrangement. The mechanical behavior of single spheres subjected to compression is first examined. Then the properties of packings of a randomized collection of deformable spheres in a box with a moving lid are tested. The transition to a state where the packing withstands finite stresses before yielding is observed. Starting from random packed states, the power law dependence of the normal force versus packing fraction or strain at different velocities is quantified. Furthermore, a compression-decompression sequence at low velocities resulted in rearrangements of the spheres. At larger packing fractions, a saturation of the mean coordination number took place, indicating the deformation and faceting of the spheres.

Spatial distribution functions of random packed granular spheres obtained by direct particle imaging

We measure the two-point density correlations and Voronoi cell distributions of cyclically sheared granular spheres obtained with a fluorescence technique and compare them with random packing of frictionless spheres. We find that the radial distribution function g(r) is captured by the Percus-Yevick equation for initial volume fraction ϕ=0.59. However, small but systematic deviations are observed because of the splitting of the second peak as ϕ is increased toward random close packing. The distribution of the Voronoi free volumes deviates from postulated Γ distributions, and the orientational order metric Q6 shows local order but no long range order. Overall, these measures show significant similarity of random packing of granular and frictionless spheres, but some systematic differences as well.

Physical test of a particle simulation model in a sheared granular system

We report a detailed comparison of a slow gravity-driven sheared granular flow with a discrete-element simulation performed in the same geometry. In the experiments, grains flow inside a silo with a rectangular cross section and are sheared by a rough boundary on one side and smooth boundaries on the other sides. Individual grain position and motion are measured using a particle index-matching imaging technique where a fluorescent dye is added to the interstitial liquid which has the same refractive index as the glass beads. The simulations use a Cundall-Strack contact model between the grains using contact parameters that have been used in many other previous studies and ignore the hydrodynamic effects of the interstitial liquid. Computations are performed to understand the effect of particle coefficient of friction, elasticity, contact model, and polydispersity on mean flow properties. We then perform a detailed comparison of the particle fluctuation properties as measured by the displacement probability distribution function and the mean square displacement. All in all, our study suggests a high level of quantitative agreement between the simulations and experiments.

Long-time tails and cage effect in driven granular fluids

We study the velocity autocorrelation function of a driven granular fluid in the stationary state in three dimensions. As the critical volume fraction of the glass transition in the corresponding elastic system is approached, we observe pronounced cage effects in the velocity autocorrelation function as well as a strong decrease of the diffusion constant, depending on the inelasticity. At moderate densities the velocity autocorrelation function is shown to decay algebraically in time, like t(-3/2), if momentum is conserved locally, and like t(-1), if momentum is not conserved by the driving. A simple scaling argument supports the observed long-time tails.

Spatial correlations in sheared isothermal liquids: from elastic particles to granular particles

Spatial correlations in sheared isothermal liquids for both elastic and granular cases are theoretically investigated. Using the generalized fluctuating hydrodynamics, correlation functions for both the microscopic scale and the macroscopic scale are obtained. We find the existence of long-range correlations obeying power laws. The validity of our theoretical predictions has been verified from molecular-dynamics simulation.

Transition from rolling to jamming in thin granular layers

We study the granular jamming transition for sheared layers of spherical beads ranging in thickness from 1 to 3 times the grain diameter d. As the layer thickness increases slightly above d, the measured friction jumps discontinuously from 0.02 to >0.1, marking the transition from rolling to jamming. Above a critical layer thickness for jamming, the effective granular pressure displays a power law increase with thickness. For thin layers, friction and P increases as the packing fraction decreases near the jamming transition, in contrast to expectations for bulk granular matter.

Correlated motion in the bulk of dense granular flows

Numerical simulations of two-dimensional stationary dense granular flows are performed. We check that the system obeys the h_{stop} phenomenology. Focusing on the spatial correlations of the instantaneous velocity fluctuations of the grains, we give evidence of the existence of correlated motion over several grain diameters in the bulk of the flow. Investigating the role of contact friction and restitution, we show that the associated typical length scale lambda is essentially independent of the grain properties. Moreover, we show that lambda is not controlled by the packing compacity. However, in agreement with previous experimental work, we observe that the correlation length decreases with the shear rate. Computing the flows inertia number I , we show a first-order dependence of lambda on I .