Stress transmission in wet granular materials

Dynamic compaction of cohesive granular materials: scaling behavior and bonding structures

The compaction of cohesive granular materials is a common operation in powder-based manufacture of many products. However, the influence of particle-scale parameters such as bond strength on the packing structure and the general scaling of the compaction process are still poorly understood. We use particle dynamics simulations to analyze jammed configurations obtained by dynamic compaction of sticky particles under a fixed compressive pressure for a broad range of system parameter values. We show that relative porosity, representing the relative importance of porosity with respect to its minimum and maximum values, is a unique function of a modified cohesion number that combines adhesion force, confining pressure, and particle size, as well as contact stiffness, which is often assumed to be ineffective but is shown here to play an essential role in compaction. An asymmetric sigmoidal form based on two power laws provides an excellent fit to the data. The statistical properties of the bond network reveal self-balanced force structures and an exponential fall-off of the number of both tensile and compressive forces. Remarkably, the properties of the bond network depend on the cohesion number rather than the modified cohesion number, implying that similar bond network characteristics are compatible with a broad range of porosities mainly due to the effect of contact stiffness. We also discuss the origins of data points escaping the general scaling of porosity and show that they reflect either finite system size or rigid confining walls.

Collapse dynamics and deposition morphology of low-viscocohesive granular columns on a rough horizontal surface

Using the three-dimensional discrete element method, we numerically investigate the collapse dynamics and deposition morphology of low-viscocohesive granular columns on a rough-horizontal plane by systematically varying a broad range of values of the initial column aspect ratio, cohesive stress, and liquid viscosity. The results show that the kinetic energy, half runout time, and runout distance increase with increasing the initial column aspect ratio but decrease with increasing the cohesive and viscous effects of the binding liquid, while the toe angle and deposit height decrease with increasing the aspect ratio and increase with increasing cohesive stress and liquid viscosity. Remarkably, by defining a dimensionless scaling number that incorporates the Bond number and initial column aspect ratio, this allows us to nicely describe the kinetic energy, half runout time, deposition height, runout distance, and toe angle. These unified descriptions may provide insights into the physical properties of the collapse dynamics and deposition morphology of low-viscocohesive granular columns, leading to good explanations of the complex properties of natural disaster events.

The combined effect of cohesion and finite size on the collapse of wet granular columns

The collapse of low-saturation liquid-containing granular materials is prevalent in nature and industrial processes, and understanding the associated transient dynamics is extremely important for exploring such complex flow processes. In this paper, the collapse of a finite-size wet granular column is systematically studied and the determinants affecting its dynamics are analyzed based on the discrete element model for wet particles and the corresponding small-scale experiments. With the aid of parametric analysis, the dimensionless cohesion parameter containing the system size and grain-scale bond number is proposed, and its relevance in characterizing column stability and collapse dynamics of wet granular materials is further confirmed. For the collapse of wet granular columns with a fixed aspect ratio, the initial height contained in the cohesion parameter is verified to be a manifestation of the finite size effect, which is present in a wet granular collapse and is coupled with the cohesive effect. Such a coupling effect is taken into account in our proposed scaling laws that can be applied to uniformly describe the deposit morphology of wet granular columns after collapse.

Reorientation of Suspended Ceramic Particles in Robocasted Green Filaments during Drying

This work considers the fabrication of ceramic parts with the help of an additive manufacturing process, robocasting, in which a paste with suspended particles is robotically extruded. Within the final part, the material properties depend on the orientation of the particles. A prediction of the particle orientation is challenging as the part usually undergoes multiple processing steps with varying contributions to the orientation. As the main contribution to the final particle orientation arises from the extrusion process, many corresponding prediction models have been suggested. Robocasting involves, however, further processing steps that are less studied as they have a smaller influence on the orientation. One of the processing steps is drying by natural convection, which follows directly after the extrusion process. A quantification of the reorientation that occurs during drying is mostly unknown and usually neglected in the models. Therefore, we studied the amount of reorientation of suspended particles in robocasted green filaments during drying in detail. For our study, we applied the discrete element method, as it meets various requirements: The exact particle geometry can be resolved precisely; particle–particle interactions can be described; the paste composition is reproduced exactly; the initial particle orientation can be set in accordance with the prediction from the analytical models for the extrusion part; macroscopic force laws exist to represent capillary forces due to the remaining fluid phase that remains during drying. From our study, we concluded that the magnitude of particle reorientation during drying is small compared to the orientation occurring during the extrusion process itself. Consequently, reorientation during drying might further be neglected within analytical orientation prediction models.

