Experimental insights into deformation dynamics and intermittency in rapid granular shear flows

This work concerns the experimental study of rapid granular shear flows in annular Couette geometry. The flow is induced by continuous driving of the horizontal plate at the top of the granular bed in an annulus. The compressive pressure, driving torque, instantaneous bed height and rotational speed of the shearing plate are measured. Moreover, local stress fluctuations are measured in a medium made of steel spheres 2 and 3 mm in diameter. Both monodisperse packing and bidisperse packing are investigated to reveal the influence of size diversity in intermittent features of granular materials. Experiments are conducted in an annulus that can contain up to 15 kg of spherical steel balls. The shearing granular medium takes place via the rotation of the upper plate which compresses the material loaded inside the annulus. Fluctuations of compressive force are locally measured at the bottom of the annulus using a piezoelectric sensor. Rapid shear flow experiments are pursued at different compressive forces and shear rates and the sensitivity of fluctuations are then investigated by different means through monodisperse and bidisperse packings. Another important feature of rapid granular shear flows is the formation of ordered structures upon shearing. It requires a certain range for the amount of granular material (uniform size distribution) loaded in the system in order to obtain stable flows. This is studied more deeply in this thesis. The results of the current work bring some new insights into deformation dynamics and intermittency in rapid granular shear flows. The experimental apparatus is modified in comparison to earlier investigations. The measurements produce data for various quantities continuously sampled from the start of shearing to the end. Static failure and dynamic shearing ofa granular medium is investigated. The results of this work revealed some important features of failure dynamics and structure formation in the system. Furthermore, some computer simulations are performed in a 2D annulus to examine the nature of kinetic energy dissipation. It is found that turbulent flow models can statistically represent rapid granular flows with high accuracy. In addition to academic outcomes and scientific publications our results have a number of technological applications associated with grinding, mining and massive grain storages.

Mass transfer and particleinteractions in granular systems and suspension flows

The purpose of this study was to investigate some important features of granular flows and suspension flows by computational simulation methods. Granular materials have been considered as an independent state ofmatter because of their complex behaviors. They sometimes behave like a solid, sometimes like a fluid, and sometimes can contain both phases in equilibrium. The computer simulation of dense shear granular flows of monodisperse, spherical particles shows that the collisional model of contacts yields the coexistence of solid and fluid phases while the frictional model represents a uniform flow of fluid phase. However, a comparison between the stress signals from the simulations and experiments revealed that the collisional model would result a proper match with the experimental evidences. Although the effect of gravity is found to beimportant in sedimentation of solid part, the stick-slip behavior associated with the collisional model looks more similar to that of experiments. The mathematical formulations based on the kinetic theory have been derived for the moderatesolid volume fractions with the assumption of the homogeneity of flow. In orderto make some simulations which can provide such an ideal flow, the simulation of unbounded granular shear flows was performed. Therefore, the homogeneous flow properties could be achieved in the moderate solid volume fractions. A new algorithm, namely the nonequilibrium approach was introduced to show the features of self-diffusion in the granular flows. Using this algorithm a one way flow can beextracted from the entire flow, which not only provides a straightforward calculation of self-diffusion coefficient but also can qualitatively determine the deviation of self-diffusion from the linear law at some regions nearby the wall inbounded flows. Anyhow, the average lateral self-diffusion coefficient, which was calculated by the aforementioned method, showed a desirable agreement with thepredictions of kinetic theory formulation. In the continuation of computer simulation of shear granular flows, some numerical and theoretical investigations were carried out on mass transfer and particle interactions in particulate flows. In this context, the boundary element method and its combination with the spectral method using the special capabilities of wavelets have been introduced as theefficient numerical methods to solve the governing equations of mass transfer in particulate flows. A theoretical formulation of fluid dispersivity in suspension flows revealed that the fluid dispersivity depends upon the fluid properties and particle parameters as well as the fluid-particle and particle-particle interactions.

Parallel 3D Simulation of a Fault Gouge using the Lattice Solid Model

