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.     

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.     

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.     

3D Numerical Modeling of Flow and Sediment Transport in Open Channels

A 3D numerical model for calculating flow and sediment transport in open channels is presented. The flow is calculated by solving the full Reynolds-averaged Navier-Stokes equations with the k ? ε turbulence model. Special free-surface and roughness treatments are introduced for open-channel flow; in particular the water level is determined from a 2D Poisson equation derived from 2D depth-averaged momentum equations. Suspended-load transport is simulated through the general convection-diffusion equation with an empirical settling-velocity term. This equation and the flow equations are solved numerically with a finite-volume method on an adaptive, nonstaggered grid. Bed-load transport is simulated with a nonequilibrium method and the bed deformation is obtained from an overall mass-balance equation. The suspended-load model is tested for channel flow situations with net entrainment from a loose bed and with net deposition, and the full 3D total-load model is validated by calculating the flow and sediment transport in a 180° channel bend with movable bed. In all cases, the agreement with measurements is generally good.     

Numerical Modeling of Bed Evolution in Channel Bends

A two-dimensional numerical model is developed to predict the time variation of bed deformation in alluvial channel bends. In this model, the depth-averaged unsteady water flow equations along with the sediment continuity equation are solved by using the Beam and Warming alternating-direction implicit scheme. Unlike the present models based on Cartesian or cylindrical coordinate systems and steady flow equations, a body-fitted coordinate system and unsteady flow equations are used so that unsteady effects and natural channels may be modeled accurately. The effective stresses associated with the flow equations are modeled by using a constant eddy-viscosity approach. This study is restricted to beds of uniform particles, i.e., armoring and grain-sorting effects are neglected. To verify the model, the computed results are compared with the data measured in 140° and 180° curved laboratory flumes with straight reaches up- and downstream of the bend. The model predictions agree better with the measured data than those obtained by previous numerical models. The model is used to investigate the process of evolution and stability of bed deformation in circular bends.     

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.     

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.     

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.     

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.     

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.     

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.     

Modeling Bed Changes in Meandering Rivers Using Triangular Finite Elements

A two-dimensional depth-averaged model was used for the simulation of scour and deposition in sand-bed meandering channels with fixed banks. The model employs unstructured meshes based on triangular elements and incorporates the effects of curvature-induced helical flow and transverse bed slope in the direction of bed-load sediment transport. The model was tested using experimental data from a well-known laboratory curved channel and a full scale meandering river. The numerical results agreed well with observed data, demonstrating that the model can reproduce the main features of bed profiles along meandering rivers, such as the formation of point bars and pools.     

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.     

Two-Dimensional Depth-Averaged Flow Modeling with an Unstructured Hybrid Mesh

An unstructured hybrid mesh numerical method is developed to simulate open channel flows. The method is applicable to arbitrarily shaped mesh cells and offers a framework to unify many mesh topologies into a single formulation. A finite-volume discretization is applied to the two-dimensional depth-averaged equations such that mass conservation is satisfied both locally and globally. An automatic wetting-drying procedure is incorporated in conjunction with a segregated solution procedure that chooses the water surface elevation as the main variable. The method is applicable to both steady and unsteady flows and covers the entire flow range: subcritical, transcritical, and supercritical. The proposed numerical method is well suited to natural river flows with a combination of main channels, side channels, bars, floodplains, and in-stream structures. Technical details of the method are presented, verification studies are performed using a number of simple flows, and a practical natural river is modeled to illustrate issues of calibration and validation.     

Depth-Averaged Model of Open-Channel Flows over an Arbitrary 3D Surface and Its Applications to Analysis of Water Surface Profile

A new set of depth-averaged equations is introduced to study the flow over an arbitrary three-dimensional (3D) surface. These equations are derived based on a generalized curvilinear coordinate system attached to the 3D bed surface, therefore it allows us to include the effect of centrifugal force due to the bottom curvature. These general equations make it possible to analyze flows over complex terrain without the limitation of mild slope assumption used in conventional depth-averaged models. This new model is then applied to calculate the water surface profiles of (1) flow over a cylindrical surface; (2) flow over a circular surface; and (3) flow with an air-core vortex at a vertical intake. A simple hydraulic experiment is conducted in the laboratory to observe the water surface profile of flow over a circular surface. The results obtained from the model are in good agreement with experimental measurements and calculation by an empirical formula. Consequently, it demonstrates the applicability of the model in cases of flow over a highly curved bottom.     

Stream Temperature Dynamics in Upland Agricultural Watersheds

A numerical model to compute the free-surface flow hydrodynamics and stream temperature dynamics by solving the depth-averaged, 1D unsteady flow and heat transport equations is presented. The hydrodynamics model considers the effects of arbitrary stream geometry, variable slopes, variable flow regimes, and unsteady boundary conditions. The thermal transport model accounts for the effects of solar radiation, air temperature, relative humidity, cloud cover, wind speed, heat conduction between water and streambed, subsurface flow, and shading by riparian vegetation. The model is verified with measurements in a stream in an upland agricultural watershed located in Indiana. Diurnal variations in the streamflow and stream temperatures are highly transient. The proposed model predicted well the streamflow and stream temperatures that were measured every 15 min over 25 days. The results of this study demonstrate that the solar (shortwave) radiation and subsurface inflow are the most significant contributors to the stream heat budget.     

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.     

Evolution of Scour Depth at Circular Bridge Piers

Experiments of bridge pier scour are carried out under steady and unsteady clear-water scour conditions with uniform and nonuniform sediments. Around the pier nose, the sediment size variation of surface bed materials is investigated, and a regressed formula is obtained for estimating the mixing layer thickness in terms of median sediment size and geometric standard deviation of grain size distribution. A method based on the mixing layer concept is developed for calculating the equilibrium scour depth in nonuniform sediment. Based on the experimental data of scour rate, a model simulating the scour-depth evolution under steady flow in nonuniform sediment is presented. By analyzing experimental data, a scheme is proposed for computing the scour-depth evolution under unsteady 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.