Numerical Model for Sediment Transport and Bed Degradation in the Yangtze River Channel Downstream of Three Gorges Reservoir

This paper presents a two-dimensional morphological model for unsteady flow and both suspended-load and bed-load transport of multiple grain size to simulate transport of graded sediments downstream from the Three Gorges Reservoir. The model system includes a hydrodynamic module and a sediment module. The hydrodynamic module is based on the depth-averaged shallow water equations in orthogonal curvilinear coordinates. The sediment module describing nonuniform sediment transport is developed to include nonequilibrium transport processes, bed deformation, and bed material sorting. The model was calibrated using field observations through application to a 63-km-long alluvial river channel on the middle Yangtze River in China. A total of 16 size groups and a loose layer method of three sublayers were considered for the transport of the nonuniform bed materials in a long-term simulation. Predictions are compared with preliminary results of field observations and factors affecting the reliability of the simulated results are discussed. The results may be helpful to the development of more accurate simulation models in the future.     

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.     

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.     

Predicting the Morphodynamic Response of Silt-Laden Rivers to Water and Sediment Release from Reservoirs: Lower Yellow River, China

The high sediment load of the Yellow River results in rapid infilling of its reservoirs when sediment is not regularly flushed. Simultaneously, the downstream reaches of the Yellow River experience extremely high siltation rates, which are reduced when sediment is retained in its reservoirs. To minimize siltation in the reservoirs and the downstream river bed, water and sediment are released from the reservoir in a controlled way through flushing experiments. In this paper, we analyze the effect of such a flushing event on the downstream river bed through data analysis and numerical modeling. Sedimentation may be minimized by relating the amount of sediment released from the reservoir to the sediment available for release through operational monitoring and by releasing relatively clear water after turbid water. Despite this flushing of sediment, the reservoir will eventually fill up, and more sediment released again into the lower Yellow River. The change in discharge magnitude and frequency brought about by the reservoir will then probably result in increased siltation rates in the lower Yellow River compared to the predam situation.     

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.     

Multiple Linear Regression Model for Total Bed Material Load Prediction

A new total bed material load equation that is applicable for rivers in Malaysia was developed using multiple linear regression analyses. A total of 346 hydraulic and sediment data were collected from nine natural and channelized rivers having diverse catchment characteristics in Malaysia. The governing parameters were carefully selected based on literature survey and field experiments, examined and grouped into five categories namely mobility, transport, sediment, shape, and flow resistance parameters. The most influential parameters from each group were selected by using all possible regression model method. The suitable model selection criteria namely the R-square, adjusted R-square, mean square error, and Mallow’s Cp statistics were employed. The accuracy of the derived model is determined using the discrepancy ratio, which is a ratio of the calculated values to the measured values. The best performing models that give the highest percentage of prediction from the validation data were chosen. In general, the newly derived model is best suited for rivers with uniform sediment size distribution with a d50 value within the range of 0.37–4.0 mm and performs better than the commonly used Graf, Yang, and Ackers–White total bed material load equations.     

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.     

Oxygen Demand by a Sediment Bed of Finite Length

A model of sedimentary oxygen demand (SOD) for a sediment bed of finite length is presented. The responses of diffusive oxygen transfer in turbulent flow above the sediment surface and of microbial activity inside the sediment to a developing diffusive boundary layer are modeled numerically. The developing diffusive boundary layer above the sediment/water interface is modeled based on shear velocity and turbulent boundary layer concepts, and dissolved oxygen (DO) uptake inside the sediment is modeled as a function of the microbial growth rate. The model predicts that the diffusive boundary layer above the sediment/water interface thickens in flow direction, and that DO penetration depth into the sediment is practically constant over the length of the sediment bed. The effect of the developing diffusive boundary layer on SOD is minor, except at very low shear/flow velocities (shear velocity U*<0.01?cm/s) and/or high microbial density inside the sediment. The average SOD over the sediment bed therefore varies only slightly with its length. SOD varies somewhat in flow direction, i.e., SOD is largest near the leading edge (x = 0), decreases with distance, and finally, approaches a nearly constant value for fully developed boundary layer. Including microbial activity in the sediment makes the change of SOD in flow direction much smaller than is predicted by a pure vertical diffusive flux model. The diffusive boundary layer is nearly fully developed at a dimensionless distance x+ = 10,000, regardless of microbial activity inside the sediment. Longer sediment beds are required to eliminate the small leading edge effect on any measured average SOD value. SOD depends strongly on the diffusion coefficient of DO inside the sediment bed. This effect becomes more significant as shear/flow velocity is increased. Overall, SOD is found to be controlled principally by shear velocity of the water flowing above the sediment/water interface, microbial activity inside the sediment, and diffusion of DO inside the sediment. The length of the sediment bed is of lesser influence.     

