Fluctuations of Suspended Sediment Concentration and Turbulent Sediment Fluxes in an Open-Channel Flow

A complex problem of turbulent-sediment interactions in an open-channel flow is approached experimentally, using specially designed field experiments in an irrigation canal. The experimental design included synchronous measurements of instantaneous three-dimensional (3D) velocities and suspended sediment concentration using acoustic Doppler velocimeters (ADV) and a water sampling system. Various statistical measures of sediment concentration fluctuations, turbulent sediment fluxes, and diffusion coefficients for fluid momentum and sediment are considered. Statistics, fractal behavior, and contributions of bursting events to vertical fluxes of fluid momentum and sediment are evaluated using quadrant analysis. It has been found that both turbulence and sediment events are organized in fractal clusters which introduce additional characteristic time and spatial scales into the problem and should be further explored. It is also shown that Barenblatt’s theory of sediment-laden flows appears to be a good approximation of experimental data.     

Analysis of Suspended Sediment Transport in Open-Channel Flows: Kinetic-Model-Based Simulation

This paper presents a comprehensive analysis of suspended sediment transport in open channels under various flow conditions through a kinetic-model-based simulation. The kinetic model, accounting for both sediment-turbulence and sediment-sediment interactions, successfully represents experimentally observed diffusion and transport characteristics of suspended sediments with different densities and sizes. Without tuning any model coefficients, the nonmonotonic concentration distribution and the noticeable lag velocity with a negative value close to the wall are reasonably reproduced. Examination of flow conditions typical of suspension dominated rivers shows that the conventional method may overestimate or underestimate unit suspended-sediment discharge, depending on the Rouse number, sediment size, as well as shear velocity. The error may be less than 20% for dp<0.5?mm and might exceed 60% for dp>1.0?mm under typical flow conditions where shear velocity ranges from 1.0?to?12.5?cm/s and flow depth ranges from 1.0?to?5.0?m.     

Vertical Mixing in the Fully Developed Turbulent Layer of Sediment-Laden Open-Channel Flow

Eulerian equations for the vertical flux and momentum of suspended particles in dilute sediment-laden open-channel flow in equilibrium have been derived using the two-fluid approach. Reynolds averaging has been applied in order to allow validation of individual terms with experimental data. Consideration of the various terms of the vertical momentum balance with experimental flume data indicates that the drag force has to be separated into a mean and a turbulent contribution with different timescales, respectively, the (gravitational) particle timescale and the integral turbulence timescale. The resulting formulation further provides a new theoretical closure for the turbulent Schmidt number.     

Modeling Nonuniform Suspended Sediment Transport in Alluvial Rivers

Problems and difficulties in modeling sediment transport in alluvial rivers arise when one uses the theory of equilibrium transport of uniform sediment to simulate riverbed variation. A two-dimensional mathematical model for nonuniform suspended sediment transport is presented to simulate riverbed deformation. Through dividing sediment mixture into several size groups in which the sediment particles are thought to be uniform, the nonuniformity and the exchange between suspended sediment and bed material are considered. The change of concentration along the flow direction, size redistribution, and cross-sectional bed variation can then be described reasonably well by the model. In simulating the flow field with big dry-wet flats, moving boundary problems are solved very well by introducing a so-called finite-slot technique. Verification with laboratory data shows that the model has a good ability to simulate channel bed variations. Last, the model was applied to a real alluvial river system. Variables such as water level, sediment concentration, suspended sediment size distribution, and riverbed variation were obtained with encouraging results.     

Monitoring Suspended Sediment Dynamics Using MBES

This technical note describes use of a multibeam echosounder (MBES) to quantify the dynamics of suspended sediment in a large open channel. A methodology is detailed that uses the backscatter magnitude from the MBES water-column data to adjust the magnitude of sonar returns for the various sonar settings, spatially and temporally average the data to account for the random nature of acoustic backscatter from the suspended sediment, and calibrate the processed data with direct samples. A case study of flow at the confluence of the Rio Paraná and Rio Paraguay, Argentina, where there is a distinct turbidity difference along the mixing interface of the two flows, is used to demonstrate the unique capabilities of MBES to quantify sediment concentrations and dynamics within the water column.     

Diagonal Cartesian Method for the Numerical Simulation of Flow and Suspended Sediment Transport over Complex Boundaries

