Numerical stress response functions of static granular layers

We investigate the stress response function of a layer of grains, i.e. the stress profile in response to a localized overload. The shape of the profile is very sensitive to the packing arrangement, and is thus a good signature of the preparation procedure of the layer. This study has been done by the use of molecular dynamics numerical simulations. Here, for a given rain-like preparation, we present the scaling properties of the response function, and in particular the influence of the thickness of the layer, and the importance of the location of the overload and measurement points (at the boundaries, in the bulk).Comment: 6 pages, 4 figures, to appear in the proceedings of the "Traffic and Granular Flow 2003" conferenc

Comment on "Minimal size of a barchan dune"

It is now an accepted fact that the size at which dunes form from a flat sand bed as well as their `minimal size' scales on the flux saturation length. This length is by definition the relaxation length of the slowest mode toward equilibrium transport. The model presented by Parteli, Duran and Herrmann [Phys. Rev. E 75, 011301 (2007)] predicts that the saturation length decreases to zero as the inverse of the wind shear stress far from the threshold. We first show that their model is not self-consistent: even under large wind, the relaxation rate is limited by grain inertia and thus can not decrease to zero. A key argument presented by these authors comes from the discussion of the typical dune wavelength on Mars (650 m) on the basis of which they refute the scaling of the dune size with the drag length evidenced by Claudin and Andreotti [Earth Pla. Sci. Lett. 252, 30 (2006)]. They instead propose that Martian dunes, composed of large grains (500 micrometers), were formed in the past under very strong winds. We show that this saltating grain size, estimated from thermal diffusion measurements, is not reliable. Moreover, the microscopic photographs taken by the rovers on Martian aeolian bedforms show a grain size of 87 plus or minus 25 micrometers together with hematite spherules at millimetre scale. As those so-called ``blueberries'' can not be entrained by reasonable winds, we conclude that the saltating grains on Mars are the small ones, which gives a second strong argument against the model of Parteli et al.Comment: A six page comment on ``Minimal size of a barchan dune'' by Parteli, Duran and Herrmann [Phys. Rev. E 75, 011301 (2007) arXiv:0705.1778

Force chain splitting in granular materials: a mechanism for large scale pseudo-elastic behaviour

We investigate both numerically and analytically the effect of strong disorder on the large scale properties of the hyperbolic equations for stresses proposed in \protect\cite{bcc,wcc}. The physical mechanism that we model is the local splitting of the force chains (the characteristics of the hyperbolic equation) by packing defects. In analogy with the theory of light diffusion in a turbid medium, we propose a Boltzmann-like equation to describe these processes. We show that, for isotropic packings, the resulting large scale effective equations for the stresses have exactly the same structure as those of an elastic body, despite the fact that no displacement field needs to be introduced at all. Correspondingly, the response function evolves from a two peak structure at short scales to a broad hump at large scales. We find, however, that the Poisson ratio is anomalously large and incompatible with classical elasticity theory that requires the reference state to be thermodynamically stable.Comment: 7 pages, 6 figures, An incorrect definition of the Poisson ratio in dimensions not equal to 3 was amended. The conclusions are unchange

Selection of dune shapes and velocities. Part 2: A two-dimensional modelling

We present in this paper a simplification of the dune model proposed by Sauermann et al. which keeps the basic mechanisms but allows analytical and parametric studies. Two kinds of purely propagative two dimensional solutions are exhibited: dunes and domes, which, by contrast to the former, do not show avalanche slip face. Their shape and velocity can be predicted as a function of their size. We recover in particular that dune profiles are not scale invariant (small dunes are flatter than the large ones), and that the inverse of the velocity grows almost linearly with the dune size. We furthermore get the existence of a critical mass below which no stable dune exists. However, the linear stability analysis of a flat sand sheet shows that it is unstable at large wavelengths and suggests a mechanism of dune initiation.Comment: submitted to Eur. Phys. J. B, 13 pages, 17 figure

Transport relaxation time and length scales in turbulent suspensions

We show that in a turbulent flow transporting suspended sediment, the unsaturated sediment flux $q(x,t)$ can be described by a first-order relaxation equation. From a mode analysis of the advection-diffusion equation for the particle concentration, the relaxation length and time scales of the dominant mode are shown to be the deposition length $H U/V_{\rm fall}$ and deposition time $H/V_{\rm fall}$, where $H$ is the flow depth, $U$ the mean flow velocity and $V_{\rm fall}$ the sediment settling velocity. This result is expected to be particularly relevant for the case of sediment transport in slowly varying flows, where the flux is never far from saturation. Predictions are shown to be in quantitative agreement with flume experiments, for both net erosion and net deposition situations.Comment: 16 pages, 8 figures, accepted for publication in J. Fluid Mec

Directed force chain networks and stress response in static granular materials

A theory of stress fields in two-dimensional granular materials based on directed force chain networks is presented. A general equation for the densities of force chains in different directions is proposed and a complete solution is obtained for a special case in which chains lie along a discrete set of directions. The analysis and results demonstrate the necessity of including nonlinear terms in the equation. A line of nontrivial fixed point solutions is shown to govern the properties of large systems. In the vicinity of a generic fixed point, the response to a localized load shows a crossover from a single, centered peak at intermediate depths to two propagating peaks at large depths that broaden diffusively.Comment: 18 pages, 12 figures. Minor corrections to one figur