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.     

Experimental Study of Bed Load Transport through Emergent Vegetation

Vegetation is an important agent in fluvial geomorphology and sedimentary processes, through its influence on the local hydraulics that determine sediment transport. Within stands of emergent vegetation, bed shear is substantially reduced through the absorption of momentum by drag on the stems. This stimulates deposition of sediment and reduces capacity for bed load transport. The effect of emergent vegetation on hydraulic parameters (including equilibrium bed gradient, flow depth, and velocity) and on bed load transport rate has been investigated experimentally for one sediment size, stem diameter, and stem spacing. Bed load transport rate was found to be closely related to bed-shear stress, which must be estimated by partitioning total flow resistance between stem drag and bed shear.     

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.     

Two-Dimensional Nonequilibrium Noncohesive and Cohesive Sediment Transport Model

The purpose of this paper is to develop an unsteady 2D depth-averaged model for nonuniform sediment transport in alluvial channels. In this model, the orthogonal curvilinear coordinate system is adopted; the transport mechanisms of cohesive and noncohesive sediment are both embedded; the suspended load and bed load are treated separately. In addition, the processes of hydraulic sorting, armoring, and bed consolidation are also included in the model. The implicit two-step split-operator approach is used to solve the flow governing equations and the coupling approach with iterative method are used to solve the mass-conservation equation of suspended sediment, mass-conservation equation of active-layer sediment, and global mass-conservation equation for bed sediment simultaneously. Three sets of data, including suspension transport, degradation and aggradation cases for noncohesive sediment, and aggradation, degradation, and consolidation cases for cohesive sediment, have been demonstrated to show the rationality and accuracy of the model. Finally, the model is applied to evaluate the desilting efficiency for Ah Gong Diann Reservoir located in Taiwan to show its applicability.     

Equilibrium Near-Bed Concentration of Suspended Sediment

A new approach is presented for calculating the equilibrium near-bed concentration of suspended sediment in an alluvial channel flow. It is formulated from the balance between bed sediment entrainment and suspended sediment deposition across the near-bed boundary. The entrainment flux is determined making use of a turbulent bursting outer-scale-based function and the flux of deposition by the product of near-bed concentration and hindered settling velocity of sediment. A number of flume data records in the literature are analyzed to calibrate and verify the present approach. The observed near-bed concentrations for the data records are obtained by first isolating the suspended load transport rate from the observed total load transport rate using Engelund and Fredsoe's bed-load formula and then equating the suspended load transport rate to the shape integration of Dyer and Soulsby. The present approach is shown to perform satisfactorily compared to the results of data analysis. It is found that the near-bed concentration is evidently dependent on sediment particle size in addition to the Shields parameter due to skin friction. This finding seems to challenge previous relationships that simply represent the near-bed concentration as empirical functions of the purely skin-friction-related Shields parameter.     

Modeling Sediment Transport Using Depth-Averaged and Moment Equations

A 1D mathematical model to calculate bed variations in alluvial channels is presented. The model is based on the depth-averaged and moment equations for unsteady flow and sediment transport in open channels. Particularly, the moment equation for suspended sediment transport is originally derived by the assumption of a simple vertical distribution for suspended sediment concentration. By introducing sediment-carrying capacity, suspended sediment concentration can be solved directly from sediment transport and its moment equations. Differential equations are then solved by using the control-volume formulation, which has been proven to have good convergence. Numerical experiments are performed to test the sensitivity of the calibrated coefficients α and k in the modeling of the bed deposition and erosion. Finally, the computed results are compared with available experimental data obtained in laboratory flumes. Comparisons of this model with HEC-6 and other numerical models are also presented. Good agreement is found in the comparisons.     

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.     

Numerical Modeling of Bed Deformation in Laboratory Channels