There is increasing demand for simulation tools of flow and suspended sediment transport over complex boundaries in hydraulic engineering. The diagonal Cartesian method, which approximates complex boundaries using both Cartesian grid lines and diagonal lines segments, is presented in the paper to simulate the complex boundaries of two-dimensional shallow-water turbulence equations and nonequilibrium suspended sediment transport equation. The method, which utilizes cell-centered nodes on a nonstaggered grid, uses boundary velocity information at the wall boundary to avoid the specification of water level. An enlarged finite-difference method is introduced for momentum and suspended sediment equations on the complex boundary. This paper describes an application of the diagonal Cartesian method to calculate the tidal current and suspended sediment concentration of Quanzhou Bay in the Fujian province of China. The results show that the method predicts the flow and suspended sediment concentration well, and the calculations agree well with the measurement.     

Reynolds Stress Anisotropy in Open-Channel Flow

This paper reports the results of an experimental study characterizing turbulence and turbulence anisotropy in smooth and rough shallow open-channel flows. The rough bed consists of a train of two-dimensional transverse square ribs with a ratio of the roughness height (k) to the total depth of flow (d) equal to 0.10. Three rib separations (p/k) of 4.5, 9, and 18 were examined. Here, p is the pitch between consecutive roughness elements and was varied to reproduce the classical condition of d- and k-type roughness. For each case, two-component velocity measurements were obtained using a laser Doppler velocimetry system at two locations for p/k = 4.5 and 9: on the top of the rib and above the cavity, and an additional location for p/k = 18. The measurements allow examination of the local variations of the higher-order turbulent moments, stress ratios as well as turbulence anisotropy. Large variations of the turbulence intensities, Reynolds shear stress, turbulent kinetic energy and turbulence production are found for y1<3k. In this region, the flow is more directly influenced by the shear layers from the preceding ribs. The higher-order moments appear to be similar for all rough surfaces beyond y1 ≈ 7k. In the outer layer (y1>3k), all higher-order turbulent moments for the k-type roughness show a substantial increase due to the complex interactions between the roughness and the remnants overlying shear layers shed from succeeding ribs. Analysis of the components of the Reynolds stress anisotropy tensor shows that at p/k = 18, the flow at y1<5k tends to be more isotropic which implies that for this particular case, the effect of the roughness density could also be important. On the smooth bed, at the shallower depths, the correlation coefficient near the free surface increases and turbulence tends to become less anisotropic.     

Unified View of Sediment Transport by Currents and Waves. II: Suspended Transport

The problem of suspended sediment transport in river and coastal flows is addressed. High-quality field data of river and coastal flows have been selected and clustered into four particle size classes (60–100, 100–200, 200–400, and 400–600?μm). The suspended sand transport is found to be strongly dependent on particle size and on current velocity. The suspended sand transport in the coastal zone is found to be strongly dependent on the relative wave height (Hs/h), particularly for current velocities in the range 0.2–0.5?m/s. The time-averaged (over the wave period) advection–diffusion equation is applied to compute the time-averaged sand concentration profile for combined current and wave conditions. Flocculation, hindered settling, and stratification effects are included by fairly simple expressions. The bed-shear stress is based on a new bed roughness predictor. The reference concentration function has been recalibrated using laboratory and field data for combined steady and oscillatory flow. The computed transport rates show reasonably good agreement (within a factor of 2) with measured values for velocities in the range of 0.6–1.8?m/s and sediments in the range of 60–600?μm. The proposed method underpredicts in the low-velocity range (<0.6?m/s). A new simplified transport formula is presented, which can be used to obtain a quick estimate of suspended transport. The modeling of wash load transport in river flow based on the energy concept of Bagnold shows that an extremely large amount of very fine sediment (clay and very fine silt) can be transported by the flow.     

Bed-Load Effects on Hydrodynamics of Rough-Bed Open-Channel Flows

The extent to which turbulent structure is affected by bed-load transport is investigated experimentally using a nonporous fixed planar bed comprising mixed-sized granular sediment with a d50 of 1.95?mm. Three different sizes of sediment (d50 = 0.77, 1.99, and 3.96?mm) were fed into the flow at two different rates (0.003 and 0.006?kg/m/s), and subsequently transported as bed load. Particle image velocimetry (PIV) was used to determine the turbulence characteristics over the fixed bed during clear water and sediment feed cases. Mean longitudinal flow velocities at any given depth were lower than their clear water counterparts for all but one of the mobile sediment cases. The exception was with the transport of fine grains at the higher feed rate. In this case, longitudinal mean flow velocities increased compared to the clear water condition. The coarse grains tended to augment bed roughness, but fine grains saturated the troughs and interstices in the bed topography, effectively causing the influence of bed irregularities to be smoothed. The PIV technique permitted examination of both temporal and spatial fluctuations in flow variables: therefore many results are presented in terms of double-averaged quantities (in temporal and spatial domains). In particular, the form-induced stress, which arises from spatially averaging the Reynolds averaged Navier–Stokes equations and is analogous to the Reynolds turbulent stress, contributed between 15 and 35% of the total measured shear stress in the roughness layer. Flow around protrusive roughness elements produced a significant proportion of the turbulent kinetic energy shear production, suggesting that this process is highly intermittent near rough beds.     

