An analysis of the characteristics of rough bed turbulent shear stresses in an open channel

Entrainment of sediment particles from channel beds into the channel flow is influenced by the characteristics of the flow turbulence which produces stochastic shear stress fluctuations at the bed. Recent studies of the structure of turbulent flow has recognized the importance of bursting processes as important mechanisms for the transfer of momentum into the laminar boundary layer. Of these processes, the sweep event has been recognized as the most important bursting event for entrainment of sediment particles as it imposes forces in the direction of the flow resulting in movement of particles by rolling, sliding and occasionally saltating. Similarly, the ejection event has been recognized as important for sediment transport since these events maintain the sediment particles in suspension. In this study, the characteristics of bursting processes and, in particular, the sweep event were investigated in a flume with a rough bed. The instantaneous velocity fluctuations of the flow were measured in two-dimensions using a small electromagnetic velocity meter and the turbulent shear stresses were determined from these velocity fluctuations. It was found that the shear stress applied to the sediment particles on the bed resulting from sweep events depends on the magnitude of the turbulent shear stress and its probability distribution. A statistical analysis of the experimental data was undertaken and it was found necessary to apply a Box-Cox transformation to transform the data into a normally distributed sample. This enabled determination of the mean shear stress, angle of action and standard error of estimate for sweep and ejection events. These instantaneous shear stresses were found to be greater than the mean flow shear stress and for the sweep event to be approximately 40 percent greater near the channel bed. Results from this analysis suggest that the critical shear stress determined from Shield's diagram is not sufficient to predict the initiation of motion due to its use of the temporal mean shear stress. It is suggested that initiation of particle motion, but not continuous motion, can occur earlier than suggested by Shield's diagram due to the higher shear stresses imposed on the particles by the stochastic shear stresses resulting from turbulence within the flow.     

Grain-size distribution in suspension over a sand-gravel bed in open channel flow

Grain-size distributions of suspended load over a sand-gravel bed at two different flow velocities were studied in a laboratory flume.The experiments had been performed to study the influence of flow velocity and suspension height on grain-size distribution in suspension over a sand-gravel bed.The experimental findings show that with an increase of flow velocity,the grain-size distribution of suspended load changed from a skewed form to a bimodal one at higher suspension heights.This study focuses on the determination of the parameter β_n which is the ratio of the sediment diffusion coefficient to the momentum diffusion coefficient of n th grain-size.A new relationship has been proposed involving β_n,the normalizing settling velocity of sediment particles and suspension height,which is applicable for widest range of normalizing settling velocity available in literature so far.A similar parameter β for calculating total suspension concentration is also developed.The classical Rouse equation is modified with β_n and β and used to compute grain-size distribution and total concentration in suspension,respectively.The computed values have shown good agreement with the measured values of experimental data.     

Hindered settling with an apparent particle diameter concept

In this work a new expression has been developed to predict the settling velocity of a sediment particle which is dispersed in a sediment-fluid mixture during a turbulent flow. A concept of apparent particle diameter has been introduced and is defined by the diameter of the spherical volume in which the particle can move randomly after collision with other particles in suspension. The effect of suspension concentration is studied on the mass density of the sediment-fluid mixture. It has been shown that the settling velocity of sediment particle in a sediment-fluid mixture is a function of different characteristics of the sediment particle such as settling velocity in clear fluid, suspension concentration, relative mass density and Reynolds number. The model has shown good agreement when compared with previously published experimental data and it’s prediction accuracy is superior than the other existing models.     

Sediment transport in ice-covered channels

The existence of ice cover has important effects on sediment transport and channel morphology for rivers in areas with an annual occurrence of an ice season. The interaction of sediment transport and s...     

Flow analysis in a channel with flexible vegetation using double-averaging method

