Elongated particles discharged with a conveyor belt in a two-dimensional silo.

The flow of elliptical particles out of a two-dimensional silo when extracted with a conveyor belt is analyzed experimentally. The conveyor belt-placed directly below the silo outlet-reduces the flow rate, increases the size of the stagnant zone, and it has a very strong influence on the relative velocity fluctuations as they strongly increase everywhere in the silo with decreasing belt speed. In other words, instead of slower but smooth flow, flow reduction by belt leads to intermittent flow. Interestingly, we show that this intermittency correlates with a strong reduction of the orientational order of the particles at the orifice region. Moreover, we observe that the average orientation of the grains passing through the outlet is modified when they are extracted with the belt, a feature that becomes more evident for large orifices.

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.

Discharge of elongated grains in silos under rotational shear.

The discharge of elongated particles from a silo with rotating bottom is investigated numerically. The introduction of a slight transverse shear reduces the flow rate Q by up to 70% compared with stationary bottom, but the flow rate shows a modest increase by further increasing the external shear. Focusing on the dependency of flow rate Q on orifice diameter D, the spheres and rods show two distinct trends. For rods, in the small-aperture limit Q seems to follow an exponential trend, deviating from the classical power-law dependence. These macroscopic observations are in good agreement with our earlier experimental findings [Phys. Rev. E 103, 062905 (2021)2470-004510.1103/PhysRevE.103.062905]. With the help of the coarse-graining methodology we obtain the spatial distribution of the macroscopic density, velocity, kinetic pressure, and orientation fields. This allows us detecting a transition from funnel to mass flow pattern caused by the external shear. Additionally, averaging these fields in the region of the orifice reveals that the strong initial decrease in Q is mostly attributed to changes in the flow velocity, while the weakly increasing trend at higher rotation rates is related to increasing packing fraction. Similar analysis of the grain orientation at the orifice suggests a correlation of the flow rate magnitude with the vertical orientation and the packing fraction at the orifice with the order of the grains. Lastly, the vertical profile of mean acceleration at the center of the silo denotes that the region where the acceleration is not negligible shrinks significantly due to the strong perturbation induced by the moving wall.

Discharge of elongated grains from silo with rotating bottom.

We study the flow of elongated grains (wooden pegs of length L=20 mm with circular cross section of diameter d_{c}=6 and 8 mm) from a silo with a rotating bottom and a circular orifice of diameter D. In the small orifice range (D/d<5) clogs are mostly broken by the rotating base, and the flow is intermittent with avalanches and temporary clogs. Here d≡(3/2d_{c}^{2}L)^{1/3} is the effective grain diameter. Unlike for spherical grains, for rods the flow rate W clearly deviates from the power law dependence W∝(D-kd)^{2.5} at lower orifice sizes in the intermittent regime, where W is measured in between temporary clogs only. Instead, below about D/d<3 an exponential dependence W∝e^{κD} is detected. Here k and κ are constants of order unity. Even more importantly, rotating the silo base leads to a strong-more than 50%-decrease of the flow rate, which otherwise does not depend significantly on the value of ω in the continuous flow regime. In the intermittent regime, W(ω) appears to follow a nonmonotonic trend, although with considerable noise. A simple picture, in terms of the switching from funnel flow to mass flow and the alignment of the pegs due to rotation, is proposed to explain the observed difference between spherical and elongated grains. We also observe shear-induced orientational ordering of the pegs at the bottom such that their long axes in average are oriented at a small angle 〈θ〉≈15^{∘} to the motion of the bottom.

Interacting jammed granular systems.