Erosion of Finite Thickness Sediment Beds by Single and Multiple Circular Jets

Sediment management in reservoirs with the help of water jets has motivated this work. Erosion caused by single and multiple submerged circular turbulent wall jets on a noncohesive sediment bed of finite thickness lying on a fixed boundary was studied with the help of laboratory experiments. Different combinations of jet diameter, jet separation, and sediment thickness to jet diameter ratio were tested. Results show a relation between dimensionless parameters characterizing the steady state bed profile and the densimetric particle Froude number F0 given by the velocity at the nozzle and the effective diameter and submerged specific density of the sediment. Evolution of scour with time confirms previous studies where the erosion was found to initially grow with the logarithm of time up to a certain reference time t*. This time, made dimensionless with a time scale tc involving the volume of sediment scoured and the rate of erosion, was also related to the densimetric Froude number. A comparison with studies regarding erosion of a semiinfinite layer of sediment is also presented.     

Modeling the Erosion of Mixtures of Organic and Granular In-Sewer Sediments

The release of fine-grained organic sediments from sediment deposits can have a detrimental impact on water quality in a number of situations. This paper examines the release of such sediment in the context of the erosion of mixed organic/granular sediment in-sewer deposits. In the European Union, sewer flow quality modeling software uses equations derived from uniform granular sediment studies. Actual sewer sediments are mixtures of organic and granular material and interactions between these fractions may account for the poor performance of current models. Laboratory experiments were carried out using surrogate sewer sediment mixtures. Impaction of the bed surface by saltating granular particles increased the erosion of fine-grained organic sediments. Changes in the composition of the bed surface over the duration of a test resulted in change in the availability of fine-grained sediment. A model that attempted to simulate these mechanisms, using an empirically based correction factor to account for the impaction mechanism and an active bed layer to account for changes in the bed surface composition was developed. The limited success of the simulations indicated that such simple modeling approaches may not be appropriate for organic/granular deposits in which grain sorting occurs.     

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.     

Three-Dimensional Modeling of Sediment Transport in a Narrow 90° Channel Bend

A three-dimensional numerical model was used for calculating the velocity and bed level changes over time in a 90° bended channel. The numerical model solved the Reynolds-averaged Navier-Stokes equations in three dimensions to compute the water flow and used the finite-volume method as the discretization scheme. The k-ε model predicted the turbulence, and the SIMPLE method computed the pressure. The suspended sediment transport was calculated by solving the convection diffusion equation and the bed load transport quantity was determined with an empirical formula. The model was enhanced with relations for the movement of sediment particles on steep side slopes in river bends. Located on a transversally sloping bed, a sediment particle has a lower critical shear stress than on a flat bed. Also, the direction of its movement deviates from the direction of the shear stress near the bed. These phenomenona are considered to play an important role in the morphodynamic process in sharp channel bends. The calculated velocities as well as the bed changes over time were compared with data from a physical model study and good agreement was found.     

Multiple Time Scales of Fluvial Processes with Bed Load Sediment and Implications for Mathematical Modeling

Fluvial bed load transport is often considered to assume a capacity regime exclusively determined by local flow conditions, but its applicability in naturally occurring unsteady flows remains to be theoretically justified. In addition, mathematical river models are often decoupled, being based on simplified conservation equations and ignoring the feedback impacts of bed deformation to a certain extent. So far whether the decoupling could have considerable impacts on the fluvial processes with bed load transport remains poorly understood. This paper presents a theoretical investigation of both issues. The multiple time scales of fluvial processes with bed load sediment are evaluated to examine the applicability of bed load transport capacity and decoupled models. Numerical case studies involving active bed load transport by highly unsteady flows complement the analysis of the time scales. It is found that bed load transport can sufficiently rapidly adapt to capacity in line with local flow because sediment exchange with the bed overwhelms the advection of bed load sediment by the mean flow. The present work provides theoretical justification of the concept of bed load transport capacity in most circumstances, which is underpinned by existing observations of bed load transport by flash floods. For fluvial processes with bed load transport, the feedback impacts of bed deformation are limited; therefore, decoupled modeling is, in this sense, appropriate.     

