Insights into the rheology of cohesive granular media.

Characterization and prediction of the "flowability" of powders are of paramount importance in many industries. However, our understanding of the flow of powders like cement or flour is sparse compared to the flow of coarse, granular media like sand. The main difficulty arises because of the presence of adhesive forces between the grains, preventing smooth and continuous flows. Several tests are used in industrial contexts to probe and quantify the "flowability" of powders. However, they remain empirical and would benefit from a detailed study of the physics controlling flow dynamics. Here, we attempt to fill the gap by performing intensive discrete numerical simulations of cohesive grains flowing down an inclined plane. We show that, contrary to what is commonly perceived, the cohesive nature of the flow is not entirely controlled by the interparticle adhesion, but that stiffness and inelasticity of the grains also play a significant role. For the same adhesion, stiffer and less dissipative grains yield a less cohesive flow. This observation is rationalized by introducing the concept of a dynamic, "effective" adhesive force, a single parameter, which combines the effects of adhesion, elasticity, and dissipation. Based on this concept, a rheological description of the flow is proposed for the cohesive grains. Our results elucidate the physics controlling the flow of cohesive granular materials, which may help in designing new approaches to characterize the "flowability" of powders.

Many industrial (wet granulation, food processing, construction, etc.) and geophysical (landslides, mudflow, etc.) processes involve the flows of an assembly of cohesive grains. The cohesion between the grains has different origins. Van der Waals or electrostatic forces are responsible for cohesion in fine grains (1, 2). Liquid capillary (3–5) or solid bridges (6) between the grains give rise to cohesion in large grains. In all cases, the cohesion introduces an additional complexity to granular materials—the flows of cohesive grains are intermittent and less homogeneous (7, 8) in comparison with coarse, cohesionless grains, leading to frequent jamming of industrial units. It is, therefore, necessary to a priori characterize and quantify the capability of flow, so-called “flowability,” of a powder to yield better handling. Different methods are used in industrial contexts for this purpose (1, 9–11). A first method measures the tapped bulk density and the freely settled bulk density of a powder to define the Hausner ratio (or Carr index), which is the ratio of the two. A powder with high Hausner ratio is shown to have poor “flowability.” A second one employs a series of measurements, using the Hosokawa powder tester, comprising angle of repose, aerated bulk density, tapped bulk density, etc., to define a weighted “flowability” index, which ranges from 0 to 100. Very cohesive powders yield “flowability” indices close to zero and the free-flowing ones close to 100. Other methods estimate the macroscopic cohesion from the yield loci of a powder using shear testers (Jenike shear tester or ring shear tester) for various preconsolidation normal stresses, which are useful in understanding the arch formation in silos. All these methods, carried out in the quasistatic limit, are useful for comparing the macroscopic properties of different powders and for characterizing their plastic behavior. However, they do not provide any information about the flow dynamics. Understanding the concept of “flowability” from a physical point of view is still a challenge.

The flow dynamics of rigid, cohesionless grains, interacting solely by contact and friction, is less complex in comparison with cohesive grains, as shown by numerous experimental and numerical studies (12). Flow rules have been evidenced and constitutive laws have been proposed for different flow regimes (13, 14). In the dense flow regime, the rheology of the grains of diameter and density , sheared at a shear rate by imposing a shear stress under a confining normal stress , is well described by a coefficient of friction and a volume fraction , which depend on a single dimensionless parameter, the inertial number (15, 16). These constitutive relations are found to be unaffected by the mechanical properties of the grains—for example, stiffness and inelasticity—as long as the grains are sufficiently rigid and inelastic (14, 17–19). This rough description of the rheology has proven to be useful in describing flows in different configurations from inclines to silos (20–23). This rheological framework has been extended to the flows of cohesive grains (24–27) by using discrete numerical simulations. Different force models with different levels of realism have been used in the simulations to account for the cohesive interactions between the grains (28, 29). In the simplest approach used in rheological studies (24, 27), the adhesion is characterized by a minimum pull-off force necessary to detach two grains. The existence of this additional force scale implies that a second dimensionless number, called cohesion number , exists besides the inertial number, which compares the adhesive force between the grains with the confining normal stress. The rheology is then described by a coefficient of friction and a volume fraction , which are functions of and only (24–27, 30).