A depth-average model using a finite-volume method with boundary-fitted grids has been developed to calculate bed deformation in alluvial channels. The model system consists of an unsteady hydrodynamic module, a sediment transport module and a bed-deformation module. The hydrodynamic module is based on the two-dimensional shallow water equations. The sediment transport module is comprised of semiempirical models of suspended load and nonequilibrium bedload. The bed-deformation module is based on the mass balance for sediment. The secondary flow transport effects are taken into account by adjusting the dimensionless diffusivity coefficient in the depth-average version of the k–ε turbulence model. A quasi-three-dimensional flow approach is used to simulate the effect of secondary flows due to channel curvature on bed-load transport. The effects of bed slope on the rate and direction of bed-load transport are also taken into account. The developed model has been validated by computing the scour hole and the deposition dune produced by a jet discharged into a shallow pool with movable bed. Two further applications of the model are presented in which the bed deformation is calculated in curved alluvial channels under steady- and unsteady-flow conditions. The predictions are compared with data from laboratory measurements. Generally good agreement is obtained.     

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.     

3D Numerical Modeling of Flow and Sediment Transport in Laboratory Channel Bends

The development of a fully three-dimensional finite volume morphodynamic model, for simulating fluid and sediment transport in curved open channels with rigid walls, is described. For flow field simulation, the Reynolds-averaged Navier–Stokes equations are solved numerically, without reliance on the assumption of hydrostatic pressure distribution, in a curvilinear nonorthogonal coordinate system. Turbulence closure is provided by either a low-Reynolds number k?ω turbulence model or the standard k?ε turbulence model, both of which apply a Boussinesq eddy viscosity. The sediment concentration distribution is obtained using the convection-diffusion equation and the sediment continuity equation is applied to calculate channel bed evolution, based on consideration of both bed load and suspended sediment load. The governing equations are solved in a collocated grid system. Experimental data obtained from a laboratory study of flow in an S-shaped channel are utilized to check the accuracy of the model’s hydrodynamic computations. Also, data from a different laboratory study, of equilibrium bed morphology associated with flow through 90° and 135° channel bends, are used to validate the model’s simulated bed evolution. The numerically-modeled fluid and sediment transportation show generally good agreement with the measured data. The calculated results with both turbulence models show that the low-Reynolds k?ω model better predicts flow and sediment transport through channel bends than the standard k?ε model.     

3D Numerical Simulation for Water Flows and Sediment Deposition in Dam Areas of the Three Gorges Project

This paper presents a three-dimensional (3D) mathematical model for suspended load transport in turbulent flows. Based on the stochastic theory of turbulent flow proposed by Dou, numerical schemes of Reynolds stresses for anisotropic turbulent flows are obtained. Instead of a logarithmic law, a specific wall function is used to describe the velocity profile close to wall boundaries. The equations for two-dimensional suspended load motion and sorting of bed material have been improved for a 3D case. Numerical results are in good agreement with the measured data of the Gezhouba Project. The present method has been employed to simulate sediment erosion and deposition in the vicinity of the Three Gorges Dam. The size distribution of the deposits and bed material, and flow and sediment concentration at different times and elevations, are predicted. The results agree well with the observations in physical experiments. Thus, a new method is established for 3D simulation of sediment motion in the vicinity of dams.     

Calculation of Bed Changes in Mountain Streams

Simulation of flow and sediment transport in mountain streams is complicated by the presence of high gradients, abrupt changes in geometry, variations in regime of flow, and large roughness elements. Most of the numerical models to predict aggradation and degradation in alluvial channels have been developed for low-gradient rivers. This paper is devoted to the development of a numerical model to calculate bed elevation and grain size distribution changes in mountain streams where the maximum bed material size is in the range of boulders. An attempt is made to validate the model by using observed field data collected upstream from a small retention dam in a Venezuelan stream. After calibration of the sediment transport equation, reasonable agreement is obtained for the variations in the grain size distribution of the bed-surface material. An additional application is presented in the Cocorotico River, a small mountain stream located in the northwest region of Venezuela, which illustrates the adaptability of the model to handle a case of coarsest-bed-material removal from the active channel and to simulate the armoring process.     

Unified View of Sediment Transport by Currents and Waves. III: Graded Beds

The problem of suspended load and bed load transport in river and coastal flows over graded beds is addressed. Two effects are important: the degree of exposure of the sediment particles of unequal size within a mixture (hiding of smaller particles resting or moving between the larger particles) and the nonlinear dependence of transport on particle diameter. The former effect can be modeled by modifying the critical bed-shear stress through a correction factor and by modifying the effective grain roughness through another correction factor. The modeling of the effective bed-shear stress parameter is studied by using various alternative methods. Based on comparison with suspended load and bed load transport data for graded beds in steady and oscillatory flow, the most promising method is selected. The proposed prediction method is found to work well for the fine sand bed range as well as the coarse sand-gravel bed range.     

