Rheology of dense vibrated granular flows: non-monotonic response controlled by granular temperature

Using numerical simulations of a stress-imposed vane rheometer, this study reveals that the rheology of dense vibrated granular flows is governed by the balance between grain-scale agitation energy and confinement, leading to a non-monotonic viscosity response to vibration frequency that reconciles previously observed trends in friction and viscosity weakening.

The Gist

Imagine you have a jar filled with sand, marbles, or sugar. If you try to stir it with a spoon, it feels hard and solid, right? That's because the particles are packed tight, locking each other in place like a crowded dance floor where no one can move.

Now, imagine you start shaking the jar. Suddenly, the sand starts to flow like water. This is the basic idea behind granular materials: they can act like solids or liquids depending on how you treat them.

This paper is about a specific experiment: What happens when you shake a jar of sand while trying to stir it? The researchers wanted to understand the "rules of the road" for this shaking, flowing sand. They used powerful computer simulations to act as a virtual lab, testing how different shaking speeds and strengths change how "thick" or "runny" the sand feels.

Here is the story of their findings, broken down into simple concepts:

1. The Setup: The Shaking Spoon

Imagine a spoon (called a "vane") stuck in the middle of a jar of sand. You apply a constant twisting force to the spoon, trying to make it spin.

• Without shaking: The sand is too tight. The spoon barely moves. The sand feels like a solid rock.
• With shaking: You start vibrating the jar up and down. The sand loosens up, and the spoon spins faster. The sand feels more like a liquid.

The researchers asked: If we change how hard we shake (amplitude) and how fast we shake (frequency), how does the "thickness" (viscosity) of the sand change?

2. The Surprise: The "Goldilocks" Frequency

They found a very interesting pattern.

• Shaking harder (larger amplitude): The sand always gets thinner and flows easier. This makes sense; more energy = more movement.
• Shaking faster (frequency): This is where it gets weird.
  • At slow speeds, shaking helps the sand flow.
  • At medium speeds, the sand flows the best. It becomes super runny.
  • But if you shake it too fast, the sand suddenly gets thicker again! It becomes harder to stir.

The Analogy: Think of a crowded dance floor.

• If you gently nudge people (slow shaking), they don't move much.
• If you shake the floor at just the right rhythm (medium frequency), everyone starts dancing and moving freely. The crowd flows!
• If you shake the floor too fast (high frequency), everyone gets jostled so violently that they bump into each other constantly, getting stuck in a chaotic pile. The crowd freezes up again.

3. The Secret Ingredient: "Granular Temperature"

Why does this happen? The paper introduces a concept called Granular Temperature.

• In normal physics, temperature is how fast atoms are jiggling.
• In sand, "Granular Temperature" is how much the individual grains are jiggling and bouncing around.

The researchers discovered that the sand's flow depends on a tug-of-war between two things:

1. The Energy Injection: The shaking tries to make the grains jiggle (raising the temperature).
2. The Energy Loss: The grains are bumpy and sticky. When they collide, they lose energy (like a bouncing ball that eventually stops).

The Sweet Spot: At medium frequencies, the shaking is perfectly timed to keep the grains jiggling without them losing too much energy in collisions. The "temperature" is high, and the sand flows.
The Trap: At very high frequencies, the grains are shaking so fast that they are constantly crashing into each other. They lose all their energy in these collisions. The "temperature" actually drops, the grains lock up, and the sand becomes thick again.

4. The Pressure Factor: The Heavy Lid

The researchers also put a heavy lid on the jar to press down on the sand (confining pressure).

• Heavy Lid: Makes the sand harder to stir (thicker).
• Light Lid: Makes it easier.

They found a simple rule: The "runniness" of the sand depends on the ratio of Shaking Power vs. Lid Weight.

• If you shake hard enough to overcome the heavy lid, the sand flows.
• If the lid is too heavy for the shake, the sand stays solid.

5. The Big Picture: A Unified Rule

The most important takeaway is that the researchers found a single "master key" to explain all these behaviors.

They realized that the sand's behavior isn't just about how fast you shake or how heavy the lid is. It's about the balance of energy.

• Energy of the Jiggle: How much the grains are bouncing (Granular Temperature).
• Energy of the Cage: How hard it is to break free from the neighbors holding you down (Pressure).

If the Jiggle Energy is strong enough to break the Cage, the sand flows. If the Cage is too strong, or if the shaking is so fast that the grains crash into each other and lose their energy, the sand stays solid.

Why Does This Matter?

This isn't just about sand in a jar. This helps us understand:

• Earthquakes: When the ground shakes, soil can turn from solid ground into a flowing liquid (liquefaction), causing landslides. Understanding the "frequency" helps predict when this happens.
• Industry: Factories use vibrating machines to move powders (like flour, cement, or medicine). Knowing the right speed prevents the powder from getting stuck or flowing too fast.

In a nutshell: The paper tells us that shaking sand is a delicate dance. Shake it too little, and it's solid. Shake it just right, and it flows like water. Shake it too frantically, and it jams up again. The secret to controlling it is understanding how much energy the grains have to break free from their neighbors.

Technical Summary

1. Problem Statement

Granular materials exhibit complex rheological behaviors, transitioning between solid-like and fluid-like states depending on external forcing. While the rheology of non-vibrated granular flows is well-characterized (e.g., by the $\mu(I)$ rheology), the behavior of dense granular materials under vertical vibration remains less understood, particularly regarding the interplay between vibration parameters (amplitude $A$, frequency $f$) and confining pressure ($p$).

Previous studies have identified conflicting or incomplete trends:

