Journal of Hydraulic Engineering

The potential for hydraulic transients in the West Area Combined Sewer Overflow Tunnel System of the City of Atlanta was investigated by numerical modeling. Different design alternatives were evaluated to mitigate pressure oscillations, backflows, and flooding during storm events. A dynamic transient model was developed to simulate the transition from gravity flow to pressurized flow that would occur during a storm event. The model is able to simulate pressurization with pipe-filling bore or gradual-flow regime transition. Numerical results compare well with experimental data on single pipe flows. Following a brief description of the model, the paper discusses the pressure transients simulated during the filling of the tunnels for both the original and improved designs.

Three analytic stage-discharge formulas, suitable for hand calculation, have been derived for prismatic open channels. The three types of prismatic geometry investigated are: simple rectangular, symmetric, and asymmetric rectangular compound channels. The formulas include three key parameters that govern the local friction factor, eddy viscosity, and secondary flow. The discharge results and the corresponding analytic depth-averaged velocity and bed shear stress distributions show good agreement with the experimental data, even when constant parameter values are assumed. The influence of these three parameters on the stage-discharge relationships is explored.

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-\omega$ turbulence model or the standard $k-\epsilon$ 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-\omega$ model better predicts flow and sediment transport through channel bends than the standard $k-\epsilon$ model.

Dividing flows in open channels are commonly encountered in hydraulic engineering systems. They are inherently three-dimensional (3D) in character. Past experimental studies were mostly limited to the collection of test data on the assumption that the flow was 1D or 2D. In the present experimental study, the flow is treated as 3D and test results are obtained for the flow characteristics of dividing flows in a $90°$ , sharp-edged, rectangular open-channel junction formed by channels of equal width. Depth measurements are made using point gauges, while velocity measurements are obtained using a Dantec laser Doppler anemometer over grids defined throughout the junction region. A 3D turbulence model is also developed to investigate the dividing open-channel flow characteristics. The predicted flow characteristics are validated using experimental data. Following proper model validation, the numerical model developed can yield design data pertaining to flow characteristics for different discharge and area ratios for other dividing flow configurations encountered in engineering practice. Energy and momentum coefficients based on the present 3D model yield more realistic energy losses and momentum transfers for dividing flow configurations. Data related to secondary flows provide information vital to bank stability, if the branch channel sides are erodible.

Using a depth-averaged model to obtain the velocity and pressure distributions in the vertical direction is difficult. A multilayer model is an option that can be used to improve on the depth-averaged model. However, the unknown flow depth needs to be predicted first and then divided into layers as an input for the multilayer model. An improved multilayer model is proposed here by introducing an implicit layer dividing interfaces that are associated with the flow velocity and pressure distribution. The formulation of interfaces also applies to boundary faces: Free surface and channel beds. Therefore, each flow layer behaves like that in the classical depth-averaged model. Subsequently the governing equations are also simplified due to the vanishing terms related to interfacial flow exchanges. This improved model has been satisfactorily applied to steady flow simulations in three cases: Flow over a slope transition from mild to steep, from steep to mild, and over a trapezoidal weir. The results demonstrate the efficiency and validity of the proposed models to simulate open channel flows with bed slope changes.

The transport and mixing of contaminants in conduits is governed by advection, dispersion, and decay. Several models are available to trace the transport of such constituents and most assume that the principal mechanisms for transport are advection and reaction only. However in pipes where low velocities prevail, longitudinal dispersion is significant and models that neglect the dispersion effects fail to properly simulate the observed concentrations in low velocity pipes. This work presents a method for simulating the advection-dispersion-reaction process of constituent transport in water networks. A Eulerian–Lagrangian method is employed whereby the dispersion term in the governing equation is approximated using finite differences and the resulting first-order partial differential equation is then integrated using the method of characteristics. Analytical solutions of the transport equation are also derived to quantify the effect of neglecting dispersion at pipe junctions and to assess the accuracy of the proposed method. The Eulerian-Lagrangian method is tested on benchmark networks and on the field study at the Cherry Hill/Brushy Plains network. Results show that the model developed is capable of simulating transport with equal accuracy for low and high velocity flows with and without significant dispersion effects. It also performs better than other models because of the nonuniform grid distribution and the interpolation schemes used.

Attempts to include river bank erosion predictors within morphological models are becoming increasingly common, but uncertainty surrounds the procedures used to couple bed deformation and bank erosion submodels in a way that maintains the mass continuity of eroded bank sediment. Herein we present a coupling procedure that comprises two discrete elements. First, immediately following bank failure, the slumped debris comes to rest at the bank toe as a planar surface inclined at an “angle of repose.” Second, we use a fractional transport model to simulate the subsequent erosion and deposition of the failed bank material debris. The method is demonstrated with an example in which an equilibrium bed topography model is combined with a river bank erosion model to predict the morphological response of a river bend to a large flow event.

A numerical solution to the finite difference of two-dimensional (2D) depth-averaged equations on nonstaggered grid points is proposed in this technical note. Following a locally one-dimensional procedure, the basic equations are split into a pair of one-dimensional equations. Therefore, the solution of a 2D problem is reduced to the solution of a sequence of two one-dimensional problems. The discretization of the split one-dimensional equations is obtained with the use of the original Preissmann operator. Using Fourier’s classic linear analysis, stability, dissipation and dispersion with frictional resistance are investigated for the variations of the Courant number and weighting time factor.

On the basis of theoretical considerations, the technical note presents two practical formulas for the dimensioning of air valves when filling a pipe with water. One is to be used for designing air valves on the basis of the maximum allowed water hammer overpressures; the other when the maximum in pipe water velocity is set. The reliability of these formulas was tested with a numerical model based on the same hypothesis, which was in turn verified with experimental tests.