Effects of size polydispersity on segregation of spherical particles in rotating drum

To get insight into the segregation process of a polydisperse granular materials flow, we numerically investigated the migration process of particles in a rotating drum operating in the rolling regime by means of the discrete element method. Particle migration is analyzed through the variation of the proportion of particles in different zones where the flow property is characterized. The proportion of particles in different zones of the drum shows to increase in the center of the flow radially and axially where a higher concentration of small particles is observed, while its decreases in other zones with a higher concentration of larger particles. Interestingly, we find that the migration process of particles leads to radial and axial segregation which is caused by a combination between the exerted fluctuation forces on particles and its surrounding pressure gradient.

Scaling behavior of the tensile strength of viscocohesive granular aggregates

We numerically analyze the tensile strength of a single wet agglomerate modeled as a viscocohesive aggregate impacting a flat surface by using the discrete-element simulations. The viscocohesive agglomerate composed of primary spherical particles with the inclusion of the interstitial liquid in the form of the capillary bridges characterized by the cohesive and viscous forces between particles is extracted from a cuboidal sample of granular materials by applying a spherical probe. The tensile strength is measured from the impact test of a wet agglomerate by systematically varying different values of the surface tension of the interstitial liquid, the liquid viscosity, and the impact speed. We show that the tensile stress increases immediately when the collision occurs between the agglomerate and the flat surface. The peak of the tensile stress obtained after the collision, then decreases smoothly with increasing the particle movement. The maximum tensile stress is defined to be the tensile strength of such agglomerate. It is remarkable that the normalized tensile strength of such agglomerate can be well described as a function of a dimensionless impact number that incorporates the capillary number and Stokes number (calculated from the surface tension and the viscosity of the liquid and the impact rate of the agglomerate), thus providing the confirmation for the unified representation of the liquid properties and the impact rate of wet granular media.

Anisotropic nature of the capillary stress tensor

The paper describes a micromechanical approach that explores the anisotropic nature of the capillary stress tensor and its evolution in pendular granular materials via Discrete Element Modeling (DEM) simulations. Dimensionless parameters are used to address the conditions under which the contribution of capillarity (or cohesive interparticle forces) to the stress transmission within a Representative Elementary Volume (REV) is expected to be considerable. From a series of suction-controlled conventional triaxial tests, numerical results show that the significance of the capillary stress and the relative magnitude of its mean to deviatoric components is directly connected to the characteristic particle size and applied stress. In addition, it is shown that the anisotropic character of the capillary stress tensor intensifies with increasing suction. Furthermore, a simple shear test is conducted at constant mean stress to reveal the development of deviatoric capillary stresses in the absence of any change in mean stress, which cannot be captured by the commonly used Bishop’s stress expression.

The role of inter-particle friction on rheology and texture of wet granular flows

In order to get insight into the rheology and texture of rough unsaturated granular flows, we study the effects of the inter-particle friction coefficient on the macroscopic attributes and the texture variables of steady-state shearing flow of wet granular materials by relying on three-dimensional (3D) particle dynamics simulations. The macroscopic attributes are characterized by the macroscopic friction coefficient, macroscopic cohesion, and the packing fraction. The microstructural variables are characterized by the fabric and force anisotropies, the coordination number, and the stress transmission ratio. We show that the macroscopic observables behave as a function of the inertial number as a dry case for different values of the inter-particle friction coefficient. In particular, the macroscopic friction coefficient increases and the packing fraction decreases rapidly for small values of the friction coefficient, then they almost reach plateaus for higher values of the friction coefficient. Interestingly, all the macroscopic observables nicely behave as a function of the small values of the friction coefficient. Similarly, we also observe these characteristics for the fabric and force anisotropies and the coordination number as well as the stress transmission ratio which reflects the intermediate relationship between the microstructure and the mechanical behavior of such flows.

Stress transmission in cemented bidisperse granular materials

