Sand settling through bedform-generated turbulence in rivers

Fluvial bedforms generate a turbulent wake that can impact suspended-sediment settling in the passing flow. This impact has implications for local suspended-sediment transport, bedform stability, and channel evolution; however, it is typically not well-considered in geomorphologic models. Our study uses a three-dimensional OpenFOAM hydrodynamic and particle-tracking model to investigate how turbulence generated from bedforms and the channel bed influences medium sand-sized particle settling, in terms of the distribution of suspended particles within the flow field and particle-settling velocities. The model resolved the effect of an engineered bedform, which altered the flow field in a manner similar to a natural dune. The modelling scenarios alternated bed morphology and the simulation of turbulence, using detached eddy simulation (DES), to differentiate the influence of bedform-generated turbulence relative to that of turbulence generated from the channel bed. The bedform generated a turbulent wake that was composed of eddies with significant anisotropic properties. The eddies and, to a lesser degree, turbulence arising from velocity shear at the bed substantially reduced settling velocities relative to the settling velocities predicted in the absence of turbulence. The eddies tended to advect sediment particles in their primary direction, diffuse particles throughout the flow column, and reduced settling likely due to production of a positively skewed vertical-velocity fluctuation distribution. Study results suggest that the bedform wake has a significant impact on particle-settling behaviour (up to a 50% reduction in settling velocity) at a scale capable of modulating local suspended transport rates and bedform dynamics. © 2020 John Wiley & Sons, Ltd.     

Significance of particle interaction to the modelling of cohesive sediment transport in rivers

Aggregation processes of fine sediments have rarely been integrated in numerical simulations of cohesive sediment transport in riverine systems. These processes, however, can significantly alter the hydrodynamic characteristics of suspended particulate matter (SPM), modifying the particle settling velocity, which is one of the most important parameters in modelling suspended sediment dynamics. The present paper presents data from field measurements and an approach to integrate particle aggregation in a hydrodynamic sediment transport model. The aggregation term used represents the interaction of multiple sediment classes (fractions) with corresponding multiple deposition behaviour. The k–ε–turbulence model was used to calculate the coefficient of vertical turbulent mixing needed for the two‐dimensional vertical‐plane simulations. The model has been applied to transport and deposition of tracer particles and natural SPM in a lake‐outlet lowland river (Spree River, Germany). The results of simulations were evaluated by comparison with field data obtained for two levels of river discharge. Experimental data for both discharge levels showed that under the prevailing uniform hydraulic conditions along the river reach, the settling velocity distribution did not change significantly downstream, whereas the amount of SPM declined. It was also shown that higher flow velocities (higher fluid shear) resulted in higher proportions of fast settling SPM fractions. We conclude that in accordance with the respective prevailing turbulence structures, typical aggregation mechanisms occur that continuously generate similar distribution patterns, including particles that settle toward the river bed and thus mainly contribute to the observed decline in the total SPM concentration. In order to determine time‐scales of aggregation and related mass fluxes between the settling velocity fractions, results of model simulations were fitted to experimental data for total SPM concentration and of settling velocity frequency distributions. The comparison with simulations for the case of non‐interacting fractions clearly demonstrated the practical significance of particle interaction for a more realistic modelling of cohesive sediment and contaminant transport. Copyright © 2004 John Wiley & Sons, Ltd.     

The influence of asymmetries in overlying stratification on near-bed turbulence and sediment suspension in a partially mixed estuary

Data collected from the York River estuary demonstrate the importance of asymmetries in stratification to the suspension and transport of fine sediment. Observations collected during two 24-h deployments reveal greater concentrations of total suspended solids during the flood phase of the tide despite nearly symmetric near-bed tidal current magnitude. In both cases, tidally averaged net up-estuary sediment transport near the bed was clearly observed despite the fact that tidally averaged residual near-bed currents were near zero. Tidal straining of the along-channel salinity gradient resulted in a stronger pycnocline lower in the water column during the ebb phase of the tide and appeared to limit sediment suspension. Indirect measurements suggest that the lower, more intense, pycnocline on the ebb acted as a barrier, limiting turbulent length scales and reducing eddy diffusivity well below the pycnocline, even though the lower water column was locally well mixed. In order to more conclusively link changes in stratification to properties of near-bed eddy viscosity and diffusivity, longer duration tripod and mooring data from an additional experiment are examined, that included direct measurement of turbulent velocities. These additional data demonstrate how slight increases in stratification can limit vertical mixing near the bed and impact the structure of the eddy viscosity below the pycnocline. We present evidence that the overlying pycnocline can remotely constrain the vertical turbulent length scale of the underlying flow, limiting sediment resuspension. As a result, the relatively small changes in stratification caused by tidal straining of the pycnocline allow sediment to be resuspended higher in the water column during the flood phase of the tide, resulting in preferential up-estuary transport of sediment.Responsible Editor: Iris Grabemann     