Despite the insight gained from 2-D particle models, and given that the dynamics of crustal faults occur in 3-D space, the question remains, how do the 3-D fault gouge dynamics differ from those in 2-D? Traditionally, 2-D modeling has been preferred over 3-D simulations because of the computational cost of solving 3-D problems. However, modern high performance computing architectures, combined with a parallel implementation of the Lattice Solid Model (LSM), provide the opportunity to explore 3-D fault micro-mechanics and to advance understanding of effective constitutive relations of fault gouge layers. In this paper, macroscopic friction values from 2-D and 3-D LSM simulations, performed on an SGI Altix 3700 super-cluster, are compared. Two rectangular elastic blocks of bonded particles, with a rough fault plane and separated by a region of randomly sized non-bonded gouge particles, are sheared in opposite directions by normally-loaded driving plates. The results demonstrate that the gouge particles in the 3-D models undergo significant out-of-plane motion during shear. The 3-D models also exhibit a higher mean macroscopic friction than the 2-D models for varying values of interparticle friction. 2-D LSM gouge models have previously been shown to exhibit accelerating energy release in simulated earthquake cycles, supporting the Critical Point hypothesis. The 3-D models are shown to also display accelerating energy release, and good fits of power law time-to-failure functions to the cumulative energy release are obtained.

The Dynamics of Force Fluctuations in Shear Granular Flows

Tämän työn tavoitteena oli tutkia rakeisen materiaalin kinematiikkaa ja rakentaa koelaitteisto rakeisen materiaalin leikkausjännitysvirtauksien tutkimiseen. Kokeellisessa osassa on keskitytty sisäisiin voimaheilahteluihin ja niiden ymmärtämiseen. Teoriaosassa on käyty läpi rakeisen materiaalin yleisiä ominaisuuksia ja lisäksi on esitetty kaksi eri tapaa mallintaa fysikaalisien ominaisuuksien heilahteluja rakeisessa materiaalissa. Nämä kaksi esitettyä mallinnusmenetelmää ovat skalaarinen q-malli ja simulointi. Skalaarinen q-malli määrittelee jokaiseen yksittäiseen rakeeseen kohdistuvan jännityksen, rakeen ollessa osa 2- tai 3-dimensionaalista asetelmaa. Tämän mallin perusidea on kuvata jännityksien epähomogeenisuutta, joka johtuu rakeiden satunnaisasettelusta. Simulointimallinnus perustuu event-driven algoritmiin, missä systeemin dynamiikkaa kuvataan yksittäisillä partikkelien törmäyksillä. Törmäyksien vaiheet ratkaistiin käyttämällä liikemääräyhtälöitä ja restituution määritelmää. Teoriaosuudessa käytiin vielä pieniltä osin läpi syitä jännitysheilahteluihin ja rakeisen materiaalin lukkiintumiseen. Tutkimuslaitteistolla tutkittiin rakeisen materiaalin käyttäytymistä rengasmaisessa leikkausjännitysvirtauksessa. Tutkimusosuuden päätavoitteena oli mitata partikkelien kosketuksista ja törmäyksistä johtuvia hetkellisiä voimaheilahteluja rengastilavuuden pohjalta. Rakeisena materiaalina tutkimuksessa käytettiin teräskuulia. Jännityssignaali ajan funktiona osoittaa suurta heilahtelua, joka voi olla jopa kertalukua keskiarvosta suurempaa. Tällainen suuren amplitudin omaava heilahtelu on merkittävä haittapuoli yleisesti rakeisissa materiaaleissa käytettyjen jatkuvuusmallien kanssa. Tällainen heilahtelu tekee käytetyt jatkuvuusmallit epäpäteviksi. Yleisellä tasolla jännityksien todennäköisyysjakauma on yhtäpitävä skalaarisen q-mallin tuloksien kanssa. Molemmissa tapauksissa todennäköisyysjakaumalla on eksponentiaalinen muoto.

The behaviour of granular material in pure shear, direct shear and simple shear

In biaxial compression tests, the stress calculations based on boundary information underestimate the principal stresses leading to a significant overestimation of the shear strength. In direct shear tests, the shear strain becomes highly concentrated in the mid-plane of the sample during the test. Although the stress distribution within the specimen is heterogeneous, the evolution of the stress ratio inside the shear band is similar to that inferred from the boundary force calculations. It is also demonstrated that the dilatancy in the shear band significantly exceeds that implied from the boundary displacements. In simple shear tests, the stresses acting on the wall boundaries do not reflect the internal state of stress but merely provide information about the average mobilised wall friction. It is demonstrated that the results are sensitive to the initial stress state defined by K0 = sh/sv. For all cases, non-coaxiality of the principal stress and strain-rate directions is examined and the corresponding flow rule is identified. Periodic cell simulations have been used to examine biaxial compression for a wide range of initial packing densities. Both constant volume and constant mean stress tests have been simulated. The characteristic behaviour at both the macroscopic and microscopic scales is determined by whether or not the system percolates (enduring connectivity is established in all directions). The transition from non-percolating to percolating systems is characterised by transitional behaviour of internal variables and corresponds to an elastic percolation threshold, which correlates well with the establishment of a mechanical coordination number of ca. 3.0. Strong correlations are found between macroscopic and internal variables at the critical state.

