Numerical Simulation of Contraction Scour in an Open Laboratory Channel

A finite-volume computer code developed at the Institute for Hydromechanics, University of Karlsruhe, has been used to calculate the flow and sediment transport in a laboratory channel with constriction and movable bed. The flow is calculated by solving the fully three dimensional Reynolds-averaged Navier-Stokes equations with k?ε turbulence model. The bed deformation is obtained from an overall mass-balance equation for sediment transport and the bed-load transport is simulated with a nonequilibrium model. The calculated results for flow and scour development in the laboratory channel are compared with experimental measurements. The sensitivity of the simulated results to the nonequilibrian adaptation-length parameter in the nonequilibrium bed-load transport model is investigated systematically, which represents the main contribution of this paper.     

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.     

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.     

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.     

One-Dimensional Numerical Model for Nonuniform Sediment Transport under Unsteady Flows in Channel Networks

In this study, the proposed one-dimensional model simulates the nonequilibrium transport of nonuniform total load under unsteady flow conditions in dendritic channel networks with hydraulic structures. The equations of sediment transport, bed changes, and bed-material sorting are solved in a coupling procedure with a direct solution technique, while still decoupled from the flow model. This coupled model for sediment calculation is more stable and less likely to produce negative values for bed-material gradation than the traditional fully decoupled model. The sediment transport capacity is calculated by one of four formulas, which have taken into consideration the hiding and exposure mechanism of nonuniform sediment transport. The fluvial erosion at bank toes and the mass failure of banks are simulated to complement the modeling of bed morphological changes in channels. The tests in several cases show that the present model is capable of predicting sediment transport, bed changes, and bed-material sorting in various situations, with reasonable accuracy and reliability.     

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.     

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.     

Integrated Two-Dimensional Surface and Three-Dimensional Subsurface Contaminant Transport Model Considering Soil Erosion and Sorption

To investigate the complex hydrological, morphodynamic, and environmental processes in watersheds, a physically-based integrated two-dimensional (2D) surface and three-dimensional (3D) subsurface model for flow, soil erosion and transport, and contaminant transport in the surface-subsurface system is presented in this paper. The model simulates the rainfall-induced surface flow by solving the depth-averaged 2D diffusion wave equation and the variably-saturated subsurface flow by solving the 3D mixed-form Richards equation. The surface and subsurface flow equations are coupled using the continuity conditions of pressure and exchange flux at the ground surface. The model uses the concept of nonequilibrium in the depth-averaged 2D simulation of nonuniform total-load sediment transport in upland fields, considering detachments by rainsplash and hydraulic erosion driven by surface flow. The integrated 2D surface and 3D subsurface contaminant transport model takes into account the contaminant changes due to sediment sorption and desorption, as well as exchanges between surface and subsurface domains due to infiltration, diffusion, and bed change. The model applies the same set of surface equations of flow, sediment, and contaminant transport for describing both upland areas and streams, so that no special treatments are required at their interface. The established model has been evaluated by comparisons with published experimental, numerical, and analytical data and then applied in an agricultural watershed. The model is suitable for wetland areas and agricultural watersheds in which streams are not very narrow and deep, and meanwhile a relatively fine mesh that can distinguish the streams is preferred.     

Applicability of Sediment Transport Capacity Formulas to Dam-Break Flows over Movable Beds

In this paper, we investigate the extent to which well-known sediment transport capacity formulas can be used in one-dimensional (1D) numerical modeling of dam-break waves over movable beds. The 1D model considered here is a one-layer model based on the shallow-water equations, a bed update (Exner) equation, a space-lag equation for the nonequilibrium sediment transport and an empirical formula calculating the sediment transport capacity of the flow. The model incorporates a variety of sediment transport capacity formulas proposed by Meyer-Peter and Müller, Bagnold, Engelund and Hansen, Ackers and White, Smart and Jaeggi, van Rijn, Rickenmann, Cheng, Abrahams and Camenen, and Larson. We examine the performance of each formula by simulating four idealized laboratory cases on dam-break waves over sandy beds. Comparisons between numerical results and measurements show that for each case better predictions are obtained using a particular formula, but overall, formulas proposed by Meyer-Peter and Müller (with the factor 8 being replaced by 12), Smart and J?ggi, Cheng, Abrahams and Camenen, and Larson rank as the best predictors for the entire range of conditions studied here. Moreover, results show that in the cases where a bed step exists, implementing a mass failure mechanism in the numerical modeling plays an important role in reproducing the bed and water profiles.     

Experimental and Numerical Investigation of Bed-Load Transport under Unsteady Flows