Parameterization of near-bed processes under collinear wave and current flows from a two-phase sheet flow model

Sediment transport models require appropriate representation of near-bed processes. We aim here to explore the parameterizations of bed shear stress, bed load transport rate and near-bed sediment erosion rate under the sheet flow regime. To that end, we employ a one-dimensional two-phase sheet flow model which is able to resolve the intrawave boundary layer and sediment dynamics at a length scale on the order of the sediment grain. We have conducted 79 numerical simulations to cover a range of collinear wave and current conditions and sediment diameters in the range 210–460 μm$\mathrm{\mu }\mathrm{m}$. The numerical results confirm that the intrawave bed shear stress leads the free stream velocity, and we assess an explicit expression relating the phase lead to the maximum velocity, wave period and bed roughness. The numerical sheet flow model is also used to provide estimates for the bed load transport rate and to inspect the near-bed sediment erosion. A common bed load transport rate formulation and two typical reference concentration approaches are assessed. A dependence of the bed load transport rate on the sediment grain diameter is observed and parameterized. Finally, the intrawave near-bed vertical sediment flux is further investigated and related to the time derivative of the bed shear stress.     

A two-dimensional model for suspended sediment transport in the southern branch of the Rhine–Meuse estuary,The Netherlands

For the southern branch of the Rhine–Meuse estuary, The Netherlands, a two-dimensional horizontal suspended sediment transport model was constructed in order to evaluate the complicated water quality management of the area. The data needed to calibrate the model were collected during a special field survey at high river runoff utilizing a number of techniques: (1) turbidity probes were used to obtain suspended sediment concentration profiles; (2) air-borne remote sensing video recordings were applied in order to obtain information concerning the spatial distribution of the suspended sediment concentration; (3) an acoustic probe (ISAC) was used to measure cohesive bed density profiles and (4) an in situ underwater video camera (VIS) was deployed to collect video recordings of the suspended sediment. These VIS data were finally processed to fall velocity and diameter distributions and were mainly used to improve insight into the relevant transport processes, indicating significant erosion of sand from the upstream Rhine branch. For quantitative calibration of the model, the data from the turbidity profiles were used. Sedimentation and erosion were modelled according to Krone and Partheniades. The model results showed a good overall fit to the measurements, with a mean absolute error of 18 per cent (standard fault = 1 per cent), corresponding to concentrations of about 0·020 (upstream) to 0·005 kg m⁻³ (downstream). The overall correlation between observed and simulated suspended sediment concentrations was 0·85. The remote sensing video recordings were used for a qualitative calibration of the model. The distribution pattern of the suspended sediment on these photos was reproduced quite well by the model. However, a more accurate calibration technique is needed to enable the use of aerial remote sensing as a quantitative calibration method. Copyright © 1999 John Wiley & Sons, Ltd.     

Theoretical/numerical model for the transport of non-uniform suspended sediment in open channels

