Journal of Hydraulic Engineering

Local clear-water scour tests were performed with three different diameter circular piles (0.114, 0.305, and $0.914m$), three different uniform cohesionless sediment diameters (0.22, 0.80, and $2.90\mathrm{mm}$) and a range of water depths and flow velocities. The tests were performed in the $6.1m$ wide, $6.4m$ deep, and $38.4m$ long flume at the United States Geological Survey Conte Research Center in Turners Falls, Mass. These tests extend local scour data obtained in controlled experiments to prototype size piles and ratios of pile diameter to sediment diameter to 4,155. Supply water for this flow through flume was supplied by a hydroelectric power plant reservoir and the concentration of suspended fine sediment (wash load) could not be controlled. Equilibrium scour depths were found to depend on the wash load concentration.

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.

An upwind conservative scheme with a weighted average water-surface-gradient approach is proposed to compute one-dimensional open channel flows. The numerical scheme is based on the control volume method. The intercell flux is computed by the one-sided upwind method. The water surface gradient is evaluated by the weighted average of both upwind and downwind gradients. The scheme is tested with various examples, including dam-break problems in channels with rectangular and triangular cross-sections, hydraulic jump, partial dam-break problem, overtopping flow, a steady flow over bump with hydraulic jump, and a dam-break flood case in a natural river valley. Comparisons between numerical and exact solutions or experimental data demonstrated that the proposed scheme is capable of accurately reproducing various open channel flows, including subcritical, supercritical, and transcritical flows. The scheme is inherently robust, stable, and monotone. The scheme does not require any special treatment, such as artificial viscosity or front tracking technique, to capture steep gradients or discontinuities in the solution.

Mechanism of sediment transport is composed of complicated interactions between turbulent flow, particle motion, and bed configurations. Of particular significance is the interaction between turbulence and particle motion, although turbulence measurements of particle-laden two phase flow have been a problem for a long time, especially in the near-wall region. In this study, simultaneous measurements of both the particles and fluid (water) were conducted in particle-laden two phase open channel flows by means of a discriminator particle-tracking velocimetry. The mean velocity and turbulence characteristics for fluid and particles each were examined in comparison with those in clear-water (particle-free) flow, together with previous existing data measured by laser Doppler anemometer and phase Doppler anemometer. The relative velocity and the turbulence modulation, which are the most important topics in two phase-flow approach, were revealed by varying the particle diameter and specific density. The fluid-sweeps are more contributory to the motion of particles than the fluid ejections in the near-wall region. In turn, the particle-sweeps transport the high momentum to the carrier fluid and enhance the turbulence intensities of fluid.

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

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.

This paper presents a rapid forecast model for simulating hyperconcentrated sediment-laden floods in the Lower Yellow River. The model is a hybrid of a conventional one-dimensional mathematical model for unsteady sediment-laden flow and an artificial neural networks model for encapsulation of numerical results. The former provides detailed river flood routing information under typical scenarios, whereas the latter extracts modeling outputs from the former and establishes a station-specific model for efficient flood forecasting. Three typical floods that occurred in the Lower Yellow River in 1977, 1982, and 1996 are simulated. Not only the hybrid model predictions are found to be in close agreement with measured data, but also the computational speed is significantly enhanced. It is found that sediment transport is of significance with regard to the flooding behavior of hyperconcentrated flows. Therefore, the model presented herein is of particular use for rivers with high sediment concentration.

The analysis of water distribution networks has been and will continue to be a core component of civil engineering water resources curricula. Since its introduction in 1936, the Hardy Cross method has been used in virtually every water resources engineering text to introduce students to network analysis. The technique gained widespread popularity primarily because it is amenable to manual calculation techniques. However, the same subtle elegance that facilitates manual calculations often obscures the primary engineering and physical principles of water distribution systems relative to the nuances of algorithm implementation. Herein, the authors illustrate the application of commonly available spreadsheet software (MicroSoft Excel) to more concisely and effectively solve typical undergraduate network distribution problems using linear theory. Application development is much more efficient and straightforward than the corresponding Hardy Cross implementation enabling students to concentrate upon the engineering system and relevant design issues. The technique presented utilizes commonly available technology and is presented as a supplement to alternatives discussed in recent literature.

Several models for the vertical distribution of the double-averaged (in time and in the plane parallel to the mean bed) longitudinal velocity in the flow region between roughness troughs and roughness tops are suggested. We found that the same model for velocity distribution may be applicable to a range of flow conditions and roughness types, which share some common features. The suggested models for velocity distribution in the near-bed region are: (1) Constant velocity; (2) exponential velocity distribution; and (3) linear velocity distribution. The measured velocity distributions may be approximated by a single model or by a combination of models depending on roughness geometry and flow conditions. The validity of these models for velocity distribution is supported by laboratory data.

Improved relationships for the flocculation factor and floc settling velocity are proposed. They differ from the existing formulas by accounting for the effect of floc porosity on the hydrodynamic drag and for the fractal properties of the floc projection on the plane. These relationships predict that the fractal dimension $d$ influences the floc settling velocity $W_s$ only if $d⩾2$. A flocculation effect on optical backscatter is also considered, and a fractal-based relationship for quantifying this effect is proposed. This relationship serves as a basis for a simple optical backscatter sensor (OBS) based technique for estimating some fractal-related parameters of sediment flocs involved in the proposed relationships for the settling velocity and flocculation factor. The technique is tested using data from three field sites. The potential effect of particle size on backscatter efficiency is highlighted, and an approach for accounting for this effect is suggested.