New stress and velocity fields for highly frictional granular materials

The idealised theory for the quasi-static flow of granular materials which satisfy the Coulomb-Mohr hypothesis is considered. This theory arises in the limit that the angle of internal friction approaches $\pi/2$, and accordingly these materials may be referred to as being `highly frictional'. In this limit, the stress field for both two-dimensional and axially symmetric flows may be formulated in terms of a single nonlinear second order partial differential equation for the stress angle. To obtain an accompanying velocity field, a flow rule must be employed. Assuming the non-dilatant double-shearing flow rule, a further partial differential equation may be derived in each case, this time for the streamfunction. Using Lie symmetry methods, a complete set of group-invariant solutions is derived for both systems, and through this process new exact solutions are constructed. Only a limited number of exact solutions for gravity driven granular flows are known, so these results are potentially important in many practical applications. The problem of mass flow through a two-dimensional wedge hopper is examined as an illustration.

Lie group symmetry analysis for granular media stress equations

The Airy stress function, although frequently employed in classical linear elasticity, does not receive similar usage for granular media problems. For plane strain quasi-static deformations of a cohesionless Coulomb���Mohr granular solid, a single nonlinear partial differential equation is formulated for the Airy stress function by combining the equilibrium equations with the yield condition. This has certain advantages from the usual approach, in which two stress invariants and a stress angle are introduced, and a system of two partial differential equations is needed to describe the flow. In the present study, the symmetry analysis of differential equations is utilised for our single partial differential equation, and by computing an optimal system of one-dimensional Lie algebras, a complete set of group-invariant solutions is derived. By this it is meant that any group-invariant solution of the governing partial differential equation (provided it can be derived via the classical symmetries method) may be obtained as a member of this set by a suitable group transformation. For general values of the parameters (angle of internal friction and gravity g) it is found there are three distinct classes of solutions which correspond to granular flows considered previously in the literature. For the two limiting cases of high angle of internal friction and zero gravity, the governing partial differential equation admit larger families of Lie point symmetries, and from these symmetries, further solutions are derived, many of which are new. Furthermore, the majority of these solutions are exact, which is rare for granular flow, especially in the case of gravity driven flows.

Perturbation solutions for flow through symmetrical hoppers with inserts and asymmetrical wedge hoppers

Under certain circumstances, an industrial hopper which operates under the "funnel-flow" regime can be converted to the "mass-flow" regime with the addition of a flow-corrective insert. This paper is concerned with calculating granular flow patterns near the outlet of hoppers that incorporate a particular type of insert, the cone-in-cone insert. The flow is considered to be quasi-static, and governed by the Coulomb-Mohr yield condition together with the non-dilatant double-shearing theory. In two dimensions, the hoppers are wedge-shaped, and as such the formulation for the wedge-in-wedge hopper also includes the case of asymmetrical hoppers. A perturbation approach, valid for high angles of internal friction, is used for both two-dimensional and axially symmetric flows, with analytic results possible for both leading order and correction terms. This perturbation scheme is compared with numerical solutions to the governing equations, and is shown to work very well for angles of internal friction in excess of 45 degree.

Anatomy of a submarine pyroclastic flow and associated turbidity current : July 2003 dome collapse, Soufri��re Hills volcano, Montserrat, West Indies

The 12 to 13 July 2003 andesite lava dome collapse at the Soufri��re Hills volcano, Montserrat, provides the first opportunity to document comprehensively both the sub-aerial and submarine sequence of events for an eruption. Numerous pyroclastic flows entered the ocean during the collapse, depositing approximately 90% of the total material into the submarine environment. During peak collapse conditions, as the main flow penetrated the air���ocean interface, phreatic explosions were observed and a surge cloud decoupled from the main flow body to travel 2 to 3 km over the ocean surface before settling. The bulk of the flow was submerged and rapidly mixed with sea water forming a water-saturated mass flow. Efficient sorting and physical differentiation occurred within the flow before initial deposition at 500 m water depth. The coarsest components (���60% of the total volume) were deposited proximally from a dense granular flow, while the finer components (���40%) were efficiently elutriated into the overlying part of the flow, which evolved into a far-reaching turbidity current.

Rheology of dense sheared granular flows

Rapid granular flows are defined as flows in which the time scales for the particle interactions are small compared to the inverse of the strain rate, so that the particle interactions can be treated as instantaneous collisions. We first show, using Discrete Element simulations, that even very dense flows of sand or glass beads with volume fraction between 0.5 and 0.6 are rapid granular flows. Since collisions are instantaneous, a kinetic theory approach for the constitutive relations is most appropriate, and we present kinetic theory results for different microscopic models for particle interaction. The significant difference between granular flows and normal fluids is that energy is not conserved in a granular flow. The differences in the hydrodynamic modes caused by the non-conserved nature of energy are discussed. Going beyond the Boltzmann equation, the effect of correlations is studied using the ring kinetic approximation, and it is shown that the divergences in the viscometric coefficients, which are present for elastic fluids, are not present for granular flows because energy is not conserved. The hydrodynamic model is applied to the flow down an inclined plane. Since energy is not a conserved variable, the hydrodynamic fields in the bulk of a granular flow are obtained from the mass and momentum conservation equations alone. Energy becomes a relevant variable only in thin 'boundary layers' at the boundaries of the flow where there is a balance between the rates of conduction and dissipation. We show that such a hydrodynamic model can predict the salient features of a chute flow, including the flow initiation when the angle of inclination is increased above the 'friction angle', the striking lack of observable variation of the volume fraction with height, the observation of a steady flow only for certain restitution coefficients, and the density variations in the boundary layers.

Dense shallow granular flows

Simplified equations are derived for a granular flow in the `dense' limit where the volume fraction is close to that for dynamical arrest, and the `shallow' limit where the stream-wise length for flow development (L) is large compared with the cross-stream height (h). The mass and diameter of the particles are set equal to 1 in the analysis without loss of generality. In the dense limit, the equations are simplified by taking advantage of the power-law divergence of the pair distribution function chi proportional to (phi(ad) - phi)(-alpha), and a faster divergence of the derivativ rho(d chi/d rho) similar to (d chi/d phi), where rho and phi are the density and volume fraction, and phi(ad) is the volume fraction for arrested dynamics. When the height h is much larger than the conduction length, the energy equation reduces to an algebraic balance between the rates of production and dissipation of energy, and the stress is proportional to the square of the strain rate (Bagnold law). In the shallow limit, the stress reduces to a simplified Bagnold stress, where all components of the stress are proportional to (partial derivative u(x)/partial derivative y)(2), which is the cross-stream (y) derivative of the stream-wise (x) velocity. In the simplified equations for dense shallow flows, the inertial terms are neglected in the y momentum equation in the shallow limit because the are O(h/L) smaller than the divergence of the stress. The resulting model contains two equations, a mass conservation equations which reduces to a solenoidal condition on the velocity in the incompressible limit, and a stream-wise momentum equation which contains just one parameter B which is a combination of the Bagnold coefficients and their derivatives with respect to volume fraction. The leading-order dense shallow flow equations, as well as the first correction due to density variations, are analysed for two representative flows. The first is the development from a plug flow to a fully developed Bagnold profile for the flow down an inclined plane. The analysis shows that the flow development length is ((rho) over barh(3)/B) , where (rho) over bar is the mean density, and this length is numerically estimated from previous simulation results. The second example is the development of the boundary layer at the base of the flow when a plug flow (with a slip condition at the base) encounters a rough base, in the limit where the momentum boundary layer thickness is small compared with the flow height. Analytical solutions can be found only when the stream-wise velocity far from the surface varies as x(F), where x is the stream-wise distance from the start of the rough base and F is an exponent. The boundary layer thickness increases as (l(2)x)(1/3) for all values of F, where the length scale l = root 2B/(rho) over bar. The analysis reveals important differences between granular flows and the flows of Newtonian fluids. The Reynolds number (ratio of inertial and viscous terms) turns out to depend only on the layer height and Bagnold coefficients, and is independent of the flow velocity, because both the inertial terms in the conservation equations and the divergence of the stress depend on the square of the velocity/velocity gradients. The compressibility number (ratio of the variation in volume fraction and mean volume fraction) is independent of the flow velocity and layer height, and depends only on the volume fraction and Bagnold coefficients.

