The clustering of debris floating on liquid interfaces such as water surfaces is a complex phenomenon that finds its applications in numerous examples from industrial processing and environmental systems. The recent paper by Shin & Coletti (J. Fluid Mech., vol. 984, 2024, R7) presents an experimental campaign investigating the three effects of turbulence, particle interactions and interfacial effects, to elucidate how the three force scales drag, capillary forces and lubrication give rise to three distinct regimes of clustering in dense suspensions. The study, hence, provides a useful systematic to categorize the clustering mechanisms. As an important finding, it is shown that, depending on volume fraction and non-dimensional turbulent shear, particles either tend to cluster into aggregate sizes larger than the Kolmogorov scale or can break into pieces that are as small as the primary particle size.

Cohesive particulate flows play an important role in environmental fluid dynamics, as well as in a wide variety of civil and process engineering applications. However, the scaling laws, constitutive equations and continuum field descriptions governing such flows are currently less well understood than for their non-cohesive counterparts. Grain-resolved simulations represent an essential tool for addressing this shortcoming, along with theoretical investigations, laboratory experiments and field studies. Here we provide a tutorial introduction to simulations of fine-grained sediments in viscous fluids, along with an overview of some representative insights that this approach has yielded to date. After a brief review of the key physical concepts governing van der Waals forces as the main cohesive effect for subaqueous sediment suspensions, we discuss their incorporation into particle-resolved simulations based on the immersed boundary method. We subsequently describe simulations of cohesive particles in several model turbulent flows, which demonstrate the emergence of a statistical equilibrium between the growth and break-up of aggregates. As a next step, we review the influence of cohesive forces on the settling behaviour of dense suspensions, before moving on to submerged granular collapses. Throughout the article, we highlight open research questions in the field of cohesive particulate flows whose investigation may benefit from grain-resolved simulations.

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.

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.

We investigate the dynamics of cohesive particles in homogeneous isotropic turbulence, based on one-way coupled simulations that include Stokes drag, lubrication, cohesive and direct contact forces. We observe a transient flocculation phase, followed by a statistically steady equilibrium phase. We analyse the temporal evolution of floc size and shape due to aggregation, breakage and deformation. Larger turbulent shear and weaker cohesive forces yield smaller elongated flocs. Flocculation proceeds most rapidly when the fluid and particle time scales are balanced and a suitably defined Stokes number is $O(1)$. During the transient stage, cohesive forces of intermediate strength produce flocs of the largest size, as they are strong enough to cause aggregation, but not so strong as to pull the floc into a compact shape. Small Stokes numbers and weak turbulence delay the onset of the equilibrium stage. During equilibrium, stronger cohesive forces yield flocs of larger size. The equilibrium floc size distribution exhibits a preferred size that depends on the cohesive number. We observe that flocs are generally elongated by turbulent stresses before breakage. Flocs of size close to the Kolmogorov length scale preferentially align themselves with the intermediate strain direction and the vorticity vector. Flocs of smaller size tend to align themselves with the extensional strain direction. More generally, flocs are aligned with the strongest Lagrangian stretching direction. The Kolmogorov scale is seen to limit floc growth. We propose a new flocculation model with a variable fractal dimension that predicts the temporal evolution of the floc size and shape.