[Paper Review] Slow Drag in 2D Granular Media

Junfei Geng, Robert Behringer|arXiv (Cornell University)|Jun 15, 2004

Granular flow and fluidized beds5 references3 citations

TL;DR

This study investigates slow drag forces in 2D granular media composed of bidisperse disks, showing that the mean drag force increases as a power-law (exponent ~1.5) with reduced packing fraction above a critical value, while exhibiting weak logarithmic dependence on drag velocity. The system displays strong statistical invariance, with force, power spectrum, and avalanche distributions collapsing under appropriate scaling, and a simple failure model with stochastic force chain rupture reproduces key experimental observations including power-law spectra (α = -2), exponential avalanche distributions, and exponential force tail decay.

ABSTRACT

We study the drag force experienced by an object slowly moving at constant velocity through a 2D granular material consisting of bidisperse disks. The drag force is dominated by force chain structures in the bulk of the system, thus showing strong fluctuations. We consider the effect of three important control parameters for the system: the packing fraction, the drag velocity and the size of the tracer particle. We find that the mean drag force increases as a power-law (exponent of 1.5) in the reduced packing fraction, $(γ- γ_c)/γ_c$, as $γ$ passes through a critical packing fraction, $γ_c$. By comparison, the mean drag grows slowly (basically logarithmic) with the drag velocity, showing a weak rate-dependence. However, the system nevertheless exhibits strong statistical invariance in the sense that many physical quantities collapse onto a single curve under appropriate scaling. We also show that the system can be understood using simple failure models, which reproduce many experimental observations. These experimental data and simulations indicate that fluctuations in the drag force seem to be associated with the force chain formation and breaking in the system. Moreover, our simulations suggest that the logarithmic increase of the mean drag force with rate can be accounted for if slow relaxation of the force chain networks is included.

Motivation & Objective

• To understand the origin and scaling of drag forces in dense 2D granular media during slow, constant-velocity motion.
• To investigate the influence of three key control parameters: packing fraction, drag velocity, and tracer particle size on the mean and fluctuating components of the drag force.
• To determine whether statistical invariance (scaling collapse) holds for force distributions, power spectra, and avalanche statistics under appropriate rescaling.
• To test whether a simple failure model based on stochastic force chain rupture can reproduce experimental observations, including power-law spectra and exponential distributions.

Proposed method

• Experiments using bidisperse disks in a 2D granular system to measure drag force on a tracer particle moving at constant velocity.
• Systematic variation of packing fraction, drag velocity, and tracer particle diameter to probe their effects on mean drag and fluctuations.
• Use of photoelastic techniques to visualize force chain structures and their dynamics during particle motion.
• Statistical analysis of force time series, including power spectral density, avalanche size and duration distributions, and force distribution P(f).
• Development of a stochastic failure model where force chains are modeled as springs with failure thresholds drawn from a distribution.
• Incorporation of slow relaxation in the failure model to account for logarithmic velocity dependence in mean drag force.

Experimental results

Research questions

• RQ1 How does the mean drag force scale with packing fraction above the critical packing fraction γc?
• RQ2 What is the dependence of the mean drag force on the drag velocity, and why is it weakly logarithmic?
• RQ3 How do the size and shape of the tracer particle affect the drag force, especially when comparable to surrounding grains?
• RQ4 Do force distributions, power spectra, and avalanche statistics collapse under appropriate scaling, indicating statistical invariance?
• RQ5 Can a simple failure model with stochastic force chain rupture reproduce the observed power-law spectra, exponential avalanche distributions, and exponential force tail decay?

Key findings

• The mean drag force increases as a power-law with exponent ~1.5 in the reduced packing fraction (γ - γc)/γc above the critical packing fraction γc.
• The mean drag force increases logarithmically with drag velocity, indicating weak rate dependence, which is explained by including slow relaxation in the failure model.
• The force distribution P(f) collapses when scaled by the mean force, showing exponential decay at large forces.
• The power spectrum P(ω) collapses when scaled by the drag velocity, with a high-frequency power-law decay characterized by exponent α = -2.
• Avalanche size and duration distributions collapse when scaled by their respective mean values, indicating exponential distributions.
• The failure model successfully reproduces key experimental features: α = -2 power-law spectrum, exponential avalanche size and duration distributions, and exponential force tail decay.