Sediment-Laden Flow in Open-Channels under Noncapacity and Capacity Conditions

To determine the suspended load component of sediment transport in open-channel flow, the vertical distribution of the concentration of suspended particles is of importance. It is usual to determine this distribution by solving the diffusion-convection equation under appropriate boundary conditions. The exponent in the resulting equation is the Rouse number, defined as z′ = vss/βu*. The -value has been the subject of much research. In natural alluvial channels the sediment-laden flow is usually in capacity (saturation) condition, implying that the flow will charge (saturate) itself with particles available in the bed load and∕or on the bed itself. However, simulation of sediment-laden flow in a laboratory flume is achieved typically by externally adding particles to the flow. Consequently, it is not certain that the flow was in capacity condition. The resulting -values are often values for noncapacity conditions. They should not be used for natural alluvial channels, because they are misleading. Reported herein are experiments performed in the laboratory under noncapacity, as well as under capacity, conditions. This study focuses on the experimental determination of the -value, which incorporates the ratio of the sediment flux and the momentum flux , as well as the velocity and concentration profiles, ?ū∕?y and ?/?y, respectively. For experiments with small particles, d50 = 0.135 mm, the -values at capacity condition are smaller than unity; at noncapacity condition, the -values are usually larger than the ones at capacity condition, but all are still smaller than unity.     

Longitudinal Dispersion in Ice-Covered Rivers

Longitudinal dispersion is the spreading of suspended or dissolved substances caused by the combined action of differential advection and lateral mixing. Tracer tests in many rivers under open-water conditions have indicated that the temporal spread first increases linearly with downstream distance (linear range). It may ultimately grow in proportion to the square root of distance (Fickian range), very far below the point of release. Test data have been obtained in ice-covered rivers, starting in 1975. They show linear-range behavior even though the test reaches were as long as hundreds of kilometers. As in the case of open-water conditions, the rates of spread under an ice cover are related to the friction factor of the flow, being generally higher than the open-water values for equal friction factors. Based on early results in pipe flow, many water quality models assume Fickian dispersion, contrary to the field evidence. A more correct quantification of the dispersion process is outlined and implications to the method of computation are discussed.     

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.     

Hydraulic Interpretation of Cross-Stream Variations in Bed-Load Transport

The hydraulic control of bed-load transport rates in Nahal Yatir and Nahal Eshtemoa, two coarse-grained ephemeral channels in the semiarid northern Negev, Israel, provides a rare opportunity to infer the spanwise variation in bed-shear stress from an analysis of cross-stream variations in bed-load transport rate. Automatic sediment transport monitoring stations were used to obtain synchronous measurements of bed-load discharge at a number of locations across the widths of two straight channel reaches. In both streams, channel-average bed-load fluxes demonstrated a common and well-defined response to changing channel-average shear stress and approximated the transporting capacity of the flow over much of the range of monitored discharges. However, transport rates measured at the channel margins are only half those at the channel centerline, and, at high discharges, a marked asymmetry in the pattern of bed-load transport develops across the central section of the widest channel. This variation in bed-load discharge over the two channel cross sections is thought to reflect lateral variations in shear stress induced by sidewall drag and, more tentatively, the generation and disposition of cellular secondary currents. But no systematic relation is found for the ratios of sediment fluxes at off-center sampling locations and those recorded at the channel center, even though the off-center locations are thought to move into and out of the region affected by sidewall drag as aspect ratio of the flow decreases and increases with changing water-stage. The results suggest that it is difficult to generalize about the changing influence of the sidewall on local shear and bed load as aspect ratio changes during the course of a flood.     

Effect of Sand Movement on a Cohesive Substrate

Flume experiments investigated the effect of mobile sand on the erosion of cohesive beds. The fluid-induced stress alone was not enough to cause erosion, and sand motion as bed load was needed. Erosion rates and suspended sediment concentration were found to increase with increasing sand transport and to decrease with increasing median grain size. The erosion rate was found to be at a maximum during saltation, intermediate during creep, and lowest during suspension.