In this paper, we address the transport of multi-disperse suspended sediment mixtures in open channels, via the use of the two-fluid model. To that end, we extend previously developed frameworks for the dilute and non-dilute transport of suspended sediment. Within the scope of the Reynolds-averaged Navier-Stokes (RANS) equations, these modeling frameworks comprise mass and momentum equations for both phases (water and sediment). Here, we calculate the distribution of total volumetric concentration of sediment using two approaches: (1) by considering the mixture as represented by a single size; we call this approach Partial two-fluid model for uniform sediments (PTFMU); and (2) by combining the volumetric concentration of the sediment corresponding to several particle size classes; we call this approach Partial two-fluid model for non-uniform sediments (PTFMNU). In the second approach, we propose a methodology for the computation of the overall velocity of the disperse phase as a function of the velocities of each size class. k-ε type closures to account for the turbulence in the carrier phase (water) are applied. We also consider the coupling between the two phases through the drag force. Velocities of the carrier and disperse phases, and concentrations for each sediment class size are numerically solved by integrating the differential equations over control volumes. In order to validate our models, we compare numerical results to experimental data of Einstein and Chien [H.A. Einstein, N. Chien, Effects of heavy sediment concentration near the bed on velocity and sediment distribution, MRD sediment series report, University of California, Berkley, 1955] and Taggart et al. [W.C. Taggart, C.A. Yermoli, S. Montes, A. Ippen, Effects of sediment size and gradation on concentration profiles for turbulent flow, Massachusetts Institute of Technology, 1972]. Results of mean velocity of the carrier phase are in close agreement with the experimental data. For the prediction of sediment concentrations, we observe that there is a difference in the results using the two approaches mentioned above. We additionally obtain values of the Schmidt number needed to improve the agreement between predictions of the distribution of suspended sediment and the experimental data, and discuss the effect of sediment size and increasing sediment concentration on the values of the Schmidt number.     

What is suspended sediment?

Suspended sediment is conventionally regarded as that sediment transported by a fluid that it is fine enough for turbulent eddies to outweigh settling of the particles through the fluid. Early work in the fluvial field attributed suspension to turbulence, and led to the notion of a critical threshold for maintaining sediment in suspension. However, research on both turbulence structures and the interactions between suspended sediment and bedforms in rivers has shown a more complex story and, although there appear to have been no studies of the impact of bedforms on aeolian suspended sediment concentrations, turbulent flow structures and transport rates of saltating particles have been shown to be affected. This research indicates that suspended sediment neither travels with the same velocity as the flow in which it is suspended, nor is it likely to remain in suspension in perpetuity, even under conditions of steady flow or in unsteady flow the where dimensionless critical threshold is permanently exceeded. Rather, like bedload, it travels in a series of hops, and is repeatedly deposited on the bed where it remains until it is re‐entrained. Is there, therefore, a qualitative difference between suspended and saltating sediment, or is it just a quantitative difference in the size of the jump length and the frequency of re‐entrainment? It is our contention that the distinction of suspension as a separate class of sediment transport is both arbitrary and an unhelpful anthropocentric artefact. If we recognize that sediment transport is a continuum and applies to any fluid medium rather than split into different “processes” based on arbitrary thresholds and fluids, then recognizing the continuity will enable development of an holistic approach sediment transport, and thus sediment‐transport models that are likely to be viable across a wider range of conditions than hitherto. Copyright © 2015 John Wiley & Sons, Ltd.     

Modeling river dune development and dune transition to upper stage plane bed

Large asymmetric bedforms known as dunes commonly dominate the bed of sand rivers. Due to the turbulence generation over their stoss and lee sides, dunes are of central importance in predicting hydraulic roughness and water levels. During floods in steep alluvial rivers, dunes are observed to grow rapidly as flow strength increases, undergoing an unstable transition regime, after which they are washed out in what is called upper stage plane bed. This transition of dunes to upper stage plane bed is associated with high transport of bed sediment in suspension and large decrease in bedform roughness. In the present study, we aim to improve the prediction of dune development and dune transition to upper stage plane bed by introducing the transport of suspended sediment in an existing dune evolution model. In addition, flume experiments are carried out to investigate dune development under bed load and suspended load dominated transport regimes, and to get insight in the time scales related to the transition of dunes to upper stage plane bed. Simulations with the extended model including the transport of suspended sediment show significant improvement in the prediction of equilibrium dune parameters (e.g. dune height, dune length, dune steepness, dune migration rate, dune lee side slope) both under bed load dominant and suspended load dominant transport regimes. The chosen modeling approach also allows us to model the transition of dunes to upper stage plane bed which was not possible with the original dune evolution model. The extended model predicts change in the dune shapes as was observed in the flume experiments with decreasing dune heights and dune lee slopes. Furthermore, the time scale of dune transition to upper stage plane bed was quite well predicted by the extended model. Copyright © 2015 John Wiley & Sons, Ltd.     

