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.     

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.     

Modeling In-Sewer Deposit Erosion to Predict Sewer Flow Quality

High levels of suspended solids are typically observed during the initial part of storms. Field evidence suggests that these suspended solids derive from the erosion of in-sewer sediment beds accumulated during dry and previous wet weather periods. Suspended sediment transport rate models within existing sewer network modeling tools have utilized inappropriate transport rate relationships developed mainly in fluvial environments. A process model that can account for the erosion of fine-grained highly organic in-sewer sediment deposits has been formulated. Values of parameters describing the increase in deposit strength with depth are required. These values are obtained using a genetic algorithm based calibration routine that ensures model simulations of suspended sediment concentrations that correspond to field data collected in a discrete length of sewer in Paris under known hydraulic event conditions. These results demonstrate the applicability of this modeling approach in simulating the magnitude and temporal distribution of suspended in-sewer sediment eroded by time varying flow. Further work is developing techniques to enable the application of this type of model at the network level.     

Rate of Deposition of Fine Sediment from Suspension

Standard depth-integrated models of sediment dynamics predict that concentrations of suspended fine sediment should decay at a characteristic rate that is controlled by the particle settling velocity and the depth of the water. In contrast, a model which resolves the processes of settling and dispersion in the water column has suggested that this decay rate should be independent of the settling velocity, and is controlled by dispersive processes in the water column. By revisiting the problem of sediment dispersion and settling following a point release of material, we resolve this discrepancy and confirm that depth-integrated models capture the correct physical behavior.     

Sediment Resuspension and Hydrodynamics in Lake Okeechobee during the Late Summer

In this study, in order to better understand the mechanisms affecting sediment resuspension, extensive data sets were collected in September and October of 2002 including wind velocity, wave, current velocity, water temperature, total suspended solid, suspended sediment concentration, and total phosphorus. Analyses of these data indicate that waves are the dominant factors in sediment resuspension while wind-induced currents are the primary forces to transport suspended sediments. Surface and bottom currents frequently flow in opposite directions, forming a stratified water column. A time lag exists between currents near lake bottom and wind forcing at the surface. The diurnal thermal stratification occurs in the deep region of the lake. A time lag is also found between suspended sediment concentration and wind speed. The study provides valuable storm-event data and mechanism analyses, which can improve our understanding of the lake circulation, wave dynamics, and their impact on sediment resuspension and vertical mixing in Lake Okeechobee. The data resulting from this study will be used to validate the Lake Okeechobee Environment Model which is used to predict the impacts of different management scenarios on lake activities.     

Depth-Averaged Two-Dimensional Numerical Modeling of Unsteady Flow and Nonuniform Sediment Transport in Open Channels

A depth-averaged two-dimensional (2D) numerical model for unsteady flow and nonuniform sediment transport in open channels is established using the finite volume method on a nonstaggered, curvilinear grid. The 2D shallow water equations are solved by the SIMPLE(C) algorithms with the Rhie and Chow’s momentum interpolation technique. The proposed sediment transport model adopts a nonequilibrium approach for nonuniform total-load sediment transport. The bed load and suspended load are calculated separately or jointly according to sediment transport mode. The sediment transport capacity is determined by four formulas which are capable of accounting for the hiding and exposure effects among different size classes. An empirical formula is proposed to consider the effects of the gravity on the sediment transport capacity and the bed-load movement direction in channels with steep slopes. Flow and sediment transport are simulated in a decoupled manner, but the sediment module adopts a coupling procedure for the computations of sediment transport, bed change, and bed material sorting. The model has been tested against several experimental and field cases, showing good agreement between the simulated results and measured data.     

Aspect Ratio to Maximize Sediment Transport in Rigid Bank Channels

The continuity equation, Manning’s equation, Einstein’s wall correction procedure and sediment transport equations are combined to indicate channel aspect ratios which maximize sediment transport for a given water discharge in rigid-bank trapezoidal and rectangular channels with fixed slope. Higher aspect ratios are required to maximize sediment transport for channels conveying bed load than for those with a dominant suspended load. A total load equation predicts optimum aspect ratios lying in between those for bed load and suspended load channels. The equations imply that the optimum aspect ratio increases markedly as the channel bank to channel bed roughness ratio increases. The resulting optimum ratios are smaller than the aspect ratios of many natural rivers.     

Hydrodynamic Forces Generated on a Spherical Sediment Particle during Entrainment

