Search

Abstract:

Abstract:

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.

We present a systematic simulation campaign to investigate the pairwise interaction of two mobile, monodisperse particles submerged in a viscous fluid and subjected to monochromatic oscillating flows. To this end, we employ the immersed boundary method to geometrically resolve the flow around the two particles in a non-inertial reference frame. We neglect gravity to focus on fluid–particle interactions associated with particle inertia and consider particles of three different density ratios aligned along the axis of oscillation. We systematically vary the initial particle distance and the frequency based on which the particles show either attractive or repulsive behaviour by approaching or moving away from each other, respectively. This behaviour is consistently confirmed for the three density ratios investigated, although particle inertia dictates the overall magnitude of the particle dynamics. Based on this, threshold conditions for the transition from attraction to repulsion are introduced that obey the same power law for all density ratios investigated. We furthermore analyse the flow patterns by suitable averaging and decomposition of the flow fields and find competing effects of the vorticity induced by the fluid–particle interactions. Based on these flow patterns, we derive a circulation-based criterion that provides a quantitative measure to categorize the different cases. It is shown that such a criterion provides a consistent measure to distinguish the attractive and repulsive arrangements. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Many aquatic environments contain cohesive sediments that flocculate and create flocs with a wide range of sizes. The Population Balance Equation (PBE) flocculation model is designed to predict the time-dependent floc size distribution and should be more complete than models based on median floc size. However, a PBE flocculation model includes many empirical parameters to represent important physical, chemical, and biological processes. We report a systematic investigation of key model parameters of the open-source PBE-based size class flocculation model FLOCMOD (Verney, Lafite, Claude Brun-Cottan and Le Hir, 2011) using the measured temporal floc size statistics reported by Keyvani and Strom (2014) at a constant turbulent shear rate S. Results show that the median floc size d 50 , in terms of both the equilibrium floc size and the initial floc growth, is insufficient to constrain the model parameters. A comprehensive error analysis shows that the model is capable of predicting three floc size statistics d 16 , d 50 and d 84 , which also reveals a clear trend that the best calibrated fragmentation rate (inverse of floc yield strength) is proportional to the floc size statistics considered. Motivated by this finding, the importance of floc yield strength is demonstrated in the predicted temporal evolution of floc size by modeling the floc yield strength as microflocs and macroflocs giving two corresponding fragmentation rates. The model shows a significantly improved agreement in matching the measured floc size statistics., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper, (Copyright © 2023 Elsevier Ltd. All rights reserved.)