The paper addresses the problem of the resistance due to vegetation in an open channel flow, characterized by partially and fully submerged vegetation formed by colonies of bushes. The flow is characterized by significant spatial variations of velocity between vertical profiles that make the traditional approach based on time averaging of turbulent fluctuations inconvenient. A more useful procedure, based on time and spatial averaging (Double-Averaging Method) is applied for the flow field analysis and characterization. The vertical distribution of mean velocity and turbulent stresses at different spatial locations has been measured with a 3D Acoustic Doppler Velocimeter (ADV) for two different vegetation densities where fully submerged real bushes (salix pentandra) have been used. Velocity measurements were completed together with the measurements of drag exerted on the flow by bushes at different flow depths. The analysis of velocity measurements allows depicting the fundamental characteristics of both the mean flow field and turbulence. The experimental data show that the contribution of form-induced stresses to the momentum balance cannot be neglected. The mean velocity profiles and the spatially averaged turbulent intensity profiles allow inferring that the vegetation density is a driving parameter for the development of a mixing layer at the canopy top in the case of submerged vegetation. Moreover, the net upward turbulent momentum flux, evaluated with the methodology proposed by Lu and Willmarth (1973), appears to be damped for increased vegetation density; this finding can rationally explain the reduction of the suspended sediment transport capacity typically observed in free surface flows over a vegetated bed.     

Effects of vegetation on flow and sediment transport: comparative analyses and validation of predicting models

The presence of vegetation modifies flow and sediment transport in alluvial channels and hence the morphological evolution of river systems. Plants increase the local roughness, modify flow patterns and provide additional drag, decreasing the bed‐shear stress and enhancing local sediment deposition. For this, it is important to take into account the presence of vegetation in morphodynamic modelling. Models describing the effects of vegetation on water flow and sediment transport already exist, but comparative analyses and validations on extensive datasets are still lacking. In order to provide practical information for modelling purposes, we analysed the performance of a large number of models on flow resistance, vegetation drag, vertical velocity profiles and bed‐shear stresses in vegetated channels. Their assessments and applicability ranges are derived by comparing their predictions with measured values from a large dataset for different types of submerged and emergent vegetation gathered from the literature. The work includes assessing the performance of the sediment transport capacity formulae of Engelund and Hansen and van Rijn in the case of vegetated beds, as well as the value of the drag coefficient to be used for different types of vegetation and hydraulic conditions. The results provide a unique comparative overview of existing models for the assessment of the effects of vegetation on morphodynamics, highlighting their performances and applicability ranges. Copyright © 2014 John Wiley & Sons, Ltd.     

Vertical concentration profile of nonuniform sediment

A new formula for the concentration profile of nonuniform sediment is derived using the diffusion differential equation with equilibrium bed concentration.The interaction coefficient for nonuniform sediment particles as a function of both relative diameter and geometric standard deviation of nonuniform sediment is taken into account in the settling velocity.The diffusion coefficient is obtained from the logarithmic velocity profile.This new formula possesses several advantages as follows:(1)avoids theoretical defect of the Rouse formula,which states that the sediment concentration is infinite at the bottom and zero at the water surface.(2)suitable for an arbitrary fraction of nonuniform sediment,and(3)easy to apply with a simple form.The formula yields less difference among concentration profiles for various particles than that for uniform sediment and reveals the lawfor fractional concentration distribution of nonuniform sediment.The calculated concentration agrees well with the measured data from the Yangtze River.     

Turbulence and suspended sediment processes in the Garonne River tidal bore in November 2016

A tidal bore is a water discontinuity at the leading edge of a ood tide wave in estuaries with a large tidal range and funneling topography. New measurements were done in the Garonne River tidal bore on 14 15 November 2016, at a site previously investigated between 2010 and 2015. The data focused on long, continuous, high-frequency records of instantaneous velocity and suspended sediment con- centration (SSC) estimate for several hours during the late ebb, tidal bore passage and ood tide. The bore passage drastically modi ed the ow eld, with very intense turbulent and sediment mixing. This was evidenced with large and rapid uctuations of both velocity and Reynolds stress, as well as large SSCs during the ood tide. Granulometry data indicated larger grain sizes of suspended sediment in water samples compared to sediment bed material, with a broader distribution, shortly after the tidal bore. The tidal bore induced a sudden suspended sediment ux reversal and a large increase in suspended sedi- ment ux magnitude. The time-variations of turbulent velocity and suspended sediment properties indicated large uctuations throughout the entire data set. The ratio of integral time scales of SSC to velocity in the x-direction was on average TE,SSC/TE,x 0.16 during the late ebb tide, compared to TE,SSC/ TE,x 0.09 during the late ood tide. The results imply different time scales between turbulent velocities and suspended sediment concentrations.     

