Stresses in Smooth Flows of Dense Granular Media

The form of the stress tensor is investigated in smooth, dense granular flows which are generated in split-bottom shear geometries. We find that, within a fluctuation fluidized spatial region, the form of the stress tensor is directly dictated by the flow field: The stress and strain-rate tensors are co-linear. The effective friction, defined as the ratio between shear and normal stresses acting on a shearing plane, is found not to be constant but to vary throughout the flowing zone. This variation can not be explained by inertial effects, but appears to be set by the local geometry of the flow field. This is in agreement with a recent prediction, but in contrast with most models for slow grain flows, and points to there being a subtle mechanism that selects the flow profiles.Comment: 5 pages, 4 figure

Continuum approach to wide shear zones in quasi-static granular matter

Slow and dense granular flows often exhibit narrow shear bands, making them ill-suited for a continuum description. However, smooth granular flows have been shown to occur in specific geometries such as linear shear in the absence of gravity, slow inclined plane flows and, recently, flows in split-bottom Couette geometries. The wide shear regions in these systems should be amenable to a continuum description, and the theoretical challenge lies in finding constitutive relations between the internal stresses and the flow field. We propose a set of testable constitutive assumptions, including rate-independence, and investigate the additional restrictions on the constitutive relations imposed by the flow geometries. The wide shear layers in the highly symmetric linear shear and inclined plane flows are consistent with the simple constitutive assumption that, in analogy with solid friction, the effective-friction coefficient (ratio between shear and normal stresses) is a constant. However, this standard picture of granular flows is shown to be inconsistent with flows in the less symmetric split-bottom geometry - here the effective friction coefficient must vary throughout the shear zone, or else the shear zone localizes. We suggest that a subtle dependence of the effective-friction coefficient on the orientation of the sliding layers with respect to the bulk force is crucial for the understanding of slow granular flows.Comment: 11 pages, 7 figure

The tail of the contact force distribution in static granular materials

We numerically study the distribution P(f) of contact forces in frictionless bead packs, by averaging over the ensemble of all possible force network configurations. We resort to umbrella sampling to resolve the asymptotic decay of P(f) for large f, and determine P(f) down to values of order 10^{-45} for ordered and disordered systems in two and three dimensions. Our findings unambiguously show that, in the ensemble approach, the force distributions decay much faster than exponentially: P(f) ~ exp(-f^{\alpha}), with alpha \approx 2.0 for 2D systems, and alpha \approx 1.7 for 3D systems.Comment: 4 pages, 4 figures, submitted to Phys. Rev.

Flow in linearly sheared two dimensional foams: from bubble to bulk scale

We probe the flow of two dimensional foams, consisting of a monolayer of bubbles sandwiched between a liquid bath and glass plate, as a function of driving rate, packing fraction and degree of disorder. First, we find that bidisperse, disordered foams exhibit strongly rate dependent and inhomogeneous (shear banded) velocity profiles, while monodisperse, ordered foams are also shear banded, but essentially rate independent. Second, we introduce a simple model based on balancing the averaged drag forces between the bubbles and the top plate and the averaged bubble-bubble drag forces. This model captures the observed rate dependent flows, and the rate independent flows. Third, we perform independent rheological measurements, both for ordered and disordered systems, and find these to be fully consistent with the scaling forms of the drag forces assumed in the simple model, and we see that disorder modifies the scaling. Fourth, we vary the packing fraction $\phi$ of the foam over a substantial range, and find that the flow profiles become increasingly shear banded when the foam is made wetter. Surprisingly, our model describes flow profiles and rate dependence over the whole range of packing fractions with the same power law exponents -- only a dimensionless number $k$ which measures the ratio of the pre-factors of the viscous drag laws is seen to vary with packing fraction. We find that $k \sim (\phi-\phi_c)^{-1}$, where $\phi_c \approx 0.84$, corresponding to the 2d jamming density, and suggest that this scaling follows from the geometry of the deformed facets between bubbles in contact. Overall, our work suggests a route to rationalize aspects of the ubiquitous Herschel-Bulkley (power law) rheology observed in a wide range of disordered materials.Comment: 16 pages, 14 figures, submitted to Phys. Rev. E. High quality version available at: http://www.physics.leidenuniv.nl/sections/cm/gr