We analyze stress distributions in a two-dimensional bidisperse cemented granular packing for a broad range of the values of particle-size ratio, the volumes of large and small particles, and the amount of cementing matrix. In such textured porous materials, the stress concentration, which controls the fracture and fragmentation of the material under tensile loading or in grinding processes, reflects not only the porosity but also the contact network of the particle phase and the resulting stress chains. By means of peridynamic simulations under tensile loading, we show how both the texture and stress distribution depend on size ratio, volume ratio, and the amount of the cementing matrix. In particular, the volume fraction of the class of small particles plays a key role in homogenizing stresses across the system by reducing porosity. Interestingly, the texture controls not only the porosity but also the distribution of pores inside the system with its statistical variability, found to be strongly correlated with the homogeneity of stresses inside the large particles. The most homogeneous stress distribution occurs for the largest size ratio and largest volume fraction of small particles, corresponding to the lowest pore size dispersion and the cushioning effect of small particles and its similar role to the binding matrix for stress redistribution across the packing. At higher porosity, the tensile stresses above the mean stress fall off exponentially in all phases with an exponent that strongly depends on the texture. The exponential part broadens with decreasing matrix volume fraction and particle-size ratio. These correlations reveal the strong interplay between size polydispersity and the cohesive action of the binding matrix for stress distribution, which is significant for the behavior of textured materials in grinding operations.

Additive rheology of complex granular flows

Granular flows are omnipresent in nature and industrial processes, but their rheological properties such as apparent friction and packing fraction are still elusive when inertial, cohesive and viscous interactions occur between particles in addition to frictional and elastic forces. Here we report on extensive particle dynamics simulations of such complex flows for a model granular system composed of perfectly rigid particles. We show that, when the apparent friction and packing fraction are normalized by their cohesion-dependent quasistatic values, they are governed by a single dimensionless number that, by virtue of stress additivity, accounts for all interactions. We also find that this dimensionless parameter, as a generalized inertial number, describes the texture variables such as the bond network connectivity and anisotropy. Encompassing various stress sources, this unified framework considerably simplifies and extends the modeling scope for granular dynamics, with potential applications to powder technology and natural flows.

Granular materials are abundant in nature, but we haven’t fully understood their rheological properties as complex interactions between particles are involved. Here, Vo et al. show that granular flows can be described by a generalized dimensionless number based on stress additivity.

Evolution of wet agglomerates inside inertial shear flow of dry granular materials

We use particle dynamics simulations to investigate the evolution of a wet agglomerate inside homogeneous shear flows of dry particles. The agglomerate is modeled by introducing approximate analytical expressions of capillary and viscous forces between particles in addition to frictional contacts. During shear flow, the agglomerate may elongate, break, or be eroded by loss of its capillary bonds and primary particles. By systematically varying the shear rate and surface tension of the binding liquid, we characterize the rates of these dispersion modes. All the rates increase with increasing inertial number of the flow and decreasing cohesion index of the agglomerate. We show that the data points for each mode collapse on a master curve for a dimensionless scaling parameter that combines the inertial number and the cohesion index. The erosion rate vanishes below a cutoff value of the scaling parameter. This leads to a power-law borderline between the vanishing erosion states and erosion states in the phase space defined by the inertial number and the cohesion index.

Agglomeration of wet particles in dense granular flows

In order to get insight into the wet agglomeration process, we numerically investigate the growth of a single granule inside a dense flow of an initially homogeneous distribution of wet and dry particles. The simulations are performed by means of the discrete element method and the binding liquid is assumed to be transported by the wet particles, which interact via capillary and viscous force laws. The granule size is found to be an exponential function of time, reflecting the conservation of the amount of liquid and the decrease of the number of available wet particles inside the flow during agglomeration. We analyze this behavior in terms of the accretion and erosion rates of wet particles for a range of different values of material parameters such as mean particle size, size polydispersity, friction coefficient and liquid viscosity. In particular, we propose a phase diagram of the granule growth as a function of the mean primary particle diameter and particle size span, which separates the parametric domain in which the granule grows from the domain in which the granule does not survive.

Force Transmission Modes of Non-Cohesive and Cohesive Materials at the Critical State

