Experimental Studies on the Mechanics of Cohesive Frictional Granular Media

The mechanical behaviour of cohesive-frictional granular materials is a combination of the strength pervading as intergranular friction (represented as an angle of internal friction - Phi), and the cohesion (C) between these particles. Most behavioral or constitutive models of this class of granular materials comprise of a cohesion and frictional component with no regard to the length scale i.e. from the micro structural models through the continuum models. An experimental study has been made on a model granular material, viz. angular sand with different weights of binding agents (varying degrees of cohesion) at multiple length scales to physically map this phenomenon. Cylindrical specimen of various diameters - 10, 20, 38, 100, 150 mm (and with an aspect ratio of 2) are reconstituted with 2, 4 and 8% by weight of a binding agent. The magnitude of this cohesion is analyzed using uniaxial compression tests and it is assumed to correspond to the peak in the normalized stress-strain plot. Increase in the cohesive strength of the material is seen with increasing size of the specimen. A possibility of ``entanglement'' occurring in larger specimens is proposed as a possible reason for deviation from a continuum framework.

The effect of base roughness on the development of a dense granular flow down an inclined plane

The development of the flow of a granular material down an inclined plane starting from rest is studied as a function of the base roughness. In the simulations, the particles are rough frictional spheres interacting via the Hertz contact law. The rough base is made of a random configuration of fixed spheres with diameter different from the flowing particles, and the base roughness is decreased by decreasing the diameter of the base particles. The transition from an ordered to a disordered flowing state at a critical value of the base particle diameter, first reported by Kumaran and Maheshwari Phys. Fluids 24, 053302 (2012)] for particles with the linear contact model, is observed for the Hertzian contact model as well. The flow development for the ordered and disordered flows is very different. During the development of the disordered flow for the rougher base, there is shearing throughout the height. During the development of the ordered flow for the smoother base, there is a shear layer at the bottom and a plug region with no internal shearing above. In the shear layer, the particles are layered and hexagonally ordered in the plane parallel to the base, and the velocity profile is well approximated by Bagnold law. The flow develops in two phases. In the first phase, the thickness of the shear layer and the maximum velocity increase linearly in time till the shear front reaches the top. In the second phase, after the shear layer encompasses the entire flow, there is a much slower increase in the maximum velocity until the steady state is reached. (C) 2013 AIP Publishing LLC.

Controlling wave propagation through nonlinear engineered granular systems

We study the fundamental dynamic behavior of a special class of ordered granular systems in order to design new, structured materials with unique physical properties. The dynamic properties of granular systems are dictated by the nonlinear, Hertzian, potential in compression and zero tensile strength resulting from the discrete material structure. Engineering the underlying particle arrangement of granular systems allows for unique dynamic properties, not observed in natural, disordered granular media. While extensive studies on 1D granular crystals have suggested their usefulness for a variety of engineering applications, considerably less attention has been given to higher-dimensional systems. The extension of these studies in higher dimensions could enable the discovery of richer physical phenomena not possible in 1D, such as spatial redirection and anisotropic energy trapping. We present experiments, numerical simulation (based on a discrete particle model), and in some cases theoretical predictions for several engineered granular systems, studying the effects of particle arrangement on the highly nonlinear transient wave propagation to develop means for controlling the wave propagation pathways. The first component of this thesis studies the stress wave propagation resulting from a localized impulsive loading for three different 2D particle lattice structures: square, centered square, and hexagonal granular crystals. By varying the lattice structure, we observe a wide range of properties for the propagating stress waves: quasi-1D solitary wave propagation, fully 2D wave propagation with tunable wave front shapes, and 2D pulsed wave propagation. Additionally the effects of weak disorder, inevitably present in real granular systems, are investigated. The second half of this thesis studies the solitary wave propagation through 2D and 3D ordered networks of granular chains, reducing the effective density compared to granular crystals by selectively placing wave guiding chains to control the acoustic wave transmission. The rapid wave front amplitude decay exhibited by these granular networks makes them highly attractive for impact mitigation applications. The agreement between experiments, numerical simulations, and applicable theoretical predictions validates the wave guiding capabilities of these engineered granular crystals and networks and opens a wide range of possibilities for the realization of increasingly complex granular material design.

Analysis and modeling of heaping behavior of granular mixtures within a computational mechanics framework

