Treatment of Natural Geometry in Finite Volume River Flow Computations

A method is proposed for the treatment of irregular bathymetry in one-dimensional finite volume computations of open-channel flow. The strategy adopted is based on a reformulation of the Saint-Venant equations. In contrast with the usual treatment of topography effects as source terms, the method accounts for slope and nonprismaticity by modifying the momentum flux. This makes it possible to precisely balance the hydrostatic pressure contributions associated with variations in valley geometry. The characteristic method is applied to the revised equations, yielding topographic corrections to the numerical fluxes of an upwind scheme. Further adaptations endow the scheme with an ability to capture transcritical sections and wetting fronts in channels of abrupt topography. To test the approach, the scheme is first applied to idealized benchmark problems. The method is then used to route a severe flood through a complex river system: the Tanshui in Northern Taiwan. Computational results compare favorably with gauge records. Discrepancies in water stage represent no more than a fraction of the magnitude of typical bathymetry variations.     

Modeling of One-Dimensional, Unsteady, Free-Surface, and Pressurized Flows

The utilization of mathematical models in hydraulic engineering for the analysis of one-dimensional, unsteady free-surface, and pressurized flows is discussed, with an emphasis when the models performed well as well as when they did not. For illustration purposes, the applications to a number of real-life projects are presented, outlining limitations, successes, and failures.     

Dynamic Wave Study of Flow in Tidal Channel System of San Juan River

In this work the complete equations of one-dimensional unsteady flow in open channels in integral form, and compatibility equations at the junctions of a channel network, are solved numerically. Analytical integration in space is used between each pair of consecutive irregular sections of a channel, and the nonprismatic term is expressed in terms of uncoupled functions of the geometry at the sections. The linearized system of equations for each time interval is solved by an elimination method based on a double-sweep algorithm. The model is applied to the estuary of the San Juan River in Venezuela, where oscillating currents by effect of semidiurnal tides take place and the amplitude of the wave at the mouth is amplified toward the inland direction. Alternating drying and filling is simulated by means of slight modifications in the bed geometry of upper river sections. Measured water elevation and flow rates available at two stations are used to calibrate the model, and a very accurate adjustment of the tidal levels observed in the river is obtained.     

Conservative Scheme for Numerical Modeling of Flow in Natural Geometry

A numerical model is proposed to compute one-dimensional open channel flows in natural streams involving steep, nonrectangular, and nonprismatic channels and including subcritical, supercritical, and transcritical flows. The Saint-Venant equations, written in a conservative form, are solved by employing a predictor-corrector finite volume method. A recently proposed reformulation of the source terms related to the channel topography allows the mass and momentum fluxes to be precisely balanced. Conceptually and algorithmically simple, the present model requires neither the solution of the Riemann problem at each cell interface nor any special additional correction to capture discontinuities in the solution such as artificial viscosity or shock-capturing techniques. The resulting scheme has been extensively tested under steady and unsteady flow conditions by reproducing various open channel geometries, both ideal and real, with nonuniform grids and without any interpolation of topographic survey data. The proposed model provides a versatile, stable, and robust tool for simulating transcritical sections and conserving mass.     

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.     

One-Dimensional Numerical Model for Nonuniform Sediment Transport under Unsteady Flows in Channel Networks

In this study, the proposed one-dimensional model simulates the nonequilibrium transport of nonuniform total load under unsteady flow conditions in dendritic channel networks with hydraulic structures. The equations of sediment transport, bed changes, and bed-material sorting are solved in a coupling procedure with a direct solution technique, while still decoupled from the flow model. This coupled model for sediment calculation is more stable and less likely to produce negative values for bed-material gradation than the traditional fully decoupled model. The sediment transport capacity is calculated by one of four formulas, which have taken into consideration the hiding and exposure mechanism of nonuniform sediment transport. The fluvial erosion at bank toes and the mass failure of banks are simulated to complement the modeling of bed morphological changes in channels. The tests in several cases show that the present model is capable of predicting sediment transport, bed changes, and bed-material sorting in various situations, with reasonable accuracy and reliability.     

URANS Computations of Shallow Grid Turbulence

This paper describes the unsteady Reynolds-averaged Navier–Stokes (URANS) computations of a quasi-two-dimensional (2D) grid turbulence in shallow open-channel flows, generated downstream of multiple piers aligned at regular intervals over the channel width. In shallow open-channel flows, the vertical confinement of the flow generally suppresses the three dimensionality and attains two-dimensional features with up-cascading of turbulent kinetic energy from small-scale toward large-scale structures. In this study, 2D depth averaged and 3D Reynolds-averaged equations with linear and nonlinear URANS turbulence models are applied to a shallow open-channel flow downstream of multiple piers and numerical results are discussed through a comparison with the experimental results performed by Uijttewaal and Jirka in 2003. We employed 0-equation models and k-ε models for the 2D and 3D computations, respectively. In 2D computations, vortices downstream of the grid occurred synchronously in the computation with both the linear and nonlinear 0-equation models. In the 3D computations, vortex merging and up-cascading of the kinetic energy were captured when artificial disturbance is added at the inlet. The measured decay of the turbulent kinetic energy in the streamwise direction, with a slope of ?1.3, was well captured by computation with the 3D models with inlet disturbance. The flow sensitivity on the inlet disturbance was rather small in the wide range of the disturbance ratios.     