In this article, we examine the flow dynamics of cohesive grains down a rough inclined plane using intensive discrete numerical simulations to gain physical insights about the “flowability” of powders. The chosen configuration has inhomogeneous stress distributions and, hence, turns out to be very rich to explore the rheology of cohesive granular materials. We perform a detailed parametric study of the flow to reveal that the flow is significantly affected by the stiffness and the inelasticity of the grains, unlike in the case of cohesionless granular media. We introduce the concept of a dynamic “effective” adhesive force to take into account the effect of the material properties, along with the interparticle adhesion, which is shown to control entirely the flow dynamics. We then define an “effective” cohesion number based on this force, which replaces the cohesion number defined above to form a pair of constitutive relations. The rheology is shown to be well described in this framework.

Results

Simulation of the Flow of Cohesive Grains Down a Rough Inclined Plane.

We investigate the flow of frictional, inelastic, cohesive grains down a rough inclined plane (Fig. 1A) using three-dimensional (3D) discrete element method simulations; an in-house code is used, which is validated by comparing the mean velocity and volume fraction profiles of (cohesionless) monodisperse spheres with those in the study of ref. 17. The grains are spherical and have an average size , with a polydispersity of 20% and an average mass . The rough base (shown in red in Fig. 1) comprises a packed bed of the same grains of height 1.8. The simulation box has length , width , and height . Periodic boundary conditions are applied in the and directions. The flow initiates from the collapse of a cubic array of grains under gravity over the base inclined at a high angle and is continued for a short duration. The flow then attains a steady state after the inclination is lowered and maintained at (the range of angle for a steady flow depends on the interparticle adhesion). The flow is controlled by changing the inclination and four parameters describing the particle properties: the friction coefficient (kept constant at in this study), the stiffness , the quality factor (a measure of dissipation in an interparticle collision), and the adhesion . These parameters are defined below.

Fig. 1.

Simulation of flow of cohesive grains down a rough inclined plane. (A) Snapshot of a simulation, showing the flow of particles (in brown) under gravity over a bumpy base (in red) inclined at an angle . (B) Depiction of the visco-elastic, frictional, and adhesive interactions between the particles using a spring-dashpot model with a frictional slider and an interparticle adhesion . (C) The nonviscous normal contact forces normalized by the interparticle adhesion as a function of the normalized overlap : elastic part (green), adhesive part (blue), and the sum of the two = + (red). See the main text for .

The interparticle contact forces were computed by using the Hookean spring-dashpot model with a frictional slider (Fig. 1) and a Johnson–Kendall–Roberts (JKR)-like (31), yet nonhysteretic, adhesive force model (24). The normal contact force () between particles and comprises three forces: 1) elastic force , where is the normal stiffness and is the normal overlap; 2) viscous force , where is the normal relative velocity, is the effective mass, and is the normal damping coefficient; and 3) adhesive force , which is assumed to be proportional to the area of contact; , the minimum pull-off force, is given by the minimum in the plot of the total nonviscous normal force as a function of the normal overlap (red curve) in Fig. 1. Note that the adhesive interaction model is short-range, meaning that the force vanishes when two grains are not in contact, unlike in wet capillary bridges. The tangential contact force () comprises only the elastic force , where ( in this study) is the tangential stiffness, and is the relative tangential overlap from the beginning of a contact. It is set as , where is the interparticle friction coefficient, to account for the Coulomb yield criterion. Finally, the total contact force acting on particle by particle is = + , where and are normal and tangential unit vectors, respectively. The details of the calculation of viscoelastic forces can be found elsewhere (17). For checking the generality of the results, we also have implemented two other models. One comprises the Hertzian spring-dashpot model and a Derjaguin–Muller–Toporov (DMT)-like adhesive force model (32) (), where the (nonhysteretic) adhesive force is assumed to be independent of the area of contact and is constant. The other one comprises the Hookean spring-dashpot model and a DMT-like, yet hysteretic (), force model, where the grains experience a constant attractive force only during the detachment of a contact, quite similar to a capillary bridge model, but without introducing a finite distance for the detachment.