Multi-scale modelling of granular avalanches

Avalanches, debris flows, and landslides are geophysical hazards, which involve rapid mass movement of granular solids, water and air as a single-phase system. The dynamics of a granular flow involve at least three distinct scales: the micro-scale, meso-scale, and the macro-scale. This study aims to understand the ability of continuum models to capture the micro-mechanics of dry granular collapse. Material Point Method (MPM), a hybrid Lagrangian and Eulerian approach, with Mohr-Coulomb failure criterion is used to describe the continuum behaviour of granular column collapse, while the micromechanics is captured using Discrete Element Method (DEM) with tangential contact force model. The run-out profile predicted by the continuum simulations matches with DEM simulations for columns with small aspect ratios ('h/r' < 2), however MPM predicts larger run-out distances for columns with higher aspect ratios ('h/r' > 2). Energy evolution studies in DEM simulations reveal higher collisional dissipation in the initial free-fall regime for tall columns. The lack of a collisional energy dissipation mechanism in MPM simulations results in larger run-out distances. Micro-structural effects, such as shear band formations, were observed both in DEM and MPM simulations. A sliding flow regime is observed above the distinct passive zone at the core of the column. Velocity profiles obtained from both the scales are compared to understand the reason for a slow flow run-out mobilization in MPM simulations. �� 2013 AIP Publishing LLC.