This paper presents a continuum model of the flow of granular material during filling of a silo, using a viscoplastic constitutive relation based on the Drucker-Prager plasticity yield function. The performed simulations demonstrate the ability of the model to realistically represent complex features of granular flows during filling processes, such as heap formation and non-zero inclination angle of the bulk material-air interface. In addition, micro-mechanical parametrizations which account for particle size segregation are incorporated into the model. It is found that numerical predictions of segregation phenomena during filling of a binary granular mixture agree well with experimental results. Further numerical tests indicate the capability of the model to cope successfully with complex operations involving granular mixtures.

Granular anchors under vertical loading - axial pull

Granular anchors are a relatively new concept in ground engineering with relatively little known regarding their load–displacement behaviour, failure modes, ultimate pullout capacity and also potential applications. A granular anchor consists of three main components: a base plate; tendon and compacted granular backfill. The tendon is used to transmit the applied load to the base plate which compresses the granular material to form the anchor. A study of the load–displacement response and ultimate pullout capacity of granular anchors constructed in intact lodgement till and made ground deposits is reported in this paper. Parallel tests were also performed on cast insitu concrete anchors which are traditionally used for anchoring purposes. A new method of analysis for the determination of the ultimate pullout capacity of granular anchors is presented and verified experimentally, with the dominant mode of failure controlled by the column length to diameter ratio. Granular anchors with L/D > 7 principally failed on bulging whereas short granular anchors failed on shaft resistance, with the latter mobilising similar pullout capacities as conventional concrete anchors.

Micromechanical analysis of cohesive granular materials using the discrete element method with an adhesive elasto-plastic contact model

An adhesive elasto-plastic contact model for the discrete element method with three dimensional non-spherical particles is proposed and investigated to achieve quantitative prediction of cohesive powder flowability. Simulations have been performed for uniaxial consolidation followed by unconfined compression to failure using this model. The model has been shown to be capable of predicting the experimental flow function (unconfined compressive strength vs. the prior consolidation stress) for a limestone powder which has been selected as a reference solid in the Europe wide PARDEM research network. Contact plasticity in the model is shown to affect the flowability significantly and is thus essential for producing satisfactory computations of the behaviour of a cohesive granular material. The model predicts a linear relationship between a normalized unconfined compressive strength and the product of coordination number and solid fraction. This linear relationship is in line with the Rumpf model for the tensile strength of particulate agglomerate. Even when the contact adhesion is forced to remain constant, the increasing unconfined strength arising from stress consolidation is still predicted, which has its origin in the contact plasticity leading to microstructural evolution of the coordination number. The filled porosity is predicted to increase as the contact adhesion increases. Under confined compression, the porosity reduces more gradually for the load-dependent adhesion compared to constant adhesion. It was found that the contribution of adhesive force to the limiting friction has a significant effect on the bulk unconfined strength. The results provide new insights and propose a micromechanical based measure for characterising the strength and flowability of cohesive granular materials. 

Granular Attrition due to Rotary Valve in a Pneumatic Conveying System

The rotary valve is a widely used mechanical device in many solids-handling industrial processes. However, it may also be responsible for most of the attrition effects occurring in a typical process. In this study, the attrition effects occurring in a rotary valve operating as a stand-alone device and as part of a pneumatic conveying system were investigated. In the former case granular attrition was carried out at three different rotary valve speeds and the experimental results obtained were found to be in good agreement with the Gwyn correlation. In the latter case three typical air flow rates were used in the pneumatic conveying system. The size distribution of the attrition product obtained at the lowest air flow rate used was not adequately described by the Gwyn correlation. The attrition process and mechanisms involved were analysed and the minimum size of the attrition product obtained from both modes of operations was found to be similar.

Study on Granular Dynamics in Vertically Vibrated Beds using Tracking Technique