Using the Hookean-JKR (first) model, the dynamics of two identical () contacting particles in the absence of any external forces is given by the following equation of a nonlinear damped oscillator:In the static limit, the left-hand side and the second term on the right-hand side in Eq. are zero, and the balance between the attractive adhesive force and the repulsive elastic force (Eq. ) then yields an equilibrium overlap . The quality factor of the oscillator is estimated after linearizing the equation around as ; high indicates less collisional dissipation. In the cohesionless case, is related to the restitution coefficient as (17).

All of the equations are made dimensionless by using as the length scale, as the time scale, and as the force scale, and the equations of motion are then solved for each particle by using these interparticle contact forces and gravity. All of the steady-state flow quantities reported below are computed by using the coarse-graining procedure given in .

Not Only the Interparticle Adhesion, but Also the Stiffness and Inelasticity of the Grains Affect the Flow.

Fig. 2 shows a typical velocity profile for the cohesive grains for an intermediate value of adhesion at a given inclination angle . Unlike the flow of dry, cohesionless grains down a rough incline, where a Bagnold velocity profile corresponding to shearing throughout the entire pile is noticed (17, 33), a plug appears near the free surface for the cohesive case (25, 34), where the shear rate () is zero (Fig. 2). The shear rate then gradually increases toward the base (Fig. 2). The inertial number profile () is similar to the shear rate profile (Fig. 2) and is not uniform over the depth of the pile, unlike in the cohesionless case (17, 33). The volume fraction profile () is also not uniform (Fig. 2) (34) with a highdensity region in the plug, unlike in the cohesionless case (17, 33). The cohesive grains possess a finite yield stress, which is reached at a finite depth in the flowing layer. This, therefore, explains the plug formation near the free surface for the flow of cohesive grains.

Fig. 2.

(A–C) Steady flow profiles: mean velocity () (A), volume fraction () (B), and shear rate () and inertial number () (C) at for , , and . The free surface velocity () and the thickness of plug () are denoted in A. (D–F) Effects of different parameters: interparticle adhesion () keeping and fixed (D), particle stiffness () keeping and fixed (E), and quality factor () keeping and fixed (F), on the velocity profile at . The dashed line in D shows a Bagnold fit.

We now examine the effect of the interparticle adhesion on the velocity profile (Fig. 2) at a given inclination angle . When goes to zero, one recovers the classical Bagnold profile (shown by the fitted dashed line). The Bagnold profile is also observed (fits are not shown for clarity) for small, nonzero values of considered, , 25, and 50; however, the free-surface velocity decreases in comparison with the cohesionless case with increasing . When becomes greater than 50, a plug appears near the free surface, which grows in size with increasing until the flow comes to a halt for . Such behavior of the velocity profile with increasing has already been reported in some studies (25, 34). More surprising results appear while changing the stiffness of the particles and the dissipation through the quality factor , as shown in Fig. 2 . The free-surface velocity increases, and the thickness of the plug decreases with increasing the stiffness, keeping other parameters fixed (Fig. 2). The plug completely disappears for a sufficiently high value of stiffness. A similar observation is made while increasing the quality factor, i.e., decreasing the dissipative nature of contact (Fig. 2). These observations clearly indicate that the bulk cohesion is not solely controlled by the interparticle adhesion but is dependent on the material properties and . The sensitivity of the flow to the stiffness and dissipation in the cohesive case contrasts with the flow in the cohesionless case (14, 17–19). This observation has a crucial consequence—the characterization of the degree of bulk cohesion based on the interparticle adhesion using the granular Bond number () (35, 36) or the cohesion number (24, 25, 28) may not be sufficient. We next attempt to collate the effects of these three independent parameters into a single parameter, the dynamic “effective” adhesion .