Hybridization of a class of "second order" models of traffic flow

Despite the simultaneous progress of traffic modelling both on the macroscopic and microscopic front, recent works [E. Bourrel, J.B. Lessort, Mixing micro and macro representation of traffic flow: a hybrid model based on the LWR theory, Transport. Res. Rec. 1852 (2003) 193���200; D. Helbing, M. Treiber, Critical discussion of ���synchronized flow���, Coop. Transport. Dyn. 1 (2002) 2.1���2.24; A. Hennecke, M. Treiber, D. Helbing, Macroscopic simulations of open systems and micro���macro link, in: D. Helbing, H.J. Herrmann, M. Schreckenberg, D.E. Wolf (Eds.), Traffic and Granular Flow ���99, Springer, Berlin, 2000, pp. 383���388] highlighted that one of the most promising way to simulate efficiently traffic flow on large road networks is a clever combination of both traffic representations: the hybrid modelling. Our focus in this paper is to propose two hybrid models for which the macroscopic (resp. mesoscopic) part is based on a class of second order model [A. Aw, M. Rascle, Resurection of second order models of traffic flow?, SIAM J. Appl. Math. 60 (2000) 916���938] whereas the microscopic part is a Follow-the Leader type model [D.C. Gazis, R. Herman, R.W. Rothery, Nonlinear follow-the-leader models of traffic flow, Oper. Res. 9 (1961) 545���567; R. Herman, I. Prigogine, Kinetic Theory of Vehicular Traffic, American Elsevier, New York, 1971]. For the first hybrid model, we define precisely the translation of boundary conditions at interfaces and for the second one we explain the synchronization processes. Furthermore, through some numerical simulations we show that the waves propagation is not disturbed and the mass is accurately conserved when passing from one traffic representation to another.

Numerical and experimental investigation of structural stiffness influence on ratcheting convection cell in granular soils under cyclic loading

Lateral cyclic loaded structures in granular soils can lead to an accumulation of irreversible strains by changing their mechanical response (densification) and forming a closed convective cell in the upper layer of the bedding. In the present thesis the convective cell dimension, formation and grain migration inside this closed volume have been studied and presented in relation to structural stiffness and different loads. This relation was experimentally investigated by applying a cyclic lateral force to a scaled flexible vertical element embedded in dry granular soil. The model was monitored with a camera in order to derive the displacement field by means of the PIV technique. Modelling large soil deformation turns out to be difficult, using mesh-based methods. Consequently, a mesh-free approach (DEM) was chosen in order to investigate the granular flow with the aim of extracting interesting micromechanical information. In both the numerical and experimental analyses the effect of different loading magnitudes and different dimensions of the vertical element were considered. The main results regarded the different development, shape and dimensions of the convection cell and the surface settlements. Moreover, the Discrete Element Method has proven to give satisfactory results in the modelling of large deformation phenomena such as the ratcheting convective cell.

Submarine pyroclastic deposits formed during the 20th May 2006 dome collapse of the Soufri��re Hills Volcano, Montserrat

The 20th May 2006 lava dome collapse of the Soufri��re Hills Volcano, Montserrat, had a total non-dense rock equivalent (non-DRE) collapse volume of approximately 115 �� 10 6 m 3. The majority of this volume was deposited into the ocean. The collapse was rapid, 85% of the mobilized volume being removed in just 35 min, giving peak pyroclastic flow flux of 66 �� 10 3 m 3 s -1. Channel and levee facies on the submarine flanks of the volcano and formation of a thick, steep-sided ridge, suggest that the largest and most dense blocks were transported proximally as a high concentration granular flow. Of the submerged volume, 30% was deposited from the base of this granular flow, forming a linear, high-relief ridge that extends 7 km from shore. The remaining 70% of the submerged volume comprises the finer grain sizes, which were transported at least 40 km by turbidity currents on gradients of &lt;2��. At several localities, the May 2006 distal turbidity currents ran up 200 m of topography and eroded up to 20 cm of underlying substrate. Multiple turbidites are preserved, representing current reflection from the graben margins and deflection around topography. The high energy of the May 2006 collapse resulted in longer submarine run out than the larger (210 �� 10 6 m 3) Soufri��re Hills dome collapse in July 2003.