The study of granular material is of great interest to many researchers in both engineering and science communities. The importance of such a study derives from its complex rheological character and also its significant role in a wide range of industrial applications, such as coal, food, plastics, pharmaceutical, powder metallurgy and mineral processing. A number of recent reports have been focused on the physics of non-cohesive granular material submitted to vertical vibration in either experimental or theoretical approaches. Such a kind of system can be used to separate, mix and dry granular materials in industries. It exhibits different instability behaviour on its surface when under vertical vibration, for example, avalanching, surface fluidization and surface wave, and these phenomena have attracted particular interest of many researchers. However, its fundamental understanding of the instability mechanism is not yet well-understood. This paper is therefore to study the dynamics of granular motion in such a kind of system using Positron Emission Particle Tracking (PEPT), which allows the motion of a single tracer particle to be followed in a non-invasive way. Features of the solids motion such as cycle frequency and dispersion index were investigated via means of authors’ specially-written programmes. Regardless of the surface behaviour, particles are found to travel in rotational movement in horizontal plane. Particle cycle frequency is found to increase strongly with increasing vibration amplitude. Particle dispersion also increased strongly with vibration amplitude. Horizontal dispersion is observed to always exceed vertical dispersion.

Granular flow in the marginal ice zone

The region of sea ice near the edge of the sea ice pack is known as the marginal ice zone (MIZ), and its dynamics are complicated by ocean wave interaction with the ice cover, strong gradients in the atmosphere and ocean and variations in sea ice rheology. This paper focuses on the role of sea ice rheology in determining the dynamics of the MIZ. Here, sea ice is treated as a granular material with a composite rheology describing collisional ice floe interaction and plastic interaction. The collisional component of sea ice rheology depends upon the granular temperature, a measure of the kinetic energy of flow fluctuations. A simplified model of the MIZ is introduced consisting of the along and across momentum balance of the sea ice and the balance equation of fluctuation kinetic energy. The steady solution of these equations is found to leading order using elementary methods. This reveals a concentrated region of rapid ice flow parallel to the ice edge, which is in accordance with field observations, and previously called the ice jet. Previous explanations of the ice jet relied upon the existence of ocean currents beneath the ice cover. We show that an ice jet results as a natural consequence of the granular nature of sea ice.

Statistical mechanics approaches to granular media: between micromechanics and macromechanics

The mechanical behavior of granular materials has been traditionally approached through two theoretical and computational frameworks: macromechanics and micromechanics. Macromechanics focuses on continuum based models. In consequence it is assumed that the matter in the granular material is homogeneous and continuously distributed over its volume so that the smallest element cut from the body possesses the same physical properties as the body. In particular, it has some equivalent mechanical properties, represented by complex and non-linear constitutive relationships. Engineering problems are usually solved using computational methods such as FEM or FDM. On the other hand, micromechanics is the analysis of heterogeneous materials on the level of their individual constituents. In granular materials, if the properties of particles are known, a micromechanical approach can lead to a predictive response of the whole heterogeneous material. Two classes of numerical techniques can be differentiated: computational micromechanics, which consists on applying continuum mechanics on each of the phases of a representative volume element and then solving numerically the equations, and atomistic methods (DEM), which consist on applying rigid body dynamics together with interaction potentials to the particles. Statistical mechanics approaches arise between micro and macromechanics. It tries to state which the expected macroscopic properties of a granular system are, by starting from a micromechanical analysis of the features of the particles and the interactions. The main objective of this paper is to introduce this approach.

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.

Behavior of Granular Materials Under Triaxial Compression with Pulsating Deviator Stress, HR-131, Special Report

Behavior of granular material subjected to repeated load triaxial compression tests is characterized by a model based on rate process theory. Starting with the Arrhenius equation from chemical kinetics, the relationship of temperature, shear stress, normal stress and volume change to deformation rate is developed. The proposed model equation includes these factors as a product of exponential terms. An empirical relationship between deformation and the cube root of the number of stress applications at constant temperature and normal stress is combined with the rate equation to yield an integrated relationship of temperature, deviator stress, confining pressure and number of deviator stress applications to axial strain. The experimental program consists of 64 repeated load triaxial compression tests, 52 on untreated crushed stone and 12 on the same crushed stone material treated with 4% asphalt cement. Results were analyzed with multiple linear regression techniques and show substantial agreement with the model equations. Experimental results fit the rate equation somewhat better than the integrated equation when all variable quantities are considered. The coefficient of shear temperature gives the activation enthalpy, which is about 4.7 kilocalories/mole for untreated material and 39.4 kilocalories/mole for asphalt-treated material. This indicates the activation enthalpy is about that of the pore fluid. The proportionality coefficient of deviator stress may be used to measure flow unit volume. The volumes thus determined for untreated and asphalt-treated material are not substantially different. This may be coincidental since comparison with flow unit volumes reported by others indicates flow unit volume is related to gradation of untreated material. The flow unit volume of asphalt-treated material may relate to asphalt cement content. The proposed model equations provide a more rational basis for further studies of factors affecting deformation of granular materials under stress similar to that in pavement subjected to transient traffic loads.

