Journal of Hydraulic Engineering

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.

Predicting the effect of a pollutant spill on downstream water quality is primarily dependent on the water velocity, longitudinal mixing, and chemical/physical reactions. Of these, velocity is the most important and difficult to predict. This paper provides guidance on extrapolating travel-time information from one within bank discharge to another. In many cases, a time series of discharge (such as provided by a U.S. Geological Survey stream gauge) will provide an excellent basis for this extrapolation. Otherwise, the accuracy of a travel-time extrapolation based on a resistance equation can be greatly improved by assuming the total flow area is composed of two parts, an active and an inactive area. For 60 reaches of 12 rivers with slopes greater than about 0.0002, travel times could be predicted to within about 10% by computing the active flow area using the Manning equation with n = 0.035 and assuming a constant inactive area for each reach. The predicted travel times were not very sensitive to the assumed values of bed slope or channel width.

An analytical method is developed to determine the longitudinal dispersion coefficient in Fischer's triple integral expression for natural rivers. The method is based on the hydraulic geometry relationship for stable rivers and on the assumption that the uniform-flow formula is valid for local depth-averaged variables. For straight alluvial rivers, a new transverse profile equation for channel shape and local flow depth is derived and then the lateral distribution of the deviation of the local velocity from the cross-sectionally averaged value is determined. The suggested expression for the transverse mixing coefficient equation and the direct integration of Fischer's triple integral are employed to determine a new theoretical equation for the longitudinal dispersion coefficient. By comparing with 73 sets of field data and the equations proposed by other investigators, it is shown that the derived equation containing the improved transverse mixing coefficient predicts the longitudinal dispersion coefficient of natural rivers more accurately.

The exchange of dissolved matter between a groyne field and a main stream influences the transport and distribution of a pollutant cloud in a river. In forecasting models, groyne fields are represented as dead zones with effective properties like exchange coefficients and exchanging volume. Despite its relevance for such practical applications, little research has been done on the exchange process between a groyne field and the main stream itself. Therefore, this study is aimed at examining this exchange process and validating the dead-zone prediction model, which treats the exchange process as a first order system. A schematized physical model of a river with groynes was built in a laboratory flume. The exchange process was visualized quantitatively with dye in adjacent groyne fields. In order to couple the exchange process to the velocity field, particle tracking velocimetry measurements were performed. Two different types of exchange were observed. First, exchange takes place via the mixing layer that is formed at the river-groyne-field interface. The large eddies formed in the mixing layer are the major cause of this exchange. Second, under certain conditions, even larger eddies are shed from the upstream groyne tip. Distortions in the flow field caused by such intermittent structures cause a much larger exchange than that by the mixing layer alone. The occurrence of large shed eddies depends on the presence of a sufficiently large, stationary, secondary gyre located at the upstream corner of the groyne field. The overall exchange of matter could be characterized as a first-order process, in accordance with the dead-zone-theory. The corresponding exchange coefficients agreed reasonably well with the results of earlier experiments and the effective coefficients as found in experiments in real river flows.

Under certain tidal conditions, a saline underflow originating in the Pacific Ocean moves into Lake Ogawara, Japan. The underflow consists of a uniform saline bottom layer that is slightly warmer than the ambient and an interfacial shear layer in which the velocity and density are decreased. Within the experimental area the underflow is confined to a channel approximately 1 km wide and is essentially two-dimensional. The underflow had a bulk Richardson number, defined in terms of the mean properties, between 1 and 2. The rate of entrainment into the bottom layer was calculated using two distinct methods. The first method used the change in the maximum salinity of the underflow measured at two stations along the path of the underflow to infer the amount of ambient water entrained. The second method made direct measurements of vertical mass fluxes with a profiler. The agreement between the two methods was excellent. The measured entrainment coefficients were consistent with the derived entrainment law. The turbulent structure of the flow was mapped for a 3 h quasisteady period of the flow. Turbulence is predominantly generated on the bottom boundary and is transported vertically to the density interface, where it leads to mixing.

Three-dimensional finite-difference or finite-volume models of sinuous open channels (e.g., narrow rivers, estuaries, and reservoirs) generally require boundary-fitted grids and curvilinear flow solution. Cartesian models with square grid cells are simpler to apply, but require a larger number of cells, as the cell size is determined by cross-stream resolution. This paper presents a simplified curvilinear approach suitable for systems where the along-stream length scale is larger than the cross-stream scale. The curvilinear Navier-Stokes equations are manipulated so the left-hand side is identical to the Cartesian momentum equations. The right-hand side then consists of grid-stretching curvature terms. These terms are written as functions of a perturbation parameter, so the first-order curvilinear effects are obtained with the lowest-order perturbation terms. As the Cartesian equations' form is preserved, we can readily adapt a Cartesian model to this perturbation curvilinear approach by adding the small curvilinear terms as explicit momentum sources.

A 1D numerical model of the downslope flow and deposition of muddy subaerial and subaqueous debris flows is presented. The model incorporates the Herschel-Bulkley and bilinear rheologies of viscoplastic fluid. The more familiar Bingham model is integrated into the Herschel-Bulkley rheological model. The conservation equations of mass and momentum of single-phase laminar debris flow are layer-integrated using the slender flow approximation. They are then expressed in a Lagrangian framework and solved numerically using an explicit finite difference scheme. Starting from a given initial shape, a debris flow is allowed to collapse and propagate over a specified topography. Comparison between the model predictions and laboratory experiments shows reasonable agreement. The model is used to study the effect of the ambient fluid density, initial shape of the failed mass, and rheological model on the simulated propagation of the front and runout characteristics of muddy debris flows. It is found that initial failure shape influences the front velocity but has little bearing on the final deposit shape. In the Bingham model, the excess of shear stress above the yield strength is proportional to the strain rate to the first power. This exponent is free to vary in the Herschel-Bulkley model. When it is set at a value lower than unity, the resulting final deposits are thicker and shorter than in the case of the Bingham rheology. The final deposit resulting from the bilinear model is longer and thinner than that from the Bingham model due to the fact that the debris flow is allowed to act as a Newtonian fluid at low shear rate in the bilinear model.

Local scour at circular piers founded on clay was studied experimentally in the laboratory to compare the depth of scour in sand and in clay and to investigate the effects of the Reynolds number, Froude number, and approach flow depth on scour depth. The depths of scour in front, at the side, and in the back of the piers were measured as a function of time under steady, gradually varied flow conditions. The measured scour-depth-versus-time curves were fitted with a hyperbola to estimate the equilibrium scour depths. The extrapolated equilibrium scour depths were compared with values predicted by the Federal Highway Administration equation. The results showed that although the rates of scour were much slower in clay than in sand, equilibrium scour in clay was about the same as in sand. It was found that the shape of the scour hole correlates with the pier Reynolds number. At low Reynolds numbers, the depth of scour was about the same all around the piers. At higher Reynolds numbers, the scour holes developed mainly behind the piers with much less scour in front of the piers. It was also found that the extrapolated equilibrium scour depth correlates well with the pier Reynolds number and that the Froude number and relative water depth did not have a significant effect on the scour depth for these experimental conditions.