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.     

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 Simulation of Relatively Wide, Shallow Channels with Erodible Banks

A two-dimensional numerical model was developed to simulate relatively wide, shallow rivers with an erodible bed and banks composed of well-sorted, sandy materials. A moving boundary-fitted coordinate system was used to calculate water flow, bed change, and bank erosion. The cubic interpolated pseudoparticle method was used to calculate flow, which introduced little numerical diffusion. The sediment-transport equation for the streamline and transverse transport was used to estimate bed and bank evolution over time, while considering the secondary flow. Bank erosion was simulated when the gradient in the cross-sectional direction of the banks was steeper than the submerged angle of repose because of bed erosion near the banks. The numerical model reproduced the features of central bars well, such as bar growth, channel widening due to divergence of the flow around the bars, scour holes at the lee of the bars, and the increase of bar size with time. These features were in accordance with the observations for laboratory experiments. It also reproduced the features of braided rivers, such as the generation of new channels and abandonment of old channels, the bifurcation and confluence of channels, and the lateral migration of the channels. The model showed that the sediment discharge rate fluctuated with time, one of the dynamic features observed in braided channels.     

Development and Application of a Phosphorus Model for a Shallow Oxbow Lake

A three-dimensional numerical model was developed for simulating the phosphorus concentration in shallow lakes. In this model, the computational domain was divided into two parts: the water column and the bed sediment layer. The processes of mineralization, settling, adsorption, desorption, bed release (diffusion), growth, and death of phytoplankton were taken into account, and the concentration of organic phosphorus, phosphate, and related water quality constituents was simulated. The concentrations of adsorbed (particulate) and dissolved phosphate due to adsorption-desorption were calculated using two formulas derived based on the Langmuir equation. The release rate of phosphorus from the bed sediment layer was calculated by considering the effects of the concentration gradient across the water-sediment interface, pH, temperature, dissolved oxygen concentration, and flow conditions. The adsorption and desorption of phosphate from sediment particles, as well as its release from bed sediment, were verified using laboratory experimental data. The model was calibrated and applied to Deep Hollow Lake in the Mississippi alluvial plain. The simulated trends and magnitudes of phosphorus concentration were compared with field observations. The simulation results show that there are strong interactions between sediment-related processes and phosphorus concentration.     

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.     

Influence of Vertical Resolution and Nonequilibrium Model on Three-Dimensional Calculations of Flow and Sediment Transport

This note, using a three-dimensional model of river flow and sediment transport, examines the effect of the vertical resolution and the choice a nonequilibrium adaptation length Ls in predicting flow and sediment transport around groins in China’s Yongding River. The results show that a fine vertical grid and nonequilibrium sediment transport model provide good predictions, especially on the river bed profile with an obvious main channel and flood plain.     

Numerical Modeling of Abutment Scour with the Focus on the Incipient Motion on Sloping Beds

A three-dimensional computational fluid dynamics model is applied to predict local scour around an abutment in a rectangular laboratory flume. When modeling local scour, steep bed slopes up to the angle of repose occur. To predict the depth and the shape of the local scour correctly, the reduction of the critical shear stress due to the sloping bed must be taken into account. The focus of this study is to investigate different formulas for the threshold of noncohesive sediment motion on sloping beds. Some formulas only take the transversal angle (perpendicular to the flow direction) into account, but others also consider the longitudinal angle (streamwise direction). The numerical model solves the transient Reynolds-averaged Navier-Stokes equations in all three dimensions to compute the water flow. Sediment continuity in combination with an empirical formula is used to capture the bed load transport and the resulting bed changes. When the sloping bed exceeds the angle of repose, the bed slope is corrected with a sand-slide algorithm. The results from the numerical simulations are compared with data from physical experiments. The reduction of the bed shear stress on the sloping bed improves the results of the numerical simulation distinctly. The best results are obtained with the formulas that use both the transversal and the longitudinal angle for the reduction of the critical bed shear stress.     

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.     

Flow Heterogeneity over 3D Cluster Microform: Laboratory and Numerical Investigation

The present study examines the flow around a self-occurring cluster bed form and the use of general computation fluid dynamics methods for hydraulic and geophysical flow applications. This is accomplished through a comprehensive experimental/numerical investigation. In the laboratory, cluster bed forms are first formed from movable sediment, and laser Doppler velocimeter measurements of two-dimensional fluid velocity are then taken around a formed cluster. A three-dimensional (3D) Reynolds averaged Navier-Stokes simulation of the physical cluster and flow conditions is then conducted using near-wall, shear stress transport (SST) turbulence modeling with the inclusion of hydraulic roughness, ks (R = 31,150, ks/h = 0.1, ks+ = 274, i.e., in the fully rough regime). SST near-wall modeling is advantageous compared to the more widely used wall functions approach for flows with significant roughness and flow separation because the model equations can be integrated down to the wall. Therefore, SST near-wall modeling makes no a priori assumption that the law of the wall is valid throughout the wall region of the flow. Additionally, it has the ability to intrinsically handle boundary roughness through the boundary condition for turbulent specific dissipation at the wall, allowing for wall functions to be bypassed in accounting for roughness effects. The study shows that in the wall region surrounding the cluster, flow is 3D and quite complex, with different scales of embedded flow structures dominating the cluster wake and leading to flow heterogeneities in pressure and bed-shear stress. Results also indicate that near-wall modeling with SST compared favorably with the experimental flow data without tuning of model constants.     