Laboratory study on sediment diffusion and deposition into blind channels

Any waterway with one end closed and the other open is generally called a blind channel. The main flow tends to expand, separate, and cause circulation at the mouth of blind channels. The main flow continuously transfers momentum and sediment into the circulation region through the turbulent mixing region (TMR) between them, thus leading to a large amount of sediment deposition in the blind channels. This paper experimentally investigated the properties of the water flow and sediment diffusion in TMR, demonstrating that both water flow and sediment motion in TMR approximately coincide with a similar structure as in the free mixing layer induced by a jet. The similarity functions of flow velocity and sediment concentration are then assumed, based on observation, and the resulting calculation of these functions is substantially facilitated. For the kind of low velocity flow system of blind channels with a finite width, a simple formula for the sediment deposition rate in blind channels is established by analyzing the gradient of crosswise velocity and sediment concentration in TMR.     

On the use of horizontal acoustic Doppler profilers for continuous bed shear stress monitoring

Continuous monitoring of bed shear stress in large river systems may serve to better estimate alluvial sediment transport to the coastal ocean.Here we explore the possibility of using a horizontally deployed acoustic Doppler current profiler(ADCP) to monitor bed shear stress,applying a prescribed boundary layer model,previously used for discharge estimation.The model parameters include the local roughness length and a dip correction factor to account for sidewall effects.Both these parameters depend on river stage and on the position in the cross-section, and were estimated from shipborne ADCP data.We applied the calibrated boundary layer model to obtain bed shear stress estimates over the measuring range of the HADCP.To validate the results,co-located coupled ADCPs were used to infer bed shear stress,both from Reynolds stress profiles and from mean velocity profiles. From HADCP data collected over a period of 1.5 years,a time series of width profiles of bed shear stress was obtained for a tidal reach of the Mahakam River,East Kalimantan,Indonesia.A smaller dataset covering 25 hours was used for comparison with results from the coupled ADCPs.The bed shear stress estimates derived from Reynolds stress profiles appeared to be strongly affected by local effects causing upflow and downflow,which are not included in the boundary layer model used to derive bed shear stress with the horizontal ADCP.Bed shear stresses from the coupled ADCP are representative of a much more localized flow,while those derived with the horizontal ADCP resemble the net effect of the flow over larger scales.Bed shear stresses obtained from mean velocity profiles from the coupled ADCPs show a good agreement between the two methods,and highlight the robustness of the method to uncertainty in the estimates of the roughness length.     

Two-phase hydrodynamic and sediment transport modeling of wave-generated sheet flow

This numerical investigation was carried out to advance mechanistic understanding of sediment transport under sheet flow conditions. An Euler–Euler coupled two-phase flow model was developed to simulate fluid–sediment oscillatory sheet flow. Since the concentration of sediment particles is high in such flows, the kinematics of the fluid and sediment phases are strongly coupled. This model includes interaction forces, intergranular stresses and turbulent stress closure. Each phase was modeled via the Reynolds-Averaged Navier–Stokes equations, with interphase momentum conservation accounting for the interaction between the phases. The generation and transformation of turbulence was modeled using the two-equation k–ε$k''–''\epsilon$ turbulence model. Concentration and sediment flux profiles were compared with experimental data for sheet flow conditions considering both symmetric and asymmetric oscillatory flows. Sediment and fluid velocity variations, concentration profiles, sediment flux and turbulence parameters of wave-generated sheet flow were studied numerically with a focus on sediment transport characteristics. In all applications, the model predictions compared well with the experimental data. Unlike previous investigations in which the flow is driven by a horizontal pressure gradient, the present model solves the Navier–Stokes equations under propagating waves. The model’s ability to predict sediment transport under oscillatory sheet flow conditions underscores its potential for understanding the evolution of beach morphology.     