More than 30 years ago Edwards and co-authors proposed a model to describe the statistics of granular packings by an ensemble of equiprobable jammed states. Experimental tests of this model remained scarce so far. We introduce a simple system to analyze statistical properties of jammed granular ensembles to test Edwards theory. Identical spheres packed in a nearly two-dimensional geometrical confinement were studied in experiments and numerical simulations. When tapped, the system evolves toward a ground state, but due to incompatible domain structures it gets trapped. Analytical calculations reproduce relatively well our simulation results, which allows us to test Edwards theory on a coupled system of two subsystems with different properties. We find that the joint system can only be described by the Edwards theory if considered as a single system due to the constraints in the stresses. The results show counterintuitive effects as in the coupled system the change in the order parameter is opposite to what is expected from the change in the compactivity.

Silo discharge of mixtures of soft and rigid grains.

We study the outflow dynamics and clogging phenomena of mixtures of soft, elastic low-friction spherical grains and hard frictional spheres of similar size in a quasi-two-dimensional (2D) silo with narrow orifice at the bottom. Previous work has demonstrated the crucial influence of elasticity and friction on silo discharge. We show that the addition of small amounts, even as low as 5%, of hard grains to an ensemble of soft, low-friction grains already has significant consequences. The mixtures allow a direct comparison of the probabilities of the different types of particles to clog the orifice. We analyze these probabilities for the hard, frictional and the soft, slippery grains on the basis of their participation in the blocking arches, and compare outflow velocities and durations of non-permanent clogs for different compositions of the mixtures. Experimental results are compared with numerical simulations. The latter strongly suggest a significant influence of the inter-species particle friction.

Orientational order in dense suspensions of elliptical particles in the non-Stokesian regime.

Suspensions of neutrally buoyant elliptic particles are modeled in 2D using fully resolved simulations that provide two-way interaction between the particle and the fluid medium. Forces due to particle collisions are represented by a diffuse interface approach that allows the investigation of dense suspensions (up to 47% packing fraction). We focus on the role inertial forces play at low and high particle Reynolds numbers termed low Reynolds number and inertial regimes, respectively. The suspensions are characterized by the orientation distribution function (ODF) that reflects shear induced rotation of the particles at low Reynolds numbers, and nearly stationary (swaying) particles at high Reynolds numbers. In both cases, orientational ordering differs qualitatively from the behavior observed in the Stokesian-regime. The ODF becomes flatter with increasing packing fraction, as opposed to the sharpening previous work predicted in the Stokesian regime. The ODF at low particle concentrations differs significantly for the low Reynolds number and inertial regimes, whereas with increasing packing fraction convergence is observed. For dense suspensions, the particle-particle interactions dominate the particle motion.

Intermittent flow and transient congestions of soft spheres passing narrow orifices.

Soft, low-friction particles in silos show peculiar features during their discharge. The outflow velocity and the clogging probability both depend upon the momentary silo fill height, in sharp contrast to silos filled with hard particles. The reason is the fill-height dependence of the pressure at the orifice. We study the statistics of silo discharge of soft hydrogel spheres. The outflow is found to become increasingly fluctuating and even intermittent with decreasing orifice size, and with decreasing fill height. In orifices narrower than two particle diameters, outflow can stop completely, but in contrast to clogs formed by rigid particles, these congestions may dissolve spontaneously. We analyze such non-permanent congestions and attribute them to slow reorganization processes in the container, caused by viscoelasticity of the material.

Flow of anisometric particles in a quasi-two-dimensional hopper.

The stationary flow field in a quasi-two-dimensional hopper is investigated experimentally. The behavior of materials consisting of beads and elongated particles with different aspect ratio is compared. We show, that while the vertical velocity in the flowing region can be fitted with a Gaussian function for beads, in the case of elongated grains the flowing channel is narrower and is bordered with sharper velocity gradient. For this case, we quantify deviations from the Gaussian velocity profile. Relative velocity fluctuations are considerably larger and slower for elongated grains.

Frustrated packing in a granular system under geometrical confinement.