• Some works suggest friction weakening is controlled by the ratio $A^2/p$, implying a threshold for disrupting force chains.
• Others observe a non-monotonic dependence of viscosity and friction on vibration frequency $f$, where flowability increases to a maximum and then decreases at higher frequencies.
• There is a lack of a unified framework that explains how energy injection, dissipation, and confinement collectively determine the effective viscosity ($\eta$) across the full range of frequencies, especially in the regime where viscosity increases again at high frequencies.

2. Methodology

The authors employed Discrete Element Method (DEM) numerical simulations to replicate a realistic experimental setup known as a stress-imposed vane rheometer.

• Simulation Setup:
  • System: A monodispersed granular system ($N=2600$ particles, radius $r=2$ mm) confined in a cylinder with a conical base and a circular lid of mass $M$.
  • Forcing: The container undergoes sinusoidal vertical displacement $\Delta z(t) = A \sin(2\pi f t)$.
  • Shear Measurement: A rectangular vane immersed in the grains is subjected to a constant torque $T$. The resulting angular velocity $\Omega$ is measured to calculate the effective viscosity $\eta = T / |\langle \Omega \rangle|$.
  • Physics Model: The Hertz-Mindlin contact model was used to describe grain-grain, grain-vane, and grain-wall interactions, accounting for elastic and dissipative forces.
• Parameters Varied:
  • Vibration amplitude ($A$): $6$ to $58$ $\mu$m.
  • Vibration frequency ($f$): Varied to change the rescaled acceleration $\Gamma = (2\pi f)^2 A / g$.
  • Confining pressure ($p$): Varied by changing the lid mass.
• Key Metrics:
  • Effective viscosity ($\eta$).
  • Granular temperature ($K$): Defined as the average kinetic energy of the grains ($K \propto \langle v^2 \rangle$).
  • Energy ratio ($\bar{K}$): The ratio of grain agitation energy to the energy scale set by confinement.

3. Key Contributions

The paper provides a unified physical framework that reconciles previous observations of viscosity weakening and strengthening. The primary contributions are:

1. Identification of the Governing Variable: The authors demonstrate that the effective viscosity is not directly controlled by $A$, $f$, or $p$ individually, but by the ratio of the granular temperature ($K$) to the confining energy scale ($K_p$).
   • Defined as $\bar{K} = K / K_p$, where $K_p \propto p$.
2. Explanation of Non-Monotonicity: The paper clarifies that the non-monotonic behavior of viscosity with respect to frequency is driven by the efficiency of energy transfer.
   • At low frequencies, increasing $f$ increases energy injection, raising $K$ and lowering $\eta$ (fluidization).
   • At high frequencies, the dissipation rate (collision rate) increases faster than energy absorption, causing $K$ to drop, which in turn increases $\eta$ (thickening).
3. Reconciliation of Scaling Laws: The study validates the $A^2/p$ scaling observed in previous friction studies but limits its validity to the minimum viscosity regime (fluidized state). It shows that this scaling fails to predict behavior at high frequencies where the non-monotonic trend dominates.
4. Minimal Model: A theoretical two-block model (Appendix C) is developed to reproduce the $A^2/p$ scaling, showing it arises from non-linear elastic interactions (Hertzian contacts) rather than solely from force chain rupture.

4. Key Results

• Viscosity vs. Frequency ($\eta$ vs. $f$):

  • The effective viscosity exhibits a non-monotonic trend: it decreases as frequency increases up to a minimum ($\eta_{min}$), then increases at higher frequencies.
  • This trend correlates perfectly with the granular temperature ($K$), which also peaks at an intermediate frequency $f^*$.
  • The relationship is described by a power law: $\eta \propto K^\alpha$ with $\alpha \in [-2, -1.5]$.

• Viscosity vs. Confining Pressure ($\eta$ vs. $p$):

  • Increasing confining pressure $p$ increases the effective viscosity (stabilizing the system).
  • However, $p$ has a negligible effect on the granular temperature $K$ itself. The effect of $p$ is purely through the energy barrier required to move grains.

• The Master Curve ($\eta$ vs. $\bar{K}$):

  • When plotting $\eta$ against the dimensionless energy ratio $\bar{K} = K / (p \pi r^3)$, all data points for different amplitudes, frequencies, and pressures collapse onto a single master curve.
  • This confirms that the rheological state is determined by the competition between the agitation energy (trying to break grain cages) and the confining energy (keeping cages intact).

• Limitations of $A^2/p$:

  • While $\eta_{min}$ scales with $A^2/p$, this variable cannot predict the viscosity at high frequencies (the rising branch of the curve).
  • The $A^2/p$ scaling is a consequence of the system reaching a specific fluidized state, but it does not capture the frequency-dependent energy dissipation mechanism.

5. Significance and Implications

• Unified Framework: The paper resolves the apparent contradiction between studies focusing on $A^2/p$ (friction weakening) and those observing non-monotonic frequency dependence. It establishes that energy injection vs. dissipation is the fundamental driver.
• Granular Temperature as a Control Parameter: The work elevates "granular temperature" from a descriptive metric to a predictive control variable for rheology in vibrated systems, analogous to thermal temperature in molecular systems, though with distinct non-equilibrium dynamics.
• Practical Applications: Understanding the non-monotonic response is crucial for industrial applications (e.g., hopper flow, pharmaceutical mixing) and geophysical phenomena (e.g., earthquake-triggered landslides). It suggests that simply increasing vibration frequency does not always improve flowability; there is an optimal frequency range, beyond which the material becomes more resistant to flow.
• Theoretical Insight: The findings highlight that in non-equilibrium granular systems, the relationship between driving parameters and system response is mediated by complex energy dissipation mechanisms, preventing simple scaling laws from holding across all regimes.

In conclusion, the authors successfully demonstrate that the rheology of dense vibrated granular flows is governed by the balance between the kinetic energy of grain agitation and the stabilizing energy of confinement, providing a robust predictive model for effective viscosity across varying driving conditions.