Granular surface flows: Investigation through numerical simulations

Daniel BONAMY,

*Coll. : Mathieu RENOUF (LaMCoS, INSA de Lyon, Villeurbannes)*

* Frédéric DUBOIS, Pierre ALART (LMGC, Université Montpellier 2)
*

Top: Typical snapshot of the steady flows in the simulated 2D rotating drum. Bottom, a: Typical velocity profile . (b): Shear stress profile &Sigma:_{xz}, both measured within S. Error bars correspond to a 99% confident interval. Calling d the mean diameter of grains, ρ_{0} their density and g the gravity constant, all distances, times, velocities and stresses are given in units of d,(d/g)^{1/2}, (gd)^{1/2} and ρ_{0}gd respectively. (c) and (d): Resulting constitutive laws relating the ratios μ=Σ_{xz}/Σ_{zz} and k= Σ_{xx}/Σ_{zz} to the inertial number respectively.

Absence of thermal agitation and the dissipative nature of the interaction between the constituent macroscopic particles set granular matter apart from conventional gases, liquids or solids. They can *e.g. *flow like liquids, but under given conditions jam and resist to external shear stress without deforming, like solids - the most spectacular manifestation occurs during avalanches when a thin layer of grains starts to roll at the surface of the packing, most of the grains remaining apparently static. This solid/liquid duality makes dry granular media not easy to describe.

To find out the constitutive laws of dense granular flows, we investigated two-dimensional steady surface flow of cohesionless disks in a rotating drum through *nonsmooth contact dynamics *simulations /1/. Profiles of volume fraction, translational velocity, rms velocity, strain rate, and stress tensor **Σ** are measured at the midpoint along the length of the surface-flowing layer, where the flow is generally considered as steady and homogeneous. Analysis of these data and their interrelations suggest that the inertial number *I* (defined as the ratio between the typical time of deformation and the typical time of confinement) is the relevant dimensionless parameter to describe the transition from the quasistatic part of the packing to the flowing part. Constitutive laws relating the components of the stress tensor to the inertial number were determined /1/.

(a) Typical field of granular temperature T and (b) field of stream function Ψ (Ω=6rpm). Calling d the mean diameter of grains, ρ_{0} their density and g the gravity constant, all distances, times, velocities and stresses are given in units of d,(d/g)^{1/2}, (gd)^{1/2} and ρ_{0}gd respectively. According to the theory, all these fields can be fully described from the knowledge of two functions that relate the granular temperature, the vorticity and the stream function. We checked these predictions on the simulated rotating drum and determined these two functions. (c): the function relating the granular temperature to the stream field.

The effective friction *μ*= Σ* _{xz}*/Σ

We then considered alternative approaches inspired from the modelling of turbulent fluids. In that case, the challenge is to relate the Reynolds stresses, based on small scale fluctuations, to large scale or time averaged quantities. A new way to tackle this problem was recently suggested /4/, through the consideration of non-linear steady solution of the Euler equations, thereby disregarding any non-universal effects induced by (large scale) forcing and (small scale) dissipation. We generalized this method to the case of dense granular flow assuming that in steady states, forcing and dissipation equilibrates *locally */5/. As a result, we obtained a full characterization of the steady state through two general functions that relate granular temperature, vorticity, and stream function. We checked successfully these predictions on the steady surface flows observed in the simulated rotating drum. We are now developing a methodology to derive a “thermodynamic” of these steady states using a formalism derived from equilibrium statistical mechanics, in a spirit similar to the one applied in turbulent fluids modelling /6/.

[1] M. Renouf, D. Bonamy, F. Dubois, P. Alart, Phys. Fluids **17**, 103303 (2005)

[2] G.D.R. Midi, Eur. Phys. J. E, **14** 341 (2004)

[3] D. Bonamy, P. Mills, Europhys. Lett. **62,** 42 (2003)

[4] R. Monchaux, F. Ravelet, B. Dubrulle, A. Chiffaudel, F. Daviaud, Phys. Rev. Lett. **96** 124502 (2006)

[5] D. Bonamy, P.-H. Chavanis, F. Daviaud, B. Dubrulle, M. Renouf, Phys. Fluids, (submitted)

[6] P.-H. Chavanis, R. Sommeria, Phys. Rev. Lett. **78** 3302 (1997)

Maj : 13/04/2007 (775)

• Systèmes complexes et transition énergétique › Physique statistique et systèmes complexes Systèmes complexes et transition énergétique › Statistical physics in mechanics

• Laboratory of Physics and Chemistry of Surfaces and Interfaces • Service de Physique et Chimie des Surfaces et des Interfaces

• Laboratory of Nano-Objects and Complex Systems (LNOSC) • GMT-MSIN : Modélisation des Surfaces Interfaces et Nanostructures