Shear Strength of Granular Materials, 1967

A theory was developed to allow the separate determination of the effects of the interparticle friction and interlocking of particles on the shearing resistance and deformational behavior of granular materials. The derived parameter, angle of solid friction, is independent of the type of shear test, stress history, porosity and the level of confining pressure, and depends solely upon the nature of the particle surface. The theory was tested against published data concerning the performance of plane strain, triaxial compression and extension tests on cohesionless soils. The theory also was applied to isotropically consolidated undrained triaxial tests on three crushed limestones prepared by the authors using vibratory compaction. The authors concluded that, (1) the theory allowed the determination of solid friction between particles which was found to depend solely on the nature of the particle surface, (2) the separation of frictional and volume change components of shear strength of granular materials qualitatively corroborated the postulated mechanism of deformation (sliding and rolling of groups of particles over other similar groups with resulting dilatancy of specimen), (3) the influence of void ratio, gradation confining pressure, stress history and type of shear test on shear strength is reflected in values of the omega parameter, and (4) calculation of the coefficient of solid friction allows the establishment of the lower limit of the shear strength of a granular material.

Stress correlation function evolution in lattice solid elasto-dynamic models of shear and fracture zones and earthquake prediction

It has been argued that power-law time-to-failure fits for cumulative Benioff strain and an evolution in size-frequency statistics in the lead-up to large earthquakes are evidence that the crust behaves as a Critical Point (CP) system. If so, intermediate-term earthquake prediction is possible. However, this hypothesis has not been proven. If the crust does behave as a CP system, stress correlation lengths should grow in the lead-up to large events through the action of small to moderate ruptures and drop sharply once a large event occurs. However this evolution in stress correlation lengths cannot be observed directly. Here we show, using the lattice solid model to describe discontinuous elasto-dynamic systems subjected to shear and compression, that it is for possible correlation lengths to exhibit CP-type evolution. In the case of a granular system subjected to shear, this evolution occurs in the lead-up to the largest event and is accompanied by an increasing rate of moderate-sized events and power-law acceleration of Benioff strain release. In the case of an intact sample system subjected to compression, the evolution occurs only after a mature fracture system has developed. The results support the existence of a physical mechanism for intermediate-term earthquake forecasting and suggest this mechanism is fault-system dependent. This offers an explanation of why accelerating Benioff strain release is not observed prior to all large earthquakes. The results prove the existence of an underlying evolution in discontinuous elasto-dynamic, systems which is capable of providing a basis for forecasting catastrophic failure and earthquakes.

A conceptual model for sheet-flow drawn from rapid granular flow theories

33rd IAHR Congress: Water Engineering for a Sustainable Environment

Final Report Factors Influencing Stability of Granular Base Course Mixes

To evaluate the various factors influencing the stability of granular base course mixes, three primary goals were included in the project: (1) determination of a suitable and realistic laboratory method of compaction; (2) effect of gradation, density and mineralogy of the fines on sheara ing strength; and (3) possible improvement of the shear strength with organic and inorganic chemical stabilization additives.

Improvement of Granular Base Course Materials with Portland Cement, Special Report, HR-99

A highway base course may be defined as a layer of granular material which lies immediately below the wearing surface of a pavement and must possess high resistance to deformation in order to withstand pressures imposed by traffic. A material commonly used for base course construction is crushed limestone. Sources of limestone, acceptable for highway bases in the state of Iowa, occur almost entirely in the Pennsylvanian, Mississippian and Devonian strata. Performance records of the latter two have been quite good, while material from the Pennsylvanian stratum has failed on numerous occasions. The study reported herein is one segment of an extensive research program on compacted crushed limestone used for flexible highway base courses. The primary goals of the total study are: 1. Determination of a suitable and realistic laboratory method of compaction. 2. Effect of gradation, and mineralogy of the fines, on shearing strength. 3. Possible improvement of the shear strength with organic and inorganic chemical stabilization additives. Although the study reported herein deals primarily with the third goal, information gathered from work on the first two was required for this investigation. The primary goal of this study was the evaluation of various factors of stability of three crushed limestones when treated with small amounts of type I Portland cement. Investigation of the untreated materials has indicated that shear strength alone is not the controlling factor for stability of crushed stone bases. Thus the following observations were made in addition to shear strength parameters, to more adequately ascertain the stability of the cement treated materials: 1. Volume change during consolidation and shear testing. 2. Pore pressure during shear. The consolidated-undrained triaxial shear test was used for determination of the above factors.

