Partially fluidized shear granular flows: Continuum theory and MD simulations

The continuum theory of partially fluidized shear granular flows is tested and calibrated using two dimensional soft particle molecular dynamics simulations. The theory is based on the relaxational dynamics of the order parameter that describes the transition between static and flowing regimes of granular material. We define the order parameter as a fraction of static contacts among all contacts between particles. We also propose and verify by direct simulations the constitutive relation based on the splitting of the shear stress tensor into a``fluid part'' proportional to the strain rate tensor, and a remaining ``solid part''. The ratio of these two parts is a function of the order parameter. The rheology of the fluid component agrees well with the kinetic theory of granular fluids even in the dense regime. Based on the hysteretic bifurcation diagram for a thin shear granular layer obtained in simulations, we construct the ``free energy'' for the order parameter. The theory calibrated using numerical experiments with the thin granular layer is applied to the surface-driven stationary two dimensional granular flows in a thick granular layer under gravity.Comment: 20 pages, 19 figures, submitted to Phys. Rev.

Granular flow down a rough inclined plane: transition between thin and thick piles

The rheology of granular particles in an inclined plane geometry is studied using molecular dynamics simulations. The flow--no-flow boundary is determined for piles of varying heights over a range of inclination angles $\theta$. Three angles determine the phase diagram: $\theta_{r}$, the angle of repose, is the angle at which a flowing system comes to rest; $\theta_{m}$, the maximum angle of stability, is the inclination required to induce flow in a static system; and $\theta_{max}$ is the maximum angle for which stable, steady state flow is observed. In the stable flow region $\theta_{r}<\theta<\theta_{max}$, three flow regimes can be distinguished that depend on how close $\theta$ is to $\theta_{r}$: i) $\theta>>\theta_{r}$: Bagnold rheology, characterized by a mean particle velocity $v_{x}$ in the direction of flow that scales as $v_{x}\propto h^{3/2}$, for a pile of height $h$, ii) $\theta\gtrsim\theta_{r}$: the slow flow regime, characterized by a linear velocity profile with depth, and iii) $\theta\approx\theta_{r}$: avalanche flow characterized by a slow underlying creep motion combined with occasional free surface events and large energy fluctuations. We also probe the physics of the initiation and cessation of flow. The results are compared to several recent experimental studies on chute flows and suggest that differences between measured velocity profiles in these experiments may simply be a consequence of how far the system is from jamming.Comment: 19 pages, 14 figs, submitted to Physics of Fluid

Continuum theory of partially fluidized granular flows

A continuum theory of partially fluidized granular flows is developed. The theory is based on a combination of the equations for the flow velocity and shear stresses coupled with the order parameter equation which describes the transition between flowing and static components of the granular system. We apply this theory to several important granular problems: avalanche flow in deep and shallow inclined layers, rotating drums and shear granular flows between two plates. We carry out quantitative comparisons between the theory and experiment.Comment: 28 pages, 23 figures, submitted to Phys. Rev.

The song of the dunes as a self-synchronized instrument

Since Marco Polo (1) it has been known that some sand dunes have the peculiar ability of emitting a loud sound with a well defined frequency, sometimes for several minutes. The origin of this sustained sound has remained mysterious, partly because of its rarity in nature (2). It has been recognized that the sound is not due to the air flow around the dunes but to the motion of an avalanche (3), and not to an acoustic excitation of the grains but to their relative motion (4-7). By comparing several singing dunes and two controlled experiments, one in the laboratory and one in the field, we here demonstrate that the frequency of the sound is the frequency of the relative motion of the sand grains. The sound is produced because some moving grains synchronize their motions. The existence of a velocity threshold in both experiments further shows that this synchronization comes from an acoustic resonance within the flowing layer: if the layer is large enough it creates a resonance cavity in which grains self-synchronize.Comment: minor changes, essentially more references

Rheology of a confined granular material

We study the rheology of a granular material slowly driven in a confined geometry. The motion is characterized by a steady sliding with a resistance force increasing with the driving velocity and the surrounding relative humidity. For lower driving velocities a transition to stick-slip motion occurs, exhibiting a blocking enhancement whith decreasing velocity. We propose a model to explain this behavior pointing out the leading role of friction properties between the grains and the container's boundary.Comment: 9 pages, 3 .eps figures, submitted to PR

Scale invariance and universality of force networks in static granular matter

Force networks form the skeleton of static granular matter. They are the key ingredient to mechanical properties, such as stability, elasticity and sound transmission, which are of utmost importance for civil engineering and industrial processing. Previous studies have focused on the global structure of external forces (the boundary condition), and on the probability distribution of individual contact forces. The disordered spatial structure of the force network, however, has remained elusive so far. Here we report evidence for scale invariance of clusters of particles that interact via relatively strong forces. We analyzed granular packings generated by molecular dynamics simulations mimicking real granular matter; despite the visual variation, force networks for various values of the confining pressure and other parameters have identical scaling exponents and scaling function, and thus determine a universality class. Remarkably, the flat ensemble of force configurations--a simple generalization of equilibrium statistical mechanics--belongs to the same universality class, while some widely studied simplified models do not.Comment: 15 pages, 4 figures; to appear in Natur

Patterns and Collective Behavior in Granular Media: Theoretical Concepts

Granular materials are ubiquitous in our daily lives. While they have been a subject of intensive engineering research for centuries, in the last decade granular matter attracted significant attention of physicists. Yet despite a major efforts by many groups, the theoretical description of granular systems remains largely a plethora of different, often contradicting concepts and approaches. Authors give an overview of various theoretical models emerged in the physics of granular matter, with the focus on the onset of collective behavior and pattern formation. Their aim is two-fold: to identify general principles common for granular systems and other complex non-equilibrium systems, and to elucidate important distinctions between collective behavior in granular and continuum pattern-forming systems.Comment: Submitted to Reviews of Modern Physics. Full text with figures (2Mb pdf) avaliable at http://mti.msd.anl.gov/AransonTsimringReview/aranson_tsimring.pdf Community responce is appreciated. Comments/suggestions send to [email protected]

Partial jamming and non-locality in dense granular flows

Dense granular flows can exhibit non-local flow behaviours that cannot be predicted by local constitutive laws alone. Such behaviour is accompanied by the existence of diverging cooperativity length. Here we show that this length can be attributed to the development of transient clusters of jammed particles within the flow. By performing DEM simulation of dense granular flows, we directly measure the size of such clusters which scales with the inertial number with a power law. We then derive a general non-local relation based on kinematic compatibility for the existence of clusters in an arbitrary non-homogenous flow. The kinematic nature of this derivation suggests that non-locality should be expected in any material regardless of their local constitutive law, as long as transient clusters exist within the flow