This paper investigates the force transmission modes, mainly described by probability density distributions, in non-cohesive dry and cohesive wet granular materials by discrete element modeling. The critical state force transmission patterns are focused on with the contact model effect being analyzed. By shearing relatively dense and loose dry specimens to the critical state in the conventional triaxial loading path, it is observed that there is a unique critical state force transmission mode. There is a universe critical state force distribution pattern for both the normal contact forces and tangential contact forces. Furthermore, it is found that using either the linear Hooke or the non-linear Hertz model does not affect the universe force transmission mode, and it is only related to the grain size distribution. Wet granular materials are also simulated by incorporating a water bridge model. Dense and loose wet granular materials are tested, and the critical state behavior for the wet material is also observed. The critical state strength and void ratio of wet granular materials are higher than those of a non-cohesive material. The critical state inter-particle distribution is altered from that of a non-cohesive material with higher probability in relatively weak forces. Grains in non-cohesive materials are under compressive stresses, and their principal directions are mainly in the axial loading direction. However, for cohesive wet granular materials, some particles are in tension, and the tensile stresses are in the horizontal direction on which the confinement is applied. The additional confinement by the tensile stress explains the macro strength and dilatancy increase in wet samples.

Flow of wet granular materials: A numerical study

We simulate dense assemblies of frictional spherical grains in steady shear flow under controlled normal stress P in the presence of a small amount of an interstitial liquid, which gives rise to capillary menisci, assumed isolated (pendular regime), and attractive forces, which are hysteretic: Menisci form at contact, but do not break until grains are separated by a finite rupture distance. The system behavior depends on two dimensionless control parameters, inertial number I and reduced pressure P*=aP/(πΓ), comparing confining forces ∼a2P to meniscus tensile strength F0=πΓa, for grains of diameter a joined by menisci with surface tension Γ. We pay special attention to the quasistatic limit of slow flow and observe systematic, enduring strain localization in some of the cohesion-dominated (P*∼0.1) systems. Homogeneous steady flows are characterized by the dependence of internal friction coefficient μ* and solid fraction Φ on I and P*. We also record normal stress differences, fairly small but not negligible and increasing for decreasing P*. The system rheology is moderately sensitive to saturation within the pendular regime, but would be different in the absence of capillary hysteresis. Capillary forces have a significant effect on the macroscopic behavior of the system, up to P* values of several units, especially for longer force ranges associated with larger menisci. The concept of effective pressure may be used to predict an order of magnitude for the strong increase of μ* as P* decreases but such a crude approach is unable to account for the complex structural changes induced by capillary cohesion, with a significant decrease of Φ and different agglomeration states and anisotropic fabric. Likewise, the Mohr-Coulomb criterion for pressure-dependent critical states is, at best, an approximation valid within a restricted range of pressures, with P*≥1. At small enough P*, large clusters of interacting grains form in slow flows, in which liquid bonds survive shear strains of several units. This affects the anisotropies associated with different interactions and the shape of function μ*(I), which departs more slowly from its quasistatic limit than in cohesionless systems (possibly explaining the shear banding tendency).

Cohesive granular materials composed of nonconvex particles

The macroscopic cohesion of granular materials made up of sticky particles depends on the particle shapes. We address this issue by performing contact dynamics simulations of 2D packings of nonconvex aggregates. We find that the macroscopic cohesion is strongly dependent on the strain and stress inhomogeneities developing inside the material. The largest cohesion is obtained for nearly homogeneous deformation at the beginning of unconfined axial compression and it evolves linearly with nonconvexity. Interestingly, the aggregates in a sheared packing tend to form more contacts with fewer neighboring aggregates as the degree of nonconvexity increases. We also find that shearing leads either to an isotropic distribution of tensile contacts or to the same privileged direction as that of compressive contacts.

Tensile strength and fracture of cemented granular aggregates

Cemented granular aggregates include a broad class of geomaterials such as sedimentary rocks and some biomaterials such as the wheat endosperm. We present a 3D lattice element method for the simulation of such materials, modeled as a jammed assembly of particles bound together by a matrix partially filling the interstitial space. From extensive simulation data, we analyze the mechanical properties of aggregates subjected to tensile loading as a function of matrix volume fraction and particle-matrix adhesion. We observe a linear elastic behavior followed by a brutal failure along a fracture surface. The effective stiffness before failure increases almost linearly with the matrix volume fraction. We show that the tensile strength of the aggregates increases with both the increasing tensile strength at the particle-matrix interface and decreasing stress concentration as a function of matrix volume fraction. The proportion of broken bonds in the particle phase reveals a range of values of the particle-matrix adhesion and matrix volume fraction for which the cracks bypass the particles and hence no particle damage occurs. This limit is shown to depend on the relative toughness of the particle-matrix interface with respect to the particles.