A Scaling for the Dynamic “Effective” Adhesive Force.

We infer from Fig. 2 that decreasing the stiffness or decreasing the quality factor are equivalent to increasing the adhesion. Hence, we seek for an expression of the dynamic “effective” adhesive force aswhere and are unknown constants. The idea behind this scaling is to check if the influence of adhesion, stiffness, and dissipation can be embedded in a single parameter , acting as an “effective” adhesive force in dynamic conditions. We do the following to test the idea. We systematically carry out simulations by varying , , and (at a fixed inclination ) and then extract the free-surface velocity and the thickness of the plug from the steady velocity profile in each case. We then plot and separately for all of the simulations as a function of computed from Eq. and look for the best collapse of the data in each case by trying out different combinations of and . For each chosen combination of and , the two master curves ( and ) resulting from the data collapses are fitted by power laws (only the data corresponding to velocity profiles with a plug flow are considered), and the quality of the collapses is determined by the determination coefficients and of the fits. The final values of the constants and are chosen such that they lead to the highest (). Finally, two well-defined master curves, and , emerge (Fig. 3), which implies that the flow is indeed controlled by an “effective” adhesive force, which comprises three independent parameters: interparticle adhesion, stiffness, and dissipation. We have checked that two different sets of parameters, corresponding to the same value of , not only give the same free-surface velocity and plug thickness, but also yield identical velocity, density, and r.m.s. velocity profiles (). We notice two distinct regions in the figure, separated by a vertical dashed line at : 1) a plug-less region on the left, where decreases monotonically with increasing and yet ; and 2) a plug-full region where decreases and increases monotonically with increasing . In order to test the generality of the above result further, we have carried out additional simulations at a given angle () in the same spirit using the Hertzian-DMT model. The sensitivity of the flow to the stiffness and dissipation is recovered, and the velocity profile is again found () to be controlled by the dynamic “effective” adhesive force given by a similar equation (). The origin of the scaling will be discussed in the last section. We also perform some simulations using the hysteretic contact model. We again find the flow dynamics (at a given angle) to be dependent on the stiffness (), signifying that in this case as well, the contact parameters, along with the interparticle adhesion, determine the bulk cohesion. We examine below if this dynamic “effective” adhesive force is relevant in the description of the rheology.

Fig. 3.

The flow is controlled by the dynamic “effective” adhesive force. Variation of the free surface velocity () (A) and the thickness of the plug () (B) with the dynamic “effective” adhesive force at for different (different symbols; see Fig. 2 for the key), different (different colors), and different (different symbol pattern). The dashed line indicates the transition from a plug-less to a plug-full velocity profile.

Flow Cessation Is Controlled by the “Effective” Adhesive Force, but Flow Initiation Is Controlled by the “Actual” Adhesive Force.