Optimal packings of uniform spheres are solved problems in two and three dimensions. The main difference between them is that the two-dimensional ground state can be easily achieved by simple dynamical processes while in three dimensions, this is impossible due to the difference in the local and global optimal packings. In this paper we show experimentally and numerically that in 2 + ε dimensions, realized by a container which is in one dimension slightly wider than the spheres, the particles organize themselves in a triangular lattice, while touching either the front or rear side of the container. If these positions are denoted by up and down the packing problem can be mapped to a 1/2 spin system. At first it looks frustrated with spin-glass like configurations, but the system has a well defined ground state built up from isosceles triangles. When the system is agitated, it evolves very slowly towards the potential energy minimum through metastable states. We show that the dynamics is local and is driven by the optimization of the volumes of 7-particle configurations and by the vertical interaction between touching spheres.

Rheology of dense granular flows for elongated particles.

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

Shear-driven segregation of dry granular materials with different friction coefficients.

We report the first experimental demonstration of bulk segregation in a shear-driven dry granular mixture, where the particles only differ in their surface friction coefficients. We found that the smoother particles tend to sink to the bottom of the shear zone, while rough particles migrate to the top of the sample. This phenomenon is similar to the well known kinetic sieving in particle mixtures with size heterogeneity. In the present case the smooth particles have a higher probability to penetrate into voids created by the shearing than the rough ones. Discrete element simulations were carried out and reproduced the experimentally observed segregation patterns. Moreover, simulations performed in the absence of gravity revealed that rough particles tend to remain in the shear zone, while the smooth particles are being expelled from it. We propose a mechanism in which the smooth particles are driven towards regions of lower shear rate.

Segregation of granular mixtures in a spherical tumbler.

Segregation of polydisperse granular materials in rotating containers is a ubiquitous but still not satisfactorily understood phenomenon. This study describes axial segregation of bidisperse granular mixtures of glass beads in a spherical container, rotating about its horizontal axis. Depending on the filling fraction of the mixer and on the composition of the mixture, qualitatively different spontaneously formed patterns are observed. For technical applications, the well-localized segregated bands allow a convenient separation of individual components of the mixtures. It is particularly surprising that the initial compositions of the granular mixtures have a fundamental influence on the location of the segregated bands. This evidences a collective pattern forming mechanism. The spontaneous formation of these bands cannot simply be traced back to individual particle dynamics. Existing models for segregation in spherical mixers are critically examined and extensions are suggested.

Heaping, secondary flows and broken symmetry in flows of elongated granular particles.

In this paper we report experiments where we shear granular rods in split-bottom geometries, and find that a significant heap of height of least 40% of the filling height can form at the particle surface. We show that heaping is caused by a significant secondary flow, absent for spherical particles. Flow reversal transiently reverses the secondary flow, leading to a quick collapse and slower regeneration of the heap. We present a symmetry argument and experimental data that show that the generation of the secondary flow is driven by a misalignment of the mean particle orientation with the streamlines of the flow. This general mechanism is expected to be important in all flows of sufficiently anisometric grains.

Evolution of shear zones in granular materials.

The evolution of wide shear zones or shear bands was investigated experimentally and numerically for quasistatic dry granular flows in split bottom shear cells. We compare the behavior of materials consisting of beads, irregular grains, such as sand, and elongated particles. Shearing an initially random sample, the zone width was found to significantly decrease in the first stage of the process. The characteristic shear strain associated with this decrease is about unity and it is systematically increasing with shape anisotropy, i.e., when the grain shape changes from spherical to irregular (e.g., sand) and becomes elongated (pegs). The strongly decreasing tendency of the zone width is followed by a slight increase which is more pronounced for rodlike particles than for grains with smaller shape anisotropy (beads or irregular particles). The evolution of the zone width is connected to shear-induced packing density change and for nonspherical particles it also involves grain reorientation effects. The final zone width is significantly smaller for irregular grains than for spherical beads.

Effects of grain shape on packing and dilatancy of sheared granular materials.