Periodic Diffusional Mass Transfer near Sediment/Water Interface: Theory

Time-variable (periodic) flow over a lake bed, and the associated boundary layer development, have the potential to control or at least influence rates of mass transfer across the sediment/water interface. An analysis for instantaneous and time averaged flux of a material across the sediment/water interface for infinite supply in the water and infinite sink in the sediment is presented. The water flow above the interface is characterized by the shear velocity (U?) which is a periodic function of time with a maximum amplitude of (U?0) as may be typical of an internal seiche (internal standing wave) motion in a density stratified lake. The relationship between the shear velocity on the lake bed and the wind shear on the lake surface is illustrated for an extremely simplified two-layered lake of constant depth. For a less restrictive analysis, shear velocities on a lake bed have to be obtained either from field measurements or from a three-dimensional lake circulation model driven by atmospheric forcing including wind. Smaller and wind-sheltered lakes will have lower (U?0) and periodicities (T). The response of the diffusive boundary layer was related to the period of the periodic motion (T), Schmidt number (Sc), and shear velocity (U?). The vertical diffusive flux at the sediment/water interface was expressed by a Sherwood number (Sh), either instantaneous or time averaged. The mean Sherwood number (Shave) varies with shear velocity of the wave motion over the sediment bed, Schmidt number (Sc) and the period (T) due to the response of the diffusive boundary layer to the time variable water velocity. Effective diffusive boundary layers develop only at low shear velocities. Where they do, maximum and minimum boundary layer thickness depends on all three independent variables (T, Sc, and U?0). The diffusive boundary layer strongly affects sediment/water mass transfer, i.e., Sherwood numbers. Mass transfer averaged over a period can be substantially less than that produced by steady-state flow at the same U?0 and Sc. At Sc = 500, typical for dissolved oxygen, the mass transfer ratio can be reduced to 60% of steady state, depending on the internal wave period (T).     

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.     

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.     

Bed-Material Load Computations for Nonuniform Sediments

The nonuniformity of bed material affects the bed-material load calculations. A size gradation correction factor Kd is developed to account for the lognormal distribution of bed material. The use of Kd in conjunction with bed-material load equations originally developed for single particle sizes improves the accuracy of transport calculations for sediment mixtures. This method is applicable to laboratory flumes and natural rivers with median diameter d50 of bed material in the sand size ranges. The improvement on transport rate by Kd factor is significant for data with standard deviation σg of bed material greater than 2, while the correction is negligible for data with σg less than 1.5. Sediment in transport also follows a lognormal distribution with a median diameter d50t generally finer than the corresponding d50. As the size gradation increases, d50t becomes much finer than the corresponding value of d50. The relationship between d50t and d50 is defined as a function of σg and agrees well with field data. The previously recommended use of d35 as representative size of the bed material is found not to be generally applicable.     

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.     

Three-Dimensional Modeling of Nonuniform Sediment Transport in an S-Shaped Channel

A three-dimensional numerical model was applied to compute uniform and nonuniform sediment transport and bed deformation in an S-shaped laboratory channel located at the University of Innsbruck, where detailed measurements of the velocity field and bed elevation changes were made. The channel had two bends, a trapezoidal cross section, and a slope of S = 0.005. Gravel with a mean diameter of 4.2?mm was used as movable bed material and for sediment feeding. Wu’s formula for multiple grain sizes was compared with van Rijn’s formula using one grain size. Fairly good agreement was found between the computed and measured bed elevations for both approaches, whereas Wu’s formula could further improve the numerical results. Looking at the physics of the erosion pattern, the computed scour areas were located slightly more downstream than what was observed in the physical model. The current study also includes several parameter tests: grid distribution in vertical, lateral, and longitudinal direction; time step; number of inner iterations/time step; active sediment layer thickness; and the Shields coefficient. The variation of those parameters gave some differences in the results, but the overall pattern of bed elevation changes remained the same.