skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Grain-resolving simulations of submerged cohesive granular collapse
We investigate the submerged collapse of weakly polydisperse, loosely packed cohesive granular columns, as a function of aspect ratio and cohesive force strength, via grain-resolving direct numerical simulations. The cohesive forces act to prevent the detachment of individual particles from the main body of the collapsing column, reduce its front velocity, and yield a shorter and thicker final deposit. All of these effects can be captured accurately across a broad range of parameters by piecewise power-law relationships. The cohesive forces reduce significantly the amount of available potential energy released by the particles. For shallow columns, the particle and fluid kinetic energy decreases for stronger cohesion. For tall columns, on the other hand, moderate cohesive forces increase the maximum particle kinetic energy, since they accelerate the initial free-fall of the upper column section. Only for larger cohesive forces does the peak kinetic energy of the particles decrease. Computational particle tracking indicates that the cohesive forces reduce the mixing of particles within the collapsing column, and it identifies the regions of origin of those particles that travel the farthest. The simulations demonstrate that cohesion promotes aggregation and the formation of aggregates. Furthermore, they provide complete information on the temporally and spatially evolving network of cohesive and direct contact force bonds. While the normal contact forces are aligned primarily in the vertical direction, the cohesive bonds adjust their preferred spatial orientation throughout the collapse. They result in a net macroscopic stress that counteracts deformation and slows the spreading of the advancing particle front.  more » « less
Award ID(s):
2138583
PAR ID:
10489376
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Cambridge University Press
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
942
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Much confusion exists on whether force‐ or energy‐based descriptions of cohesive‐particle interactions are more appropriate. We hypothesize a force‐based description is appropriate when enduring‐contacts dominate and an energy‐based description when contacts are brief in nature. Specifically, momentum is transferred through force‐chains when enduring‐contacts dominate and particles need to overcome a cohesive force to induce relative motion, whereas particles experiencing brief contacts transfer momentum through collisions and must overcome cohesion‐enhanced energy losses to avoid agglomeration. This hypothesis is tested via an attempt to collapse the dimensionless, dependent variable characterizing a given system against two dimensionless numbers: A generalized bond number, BoG–ratio of maximum cohesive force to the force driving flow, and a new Agglomerate number, Ag–ratio of critical cohesive energy to the granular energy. A gamut of experimental and simulation systems (fluidized bed, hopper, etc.), and cohesion sources (van der Waals, humidity, etc.), are considered. For enduring‐contact systems, collapse occurs with BoGbut not Ag, and vice versa for brief‐contact systems, thereby providing support for the hypothesis. An apparent discrepancy with past work is resolved, and new insight into Geldart's classification is gleaned. 
    more » « less
  2. We develop a physical and computational model for performing fully coupled, grain-resolved direct numerical simulations of cohesive sediment, based on the immersed boundary method. The model distributes the cohesive forces over a thin shell surrounding each particle, thereby allowing for the spatial and temporal resolution of the cohesive forces during particle–particle interactions. The influence of the cohesive forces is captured by a single dimensionless parameter in the form of a cohesion number, which represents the ratio of cohesive and gravitational forces acting on a particle. We test and validate the cohesive force model for binary particle interactions in the drafting–kissing–tumbling (DKT) configuration. Cohesive sediment grains can remain attached to each other during the tumbling phase following the initial collision, thereby giving rise to the formation of flocs. The DKT simulations demonstrate that cohesive particle pairs settle in a preferred orientation, with particles of very different sizes preferentially aligning themselves in the vertical direction, so that the smaller particle is drafted in the wake of the larger one. This preferred orientation of cohesive particle pairs is found to remain influential for systems of higher complexity. To this end, we perform large simulations of 1261 polydisperse settling particles starting from rest. These simulations reproduce several earlier experimental observations by other authors, such as the accelerated settling of sand and silt particles due to particle bonding, the stratification of cohesive sediment deposits, and the consolidation process of the deposit. They identify three characteristic phases of the polydisperse settling process, viz. (i) initial stir-up phase with limited flocculation, (ii) enhanced settling phase characterized by increased flocculation, and (iii) consolidation phase. The simulations demonstrate that cohesive forces accelerate the overall settling process primarily because smaller grains attach to larger ones and settle in their wakes. For the present cohesive number values, we observe that settling can be accelerated by up to 29 %. We propose physically based parametrization of classical hindered settling functions introduced by earlier authors, in order to account for cohesive forces. An investigation of the energy budget shows that, even though the work of the collision forces is much smaller than that of the hydrodynamic drag forces, it can substantially modify the relevant energy conversion processes. 
    more » « less
  3. We interpret the Taylor–Green cellular vortex model in terms of the Kolmogorov length and velocity scales, in order to study the balance between aggregation and breakup of cohesive sediment in fine-scale turbulence. One-way coupled numerical simulations, which capture the effects of cohesive, lubrication and direct contact forces on the flocculation process, reproduce the non-monotonic relationship between the equilibrium floc size and shear rate observed in previous experiments. The one-way coupled results are confirmed by select two-way coupled simulations. Intermediate shear gives rise to the largest flocs, as it promotes preferential concentration of the primary particles without generating sufficiently strong turbulent stresses to break up the emerging aggregates. We find that the optimal intermediate shear rate increases for stronger cohesion and smaller particle-to-fluid density ratios, and we propose a simple model for the equilibrium floc size that agrees well with experimental data reported in the literature. 
    more » « less
  4. Acoustic levitation in air provides a containerless, gravity-free platform for investigating driven many-particle systems with nonconservative interactions and underdamped dynamics. In prior work the interactions among levitated particles were limited to attractive forces from scattered sound and repulsion from hydrodynamic microstreaming. We report on experiments in which contact cohesion provides a third type of interaction. When particle size and separation are both much smaller than the sound wavelength, this interplay of three interactions results in forces that are attractive over several particle diameters, become repulsive at close approach, and are again attractive at contact. In the presence of sound-induced athermal fluctuations that generate particle collisions, the interplay of these three forces enables the formation of particle chains with anisotropic interactions that depend on chain size and shape due to multibody effects. With the control of the kinetic pathways and the strength of the contact cohesion, different patterns can be assembled, from triangular lattices to labyrinthine patterns of chains to lacelike networks of interconnected rings. These results shed light on the multibody character of acoustic interactions and can be utilized to direct the self-assembly of particles. Published by the American Physical Society2025 
    more » « less
  5. When a colloidal suspension is dried, capillary pressure may overwhelm repulsive electrostatic forces, assembling aggregates that are out of thermal equilibrium. This poorly understood process confers cohesive strength to many geological and industrial materials. Here we observe evaporation-driven aggregation of natural and synthesized particulates, probe their stability under rewetting, and measure bonding strength using an atomic force microscope. Cohesion arises at a common length scale (∼5 μm), where interparticle attractive forces exceed particle weight. In polydisperse mixtures, smaller particles condense within shrinking capillary bridges to build stabilizing “solid bridges” among larger grains. This dynamic repeats across scales, forming remarkably strong, hierarchical clusters, whose cohesion derives from grain size rather than mineralogy. These results may help toward understanding the strength and erodibility of natural soils, and other polydisperse particulates that experience transient hydrodynamic forces. 
    more » « less