A granular material exposed to shear shows a variety of unique phenomena: Reynolds dilatancy, positional order and orientational order effects may compete in the shear zone. We study granular packing consisting of macroscopic prolate, oblate and spherical grains and compare their behaviour. X-ray tomography is used to determine the particle positions and orientations in a cylindrical split bottom shear cell. Packing densities and the arrangements of individual particles in the shear zone are evaluated. For anisometric particles, we observe the competition of two opposite effects. On the one hand, the sheared granules are dilated, on the other hand the particles reorient and align with respect to the streamlines. Even though aligned cylinders in principle may achieve higher packing densities, this alignment compensates for the effect of dilatancy only partially. The complex rearrangements lead to a depression of the surface above the well oriented region while neighbouring parts still show the effect of dilation in the form of heaps. For grains with isotropic shapes, the surface remains rather flat. Perfect monodisperse spheres crystallize in the shear zone, whereby positional order partially overcompensates dilatancy effects. However, even slight deviations from the ideal monodisperse sphere shape inhibit crystallization.

Shear-induced alignment and dynamics of elongated granular particles.

The alignment, ordering, and rotation of elongated granular particles was studied in shear flow. The time evolution of the orientation of a large number of particles was monitored in laboratory experiments by particle tracking using optical imaging and x-ray computed tomography. The experiments were complemented by discrete element simulations. The particles develop an orientational order. In the steady state the time- and ensemble-averaged direction of the main axis of the particles encloses a small angle with the streamlines. This shear alignment angle is independent of the applied shear rate, and it decreases with increasing grain aspect ratio. At the grain level the steady state is characterized by a net rotation of the particles, as dictated by the shear flow. The distribution of particle rotational velocities was measured both in the steady state and also during the initial transients. The average rotation speed of particles with their long axis perpendicular to the shear alignment angle is larger, while shear aligned particles rotate slower. The ratio of this fast/slow rotation increases with particle aspect ratio. During the initial transient starting from an unaligned initial condition, particles having an orientation just beyond the shear alignment angle rotate opposite to the direction dictated by the shear flow.

Orientational order and alignment of elongated particles induced by shear.

Shear induced alignment of elongated particles is studied experimentally and numerically. We show that shear alignment of ensembles of macroscopic particles is comparable even on a quantitative level to simple molecular systems, despite the completely different types of particle interactions. We demonstrate that for dry elongated grains the preferred orientation forms a small angle with the streamlines, independent of shear rate across three decades. For a given particle shape, this angle decreases with increasing aspect ratio of the particles. The shear-induced alignment results in a considerable reduction of the effective friction of the granular material.

High-speed imaging of traveling waves in a granular material during silo discharge.

We report experimental observations of sound waves in a granular material during resonant silo discharge called silo music. The grain motion was tracked by high-speed imaging while the resonance of the silo was detected by accelerometers and acoustic methods. The grains do not oscillate in phase at neighboring vertical locations, but information propagates upward in this system in the form of sound waves. We show that the wave velocity is not constant throughout the silo but considerably increases toward the lower end of the system, suggesting increased pressure in this region, where the flow changes from cylindrical to converging flow. In the upper part of the silo the wave velocity matches the sound velocity measured in the same material when standing (in the absence of flow). Grain oscillations show a stick-slip character only in the upper part of the silo.

Patterns in flowing sand: understanding the physics of granular flow.

Dense granular flows are often unstable and form inhomogeneous structures. Although significant advances have been recently made in understanding simple flows, instabilities of such flows are often not understood. We present experimental and numerical results that show the formation of longitudinal stripes that arise from instability of the uniform flowing state of granular media on a rough inclined plane. The form of the stripes depends critically on the mean density of the flow with a robust form of stripes at high density that consists of fast sliding pluglike regions (stripes) on top of highly agitated boiling material--a configuration reminiscent of the Leidenfrost effect when a droplet of liquid lifted by its vapor is hovering above a hot surface.