Characterizing the Rheology of Fluidized Granular Matter

How does a flowing granular material respond when you push it sideways? This question sits at the heart of granular rheology — the study of how dense collections of grains flow under stress. This work introduces a new tabletop geometry for probing the rheological properties of granular matter that is already in a flowing (fluidized) state, moving beyond simple viscosity measurements to characterize the full dynamic response of the material.

Schematic of the rotating drum apparatus with oscillatory secondary forcing

Experimental apparatus: a horizontal drum half-filled with glass beads is rotated to fluidize the grains into a steady flowing layer (primary flow). A secondary oscillatory shear is then superimposed by driving the entire drum sinusoidally along its axis. The high-speed camera above tracks both grain motion and drum displacement simultaneously, allowing the phase lag between forcing and response to be extracted precisely.

Experimental Approach

The setup fluidizes granular matter — glass beads of 1–3 mm diameter — by rotating a horizontal drum about its long axis. The grains form a steady flowing layer along the drum surface, captured at 500 fps by a high-speed camera. On top of this primary flow, a second shear is applied by oscillating the drum along its axis at frequency f. The grains respond to this secondary forcing with a measurable phase lag δφ — the angle by which the grain motion lags behind the drum oscillation. This phase lag encodes the rheological state of the flowing material, analogous to the loss angle in conventional viscoelastic measurements.

Linear Response and the GDR-Midi Model

At low forcing frequencies, the phase lag scales linearly with f — a non-Newtonian signature, since a Newtonian fluid would give δφ ∼ √f. This linear scaling is well described by the GDR-Midi model, a continuum framework for dense granular flow that treats the material as a frictional medium with a shear-rate-dependent effective viscosity. Fitting the linear regime yields the effective friction coefficient μ for the flowing layer, providing a direct rheological measurement that matches the expected value from steady-state inclined-plane theory.

Phase lag vs oscillation frequency for different particle sizes and rotation rates — Phase lag δφ as a function of oscillation frequency f for several particle sizes and rotation rates. The linear scaling at low f — distinct from the √f behavior of a Newtonian fluid — is a hallmark of the granular material's frictional rheology.

Dimensionless response curves collapsing to a single master curve — Over 800 response curves, spanning wide ranges of particle size, rotation rate, and inclination angle, collapse onto a single master curve when plotted against the dimensionless ratio t_η/t_f — the ratio of the material's viscous relaxation time to the forcing period. The GDR-Midi model prediction (solid line) matches the data in the linear regime but underpredicts viscosity at fast forcing.

Model Breakdown at Fast Forcing

At higher oscillation frequencies — when the forcing timescale becomes comparable to the material's own viscous relaxation time t_η — the GDR-Midi model breaks down. The experimental data show a smaller phase lag than predicted, implying the material flows more easily than the model expects. This is attributed to the slow reorientation of force chains and the fabric tensor: granular flows are transiently weak when their principal stress axis changes direction, and the material remains in this weakened state for a time comparable to t_η. The dimensionless ratio t_η/t_f neatly collapses all data onto a single master curve, identifying it as the key parameter governing the onset of model breakdown.

Binary Particle Mixtures

Extending the measurements to binary mixtures of 1 mm and 2 mm beads reveals a striking result: even a small fraction of fine particles can dramatically alter the rheology of the flowing layer. After the grains segregate by size (with fine particles migrating to the interior), the mixture response at low forcing resembles that of pure 1 mm beads — far softer than the 2 mm system — suggesting that a small number of fine particles in the flowing surface layer disproportionately control the bulk rheological response.

Phase lag for binary bead mixtures compared to monodisperse systems

Response curves for two binary mixtures (75% and 50% 1 mm beads by volume, mixed with 2 mm beads) compared to the pure monodisperse systems. At low frequencies the mixtures behave like pure 1 mm systems, while at higher frequencies the 50% mixture converges toward the 2 mm curve — suggesting frequency-dependent participation of fine particles in the oscillatory flow.