A first step toward exploring the rheology is to study the yield criteria of our model cohesive material, i.e., to study the stress conditions under which the flow stops (dynamic yielding) or starts (static yielding). The yield criteria are usually described by using a cohesive Mohr–Coulomb model stipulating that, on the plane of incipient failure, the shear stress is related to the normal stress by , where is the static (respectively [resp.] dynamic) friction coefficient, and is the static (resp. dynamic) cohesive stress. We first focus on the dynamic yielding by considering the data of the velocity profiles in the previous section. The dynamic yield criterion is satisfied at the intersection of the plug and sheared regions, denoted by “yield point” in Fig. 2. The location of the “yield point” shifts with changing the inclination angle. Thus, a yield locus for a given set of particle properties (, , and ) is obtained by noting a series of points ( and ), each corresponding to a “yield point” at a given angle. Fig. 4 , Inset shows different yield loci (shown by symbols; only a few are shown for clarity) for different sets of particle properties. All of the yield loci are well approximated by straight lines, which are the best fits of the Mohr–Coulomb model. The slope of each straight line gives the dynamic friction coefficient and the intercept with the y axis, the dynamic cohesive stress . We have systematically extracted and for all of the cases mentioned in Fig. 3. The dynamic friction coefficient is found to be constant and to be independent of the particle properties. On the other hand, the dynamic cohesive stress increases with increasing and decreasing or . However, they all collapse on a single master curve if plotted against the “effective” adhesion (Fig. 4), showing again that is the only variable controlling the dynamic yielding. A linear fit through the data givesWe next study the static yielding from the initiation of the flow. We initially bring a steady flow of a given pile thickness to a halt by gradually lowering the inclination angle. We then gradually increase the inclination by a step of 0. (a steady state in the kinetic energy profile is ensured after each increment) until the flow starts (indicated by a threshold of kinetic energy ). We finally note the stresses ( and ) at the bottom of the pile, where the yielding happens. The procedure is repeated by varying the initial pile thickness to get a series of points ( and ), which constitutes the yield locus for a given set of particle properties (, , ). Fig. 4 , Inset shows different yield loci (shown by symbols) for different sets of particle properties. The striking observation is that two yield loci for two different for a given and collapse well on each other. A similar observation is made for two different for a given and (only one point is shown). This result implies that the flow initiation is entirely controlled by the interparticle adhesion, which contrasts with the observation made above that the flow cessation is controlled by the dynamic “effective” adhesion. A careful inspection reveals that the static yield loci are not perfectly linear, especially for high values of , which have tails comprising data points at low values of the normal stress. However, at sufficiently high values of the normal stress (corresponding to thick piles), one recovers a linear variation of the shear stress with the normal stress, as expected from the Mohr–Coulomb model. Hence, we restrict the fitting of the model (shown by the straight lines) to the linear parts of the yield loci; for a given , the combined datasets corresponding to different and are considered during fitting. The fitting yields static friction coefficients (estimated by the slopes of the straight lines): for and for the others. Note that is slightly higher than . We also obtain the static cohesive stress from the intercept of the straight line in each case, which is shown as a function of the interparticle adhesion in Fig. 4. The variation of with is well approximated byDifferent studies (4, 37, 38) have predicted theoretically such a linear relation between the static cohesive stress and the adhesion aswhere is the average coordination number (number of contacts per particle). The contact network was assumed to be isotropic in the derivation. Using measured and and extracted , Eq. yields , with a proportionality constant of 0.76, which is a decade higher in comparison with that (0.07) extracted from the fitting. Anisotropy in the contact network in the present study may lead to this mismatch. A similar mismatch has been reported by Rognon et al. (25) for flows down inclines. The last noteworthy point, which can be inferred from Fig. 4, is that the maximum dynamic cohesive stress observed is smaller than the static cohesive stress for a given adhesion.

Fig. 4.

Dynamic and static yield criteria. (A, Inset) Yield shear stress () as a function of yield normal stress () for the dynamic yielding in different cases. The keys are the same as in Figs. 2 and 3. The lines are the fits of the Mohr–Coulomb failure equation. In the main plot, the variation of the dynamic macroscopic cohesive stress () with the “effective” adhesive force for all of the cases. The straight line is the best fit to the data. (B, Inset) Yield shear stress () as a function of yield normal stress () for the static yielding in different cases. The dotted lines are the guides to the eye, and the solid lines are the fits of the Mohr–Coulomb failure equation in the linear part of the data. The main plot shows variation of the static macroscopic cohesive stress () with the actual adhesive force . The dotted line is a guide to the eye, and the straight line is the best fit to the data.

Bulk Rheology Is Described by an “Effective” Cohesion Number.