Theoretical Analysis and Simulations of Vertically Vibrated Granular Materials

The dynamical properties ofshaken granular materials are important in many industrial applications where the shaking is used to mix, segregate and transport them. In this work asystematic, large scale simulation study has been performed to investigate the rheology of dense granular media, in the presence of gas, in a three dimensional vertical cylinder filled with glass balls. The base wall of the cylinder is subjected to sinusoidal oscillation in the vertical direction. The viscoelastic behavior of glass balls during a collision, have been studied experimentally using a modified Newton's Cradle device. By analyzing the results of the measurements, using numerical model based on finite element method, the viscous damping coefficient was determinedfor the glass balls. To obtain detailed information about the interparticle interactions in a shaker, a simplified model for collision between particles of a granular material was proposed. In order to simulate the flow of surrounding gas, a formulation of the equations for fluid flow in a porous medium including particle forces was proposed. These equations are solved with Large Eddy Simulation (LES) technique using a subgrid-model originally proposed for compressible turbulent flows. For a pentagonal prism-shaped container under vertical vibrations, the results show that oscillon type structures were formed. Oscillons are highly localized particle-like excitations of the granular layer. This self-sustaining state was named by analogy with its closest large-scale analogy, the soliton, which was first documented by J.S. Russell in 1834. The results which has been reportedbyBordbar and Zamankhan(2005b)also show that slightly revised fluctuation-dissipation theorem might apply to shaken sand, which appears to be asystem far from equilibrium and could exhibit strong spatial and temporal variations in quantities such as density and local particle velocity. In this light, hydrodynamic type continuum equations were presented for describing the deformation and flow of dense gas-particle mixtures. The constitutive equation used for the stress tensor provides an effective viscosity with a liquid-like character at low shear rates and a gaseous-like behavior at high shear rates. The numerical solutions were obtained for the aforementioned hydrodynamic equations for predicting the flow dynamics ofdense mixture of gas and particles in vertical cylindrical containers. For a heptagonal prism shaped container under vertical vibrations, the model results were found to predict bubbling behavior analogous to those observed experimentally. This bubbling behavior may be explained by the unusual gas pressure distribution found in the bed. In addition, oscillon type structures were found to be formed using a vertically vibrated, pentagonal prism shaped container in agreement with computer simulation results. These observations suggest that the pressure distribution plays a key rolein deformation and flow of dense mixtures of gas and particles under vertical vibrations. The present models provide greater insight toward the explanation of poorly understood hydrodynamic phenomena in the field of granular flows and dense gas-particle mixtures. The models can be generalized to investigate the granular material-container wall interactions which would be an issue of high interests in the industrial applications. By following this approach ideal processing conditions and powder transport can be created in industrial systems.

Granular flow in polymineralic rocks bearing sheet silicates: New evidence from natural examples

Natural deformation in carbonate mylonites bearing sheet silicates occurs via a complex interaction of granular flow and solution transfer processes and involves continuous cycles of dissolution, grain boundary diffusion, nucleation and growth. In this way, new sheet silicates (a) nucleate within voids formed by grain boundary sliding of calcite grains. (b) grow, and (c) rotate towards the shear plane. As a consequence, small mica grains show a wide range of orientations with respect to the shear plane, but moderate to large grains are subparallel both to each other and to the shear plane. Increases of average grain sizes with increasing temperature of sheet silicates in mica-rich layers is more pronounced than in mica-poor layers. In the calcitic matrix however, sheet silicates can only grow via solution-precipitation and mass transfer processes. Therefore, the observed grain size variability indicates drastic differences in mass transfer behavior between the individual layers, which might be related to differences in the fluid flux. Based on these observations, a conceptual model for the microfabric evolution in sheet silicate bearing mylonites is presented. © 2001 Elsevier Science B.V. All rights reserved.