Nearshore bar migration and sediment-induced buoyancy effects

A 1-D General Ocean Turbulence Model that includes the effects of sediment-induced stratification is shown to simulate the observed onshore and offshore migration of a nearshore sandbar. The only two free parameters of the model, the bed reference concentration and the sediment diffusivity, are taken from the literature, rather than tuned to the data used here. The model results suggest that predictions of onshore bar migration, in which wave-induced sediment transport confined to within a few centimeters of the bottom dominates, are not greatly affected by accounting for buoyancy effects. The model results also suggest that both mean flows and waves transport sediment during offshore bar migration, with different components of transport dominating at different cross-shore locations across the bar-trough bathymetry. Neglecting the effects of sediment-induced stratification results in higher model skill during the largest waves, likely because the excess turbulence production simulated by the non-stratified model is counterbalanced by neglected breaking-wave-generated turbulence. Considering both onshore and offshore migration, the model that includes sediment-induced stratification has higher skill than the model without stratification.     

A three dimensional model for water-sediment movement

A three-dimensional k-ε-Ap two-fluid turbulence model is proposed to study liquid-particle two-phase flow and bed deformation.By solving coupled liquid-phase and solid-phase governing equations in a finite-volume method,the model can calculate the movement of both water and sediment.The model was validated by water-sediment transport in a 180° channel bend with a movable bed.The validation concerns two-phase time-averaged velocities,bed deformation,water depth,depth-averaged streamwise velocity,cross-stream bed profiles,and two-phase secondary flow velocity vectors.The agreement between numerical results and experimental results was generally good.The comparisons of the numerical results of different models show that the three-dimensional k-ε-Ap two-fluid turbulence model has a relatively higher accuracy than one-fluid model.     

Lithology‐controlled evolution of stream bed sediment and basin‐scale sediment yields in adjacent mountain watersheds,Idaho, USA

The composition, grain‐size, and flux of stream sediment evolve downstream in response to variations in basin‐scale sediment delivery, channel network structure, and diminution during transport. Here, we document downstream changes in lithology and grain size within two adjacent ~300 km² catchments in the northern Rocky Mountains, USA, which drain differing mixtures of soft and resistant rock types, and where measured sediment yields differ two‐fold. We use a simple erosion–abrasion mass balance model to predict the downstream evolution of sediment flux and composition using a Monte Carlo approach constrained by measured sediment flux. Results show that the downstream evolution of the bed sediment composition is predictably related to changes in underlying geology, influencing the proportion of sediment carried as bedload or suspended load. In the Big Wood basin, particle abrasion reduces the proportion of fine‐grained sedimentary and volcanic rocks, depressing bedload in favor of suspended load. Reduced bedload transport leads to stronger bed armoring, and coarse granitic rocks are concentrated in the stream bed. By contrast, in the North Fork Big Lost basin, bedload yields are three times higher, the stream bed is less armored, and bed sediment becomes dominated by durable quartzitic sandstones. For both basins, the geology‐based mass balance model can reproduce within ~5% root‐mean‐square error the composition of the bed substrate using realistic erosion and abrasion parameters. As bed sediment evolves downstream, bedload fluxes increase and decrease as a function of the abrasion parameter and the frequency and size of tributary junctions, while suspended load increases steadily. Variable erosion and abrasion rates produce conditions of variable bed‐material transport rates that are sensitive to the distribution of lithologies and channel network structure, and, provided sufficient diversity in bedrock geology, measurements of bed sediment composition allow for an assessment of sediment source areas and yield using a simple modeling approach. Copyright © 2016 John Wiley & Sons, Ltd.     

Incompressible SPH scour model for movable bed dam break flows

