Search

Jamming transitions and the rheology of granular avalanches in fluids are investigated using experiments and numerical simulations. Simulations use the lattice-Boltzmann method coupled with the discrete element method, providing detailed stress and deformation data. Both simulations and experiments present a perfect match with each other in carefully conducted deposition experiments, validating the simulation method. We analyze transient rheological laws and jamming transitions using our recently introduced length-scale ratio $G$. $G$ serves as a unified metric for the pressure and shear rate capturing the dynamics of sheared fluid-granular systems. Two key transition points, $G_{Y}$ and $G_{0}$, categorize the material's state into solid-like, creeping, and fluid-like states. Yielding at $G_{Y}$ marks the transition from solid-like to creeping, while $G_{0}$ signifies the shift to the fluid-like state. The $\mu-G$ relationship converges towards the equilibrium $\mu_{eq}(G)$ after $G>G_0$ showing the critical point where the established rheological laws for steady states apply during transient conditions.

Granular column collapse is a simple but important problem to the granular material community, due to its links to dynamics of natural hazards, such as landslides and pyroclastic flows, and many industrial situations, as well as its potential of analyzing transient and non-local rheology of granular flows. This article proposes a new dimensionless number to describe the run-out behaviour of granular columns on inclined planes based on both previous experimental data and dimensional analysis. With the assistance of the sphero-polyhedral discrete element method (DEM), we simulate inclined granular column collapses with different initial aspect ratios, inter-particle frictions, and initial solid fractions on inclined planes with different inclination angles (2.5$^{\circ}$ - 20.0$^{\circ}$) to verify the proposed dimensional analysis. Detailed analyses are further provided for better understanding of the influence of different initial conditions and boundary conditions, and to help unify the description of the relative run-out distances of systems with different inclination angles. This work determines the similarity and unity between granular column collapses on inclined planes and those on horizontal planes, and helps investigate the transient rheological behaviour of granular flows, which has direct relevance to various natural and engineering systems.

In this letter, following an extensive experimental validation, we perform constant-volume shearing simulations of non-Brownian granular suspensions using the discrete element method coupled with the lattice Boltzmann method. We choose a wide range of solid fractions, shear rates, fluid viscosities, particle sizes, and inter-particle frictional coefficients to obtain a scaling solution for the viscous behavior of suspensions in both dilute and dense regimes. This letter demonstrates that, with a proposed dilute-dense transitional solid fraction, $\phi_d$, there exists a strong correlation between the inverse relative viscosity and the shear stress. This work incorporates both the $\phi$-dependence and the $\dot{\gamma}$-dependence of suspension viscosity in a universal framework, which provides a scaling solution for granular suspensions across dilute and dense regimes and sheds light on the dilute-dense transition mechanisms.

The mean stress tensor is essential to investigate the dynamics of granular material. In this paper, we use Hamilton's principle of least action to derive the averaged stress tensor of discrete granular assemblies subjected to hydraulic force fields, as well as rigorous conditions for a proper definition of the Representative Volume Element (RVE). The main goal behind our efforts is to upscale particle physics into a sound stress tensor for systems involving the complex interaction between grains and water. We identify the contributions from the unbalanced forces, hydraulic forces, gravity, external forces, and particle fluctuation to the mean stress tensor. In doing so, it is convenient to separate the influence of different force fields when the granular system is subjected to complex environments, e.g., subaqueous conditions. The obtained formula is then validated by triaxial test simulations of dry and saturated granular systems using the Discrete Element Method (DEM) and the Lattice-Boltzmann Method (LBM). The results show that the deduced formula can accurately calculate the stress tensor of discrete assemblies with various body-force fields. We used validated DEM-LBM simulations of submerged granular column collapses to explore the physics happening at the grain scale with this mathematical formalism and showcase its potential. We provide a new perspective based on the granular assembly scale to pursue the fluid-solid interaction. Due to the importance of stress analysis in the constitutive modelling of granular materials, this work could help to better obtain the stress-strain relationship of saturated or submerged granular systems.

In this communication, we present a new lengthscale-based rheology for dense sheared particle suspensions as they transition from inertial- to viscous-dominated. We derive a lengthscale ratio using straightforward physics-based considerations for a particle subjected to pressure and drag forces. In doing so, we demonstrate that an appropriately chosen length-scale ratio intrinsically provides a consistent relationship between normal stress and system proximity to its ''jammed'' or solid-like state, even as a system transitions between inertial and viscous states, captured by a variable Stokes number.

The behavior of granular column collapses is associated with the dynamics of geohazards, such as debris flows, landslides, and pyroclastic flows, yet its underlying physics is still not well understood. In this paper, we explore granular column collapses using the spheropolyhedral discrete element method (DEM), where the system contains two types of particles with different frictional properties. We impose three different mixing ratios and multiple different particle frictional coefficients, which lead to different run-out distances and deposition heights. Based on our previous work and a simple mixture theory, we propose a new effective initial aspect ratio for the bi-frictional granular mixture, which helps unify the description of the relative run-out distances. We analyze the kinematics of bi-frictional granular column collapses and find that deviations from classical power-law scaling in both the dimensionless terminal time and the dimensionless time when the system reaches the maximum kinetic energy may result from differences in the initial solid fraction and initial structures. To clarify the influence of initial states, we further decrease the initial solid fraction of granular column collapses, and propose a trial function to quantitatively describe its influence. Due to the utilization of a simple mixture theory of contact occurrence probability, this study can be associated with the friction-dependent rheology of granular systems and friction-induced granular segregations, and further generalized into applications with multiple species of particles in various natural and engineering mixtures.

