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: An efficient cellular flow model for cohesive particle flocculation in turbulence
We propose a one-way coupled model that tracks individual primary particles in a conceptually simple cellular flow set-up to predict flocculation in turbulence. This computationally efficient model accounts for Stokes drag, lubrication, cohesive and direct contact forces on the primary spherical particles, and allows for a systematic simulation campaign that yields the transient mean floc size as a function of the governing dimensionless parameters. The simulations reproduce the growth of the cohesive flocs with time, and the emergence of a log-normal equilibrium distribution governed by the balance of aggregation and breakage. Flocculation proceeds most rapidly when the Stokes number of the primary particles is $O(1)$ . Results from this simple computational model are consistent with experimental observations, thus allowing us to propose a new analytical flocculation model that yields improved agreement with experimental data, especially during the transient stages.  more » « less
Award ID(s):
1803380 1924532
PAR ID:
10143098
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
889
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    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. 
    more » « less
  2. 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
  3. 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
  4. Abstract Biophysical cohesion, introduced predominantly by Extracellular Polymeric Substances (EPS) during mineral flocculation in subaqueous environments, plays important role in morphodynamics, biogeochemical cycles and ecosystem processes. However, the mechanism of how EPS functioning with cohesive particles and affects settling behaviors remain poorly understood. We measure initial flocculation rate, floc size and settling velocity of mineral and artificial EPS (Xanthan gum) mixtures. Combining results from these and previous studies demonstrate coherent intensification of EPS-related flocculation compare with those of pure mineral and oil-mineral mixtures. Importantly, the presence of EPS fundamentally changes floc structure and reduces variability of settling velocity. Measured data shows that ratios of microfloc and macrofloc settling velocity for pure mineral flocs is 3.9 but greatly reduced to a lowest value of 1.6 due to biological EPS addition. The low variability of settling velocity due to EPS participation explains the seemingly inconsistent results previously observed between field and laboratory studies. 
    more » « less
  5. The problem of orthogonal penetration of a blade into the core of a vortex with non-zero axial flow was studied using a combination of scaling theory, a heuristic plug-flow model and full Navier–Stokes simulations. The particular focus of this paper was to understand the mechanics of the transient lift force that occurs during the initial penetration of the blade leading edge into the vortex core, and the relationship of this transient force to the steady-state lift force that develops due to the difference in vortex core radius over the blade surface. The three modelling approaches all lead to the conclusion that the maximum value of the lift coefficient for the transient blade penetration force is proportional to the impact parameter and inversely proportional to the axial flow parameter. This observation is used to develop a simple expression that collapses the predictions of the full Navier–Stokes simulations for lift coefficient over a wide range of parameter values. 
    more » « less