The last step toward exploring the rheology is to go beyond the yield criteria and analyze how the shear stress varies with the shear rate. The inclined plane configuration serves as a rheometer to enable us to measure the local shear rate , shear stress , normal stress , and volume fraction at every vertical position in the flow. Rognon and coworkers (24–26) proposed the following -rheology framework: and , for describing the rheology of cohesive grains, as mentioned in the introduction. However, the results obtained in the previous sections suggest that the above framework is not sufficient for the complete description of the rheology. The observation that the “effective” adhesion , not the interparticle adhesion , controls the flow dynamics, prompts to define a new cohesion number, the “effective” cohesion number . One can then test if the rheology can be described by using the new framework: and . We do the following to test the idea. We compute local , , and for all of the simulations done for different (, , and ) and at different angles . Recent studies (39–41) have emphasized the occurrence of nonlocal effects in inhomogeneous flows, such as the one investigated here. For each simulation, we try to minimize this effect on the analysis by neglecting the data points (for which ) adjacent to the “yield point,” where the effect is significant. We also discard six data points each adjacent to the free surface and the rough base to avoid boundary effects. We consider the rest and combine them for all of the simulations. The combined dataset (, , , and ) is then divided into eight subsets based on the value of (a bin size of 0.02 is chosen). For each subset, , , and are then averaged in small bins of size 0.01 in . Fig. 5 shows the variation of and with for the eight different values of . The first important observation is that the datasets are reasonably well sorted according to the value of , which indicates that the rheology is entirely controlled by two dimensionless numbers, and , validating our model. The curve (Fig. 5) shifts upward with increasing , showing that increases with increasing , particularly for small values of . Substituting given in Eq. into the Mohr–Coulomb model and using the definition of and , we show that increases linearly with increasing in the quasistatic limit (), following . However, the curves surprisingly seem to merge in the high- regime, indicating less impact of adhesion on the value of . appears to be independent of for sufficiently high values of . The curve (Fig. 5) shifts downward with increasing , showing that the volume fraction decreases with increasing . A nonmonotonic behavior of appears for intermediate values of —the volume fraction first decreases and then increases with increasing inertial number before decreasing again at high values. also appears to be independent of for sufficiently high values of . Berger et al. (27) observed similar behavior of the effective friction and volume fraction with the inertial and cohesion numbers and proposed an empirical function to capture the rheology (). The same function provides a reasonable fit of our data, as shown in .

Fig. 5.

Constitutive laws for the flow of cohesive granular media. Variation of the effective friction () (A) and the volume fraction () (B) with the inertial number () for different values of the “effective” cohesion number (). The assembled data are obtained by varying , , , and .

Discussion and Conclusion

In this work, we have examined the flow of cohesive grains down an inclined plane using discrete element method simulations. We use a simplified adhesive interaction law, characterized by a minimum pull-off force necessary to detach two grains, to model the interparticle adhesion. Despite this simple interaction model, a rich flow dynamics appears—the cohesive nature of the flow becomes sensitive to the stiffness and inelasticity of the grains besides the interparticle adhesion , which contrasts with the flow of cohesionless granular media. Although this sensitivity has never been reported for the flows down inclines, a few studies on the fluidization of cohesive powders observed it (42–44). This finding questions the use of granular Bond number or cohesion number (based on the interparticle adhesion) in quantifying the degree of bulk cohesion. We have demonstrated a way to take into account the effects of these parameters, along with the interparticle adhesion into a single parameter called the dynamic “effective” adhesion. This force appears to control the flow dynamics solely and comes out to be for the Hookean-JKR model and for the Hertzian-DMT model; the success of obtaining a single control variable eases the problem from rheological perspectives. One can note from the above equations that increasing and decreasing or increases the “effective” adhesion. Interestingly, the scalings obtained in both the cases can be recast by using the equilibrium overlap in the following form: , where for the Hookean-JKR model and for the Hertzian-DMT model.