We investigate granular column collapses with different column cross-sections and associate the cross-section shape influence with a finite-size analysis. Previous research, confined to initially circular configurations, reviewed the importance of granular column collapse studies and concluded that the run-out distance scales with the initial aspect ratio. In this work, granular columns with three different types of initial column cross-sections (square, equilateral triangular, and rectangular cross-sections) are simulated using the discrete element method (DEM). We explore how non-circular cross-sections lead to different run-out distances. Based on the previously obtained finite-size analysis, we further link the initial radius in different directions to the relative size of a column and perform the finite-size analysis to explain the cross-section influence. In the end, a universal relationship, which includes frictional properties, relative sizes, and the cross-section influence, is proposed. Our results are shown to have direct relevance to various natural and engineering systems.

Understanding the rheology of granular assemblies is important for natural and engineering systems, but the relationship between inter-particle friction (or microscopic friction) and macroscopic friction is still not well understood. In this study, using the the discrete element method (DEM) with spherical particles and realistic contact laws, we investigate the mechanics of granular systems with a wide range of inter-particle frictional coefficients and aim to establish a friction-dependent rheology for dry granular flows. The corresponding results show that increasing inter-particle friction dramatically increases the effective frictional coefficient, $\mu_{\textrm{eff}}$, while decreasing the solid fraction of the system and increasing the transitional inertial number that marks the division of quasi-static regimes and intermediate flow regimes. We further propose a new dimensionless number, $\mathcal{M}$, as a ratio between the inertial effect and frictional effect, which is similar to the effective aspect ratio in granular column collapses, and unifies the influence of inter-particle friction with the inertial number. We then establish a relationship between $\mathcal{M}$ and the dimensionless granular temperature, $\Theta$, to further universalize the influence of inter-particle frictions. Such study can broaden the application of the $\mu(I)$ rheology in natural and engineering systems and help establish a more general constitutive model for complex granular systems.

In this letter, we investigate the transient rheological behavior of immersed granular flows using both experiments of submerged granular column collapses and corresponding numerical simulations. The simulations are performed with the lattice-Boltzmann method (LBM) coupled with the discrete element method (DEM) and provide a significant amount of data of the stress and deformation conditions at different positions and times during the granular collapse. We derive a new dimensionless number $\mathcal{G}$ that can unify the rheology of transient granular flows in different regimes for all the simulation data points. $\mathcal{G}$ smoothly transforms from an inertial number into a viscous number, unifying both extremes of the rheology law. We also show the need to introduce the kinetic stresses to achieve a universal relation. The findings establish a transient constitutive framework for visco-inertial granular flows and are important for a better understanding of granular-fluid mixtures in both natural and engineering situations.

In this letter, we focus on the size effect of granular column collapses, which are potentially connected to the dynamics of complex geophysical flows, even if the link between microscopic structures of granular assemblies and their macroscopic behaviors is still not well understood. Using the sphero-polyhedral discrete element method (DEM), we show that the column radius/grain size ratio has a strong influence on the collapse behavior. A finite-size analysis, which is inspired by a phase transition around an inflection point, is performed to obtain a general scaling equation with critical exponents for run-out distances. We further link the size effect with the strong force network and formalize a correlation length scale that exponentially scales with the effective aspect ratio. Such a scaling solution shows similarities with the percolation problem of two-dimensional random networks and can be extended to other similar natural and engineering systems.

Granular column collapses result in an array of flow phenomena and deposition morphologies, the understanding of which brings insights into studying granular flows in both natural and engineering systems. Guided by experiments, we carried out computational studies with the discrete element method (DEM) to identify fundamental links between the macro-scale behavior and micro-scale properties of granular columns. A physics-based non-dimensional number combining particle and bulk properties of the column, $\alpha_{\textrm{eff}}$, was found to determinate three collapse regimes (quasi-static, inertial, and liquid-like), revealing universal trends of flow regimes and deposition morphologies under different conditions. This non-dimensional number represents physically the competing inertial and frictional effects that govern the behavior of the granular column collapse, including energy conversion and dissipation. The finding is important for understanding quantitatively the flow of granular materials and their deposits.

In this letter, instead of investigating the compaction of granular materials under tapping or cyclic shearing, we focused more on the gyratory shearing compaction, where particles are subjected to constant pressure and shear rate continuously. Such a phenomenon is crucial to the compaction of asphalt mixtures or soil in civil engineering and can be extended to other areas, such as powder processing and pharmaceutical engineering. In this study, we found that the gyratory speed or interstitial fluid viscosity has almost no impact on the compaction behavior, while the pressure plays a more important role. Additionally, it is the inertial time scale which dictates the compaction behavior under gyratory shearing in most cases, whereas the viscous time scale can also have influence in some conditions.

We investigate the link between particle-scale dynamics and bulk behaviors of thickly-coated particle-fluid flows using computational simulations. We find that, similar to dense fully-saturated slurries, the form the rheology takes in these systems can carry signatures of interparticle collisions and/or interparticle viscous dynamics that vary solid fraction. However, we find significant qualitative and quantitative differences in the transitions of the system between what might be called viscous, collisional, and visco-collisional behaviors. We show how these transitions arise from changes in the `fabric,' e.g. strong force network, and the `granular temperature,' or fluctuation energy, and suggest extension of these frameworks to better elucidate other particle-fluid behaviors.