Submarine deposition of volcaniclastic material from the 1995-2005 eruptions of Soufriere Hills volcano, Montserrat

Soufri��re Hills volcano, Montserrat, has been erupting since 1995. During the current eruption, a large part of the material produced by the volcano has been transported into the sea, modifying the morphology of the submarine flanks of the volcano. We present a unique set of swath bathymetric data collected offshore from Montserrat in 1999, 2002 and 2005. From 1999 to 2002, pyroclastic flows associated with numerous dome collapses entered the sea to produce 100 Mm3 deposit. From 2002 to 2005, the 290 Mm3 submarine deposit is mainly from the 12���13 July 2003 collapse. These data allow us to estimate that, by May 2005, at least 482 Mm3 of material had been deposited on the sea floor since 1995. We compare on-land characteristics and volumes of dome collapse events with the submarine deposits and propose a new analysis of their emplacement on the submarine flanks of the volcano. The deposition mechanism shows a slope dependence, with the maximum thickness of deposit before the break in the slope, probably because of the type of the dense granular flow involved. We conclude that from 1995 to 2005 more than 75% of the erupted volume entered the sea.

Saltation in wind blown sand

The Boltzmann equation of the sand particle velocity distribution function in wind-blown sand two-phase flow is established based on the motion equation of single particle in air. And then, the generalized balance law of particle property in single phase granular flow is extended to gas-particle two-phase flow. The velocity distribution function of particle phase is expanded into an infinite series by means of Grad's method and the Gauss distribution is used to replace Maxwell distribution. In the case of truncation at the third-order terms, a closed third-order moment dynamical equation system is constructed. The theory is further simplified according to the measurement results obtained by stroboscopic photography in wind tunnel tests.

A constitutive description of the inelastic response of ceramics

The objective of the article is to present a unified model for the dynamic mechanical response of ceramics under compressive stress states. The model incorporates three principal deformation mechanisms: (i) lattice plasticity due to dislocation glide or twinning; (ii) microcrack extension; and (iii) granular flow of densely packed comminuted particles. In addition to analytical descriptions of each mechanism, prescriptions are provided for their implementation into a finite element code as well as schemes for mechanism transitions. The utility of the code in addressing issues pertaining to deep penetration is demonstrated through a series of calculations of dynamic cavity expansion in an infinite medium. The results reveal two limiting behavioral regimes, dictated largely by the ratio of the cavity pressure p to the material yield strength ��Y. At low values of p/��Y, cavity expansion occurs by lattice plasticity and hence its rate diminishes with increasing ��Y. In contrast, at high values, expansion occurs by microcracking followed by granular plasticity and is therefore independent of ��Y. In the intermediate regime, the cavity expansion rate is governed by the interplay between microcracking and lattice plasticity. That is, when lattice plasticity is activated ahead of the expanding cavity, the stress triaxiality decreases (toward more negative values) which, in turn, reduces the propensity for microcracking and the rate of granular flow. The implications for penetration resistance to high-velocity projectiles are discussed. Finally, the constitutive model is used to simulate the quasi-static and dynamic indentation response of a typical engineering ceramic (alumina) and the results compared to experimental measurements. Some of the pertinent observations are shown to be captured by the present model whereas others require alternative approaches (such as those based on fracture mechanics) for complete characterization. �� 2011 The American Ceramic Society.

Damage development in an armor ceramic under quasi-static indentation

The objective of the present study is to assess the capabilities of a recently developed mechanism-based model for inelastic deformation and damage in structural ceramics. In addition to conventional lattice plasticity, the model accounts for microcrack growth and coalescence as well as granular flow following comminution. The assessment is made through a coupled experimental/computational study of the indentation response of a commercial armor ceramic. The experiments include examinations of subsurface damage zones along with measurements of residual surface profiles and residual near-surface stresses. Extensive finite element computations are conducted in parallel. Comparisons between experiment and simulation indicate that the most discriminating metric in the assessment is the spatial extent of subsurface damage following indentation. Residual stresses provide additional validation. In contrast, surface profiles of indents are dictated largely by lattice plasticity and thus provide minimal additional insight into the inelastic deformation resulting from microcracking or granular flow. A satisfactory level of correlation is obtained using property values that are either measured directly or estimated from physically based arguments, without undue reliance on adjustable (nonphysical) parameters. �� 2011 The American Ceramic Society.