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.     

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

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

Modeling Noncohesive Suspended Sediment Transport in Stream Channels Using an Ensemble-Averaged Conservation Equation

The governing conservation equation for the transport of noncohesive suspended sediment in erodible channels is recognized as a stochastic partial differential equation due to the uncertainties in the parameters, and a deterministic ensemble-averaged equation is developed. Variables in this one-dimensional equation are represented as averaged quantities, and their covariances are also taken into account. Lateral inflows and deposition and entrainment of sediment are incorporated in the formulation. A hypothetical test problem is constructed to examine the model behavior. Manning’s coefficient, bed slope and bottom width are taken as the primary random parameters. Results from the solution of the ensemble-averaged equation are compared to results from Monte Carlo simulations. For comparison purposes, predicted values are also obtained by solving the deterministic transport equation without the covariance terms. It is found that predictions obtained from this latter approach deviate significantly from Monte Carlo simulation results. On the other hand, the ensemble-averaged predictions compare favorably to the Monte Carlo simulation results indicating that this promising technique needs further exploration.     

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.     

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.     

Case Study: Sediment Transport in Proposed Geomorphic Channel for Napa River

A model study evaluates sediment transport in a geomorphic channel proposed for restoration and flood damage reduction of an 11-km tidally influenced reach of the Napa River located in California. The model study employs the unsteady quasi-2D hydrodynamic and sediment transport model MIKE 11, the simplified marsh plain accretion model MARSH 98, and the Rouse equation to predict annual average morphological changes of the geomorphic channel. The adopted modeling approach allows for the simulation of salient sediment transport processes in a river estuary, including lateral and vertical sorting of sediments, and local flushing of fine cohesive and noncohesive sediments during flooding. Accretion rates, particularly within the marsh plain terrace of the multistage channel, are found to be within acceptable limits for project maintenance and ecosystem restoration purposes. This enhanced 1D modeling approach may offer a viable and cost-effective alternative compared to fully 2D and 3D models, with relatively less model set-up and run-time requirements.     

Simulating Sediment Transport in a Flume with Forced Pool-Riffle Morphology: Examinations of Two One-Dimensional Numerical Models

One-dimensional numerical sediment transport models (DREAM-1 and DREAM-2) are used to simulate seven experimental runs designed to examine sediment pulse dynamics in a physical model of forced pool-riffle morphology. Comparisons with measured data indicate that DREAM-1 and -2 closely reproduce the sediment transport flux and channel bed adjustments following the introduction of fine and coarse sediment pulses, respectively. The cumulative sediment transport at the flume exit in a DREAM-1 simulation is within 10% of the measured values, and cumulative sediment transport at flume exit in a DREAM-2 simulation is within a factor of 2 of the measured values. Comparison of simulated and measured reach-averaged aggradation and degradation indicates that 84% of DREAM-1 simulation results have errors less than 3.3?mm, which is approximately 77% of the bed material geometric mean grain size or 3.7% of the average water depth. A similar reach-averaged comparison indicates that 84% of DREAM-2 simulation results have errors less than 7.0?mm, which is approximately 1.7 times the bed material geometric mean grain size or 11% of the average water depth. Simulations using measured thalweg profiles as the input for the initial model profile produced results with larger errors and unrealistic aggradation and degradation patterns, demonstrating that one-dimensional numerical sediment transport models need to be applied on a reach-averaged basis.     

Numerical Modeling of Breach Erosion of River Embankments

The process of breach erosion of river embankments depends on the interaction among flow, sediment transport, and the corresponding morphological changes. Levees often consist of noncohesive material with a wide range of grain sizes. The dam material is mainly eroded due to the transport capacity of the overtopping water. Both bed load and suspended load are of importance. For breach formation, the lateral erosion due to slope instabilities has a significant impact. A depth averaged, two-dimensional numerical model was developed to account for these processes. The sensitivity of the discharge through the breach related to different processes and material parameters was investigated and compared to experimental and field data. The results show that the most sensitive parameter of an erosion-based dike-breach simulation is the breach side-slope angle which determines the lateral erosion. The application of the described Model 2dMb to different embankment failures at the Elbe River illustrates its capability in simulating overtopping breaching.     

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.     

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.     

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.     

Experiments on Deposition Behavior of Fine Sediment in a Reservoir

The deposition behavior of fine sediment is an important phenomenon, and yet it is unclear to engineers concerned about reservoir sedimentation. Laboratory experiments were conducted to produce both quasi-homogeneous flow and a turbidity current region divided by a plunge section. Silica powder (a noncohesive sediment) and kaolin (a cohesive sediment) were used as the suspended material. Because the effective gravitational force is the primary driving force for velocity in turbidity currents, the velocity profile was closely related to the concentration profile. The deposition rate of noncohesive coarser particles exponentially decays along the flow path. Most of the coarser particles were deposited in the quasi-homogeneous flow region or within a small distance downstream of the plunge section. The plunge did not carry those coarser particles further downstream. Deposition in the region of the turbidity current was found mainly by cohesive particles. Hydraulic sorting exists in the quasi-homogeneous flow region for noncohesive coarser particles, but becomes less significant in the downstream portion with deposition rates becoming mildly decayed. For fine cohesive particles, hydraulic sorting for the deposited gradation is not significant.     

Method for Estimating Sediment Transport in Ice-Covered Channels

A method is proposed for estimating rates of sediment transport in ice-covered alluvial channels. The method extends existing, open-water procedures for estimating rates of sediment transport to conditions of ice-covered flow. A key aspect of the method is the assessment of flow resistance attributable to bed-surface drag. That assessment is used to estimate rates of bed load and suspended load, and thereby total bed-sediment transport rate. Estimation of ice-covered suspended load additionally entails an approximation whereby open-water suspended load is scaled in proportion to the ratio of a reference sediment concentration for ice-covered flow relative to that for open-water flow. The reference concentration is calculated in terms of bed-load rate and shear velocity attributed to bed-surface drag. Flume data are used to develop the method and tentatively verify it. Field verification of the method presently is hampered by the absence of field data on bed sediment transport in ice-covered channels.     

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.     

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.     

Characterization of Model Uncertainty in the Static Pile Design Formula

Level 1 reliability methods have been internationally accepted as the basis for development of the new generation of geotechnical design codes. A key requirement of this design approach is the identification and quantification of uncertainties associated with the geotechnical design under consideration. This paper presents four load test databases from South Africa for driven piles in noncohesive soils (29 tests), bored piles in noncohesive soils (33 tests), driven piles in cohesive soils (59 tests), and bored piles in cohesive soils (53 tests). The capacity model factor is defined as the ratio of the interpreted capacity (Chin-Davisson approach) and the predicted capacity (static pile design formula). The uncertainty in the capacity model factor is modeled as a lognormal random variable. The model factor statistics reported in this study are required for reliability-based ultimate limit state design. The uncertainty in the load-settlement behavior is characterized by fitting measured load-settlement data to a hyperbolic equation and then normalizing the hyperbolic curve with the interpreted capacity. The resulting exercise reduces uncertainties in a set of nonlinear continuous curves to uncertainties in two hyperbolic curve-fitting parameters. This approach is practical and grounded realistically on the load test database with minimal assumptions. The hyperbolic parameter statistics reported in this study are required for reliability-based serviceability limit state design.