Case Study: Bed Stability around the East Caisson of the Tacoma Narrows Bridge

A study on the hydraulic and sediment conditions at the Tacoma Narrows Bridge, in Washington State, was carried out to examine the stability of the bed material around the bridge caissons. Specifically, this was conducted around the east caisson, where the maximum velocities around either of the two caissons are experienced. This was performed for the peak tidal exchange event of May 27th to 28th, 2002. During this max flow event, multibeam surveyed bathymetry and three-dimensional acoustic doppler current profiler velocity data were collected around the east caisson in the course of both the flood and ebb. The surface of the bed material surrounding the east caisson was videotaped during the slack conditions following the yearly maximum flow event, and used to determine the particle size distribution and spatial arrangement of those distributions around the caisson. This was done by lowering a submersible video camera and appropriate lighting to the bottom of the Narrows, a depth of approximately 45?m. Flow and sediment observations were coupled to determine the commencement of sediment motion for different size classes of sediment. Two methods were utilized to calculate friction velocity in order to assess the stability of different bed particle size fractions during these high flow conditions. Friction velocity was first calculated from measured velocity profiles at various locations around the east caisson. The second method was based on the concept of a free stream power-law expression for depth-averaged velocity. Stability was then examined using the critical shear stress concept and captured video data of the bed. General results showed the particles ? 30?mm in diameter were in motion during the flood and ebb. The work is here presented as a case study because of the unique large-scale flow conditions that are present around the east caisson of the Tacoma Narrows Bridge.     

Three-Dimensional CFD Modeling of Self-Forming Meandering Channel

A three-dimensional CFD model was used to compute the formation of the meandering pattern in an initially straight alluvial channel. The numerical model was based on the finite volume method using an unstructured grid with dominantly hexahedral cells. The k-ε model was used to predict turbulence and the SIMPLE method was used to compute the pressure. The sediment transport was computed as bed load in addition to solving the convection-diffusion equation for suspended sediment transport. The bed changes were calculated and the grid was altered during the computation as channel erosion and deposition caused wetting and drying. The model was tested by comparing with results from physical model studies carried out at Colorado State Univ., Fort Collins, Colo. The results showed successfully the replication of many of the meander characteristics, including secondary currents, cross-sectional profiles, meander planform, meander wavelength, downstream meander migration, and chute formation.     

Simulation of Scour Process in Plunging Pool of Loose Bed-Material

The scouring process in a plunge pool of loose bed with uniform bed-materials due to a two-dimensional plane impinging jet was simulated computationally. The finite-element-based unsteady three-dimensional model, CCHE3D, with k-ε turbulence closure was employed to solve the flow field. It has long been recognized that the unsteady behavior of the turbulent jet fluctuation plays an important role in scouring and transporting sediment in the plunge pool. In order to model this phenomenon realistically, one has to consider the effects of both shear stress and the life force on sediment particles due to pressure fluctuation. The latter has been taken into account by using empirical relationships of flume data. Both of these effects have been incorporated in the nonequilibrium sediment transport model consisting of sediment pickup rate and step length adopted for the jet scour problem. The model constant relating to the fluctuating lift force was calibrated using an empirical equation to predict the quasi-equilibrium scour depth. The results simulated by the model proposed here agree reasonably well with experimental data.     

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.     

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 Modeling of Three-Dimensional Flow Field Around Circular Piers

A three-dimensional numerical model FLUENT is used to simulate the separated turbulent flow around vertical circular piers in clear water. Computations are performed using different turbulence models and results are compared with several sets of experimental data available in the literature. Despite commonly perceived weakness of the k-ε model in resolving three-dimensional (3D) open channel and geophysical flows, several variants of this turbulence model are found to have performed satisfactorily in reproducing the measured velocity profiles. However, model results obtained using the k-ε models show some discrepancy with the measured bed shear stress. The Reynolds stress model performed quite well in simulating velocity distribution on flat bed and scour hole as well as shear stress distribution on flat bed around circular piers. The study demonstrates that a robust 3D hydrodynamic model can effectively supplement experimental studies in understanding the complex flow field and the scour initiation process around piers of various size, shape, and dimension.     

Numerical Modeling of River Flow for Ecohydraulic Applications: Some Experiences with Velocity Characterization in Field and Simulated Data

Information regarding the spatial and temporal organization of river flow is required for many applications in river management, and is a fundamental requirement in ecohydraulics. As an alternative to detailed field surveys and to mesohabitat reconnaissance schemes, potential exists to deploy numerical flow simulation as an assessment and design tool. A key question is the extent to which complex hydrodynamic models are really practical in river management applications. This paper presents experiences using sediment simulation in intakes with multiblock, a three-dimensional modeling code, in conjunction with a statistical approach for classifying the spatiotemporal dynamics of flow behavior. Even in a simple configuration, the model is able to replicate well flow structures which associate with the mesohabitat concepts used in field reconnaissance techniques. The model also captures spatiotemporal dynamics in flow and depth behavior at these scales. However, because the model shows differential performance between flow stages and between differing channel (bed form) units, the smaller-scale and discharge-dependent dynamics of some zones within the channel may be less-well represented, and the implications of this for future research are noted.     

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.     

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.     

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.