Spatial variability of three‐dimensional Reynolds stresses in a developing channel bend

Experimental results of the mean flow field and turbulence characteristics for flow in a model channel bend with a mobile sand bed are presented. Acoustic Doppler velocimeters (ADVs) were used to measure the three components of instantaneous velocities at multiple cross sections in a 135° channel bend for two separate experiments at different stages of clear water scour conditions. With measurements at multiple cross sections through the bend it was possible to map the changes in both the spatial distribution of the mean velocity field and the three Reynolds shear stresses. Turbulent stresses are known to contribute to sediment transport and the three‐dimensionality inherent to flow in open channel bends presents a useful case for determining specific relations between three‐dimensional turbulence and sediment entrainment and transport. These measurements will also provide the necessary data for validating numerical simulations of turbulent flow and sediment transport. The results show that the magnitude and distribution of three‐dimensional Reynolds stresses increase through the bend, with streamwise‐cross stream and cross stream‐vertical components exceeding the maximum principal Reynolds stress through the bend. The most intriguing observation is that near‐bed maximum positive streamwise‐cross stream Reynolds stress coincides with the leading edge of the outer bank scour hole (or thalweg), while maximum cross stream‐vertical Reynolds stress (in combination with high negative streamwise‐cross stream Reynolds stress near the bend apex) coincides with the leading edge of the inner bank bar. Maximum Reynolds stress and average turbulent kinetic energy appear to be greater and more localized over the scour hole before final equilibrium scour is reached. This suggests that the turbulent energy in the flow is higher while the channel bed is developing, and both lower turbulent energy and a broader distribution of turbulent stresses near the bed are required for cessation of particle mobilization and transport. Copyright © 2010 John Wiley & Sons, Ltd.     

Issues in Eulerian–Lagrangian modeling of sediment transport under saltation regime

The saltation regime is very important for understanding the sediment transport mechanism. However,there is no consensus on a model for the saltation regime. This study answers several questions raised with respect to the Eulerian-Lagrangian modeling of sediment transport. The first question is why the previous saltation models that use different combinations of hydrodynamic forces yielded acceptable results? The second question is which shear lift model(i.e. a shear lift expression and its coefficient) is more appropriate? Another important question is which hydrodynamic forces have greater contributions to the saltation characteristics of a sediment particle? The last question is what are the contributions of the turbulence fluctuations as well as effects of using two-and three-dimensional(2 D and 3 D) models on the simulation results? In order to fairly answer these questions, a systematic study was done by considering different scenarios. The current study is the first attempt to clearly discuss these issues. A comprehensive 3 D saltation model for non-cohesive sediment was developed that includes all the hydrodynamic forces acting on the particle. The random nature of sediment transport was included using turbulent flow and bed-particle collision models. The eddy interaction model was applied to generate a3 D turbulent flow field. Bed-particle collisions were considered using the concept of a contact zone and a corresponding contact point. The validation of the model was done using the available experimental data for a wide range of sediment size(0.03 to 4.8 cm). For the first question, the results indicated that some of the hydrodynamic effects show opposing trends and some have negligible effects. With these opposing effects it is possible to adjust the coefficients of different models to achieve acceptable agreement with the same experimental data while omitting some aspects of the physics of the process. A suitable model for the shear lift force was developed by linking the lift coefficient to the drag coefficient and the contributions of the hydrodynamic forces and turbulence fluctuations as well as the consequences of using of 2 D and 3 D models were studied. The results indicate that the shear lift force and turbulent flow fluctuations are important factors for the saltation of both sand and gravel, and they cannot be ignored.     

Experimental investigation on local shear stress and turbulence intensities over a rough non-uniform bed with and without sediment using 2D particle image velocimetry