Bed Sediment Entrainment Function Based on Kinetic Theory

Bed sediment entrainment function is derived based on kinetic theory. The entrainment rate is expressed as an upward flux by integrating the particle velocity over all possible upward motions. The possibility of upward motion is determined by a sediment particle velocity distribution function, which is obtained by solving the Boltzmann equation instead of using a prior assumed distribution. External forces and turbulence intensity of the flow are shown to exert significant influences on the velocity distribution function and, in turn, the entrainment rate. Comparisons between available laboratory data and the entrainment function show that the calculated entrainment rates agree well with the observations. Applications of the entrainment function to the specification of the bottom boundary condition for convection-diffusion equation of suspended load are also presented, which show that the calculated concentration profiles are in good agreement with observations. The study also suggests that kinetic theory is a promising analytical approach for the study of sediment motion near a riverbed.     

Modeling Suspended Sediment during Construction in Great Barrier Reef World Heritage Area

A marina was constructed in the Great Barrier Reef World Heritage Area in close proximity to coral reefs that could be damaged by excess turbidity generated during construction. Since there was uncertainty about both the fate of suspended sediments and their effect on corals, initial water quality constraints were set very conservatively. In order to better understand the movement of suspended sediment during construction, a numerical model study was commissioned using three-dimensional, numerical, hydrodynamic, and Lagrangian particle tracking models. The study was successful in: (1) increasing the understanding of and reducing the uncertainty of sediment dispersal patterns under a range of common forcing conditions; (2) testing the variation in suspended sediment concentrations over sensitive areas for two different outfall locations; (3) offering evidence that a good choice in outfall locations will reduce the threat to corals; and importantly (4) presenting the results in a way that enhanced understanding by nontechnical reef managers. This final result was achieved by creating movies of sediment movement that clearly demonstrated the complex hydrodynamic processes involved with near-coastal water currents. Specific model results showed: (1) that a more seaward outfall increases effluent dispersal away from sensitive areas; (2) the highest concentrations of effluent over sensitive sites occur during no wind and neap tide conditions; and (3) prevailing southeast winds advect effluent offshore, away from sensitive sites.     

Flow in Open-Channel Embayments

This study investigates flows in a square and a rectangular embayment located on the side bank of an open channel. It is found that the main flow of the open channel induces a circulatory flow in the square and rectangular embayment and the center of the circulatory flow is shifted downstream in comparison with the center of the embayment. A solution of the shallow water equations solved using the method of variation yields results of the 95% confidence intervals within 10% of mean errors between the observed and computed nondimensional velocities.     

General Characteristics of Solutions to the Open-Channel Flow, Feedforward Control Problem

A dimensionless formulation of the open-channel flow equations was used to study the feedforward control problem for single-pool canals. Feedforward inflow schedules were computed for specified downstream demands using a gate-stroking model. The analysis was conducted for various design and operational conditions. Differences in the shape of the computed inflow hydrographs are largely related to the volume change resulting from the transient, the time needed to supply this volume, and the time needed by the inflow perturbation to travel down the canal. The gate-stroking method will fail to produce a solution or the solution will demand extreme and unrealistic inflow variations if the time needed to supply the canal volume change is much greater than the travel time of the upstream flow change. As an alternative, a simple feedforward-control flow schedule can be developed based on this volume change and a reasonable delay estimate. This volume compensating schedule can deliver the requested flow change and keep water levels reasonably close to the target under the range of conditions tested.     

Turbulence Characteristics of Open-Channel Flow with Bed Suction

The influence of bed suction on the characteristics of turbulent open channel flow is studied in a laboratory flume using a two-component laser Doppler velocimeter. The experimental results show how bed suction significantly affects the mean flow properties, turbulence levels, and Reynolds stress distributions. The data reveal the presence of a more negative vertical (downward) velocity. The results also show how the horizontal and vertical turbulence intensities and Reynolds shear stresses respond to suction. All these properties are found to reduce with increasing relative suctions: decreasing more rapidly around the bed region than that near the free surface. In the downstream direction, the flow structure in the suction zone undergoes a process of rapid readjustment within a transitional region. Beyond this region, the turbulence flow structures asymptotes toward an “equilibrium” region.     