Entangled granular media

We study the geometrically induced cohesion of ensembles of granular "u particles" that mechanically entangle through particle interpenetration. We vary the length-to-width ratio l/w of the u particles and form them into freestanding vertical columns. In a laboratory experiment, we monitor the response of the columns to sinusoidal vibration (with peak acceleration Γ). Column collapse occurs in a characteristic time τ which follows the relation τ∝exp(Γ/Δ). Δ resembles an activation energy and is maximal at intermediate l/w. A simulation reveals that optimal strength results from competition between packing and entanglement.

Bond anisotropy and cohesion of wet granular materials

We analyse the Coulomb cohesion of wet granular materials and its relationship with the distribution of capillary bonds between particles. We show that, within a harmonic representation of the bond and force states, the shear strength can be expressed as a state equation in terms of internal anisotropy parameters. This formulation involves a dependence of the shear strength on loading direction and leads to the fragile behaviour of granular materials. The contact dynamics simulations of a wet material, in which a capillary force law is prescribed, are in excellent agreement with the predictions of this model. We find that the fragile character decreases as the local adhesion is increased or the mean stress is decreased.

Characterization of interparticle forces in the packing of cohesive fine particles

We numerically investigate force structures in the packing of fine cohesive particles using the discrete element method. By changing the particle size and therefore the van der Waals force, the effect of cohesion on the normal contact force and the total normal force, which is the sum of the normal contact forces and the van der Waals forces, is analyzed. It is shown that, with decreasing particle size, the normal contact forces become more uniform and have a narrower and more symmetric distribution, while the distributions of the total normal forces widen. Spatial correlation between the interparticle forces exists for the packing of coarse noncohesive particles. As the particle size decreases, this correlation becomes weaker for the contact forces but stronger for the total normal forces. A comparison between the effective weight of particles and the internal force structure suggests that there are differences between the particle-particle and particle-wall forces. The bimodal distribution of the effective weight indicates that there may exist two phases in the packings when cohesion is present, governed by the compressive and tensile stresses.

Computer simulation of model cohesive powders: plastic consolidation, structural changes, and elasticity under isotropic loads

The quasistatic behavior of a simple two-dimensional model of a cohesive powder under isotropic loads is investigated by discrete element simulations. We ignore contact plasticity and focus on the effect of geometry and collective rearrangements on the material behavior. The loose packing states, as assembled and characterized in a previous numerical study [Gilabert, Roux, and Castellanos, Phys. Rev. E 75, 011303 (2007)], are observed, under growing confining pressure P , to undergo important structural changes, while solid fraction Phi irreversibly increases (typically, from 0.4-0.5 to 0.75-0.8). The system state goes through three stages, with different forms of the plastic consolidation curve, i.e., Phi as a function of the growing reduced pressure P;{*}=PaF_{0} , defined with adhesion force F0 and grain diameter a . In the low-confinement regime (I), the system undergoes negligible plastic compaction, and its structure is influenced by the assembling process. In regime II the material state is independent of initial conditions, and the void ratio varies linearly with lnP [i.e., Delta(1Phi)=lambdaDelta(lnP;{*}) ], as described in the engineering literature. Plasticity index lambda is reduced in the presence of a small rolling resistance (RR). In the last stage of compaction (III), Phi approaches an asymptotic, maximum solid fraction Phi_{max} , as a power law Phi_{max}-Phi proportional, variant(P;{*});{-alpha} , with alpha approximately 1 , and properties of cohesionless granular packs are gradually retrieved. Under consolidation, while the range xi of fractal density correlations decreases, force patterns reorganize from self-balanced clusters to force chains, with correlative evolutions of force distributions, and elastic moduli increase by a large amount. Plastic deformation events correspond to very small changes in the network topology, while the denser regions tend to move like rigid bodies. Elastic properties are dominated by the bending of thin junctions in loose systems. For growing RR those tend to form particle chains, the folding of which, rather than tensile ruptures, controls plastic compaction.