In this study an incompressible smoothed particle hydrodynamics (ISPH) approach coupled with the sediment erosion model is developed to investigate the sediment bed scour and grain movement under the dam break flows. Two-phase formulations are used in the ISPH numerical algorithms to examine the free surface and bed evolution profiles, in which the entrained sediments are treated as a different fluid component as compared with the water. The sediment bed erosion model is based on the concept of pick-up flow velocity and the sediment is initiated when the local flow velocity exceeds a critical value. The proposed model is used to reproduce the sediment erosion and follow-on entrainment process under an instantaneous dam break flow and the results are compared with those from the weakly compressible moving particle semi-implicit (WCMPS) method as well as the experimental data. It has been demonstrated that the two-phase ISPH model performed well with the experimental data. The study shows that the ISPH modelling approach can accurately predict the dynamic sediment scouring process without the need to use empirical sediment transport formulas.     

2D NUMERICAL SIMULATION OF FLOOD AND FLUVIAL PROCESS IN THE MEANDERING AND ISLAND-BRAIDED MIDDLE YANGTZE RIVER

The characteristics of water flow and sediment transport in a typical meandering and island-braided reach of the middle Yangtze River is investigated using a two-dimensional （2D） mathematical model. The major problems studied in the paper include the carrying capacity for suspended load, the incipient velocity and transport formula of non-uniform sediment, the thickness of the mixed layer on the riverbed, and the partitioning of bed load and suspended load. The model parameters are calibrated using extensive field data. Water surface profiles, distribution of flow velocities, riverbed deformation are verified with site measurements. The model is applied to a meandering and island-braided section of the Wakouzi-Majiazui reach in the middle Yangtze River, which is about 200 km downstream from the Three Gorges Dam, to study the training scheme of the navigation channels. The model predicts the processes of sediment deposition and fiver bed erosion, changes of flow stage and navigation conditions for the first 20 years of impoundment of the Three Gorges Project.     

Relating millimeter-scale turbulence to meter-scale subtidal erosion and accretion across the fringe of a coastal mangrove forest

Within a wave-exposed mangrove forest, novel field observations are presented, comparing millimeter-scale turbulent water velocity fluctuations with contemporaneous subtidal bed elevation changes. High-resolution velocity and bed level measurements were collected from the unvegetated mudflat, at the mangrove forest fringe, and within the forest interior over multiple tidal cycles (flood–ebb) during a 2-week period. Measurements demonstrated that the spatial variability in vegetation density is a control on sediment transport at sub-meter scales. Scour around single and dense clusters of pneumatophores was predicted by a standard hydraulic engineering equation for wave-induced scour around regular cylinders, when the cylinder diameter in the equations was replaced with the representative diameter of the dense pneumatophore clusters. Waves were dissipated as they propagated into the forest, but dissipation at infragravity periods (> 30 s) was observed to be less than dissipation at shorter periods (< 30 s), consistent with the predictions of a simple model. Cross-wavelet analysis revealed that infragravity-frequency fluctuations in the bed level were occasionally coherent with velocity, possibly indicating scour upstream of dense pneumatophore patches when infragravity waves reinforced tidal currents. Consequently, infragravity waves were a likely driver of sediment transport within the mangrove forest. Near-bed turbulent kinetic energy, estimated from the turbulent dissipation rate, was also correlated with bed level changes. Specifically, within the mangrove forest and over the unvegetated mudflat, high-energy events were associated with erosion or near-zero bed level change, whereas low-energy events were associated with accretion. In contrast, no single relationship between bed level changes and mean current velocity was applicable across both vegetated and unvegetated regions. These observations support the theory that sediment mobilization scales with turbulent energy, rather than mean velocity, a distinction that becomes important when vegetation controls the development of turbulence.     

Fully coupled mathematical modeling of turbidity currents over erodible bed

Turbidity currents may feature active sediment transport and rapid bed deformation, such as those responsible for the erosion of many submarine canyons. Yet previous mathematical models are built upon simplified governing equations and involve steady flow and weak sediment transport assumptions, which are not in complete accordance with rigorous conservation laws. It so far remains unknown if these could have considerable impacts on the evolution of turbidity currents. Here a fully coupled modeling study is presented to gain new insights into the evolution of turbidity currents. The recent analysis of the multiple time scales of subaerial sediment-laden flows over erodible bed [Cao Z, Li Y, Yue Z. Multiple time scales of alluvial rivers carrying suspended sediment and their implications for mathematical modeling. Adv Water Resour 2007;30(4):715–29] is extended to subaqueous turbidity currents to complement the fully coupled modeling. Results from numerical simulations show the ability of the present coupled model to reproduce self-accelerating turbidity currents. Comparison among the fully and partially coupled and decoupled models along with the analysis of the relative time scale of bed deformation explicitly demonstrate that fully coupled modeling is essential for refined resolution of those turbidity currents featuring active sediment transport and rapid bed deformation, and existing models based on simplified conservation laws need to be reformulated.     