The objective of this research is to study the relationship between the coherent flow structures and the hydrodynamic forces leading to entrainment of a spherical bed sediment particle for a rough bed uniform turbulent flow. Two types of experiments, namely, movable and fixed balls, were conducted using spherical roughness-element beds with particle image velocimetry to measure the instantaneous flow-velocity field. Miniature piezoelectric pressure sensors were used to capture the instantaneous pressure on the surface of the sphere. Movable ball experiments reveal the predominance of large sweep structures at the instant of entrainment. Fixed ball experiments carried out at entrainment conditions show the importance of both vertical and horizontal pressure gradients on the ball leading to entrainment. Probability distribution function plots of pressures based on quadrant analysis of velocities also reveal the higher probability of occurrence of high magnitude force induced by sweep (Q4) events.     

Measurements of Sediment Erosion and Transport with the Adjustable Shear Stress Erosion and Transport Flume

Soil and sediments play an important role in water management and water quality. Issues such as water turbidity, associated contaminants, reservoir sedimentation, undesirable erosion and scour, and aquatic habitat are all linked to sediment properties and behaviors. In situ analysis is necessary to develop an understanding of the erosion and transport of sediments. Sandia National Laboratories has recently patented the Adjustable Shear Stress Erosion and Transport (ASSET) Flume that quantifies in situ erosion of a sediment core with depth while affording simultaneous examination of transport modes (bedload versus suspended load) of the eroded material. Core erosion rates and ratios of bedload to suspended load transport of quartz sediments were studied with the ASSET Flume. The erosion and transport of a fine-grained natural cohesive sediment were also observed. Experiments using quartz sands revealed that the ratio of suspended load to bedload sediment transport is a function of grain diameter and shear stress at the sediment surface. Data collected from the ASSET Flume were used to formulate a novel empirical relation for predicting the ratio of bedload to suspended load as a function of shear stress and grain diameter for noncohesive sediments.     

One-Dimensional Modeling of Dam-Break Flow over Movable Beds

A one-dimensional model has been established to simulate the fluvial processes under dam-break flow over movable beds. The hydrodynamic model adopts the generalized shallow water equations, which consider the effects of sediment transport and bed change on the flow. The sediment model computes the nonequilibrium transport of bed load and suspended load. The effects of sediment concentration on sediment settling and entrainment are considered in determining the sediment settling velocity and transport capacity. In particular, a correction factor is proposed to modify the Van Rijn formulas of equilibrium bed-load transport rate and near-bed suspended-load concentration for the simulation of sediment transport under high-shear flow conditions. The governing equations are solved by an explicit finite-volume method with the first-order upwind scheme for intercell fluxes. The model has been tested in two experimental cases, with fairly good agreement between simulations and measurements. The sensitivities of the model results to parameters such as the sediment nonequilibrium adaptation length, Manning’s roughness coefficient and the proposed correction factor have been verified. The proposed model has also been compared to an existing model and the results indicate the new model is more reliable.     

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.     

Total Load Transport Formula for Flow in Alluvial Channels

A user-friendly total bed-material load transport formula for flow in alluvial channels under equilibrium transport conditions has been developed based on dimensional analysis. The main advantages of this formula are its ease of computation, accuracy in prediction, and the wide range of application. The total sediment discharge gt is computed directly and is linearly related to the new total load transport parameter, TT. The latter involves variables that can be easily measured in field conditions, i.e., flow depth, mean flow velocity, energy slope, median sediment size and density, and water temperature. The factor of proportionality k in the formula has been checked for a wide range of hydraulic conditions and it remains a constant equal to 12.5. Comparisons between the computed and measured total sediment discharge indicate that the predictions are good.     

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.     

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.     

Pollutant Transport and Mixing Zone Simulation of Sediment Density Currents

Prediction of water column concentrations of suspended sediment is often necessary for environmental impact assessment of point source industrial discharges. For example, in “flow lane” or “open water” disposal, suction dredges discharge large volumes of suspended sediment into shallow water disposal locations. A sediment density current mixing model is presented here as part of the D-CORMIX expert system for hydrodynamic simulation of mixing zone behavior. This density current model extends the CORMIX decision support system to simulate continuous negatively buoyant discharges with or without suspended sediment loads on a sloping bottom with loss of suspended particles by sedimentation. Sedimentation is modeled using Stokes settling for five particle size classes. Density current width and depth, trajectory, total solids, tracer concentration, dilution, and particle size concentration are predicted. In addition, location and widths of sediment deposits, accretion rates, including particle size fractions within the spoils deposit, are predicted. The model results are in good overall agreement with available field and laboratory data.     

Two-Dimensional Total Sediment Load Model Equations