The dynamic behavior of bed-load sediment transport under unsteady flow conditions is experimentally and numerically investigated. A series of experiments are conducted in a rectangular flume (18?m in length, 0.80?m in width) with various triangular and trapezoidal shaped hydrographs. The flume bed of 8?cm in height consists of scraped uniform small gravel of D50 = 4.8??mm. Analysis of the experimental results showed that bed-load transport rates followed the temporal variation of the triangular and trapezoidal hydrographs with a time lag on the average of 11 and 30?s, respectively. The experimental data were also qualitatively investigated employing the unsteady-flow parameter and total flow work index. The analysis results revealed that total yield increased exponentially with the total flow work. An original expression which is based on the net acceleration concept was proposed for the unsteadiness parameter. Analysis of the results then revealed that the total yield increased exponentially with the increase in the value of the proposed unsteadiness parameter. Further analysis of the experimental results revealed that total flow work has an inverse exponential variation relation with the lag time. A one-dimensional numerical model that employs the governing equations for the conservation of mass for water and sediment and the momentum was also developed to simulate the experimental results. The momentum equation was approximated by the diffusion wave approach, and the kinematic wave theory approach was employed to relate the bed sediment flux to the sediment concentration. The model successfully simulated measured sedimentographs. It predicted sediment yield, on the average, with errors of 7% and 15% of peak loads for the triangular and trapezoidal hydrograph experiments, respectively.     

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.     

Simulation of Flow and Mass Dispersion in Meandering Channels

This paper reports the development of an enhanced two-dimensional (2D) numerical model for the simulation of flow hydrodynamics and mass transport in meandering channels. The hydrodynamic model is based on the solution of the depth-averaged flow continuity and momentum equations where the density of flow varies with the concentration of transported mass. The governing equation for mass transport model is the depth-averaged convection and diffusion equation. The dispersion terms arisen from the integration of the product of the discrepancy between the mean and the actual vertical velocity distribution were included in the momentum equations to take into account the effect of secondary current. Two laboratory experimental cases, flow in mildly and sharply curved channels, were selected to test the hydrodynamic model. The comparison of the simulated velocity and water surface elevation with the measurements indicated that the inclusion of the dispersion terms has improved the simulation results. A laboratory experiment study of dye spreading in a sine-generated channel, in which dye was released at the inner bank, centerline, and outer bank, respectively, was chosen to verify the mass transport model. The simulated concentration field indicated that the Schmidt number can be used as a calibration parameter when dispersion is computed using a 2D approach with a simplified turbulence model.     

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.     

Exact Discontinuous Solutions of Exner’s Bed Evolution Model: Simple Theory for Sediment Bores

Determining the evolution of the bed of a river or channel due to the transport of sediment was first examined in a theoretical context by Exner in 1925. In his work, Exner presents a simplified bed evolution model derived from the conservation of fluid mass and an “erosion” equation that is commonly referred to as the sediment continuity or Exner equation. Given that Exner’s model takes the form of a nonlinear hyperbolic equation, one expects, depending on the given initial condition of the bed, the formation of discontinuities in the solution in finite time. The analytical solution provided by Exner for his model is the so-called classical or genuine solution of the initial-value problem, which is valid while the solution is continuous. In this paper, using the general theory of nonlinear hyperbolic equations, we consider generalized solutions of Exner’s classic bed evolution model thereby developing a simple theory for the formation and propagation of discontinuities in the bed or so-called sediment bores.     

Effect of Bed Armoring on Bed Topography of Channel Bends

The two-dimensional numerical model previously developed by the writers for modeling the bed variations in a channel bend with uniform sediment is upgraded to incorporate the nonuniformity of sediment particles as well as bed armoring. In this model, the two-dimensional, depth-averaged, unsteady flow equations along with the bed-load mass conservation equation are solved in a body-fitted coordinate system by using the Beam and Warming alternating-direction implicit (ADI) scheme. A one-dimensional bed surface armoring approach is extended herein for application to a two-dimensional domain. The model is applied to a 180° bend with a constant radius under unsteady flow conditions. Numerical simulations are carried out to study the effect of bed armoring on the bed deformations in channel bends. Results show that bed armoring reduces scour in channel bends.     

Effects of Lifting Force on Bed Topography and Bed-Surface Sediment Size in Channel Bend

In bed-load sediment transport, the lifting force plays an important role in reducing the friction between sediment particles and the bed surface, and it makes particle transportation by the shear force easier. Because the lifting force is related to vorticity, a three-dimensional (3D) numerical model incorporating large eddy simulations was applied to simulate the vorticity field in a channel bend. The results show that the distribution of vorticity is highly nonuniform, and it can lead to significant variations in lifting force and bed-load sediment transport per unit width in a channel bend. Relevant theories are modified on the basis of physical reasoning and then incorporated into numerical models to investigate the lifting-force effects on the bed topography and bed-surface sediment size gradation in a channel bend. With the lifting-force effects considered, it is shown that the errors in simulated bed topography can be reduced by approximately 40% and in bed-surface sediment size by 50%.     

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.