Micromechanical insight into the undrained instability of granular materials: deformation characteristics of geomaterials.

There is no agreement between experimental researchers whether the point where a granular material responds with a large change of stresses, strains or excess pore water pressure given a prescribed small input of some of the same variables defines a straight line or a curve in the stress space. This line, known as the instability line, may also vary in shape and position if the onset of instability is measured from drained or undrained triaxial tests. Failure of granular materials, which might be preceded by the onset of instability, is a subject that the geotechnical engineers have to deal with in the daily practice, and generally speaking it is associated to different phenomena observed not only in laboratory tests but also in the field. Examples of this are the liquefaction of loose sands subjected to undrained loading conditions and the diffuse instability under drained loading conditions. This research presents results of DEM simulations of undrained triaxial tests with the aim of studying the influence of stress history and relative density on the onset of instability in granular materials. Micro-mechanical analysis including the evolution of coordination numbers and fabric tensors is performed aiming to gain further insight on the particle-scale interactions that underlie the occurrence of this instability. In addition to provide a greater understanding, the results presented here may be useful as input for macro-scale constitutive models that enable the prediction of the onset of instability in boundary value problems.

"Candidatus Similichlamydia laticola", a novel Chlamydia-like agent of epitheliocystis in seven consecutive cohorts of farmed Australian barramundi, lates calcarifer (Bloch)

Six consecutively hatched cohorts and one cohort of pre-hatch eggs of farmed barramundi (Lates calcarifer) from south Australia were examined for Chlamydia-like organisms associated with epitheliocystis. To identify and characterise the bacteria, 59 gill samples and three pre-hatch egg samples were processed for histology, in situ hybridisation and 16S rRNA amplification, sequencing and comprehensive phylogenetic analysis. Cases of epitheliocystis were observed microscopically and characterised by membrane-enclosed basophilic cysts filled with a granular material that caused hypertrophy of the epithelial cells. In situ hybridisation with a Chlamydiales-specific probe lead to specific labelling of the epitheliocystis inclusions within the gill epithelium. Two distinct but closely related 16S rRNA chlamydial sequences were amplified from gill DNA across the seven cohorts, including from pre-hatch eggs. These genotype sequences were found to be novel, sharing 97.1 - 97.5% similarity to the next closest 16S rRNA sequence, Ca. Similichlamydia latridicola, from Australian striped trumpeter. Comprehensive phylogenetic analysis of these genotype sequences against representative members of the Chlamydiales order and against other epitheliocystis agents revealed these Chlamydia-like organisms to be novel and taxonomically placed them within the recently proposed genus Ca. Similichlamydia. Following Fredricks and Relman's molecular postulates and based on these observations, we propose the epitheliocystis agents of barramundi to be known as "Candidatus Similichlamydia laticola" (sp. nov.).

Effect of heat-treatment on strength of clays

The main objective of this investigation was to understand the strength development of clays below fusion or vitrification temperatures of 900°C. The other objective was to establish threshold temperatures to produce a satisfactory construction material from clayey sediments from the Western Beaufort Sea for shore protection of artificial islands with minimum expense of thermal energy. Studies were, therefore, conducted using kaolinite, bentonite, and a clayey sediment from the Beaufort Sea. Unconfined-compressive-strength tests were conducted on clay samples heat treated from 110 to 700°C. Furthermore, to understand the factors responsible for strength-development-thermogravimetric studies and pore-size analysis, using mercury porosimetry, were also conducted. A gradual increase in strength was obtained with an increase in firing temperature. However, substantial and permanent increase in strength occurred only after dehydroxylation of all the clays studied; Clay samples heated to temperatures above dehydroxylation became resistant to disintegration upon immersion in water. Results indicate that the clayey sediments from Western Beaufort Sea have to be heat treated to about 600°C to produce granular material for use as a fill or shore-protection material for artificial islands.