An unsteady total load equation is derived for use in depth-averaged sediment transport models. The equation does not require the load to be segregated a priori into bed and suspended but rather automatically switches to suspended load, bed load, or mixed load depending on a transport mode parameter consisting of local flow hydraulics. Further, the sediment transport velocity, developed from available data, is explicitly tracked, and makes the equation suitable for unsteady events of sediment movement. The equation can be applied to multiple size fractions and ensures smooth transition of sediment variables between bed load and suspended load for each size fraction. The new contributions of the current work are the consistent treatment of sediment concentration in the model equation and the empirical definition of parameters that ensure smooth transitions of sediment variables between suspended load and bed load.     

Dissolved Oxygen Demand at the Sediment-Water Interface of a Stream: Near-Bed Turbulence and Pore Water Flow Effects

A microbial dissolved oxygen (DO) uptake model was developed for a stream bed, including the effect of turbulence in the flow over the bed and pore water flow in the porous bed. The fine-grained sediment bed has hydraulic conductivities 0.01 ≤ k ≤ 1??cm/s, i.e., sediment particle diameter 0.006 ≤ ds ≤ 0.06??cm. The pore water flow is driven by pressure fluctuations at the sediment-water interface, mostly attributable to near-bed coherent motions in the turbulent boundary layer above the sediment bed. An effective mass transfer coefficient (De) coupled to a pore water flow model was used in the DO transport and DO uptake model. DO flux across the sediment-water interface and into the sediment, i.e., sedimentary oxygen demand (SOD), was related to hydraulic conductivity and microbial oxygen uptake rate in the sediment and shear velocity at the sediment-water interface. Simulated SOD values were validated against experimental data. For hydraulic conductivities of the sediment bed up to k ≈ 0.01??cm/s, the pore water flow effect on SOD was found negligible. Above this threshold, the effective mass (DO) transfer coefficient in the sediment bed (De) becomes larger as the hydraulic conductivity (k) becomes larger as the interstitial flow velocities increase; consequently, DO penetration depth increases with larger hydraulic conductivity of the sediment bed (k), and SOD increases as well. The enhancement of vertical DO transport into the sediment bed is strongest near the sediment-water interface, and rapidly diminishes with depth into the sediment layer. An increase in shear velocity at the sediment-water interface also enhances DO transfer. Shear velocity increases at the sediment-water interface will raise SOD regardless of the maximum oxidation rate if the hydraulic conductivity is above the threshold of k ≈ 1??cm/s. The relationship is nearly linear when U*<0.8??cm/s. At shear velocity U* = 1.6??cm/s, SOD for oxidation rates μ = 1000 and 2000??mg?l-1?d-1 are almost five times larger than those with no pore water flow. When pore water transport of DO is not limiting, SOD is a linear function of oxygen demand rate μ in the sediment when 0 ≤ μ ≤ 200??mg?l-1?d-1.     

Sand Transport on Steeply Sloping Plane and Rippled Beds

Experiments on sand transport have been carried out in the Sloping Sediment Duct at HR Wallingford. The aim of the experiments was to investigate sediment transport mechanisms, for sand of varying degree of grading, on sloping beds. The Sloping Sediment Duct is a steady flow, recirculating duct, capable of generating mean flow speeds of up to 1 m/s and tilting to +/?30°. Twenty-two tests with two different sediments were conducted. Both sediments had a median grain size of about 0.23 mm but different standard deviations. Bed slopes up to +/?20° were used in the experiments. The results show that bedforms have a significant effect on the transport rate. Since the bedforms, in turn, are affected significantly by the slope, the relation between transport rate and slope is not a monotonic function. Maximum suspended transport rates were attained for downslope flows at angles of about 10°. The transport rate for widely graded sediment was significantly larger than that for well-sorted sediment for almost all flows and slopes.     

Retention of Sediment and Nutrient Loads with Peak Runoff Control

The purpose of this research was to evaluate peak runoff control as a water protection method to reduce sediment and nutrient loads. Increased eutrophication of surface waters and risk of floods demands cost effective methods to reduce pollutant input and risks of flooding. With the peak runoff control it is possible to cut the main peaks and store the runoff water temporarily in ditches. The method decreases the suspended solids (SS) and nutrient loads by reducing flow velocities, and improving the settling of sediment particles. The method was tested in two heavily drained adjacent peat harvesting areas suffering considerable erosion. The peak flows were cut by 27–87%, the SS load by 61–94%, the total nitrogen (Ntot) by 45–91%, and the total phosphorus (Ptot) load by 47–88%. The peak runoff control method operated most effectively during extreme events when most of the SS load is transported. A detailed particle analysis of runoff water showed that water detention reduced the median particle size of SS load as the largest particles settle. The results clearly indicate that the peak runoff control is an effective method to control the sediment loads and peak flows from peatland drainage.