Scale Effects in Flume Experiments on Flow around a Spur Dike in Flatbed Channel

This paper presents the findings from a series of flume experiments conducted to determine scale effects in small-scale models of flow around a single spur dike (wing-dam, groyne, or abutment) placed in a channel whose bed is fixed and flat. The flow features of primary interest are flow-thalweg alignment (line of maximum streamwise velocity) around a dike, and area extent of the flow-separation region (wake) immediately downstream of the dike. Those features are of practical concern in the deployment of dikes for various channel control purposes. The scale effects of concern herein are those attributable to use small length scales together with a bed-shear stress parameter (e.g., Shields parameter) as the primary criterion for dynamic similitude. Small-scale models, especially micromodels, often are used for investigating channel-control issues. Also, the shear-stress criterion is commonly used for models of flow in loose-bed channels, whereas Froude number commonly is the primary similitude criterion for models of fixed-bed open-channel flows. The experiments show that use of a shear-stress parameter as the primary criterion for dynamic similitude influences the flow thalweg and flow separation region at a dike. It does so by distorting pressure gradients around the model dike and by affecting turbulence generated by the dike. It also is shown that, for a range of small models, thalweg alignment and extent of separation region do not scale with model length scales. These findings are important for interpreting results from small hydraulic models, especially micromodels, of flow in loose-bed channels.     

Modeling Unsteady Open-Channel Flow for Controller Design

An approximated linear model of unsteady open-channel flow is necessary to design the water-level controller for irrigation open channels. Toward this end, this paper presents the matrix approach to derive the linear model of open-channel system in analytical form mainly according to the Saint Venant equations and the backwater profile at the steady state of open channel. The hydraulic model of the check structure at the downstream end of open channel is also incorporated into the linear model. A practical example indicates that the frequency response of the open-channel system can be accurately analyzed with the linear model. The simulation results in the time domain show that the dynamic behavior of the linear model approximates to that of the nonlinear model of the open-channel system. Finally, the limitations of the linear model are discussed.     

Modeling Unsteady Flow Characteristics of Hydropeaking Operations and Their Implications on Fish Habitat

Reservoir releases associated with energy production and flood mitigation need to be reconciled with efforts to maintain healthy ecosystems in regulated rivers. Unsteady flow phenomena caused by hydropeaking operations typically affect riverbed erosion and fish displacement. A three-dimensional hydrodynamic model is used to simulate the flow characteristics during the passage of the rising limb of an observed hydropeaking event in a gravel-bed reach of Smith River, Virginia. The calculated time-dependent water surface elevations, velocities, and shear stresses are compared with field measurements. Further, comparison based on numerical simulations of this historical and a hypothetical “staggering” hydropeaking event reveals that the latter has the capability of reducing the area subject to erosion and prolonging refugia availability for juvenile brown trout. Issues related to the adoption of either a truly dynamic modeling approach or a quasi-steady methodology for simulating unsteady flows are examined through a proposed unsteadiness flow parameter. The insights obtained from this study can assist in properly accounting for the impact of hydropeaking operations on fish habitat and instream flow management.     

Analysis of Dynamic Wave Model for Unsteady Flow in an Open Channel

The models for flood propagation in an open channel are governed by Saint-Venant’s equations or by their simplified forms. Assuming the full form of hyperbolic type nonlinear expressions, the complete or dynamic wave model is obtained. Hence, after first-order linearization procedure, the dispersion relation is obtained by using the classical Fourier analysis. From this expression, the phase and group speed and the variations of the amplitude of the waves are defined and investigated. Adopting Manning’s resistance formula, the effects of the variations of the Froude number, Courant number, and friction parameter are examined in the wave number domain for progressive and regressive waves. For small and high wave numbers, the simplified kinematic and gravity wave models are recovered, respectively. Moreover, the analysis confirms, according to the Vedernikov criterion, the Froude number value corresponding to the stability condition to contrast the development of roll waves. In addition, for stable flow on the group speed versus wave number curves, the results show critical points, maximum and minimum for progressive and regressive waves, respectively.     

Analysis of the Friction Term in the One-Dimensional Shallow-Water Model

The numerical simulation of unsteady open channel flows is very commonly performed using the one-dimensional shallow-water model. Friction is one of the relevant forces included in the momentum equation. In this work, a generalization of the Gauckler-Manning friction model is proposed to improve the modeling approach in cases of dominant roughness, unsteady flow, and distorted cross-sectional shapes. The numerical stability conditions are revisited in cases of dominant friction terms and a new condition, complementary to the basic Courant-Friedrichs-Lewy condition, is proposed. Some test cases with measured data are used to validate the quality of the approaches.     

Local Time Stepping for Modeling Open Channel Flows