Couette Flow of Two-Dimensional Foams

We experimentally investigate flow of quasi two-dimensional disordered foams in Couette geometries, both for foams squeezed below a top plate and for freely floating foams. With the top-plate, the flows are strongly localized and rate dependent. For the freely floating foams the flow profiles become essentially rate-independent, the local and global rheology do not match, and in particular the foam flows in regions where the stress is below the global yield stress. We attribute this to nonlocal effects and show that the "fluidity" model recently introduced by Goyon {\em et al.} ({\em Nature}, {\bf 454} (2008)) captures the essential features of flow both with and without a top plate.Comment: 6 pages, 5 figures, revised versio

Forecasting the SST space-time variability of the Alboran Sea with genetic algorithms

We propose a nonlinear ocean forecasting technique based on a combination of genetic algorithms and empirical orthogonal function (EOF) analysis. The method is used to forecast the space-time variability of the sea surface temperature (SST) in the Alboran Sea. The genetic algorithm finds the equations that best describe the behaviour of the different temporal amplitude functions in the EOF decomposition and, therefore, enables global forecasting of the future time-variability.Comment: 15 pages, 3 figures; latex compiled with agums.st

Bounds on the shear load of cohesionless granular matter

We characterize the force state of shear-loaded granular matter by relating the macroscopic stress to statistical properties of the force network. The purely repulsive nature of the interaction between grains naturally provides an upper bound for the sustainable shear stress, which we analyze using an optimization procedure inspired by the so-called force network ensemble. We establish a relation between the maximum possible shear resistance and the friction coefficient between individual grains, and find that anisotropies of the contact network (or the fabric tensor) only have a subdominant effect. These results can be considered the hyperstatic limit of the force network ensemble and we discuss possible implications for real systems. Finally, we argue how force anisotropies can be related quantitatively to experimental measurements of the effective elastic constants.Comment: 17 pages, 6 figures. v2: slightly rearranged, introduction and discussion rewritte

Ensemble Theory for Force Networks in Hyperstatic Granular Matter

An ensemble approach for force networks in static granular packings is developed. The framework is based on the separation of packing and force scales, together with an a-priori flat measure in the force phase space under the constraints that the contact forces are repulsive and balance on every particle. In this paper we will give a general formulation of this force network ensemble, and derive the general expression for the force distribution $P(f)$. For small regular packings these probability densities are obtained in closed form, while for larger packings we present a systematic numerical analysis. Since technically the problem can be written as a non-invertible matrix problem (where the matrix is determined by the contact geometry), we study what happens if we perturb the packing matrix or replace it by a random matrix. The resulting $P(f)$'s differ significantly from those of normal packings, which touches upon the deep question of how network statistics is related to the underlying network structure. Overall, the ensemble formulation opens up a new perspective on force networks that is analytically accessible, and which may find applications beyond granular matter.Comment: 17 pages, 17 figure

Front propagation into unstable and metastable states in Smectic C* liquid crystals: linear and nonlinear marginal stability analysis

We discuss the front propagation in ferroelectric chiral smectics (SmC*) subjected to electric and magnetic fields applied parallel to smectic layers. The reversal of the electric field induces the motion of domain walls or fronts that propagate into either an unstable or a metastable state. In both regimes, the front velocity is calculated exactly. Depending on the field, the speed of a front propagating into the unstable state is given either by the so-called linear marginal stability velocity or by the nonlinear marginal stability expression. The cross-over between these two regimes can be tuned by a magnetic field. The influence of initial conditions on the velocity selection problem can also be studied in such experiments. SmC$^*$ therefore offers a unique opportunity to study different aspects of front propagation in an experimental system

On elliptic solutions of the cubic complex one-dimensional Ginzburg-Landau equation

The cubic complex one-dimensional Ginzburg-Landau equation is considered. Using the Hone's method, based on the use of the Laurent-series solutions and the residue theorem, we have proved that this equation has neither elliptic standing wave nor elliptic travelling wave solutions. This result amplifies the Hone's result, that this equation has no elliptic travelling wave solutions.Comment: LaTeX, 12 page