Vertical Dispersion of Fine and Coarse Sediments in Turbulent Open-Channel Flows

Through using a kinetic model for particles in turbulent solid–liquid flows, underlying mechanisms of sediment vertical dispersion as well as sediment diffusion coefficient are investigated. Four hydrodynamic mechanisms, namely gravitational settling, turbulent diffusion, effect of lift force, and that of sediment stress gradient, coexist in two-dimensional (2D) uniform and steady open-channel flows. The sediment diffusion coefficient consists of two independent components: one accounts for the advective transport of sediment probability density distribution function due to sediment velocity fluctuations, and the other results from sediment–eddy interactions. Predictions of the kinetic model are in good agreement with experimental data of 2D open-channel flows. In such flows, it is shown that: (1) the parameter γ (i.e., the inverse of the turbulent Schmidt number) may be greater than unity and increases toward the bed, being close to unity for fine sediments and considerably large for coarse ones; (2) effects of lift force and sediment stress gradient become significant and need to be considered below the 0.1 flow depth; and (3) large errors may arise from the traditional advection–diffusion equation when it is applied to flows with coarse sediments and/or high concentrations.     

Flow Resistance over Mobile Bed in an Open-Channel Flow

This paper deals with the underlying mechanism of flow resistance in an alluvial channel: The effects of sidewall and bed form on flow resistance, Einstein’s divided hydraulic radius approach and Engelund’s energy slope division approach are reexamined. These two approaches assume that the shear stress on a mobile bed is the summation of shear stresses caused by skin friction and bed form. Using a different approach, this paper presents a theoretical relationship between the total bed shear stress with grain and bed-form shear stresses. The contribution of sidewall on the total bed shear stress is also discussed. The writers found that the size of bed form plays a significant role for the flow resistance, and developed relevant expressions for the length of the separation zone behind the bed forms. In addition, a systematical approach has been developed to compute the flow velocity in an alluvial channel. This approach is tested and verified against 5,989 flume and field measurements. The computed and measured discharge/velocity are in good agreement and 83.0% of all data sets fall within the ±20% error band.     

Influence of Turbulence on Bed Load Sediment Transport

This paper summarizes the results of an experimental study on the influence of an external turbulence field on the bed load sediment transport in an open channel. The external turbulence was generated by (1) a horizontal pipe placed halfway through the depth h; (2) a series of grids with a clearance of about one-third of the depth from the bed, and extending over a finite length of the flume; and (3) a series of grids with a clearance in the range (0.1–1.0)h from the bed, but extending over the entire length of the flume. Two kinds of experiments were conducted: plane-bed experiments and ripple-covered-bed experiments. In the former case, the flow in the presence of the turbulence generator was adjusted so that the mean bed shear stress was the same as in the case without the turbulence generator in order to single out the effect of the external turbulence on the sediment transport. In the ripple-covered-bed case, the mean and turbulence quantities of the streamwise component of the velocity were measured, and the Shields parameter, due to skin friction, was determined. The Shields parameter, together with the RMS value of the streamwise velocity fluctuations, was correlated with the sediment transport rate. The sediment transport increases markedly with increasing turbulence level.     

Response of Velocity and Turbulence to Sudden Change of Bed Roughness in Open-Channel Flow

The experimental study shows how an open-channel flow would respond to a sudden change (from smooth to rough) in bed roughness. Using a two-dimensional acoustic Doppler velocimeter and a laser Doppler velocimeter, the velocity, turbulent intensities, and Reynolds stress profiles at different locations along a laboratory flume were measured. Additionally, the water surface profile was also measured using a capacitance-type wave height meter. The experimental data show the formation of an internal boundary layer as a result of the step change in bed roughness. The data show that this boundary layer grows much more rapidly than that formed in close-conduit flows. The results also show that the equivalent bed roughness, bed-shear stress, turbulent intensities, and Reynolds stress change gradually over a transitional region, although the bed roughness changes abruptly. The behavior is different from that observed in close-conduit flows, where an overshooting property—which describes the ability of the bed-shear stress to attain a high-peak value over the section with the larger roughness, was reported. A possible reason for the difference is the variation of the water surface profile when an open-channel flow is subjected to a sudden change in bed roughness.