One-dimensional morphodynamic model for retrogressive erosion based on a sediment entrainment theory at high flow velocity

Retrogressive erosion is a high-speed erosion process that usually occurs during the rapid release of stored water in reservoirs built on sandy rivers.Retrogressive erosion has been utilized in the practice of reservoir sedimentation control,but accurate prediction of the bed deformation process by numerical models has rarely been reported.The current study presents a one-dimensional morphodynamic model for simulating the evolution process of retrogressive erosion induced by high-velocity flows on steep slopes.The governing equations apply a Cartesian coordinate system with a vertically oriented z axis.The bed surface gradient and friction terms in the flow equations include correction factors to take account of the effects of high slope on flow movement.The net vertical sediment flux term in the sediment transport and bed deformation equations is calculated using an equation of erosion velocity.Particularly,this equation is based on an empirical relation between the sediment entrainment rate and the Shields parameter in contrast to the traditional sediment transport capacity,and the critical Shields parameter is modified by taking into account the permeability of the sediment layer and the stability of particles on a slope.The feedback of scoured sediment on the flow movement is considered by additional terms in the governing equations.Flume experiments of retrogressive erosion in literature were simulated to validate the model.The temporal variations of the longitudinal profiles of the free surface and channel bed and the sediment transport rate were well predicted.The algorithm calculating sediment entrainment in the proposed model also was validated for an experiment measuring entrainment rate from the literature.More importantly,it was found that the morphodynamic model using the sediment transport capacity equation predicts the trend of cumulative erosion contrary to the measurements,while results of the proposed model can follow a similar trend with the observed data in the retrogressive erosion process.     

Application of Hausdorff fractal derivative to the determination of the vertical sediment concentration distribution

The Rouse formula and its variants have been widely used to calculate the steady-state vertical concentration distribution for suspended sediment in steady sediment-laden flows, where the diffusive flux is assumed to be Fickian. Turbulent flow, however, exhibits fractal properties, leading to non-Fickian diffusive flux for sediment particles. To characterize non-Fickian dynamics of suspended sediment, the current study proposes a Hausdorff fractal derivative based advection-dispersion equation(H...     

Two-phase modeling of sheet-flow beneath waves and its dependence on grain size and streaming

We study erosion depth and sediment fluxes for wave-induced sheet-flow, and their dependency on grain size and streaming. Hereto, we adopt a continuous two-phase model, applied before to simulate sheet-flow of medium and coarse sized sand. To make the model applicable to a wider range of sizes including fine sand, it appears necessary to adapt the turbulence closure of the model. With an adapted formulation for grain–carrier flow turbulence interaction, good reproductions of measured erosion depth of fine, medium and coarse sized sand beds are obtained. Also concentration and velocity profiles at various phases of the wave are reproduced well by the model. Comparison of sediment flux profiles from simulations for horizontally uniform oscillatory flow as in flow tunnels and for horizontally non-uniform flow as under free surface waves, shows that especially for fine sand onshore fluxes inside the sheet-flow layer increase under influence of progressive wave effects. This includes both the current-related and the wave-related contribution to the period-averaged sheet-flow sediment flux. The simulation results are consistent with trends for fine and medium sized sediment flux profiles observed from tunnel and flume experiments. This study shows that the present two-phase model is a valuable instrument for further study and parameterization of sheet-flow layer processes.     

Scouring of a granular bed by dam break: Experimental study and numerical simulation by a VOF-LPT coupling

The first part of the research reported here consists of an experimental campaign to study the scouring of a granular bed(glass beads, sand) induced by a dam break in an open channel. Two configurations are considered: with and without cylinders. In the second part of this study, the volume of fluid method coupled with the shear stress transport turbulent model and the lagrangian particle tracking method is used to simulate the local scour processes. The four-way coupling is realized by consider...