We have not succeeded in understanding the scaling. However, the existence of an interplay between the interparticle adhesion and the mechanical properties of the grains can be evidenced, considering the dynamics of a binary collision. When two cohesive grains collide, one can show that they remain glued together if the relative kinetic energy before impact is less than a critical value () given by , where is a decreasing function of the quality factor, which depends on the interaction model. This shows that not only the adhesion, but also the mechanical properties of the grains influence a collision: The chance of a rebound is less if is large and also if () or is small. To determine the “effective” adhesive force , one may estimate the typical kinetic energy involved in a collision in a flow under a shear stress (45) and compare it to the critical value —a work which we aim to do in the near future.

We have shown that the initiation of the flow is controlled by a static cohesive stress proportional to , whereas the arrest of the flow, by a dynamic cohesive stress proportional to . The maximum dynamic cohesive stress observed is smaller than the static cohesive stress for a given interparticle adhesion. This implies that the characterization and quantification of the “flowability” made at the quasistatic limit may not be appropriate for the dynamic limit. Although we have not been able to understand the scaling of dynamic “effective” adhesive force, we have shown that the “effective” cohesion number based on this empirical force and the inertial number yield two constitutive relations: and , which give a reasonably good description of the rheology. The determination of the exact forms of these two relations is beyond the scope of this work and requires a separate study, possibly using a plane shear flow. In , we present some preliminary rheological data obtained in a normal stress-imposed shear cell (see for simulation details). The data of and for two different sets ( and ) for a given , yielding the same , collapse well on each other and also match reasonably the data from the inclined plane flow (), showing the generality of the proposed description.

The relevance of the “effective” cohesion number in defining the flow dynamics opens perspectives to analyze the behavior of cohesive granular media in other configurations, and we give one example below. We carry out additional simulations to measure the packing fraction in a pile of cohesive grains for various (, , and ). The grains are poured randomly under gravity in the same periodic box () as used above, and the measurements are taken once all of the grains settle down. Since the packing fraction and the “effective” cohesion number vary along the height of the pile, we obtain a series of points ( and ) from each simulation. Fig. 6 shows the variation of with for all of the simulations, neglecting the data points adjacent to the rough base and the free surface. We obtain a good collapse of the data, resulting in a master curve. However, this is not valid if the cohesion number based on the actual adhesion is considered instead (Fig. 6, Inset). More configurations need to be studied to be able to understand to which extent this concept of “effective” cohesion number is valid for the flow of cohesive grains, which might help in developing new approaches for the characterization of powders.

Fig. 6.

A validity check for the “effective” cohesion number. Variation of the volume fraction () with the “effective” cohesion number () in the case of gravity-driven compaction is shown. See Figs. 2 and 3 for additional legends. (Inset) Variation of the volume fraction () with the cohesion number ().

Although the concept of “effective” adhesion is based on a simplified model of adhesion, we find that it is generic for other kinds of adhesive interactions; for example, (hysteretic) capillary and electrostatic adhesion. In these cases, the particles experience an attractive force, even without a physical contact over a small separation distance, and this distance plays an equivalent role of the stiffness and influence the “effective” adhesion. The bulk cohesion, in studies (26, 46) using capillary bridge models, was noted to decrease when decreasing this separation distance (named as the “rupture distance,” beyond which a capillary bridge breaks). This outcome can be understood considering the argument of energy proposed above—the work needed to separate two bonded grains, in this case, is proportional to the rupture distance.

One last remark concerns the limit of rigid particles. The scaling reported above suggests that the “effective” adhesion will go to zero in the rigid limit, implying that a direct comparison of our results with those using contact dynamics simulations (27), in which the grains are treated as perfectly rigid, is difficult. We speculate that, in these simulations, the size of the time step could be crucial and could play a role similar to the stiffness in our soft-particle approach.

The main idea conveyed in this article that the “effective” adhesion is not controlled solely by the interparticle adhesion, but also dependent on material properties is a first step toward a better understanding of the flow of cohesive granular media, which may benefit engineering and geophysical communities to understand the long-standing issue of “flowability” of cohesive powders.

All data presented in this article are openly available in the Zenodo repository, https://zenodo.org/record/3699632#.XmJ048tKgaw.