The current work focuses on locally resolving velocities,turbulence,and shear stresses over a rough bed with locally non-uniform character.A nonporous subsurface layer and fixed interfacial sublayer of gravel and sand were water-worked to a nature-like bed form and additionally sealed in a hydraulic flume.Two-dimensional Particle Image Velocimetry(2 D-PIV) was applied in the vertical plane of the experimental flume axis.Runs with clear water and weak sediment transport were done under slightly supercritical flow to ensure sediment transport conditions without formation of considerable sediment deposits or dunes.The study design included analyzing the double-averaged flow parameters of the entire measurement domain and investigating the flow development at 14 consecutive vertical subsections.Local geometrical variabilities as well the presence of sediment were mainly reflected in the vertical velocity component.Whereas the vertical velocity decreased over the entire depth in presence of sediment transport,the streamwise velocity profile was reduced only within the interfacial sublayer.In the region with decelerating flow conditions,however,the streamwise velocity profile systematically increased along the entire depth extent.The increase in the main velocity(reduction of flow resistance)correlated with a decrease of the turbulent shear and main normal stresses.Therefore,effects of rough bed smoothening and drag force reduction were experimentally documented within the interfacial sublayer due to mobile sediment.Moreover,the current study leads to the conclusion that in nonuniform flows the maximum Reynolds stress values are a better predictor for the bed shear stress than the linearly extrapolated Reynolds stress profile.This is an important finding because,in natural flows,uniform conditions are rare.     

Characteristics of energy dissipation in hyperconcentrated flows

An equilibrium equation for the turbulence energy in of solid-liquid two-phase flow theory. The equation sediment-laden flows was derived on the basis was simplified for two-dimensional, uniform, steady and fully developed turbulent hyperconcentrated flows. An energy efficiency coefficient of suspended-load motion was obtained from the turbulence energy equation, which is defined as the ratio of the sediment suspension energy to the turbulence energy of the sediment-laden flows. Laboratory experiments were conducted to investigate the characteristics of energy dissipation in hyperconcentrated flows. A total of 115 experimental runs were carried out, comprising 70 runs with natural sediments and 45 runs with cinder powder. Effects of sediment concentration on sediment suspension energy and flow resistance were analyzed and the relation between the energy efficiency coefficient of suspended-load motion and sediment concentration was established on the basis of experimental data. Furthermore, the characteristics of energy dissipation in hyperconcentrated flows were identified and described. It was found that the high sediment concentration does not increase the energy dissipation; on the contrary, it decreases flow resistance.     

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...     

A field study of mean and turbulent flow characteristics upwind,over and downwind of barchan dunes

Field‐measured patterns of mean velocity and turbulent airflow are reported for isolated barchan dunes. Turbulence was sampled using a high frequency sonic anemometer, deriving near‐surface Reynolds shear and normal stresses. Measurements upwind of and over a crest‐brink separated barchan indicated that shear stress was sustained despite a velocity reduction at the dune toe. The mapped streamline angles and enhanced turbulent intensities suggest the effects of positive streamline curvature are responsible for this maintenance of shear stress. This field evidence supports an existing model for dune morphodynamics based on wind tunnel turbulence measurements. Downwind, the effect of different dune profiles on flow re‐attachment and recovery was apparent. With transverse incident flow, a re‐attachment length between 2·3 and 5·0h (h is dune brink height) existed for a crest‐brink separated dune and 6·5 to 8·6h for a crest‐brink coincident dune. The lee side shear layer produced elevated turbulent stresses immediately downwind of both dunes, and a decrease in turbulence with distance characterized flow recovery. Recovery of mean velocity for the crest‐brink separated dune occurred over a distance 6·5h shorter than the crest‐brink coincident form. As the application of sonic anemometers in aeolian geomorphology is relatively new, there is debate concerning the suitability of processing their data in relation to dune surface and streamline angle. This paper demonstrates the effect on Reynolds stresses of mathematically correcting data to the local streamline over varying dune slope. Where the streamline angle was closely related to the surface (windward slope), time‐averaged shear stress agreed best with previous wind tunnel findings when data were rotated along streamlines. In the close lee, however, the angle of downwardly projected (separated) flow was not aligned with the flat ground surface. Here, shear stress appeared to be underestimated by streamline correction, and corrected shear stress values were less than half of those uncorrected. Copyright © 2011 John Wiley & Sons, Ltd.