This paper presents two time accurate local time stepping (LTS) algorithms developed within aeronautics and develops the techniques for application to the Saint-Venant equations of open channel flow. The LTS strategies are implemented within an explicit finite volume framework based on using the Roe Riemann solver together with an upwind treatment for the source terms. The benefits of using an LTS approach over more traditional global time stepping methods are illustrated through a series of test cases, and a comparison is made between the two LTS algorithms. The results demonstrate how local time stepping can reduce computer run times, increase the reliability of the error control, and also increase the accuracy of the solution in certain regions.     

Role of Artificial Dissipation Scaling and Multigrid Acceleration in Numerical Solutions of the Depth-Averaged Free-Surface Flow Equations

A numerical model is developed for solving the depth-averaged, open-channel flow equations in generalized curvilinear coordinates. The equations are discretized in space in strong conservation form using a space-centered, second-order accurate finite-volume method. A nonlinear blend of first- and third-order accurate artificial dissipation terms is introduced into the discrete equations to accurately model all flow regimes. Scalar- and matrix-valued scaling of the artificial dissipation terms are considered and their effect on the accuracy of the solutions is evaluated. The discrete equations are integrated in time using a four-stage explicit Runge–Kutta method. For the steady-state computations, local time stepping, implicit residual smoothing, and multigrid acceleration are used to enhance the efficiency of the scheme. The numerical model is validated by applying it to calculate steady and unsteady open-channel flows. Extensive grid sensitivity studies are carried out and the potential of multigrid acceleration for steady depth-averaged computations is demonstrated.     

Systematic Evaluation of One-Dimensional Unsteady Friction Models in Simple Pipelines

In this paper, basic unsteady flow types and transient event types are categorized, and then unsteady friction models are tested for each type of transient event. One important feature of any unsteady friction model is its ability to correctly model frictional dissipation in unsteady flow conditions under a wide a range of possible transient event types. This is of importance to the simulation of transients in pipe networks or pipelines with various devices in which a complex series of unsteady flow types are common. Two common one-dimensional unsteady friction models are considered, namely, the constant coefficient instantaneous acceleration-based model and the convolution-based model. The modified instantaneous acceleration-based model, although an improvement, is shown to fail for certain transient event types. Additionally, numerical errors arising from the approximate implementation of the instantaneous acceleration-based model are determined, suggesting some previous good fits with experimental data are due to numerical error rather than the unsteady friction model. The convolution-based model is successful for all transient event types. Both approaches are tested against experimental data from a laboratory pipeline.     

Hydraulic Elements Chart for Pipe Flow Using New Definition of Hydraulic Radius

The Manning formula is used routinely to calculate the mean velocity of uniform flow. Although this empirical formula is effective when applied to uniform flow in simple rectangular or trapezoidal cross sections, the roughness coefficient of the formula is variable when examining flow in a pipe that is partially full. Thus, the coefficient must be altered depending on the relative depth of fluid in the pipe. As this seems to be due to the definition of the hydraulic radius, a new definition of hydraulic radius is proposed here that was used to calculate a hydraulic elements chart for flow in pipes with a constant roughness coefficient. The results of the calculations showed very good agreement with Camp’s chart. Furthermore, with adjustment of the “free-surface weight factor,” this method was also capable of expressing other hydraulic elements charts reported previously. This new definition of hydraulic radius can also be applied to flow in simple cross sections and may be developed further for use with compound channel flows in future studies.     

Upwind Conservative Scheme for the Saint Venant Equations

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.     

Flume for Teaching Spatially Varied Open-Channel Flow

Spatially varied flow in open channels is a topic that is often included in undergraduate open channel hydraulics courses. Physical and computational models are developed to enhance the presentation of spatially varied flow to engineering students at the late undergraduate or early graduate level. The physical model is inexpensive and easy to build and the computational model is easily developed using commercially available spreadsheet software. The physical model consists of a 30.48 cm nominal-diameter PVC pipe that is 6.1 m in length and has circular orifices approximately 1.40 cm in diameter drilled on 15.24 cm centers along the pipe invert. A relationship between the orifice discharge coefficient and a modification of the Froude number, as measured in the flume upstream of the orifice in question, was developed in repeated trials having varying flume slope, volumetric inflow rate, and end conditions. With this relationship, a stepwise solution to the energy equation is used to predict the water surface profile. Differences between the water surface profiles observed and predicted in repeated trials averaged approximately 2 mm.     

Use of Steady-State Pump Head-Discharge Curve for Unsteady Pipe Flow Applications

An experimental pipeline system with a multistage centrifugal pump was used to study the effect of transient operations on the hydrodynamic performance of a centrifugal pump. Several transient flow operations were considered, ranging from very mild to severe transients. The dynamic relationship of total pressure rise across the pump to the flow rate was compared with that of the steady state. Deviation between the dynamic pump head and the value given by the steady-state curve at the same instantaneous discharge was established and found to be a function of the severity of the transient. It was found that severe flow conditions could cause this deviation to exceed 30% of the steady-state value. The use of the steady-state pump head-discharge relationship in the solution of transient pipe flow by the method of characteristics (MOC) is discussed. It was found that the steady-state pump head-discharge curve was not accurate enough to support the solution of unsteady pipe flow application by the MOC.