Experiments on Flow at a 90° Open-Channel Junction

Although open-channel junctions are common in many hydraulic structures, no comprehensive data set has been compiled that describes the 3D flow field within the junction itself. This physical model study examined a 90°, sharp-edged, open-channel junction for channels of equal width. Depth measurements were made using a point gauge while velocity measurements were taken using an acoustic doppler velocimeter over a grid defined throughout the junction region. The average velocity and turbulence intensity were calculated from a time series of velocities that was recorded at each location. In addition, a 2D mapping of the water surface was performed on a 76.2 mm square grid throughout the channel junction. This paper presents the details of the experimental procedure and the general flow characteristics observed. The full data set generated during this experimental work is available for downloading on the Internet. Using a small portion of the data recorded, an evaluation of several previously proposed theories of combining flow in open-channel junctions is presented. This has revealed that the simplified mathematical model gives reasonable prediction of the experimental results. The complete data set describing combining flows at a 90° open-channel junction is presented as a resource for the validation of 3D computational fluid dynamics codes that utilize a free-surface model.     

Three-Dimensional Numerical Study of Flows in Open-Channel Junctions

An open-channel junction flow is encountered in many hydraulic structures ranging from wastewater treatment facilities to fish passage conveyance structures. An extensive number of experimental studies have been conducted but a comprehensive three-dimensional numerical study of junction flow characteristics has not been performed and reported. In this paper, a three-dimensional numerical model is developed to investigate the open-channel junction flow. The main objective is to present the validation of a three-dimensional numerical model with high-quality experimental data and compare additional simulations with classical one-dimensional water surface calculations. The three-dimensional model is first validated using the experimental data of a 90° junction flow under two flow conditions. Good agreement is obtained between the model simulation and the experimental measurements. The model is then applied to investigate the effect of the junction angle on the flow characteristics and a discussion of the results is presented.     

Simulation of Flow past an Open-Channel Floor Slot

In the past, solutions to the problem of flow past a floor slot in a rectangular open channel used to divert flow from one stream to another were obtained mainly on the basis of model tests or through the development of simplified theoretical expressions. In the present study, the free-surface turbulence model is applied to obtain the flow parameters such as pressure head distribution, velocity distribution, and water surface profile. The predictions of the proposed numerical model are validated using previous experimental data. In particular, the model predictions agree well with the test data related to flow parameters. The study indicates that the free-surface turbulence model developed is an efficient and useful tool for predicting characteristics of free surface flows such as flow past a floor slot. For flow past an open-channel floor slot, a model that is properly validated can be used to predict the flow characteristics under various flow configurations encountered in the field, without resorting to expensive experimental procedures.     

Two-Dimensional Simulation of Subcritical Flow at a Combining Junction: Luxury or Necessity?

Classically, in open-channel networks, the flow is numerically approximated by the one-dimensional Saint Venant equations coupled with a junction model. In this study, a comparison between the one-dimensional (1D) and two-dimensional (2D) numerical simulations of subcritical flow in open-channel networks is presented and completely described allowing for a full comprehension of the modeling of water flow. For the 1D, the mathematical model used is the 1D Saint Venant equations to find the solution in branches. For junction, various models based on momentum or energy conservation have been developed to relate the flow variables at the junction. These models are of empirical nature due to certain parameters given by experimental results and moreover they often present a reduced field of validity. In contrast, for the 2D simulation, the junction is discretized into triangular cells and we simply apply the 2D Saint Venant equations, which are solved by a second-order finite-volume method. In order to give an answer to the question of luxury or necessity of the 2D approach, the 1D and 2D numerical results for steady flow are compared to existing experimental data.     

3D Simulation of Combining Flows in 90° Rectangular Closed Conduits

Combining flows are encountered often in environmental engineering and hydraulic engineering. Experimental data are available to assist the engineers who need the various loss coefficients associated with combining flows in closed conduits. For the combining flows in 90° rectangular conduit junctions, the Reynolds averaged Navier–Stokes equations are applied, while using the three-dimensional k-ω model. The energy loss coefficients and the mean flow pattern are obtained and validated by experimental data. The numerical modeling is less time-consuming and less expensive to obtain the various flow parameters needed for engineering design.     

Characteristics of Flow around Open Channel 90° Bends with Vanes

Sharp open-channel bends are commonly encountered in hydraulic engineering design. Disturbances such as secondary flows and flow separation caused by the bend may persist for considerable distances in the downstream channel. A simple way of reducing these disturbances is through the insertion of vertical vanes in the bend section. A laser Doppler anemometry (LDA) unit was used to measure the three-dimensional mean and turbulent velocity components of flow in an experimental rectangular open-channel bend. Flow characteristics of the bend with no vanes are compared with those of bends having one or three vertical vanes. The size of the flow separation zone at the inner wall of the bend was determined from dye visualization data and confirmed with mean streamwise velocity data. Results show that the vertical vanes are effective in considerably reducing flow separation, intensity of secondary flows, and turbulence energy in the downstream channel. Furthermore, energy loss for bends with vanes is slightly less than for the no-vane case.     

Turbulence Measurement in Nonuniform Open-Channel Flow Using Acoustic Doppler Velocimeter (ADV)

Measurements of the mean and turbulence characteristics in nonuniform open-channel flows were carried out using a 3D acoustic Doppler velocimeter. Both accelerating and decelerating flows were investigated. Analyses based on the Reynolds equation and the continuity equation of 2D open-channel flow show that a flow will be in equilibrium if the pressure-gradient parameter β is a constant at different sections along the flow direction. The experimental data show that all the flows investigated are in equilibrium. The effect of the flow nonuniformity on the mean velocity and turbulence profiles was also examined. The study shows that (1) the log law is still valid for both accelerating and decelerating flows in the inner region. The Coles law can be used for the entire region, but the wake-strength parameter Π depends on the β-value; and (2) the turbulence intensities and the Reynolds stress decrease in accelerating flow and increase in decelerating flow, when compared with those in uniform flow. Finally, using the Reynolds equation and the continuity equation of 2D open-channel flow, the theoretical expressions for the distribution of vertical velocity and the Reynolds stress have been developed. The measured vertical velocity and the Reynolds stress profile compare well with the derived expressions.     

Dynamic Model for Subcritical Combining Flows in Channel Junctions

A one-dimensional theoretical model for subcritical flows in combining open channel junctions is developed. Typical examples of these junctions are encountered in urban water treatment plants, irrigation and drainage canals, and natural river systems. The model is based on applying the momentum principle in the streamwise direction to two control volumes in the junction together with overall mass conservation. Given the inflow discharges and the downstream depth, the proposed model solves for each of the upstream depths. The interfacial shear force between the two control volumes, the boundary friction force, and the separation zone shear force downstream of the lateral channel entrance are included. Predictions based on the proposed approach are shown to compare favorably with existing experimental data, previous theories, and conventional junction modeling approaches. The main advantages of the proposed model are that the proposed model does not assume equal upstream depths and that the dynamic treatment of the junction flow is consistent with that of the channel reaches in a network model.     

New Approach for Predicting Flow Bifurcation at Right-Angled Open-Channel Junction

An unsteady mathematical model for predicting flow divisions at a right-angled open-channel junction is presented. Existing dividing models depend on a prior knowledge of a constant flow regime. In addition, their strong nonlinearity does not guarantee compatibility with the St. Venant solutions in the context of an internal boundary condition treatment. Assuming zero crest height at the junction region, a side weir model explicitly introduced within the one-dimensional St. Venant equations is used to cope with the two-dimensional pattern of the flow. An upwind implicit numerical solver is employed to compute the new governing equations. The performance of the proposed technique in predicting super-, trans-, and subcritical flow bifurcations is illustrated by comparing with experimental data and/or theoretical predictions. In all the tests, lateral-to-upstream discharge ratios (Rq) are successfully reproduced by the present technique with a maximum error magnitude of less than 9%.     

Fluctuations of Suspended Sediment Concentration and Turbulent Sediment Fluxes in an Open-Channel Flow

A complex problem of turbulent-sediment interactions in an open-channel flow is approached experimentally, using specially designed field experiments in an irrigation canal. The experimental design included synchronous measurements of instantaneous three-dimensional (3D) velocities and suspended sediment concentration using acoustic Doppler velocimeters (ADV) and a water sampling system. Various statistical measures of sediment concentration fluctuations, turbulent sediment fluxes, and diffusion coefficients for fluid momentum and sediment are considered. Statistics, fractal behavior, and contributions of bursting events to vertical fluxes of fluid momentum and sediment are evaluated using quadrant analysis. It has been found that both turbulence and sediment events are organized in fractal clusters which introduce additional characteristic time and spatial scales into the problem and should be further explored. It is also shown that Barenblatt’s theory of sediment-laden flows appears to be a good approximation of experimental data.     

Three-Dimensional Numerical Model of Lateral-Intake Inflows

A three-dimensional (3D) numerical model for predicting steady, in the mean, turbulent flows through lateral intakes with rough walls is developed, validated, and employed in a parametric study. The method solves the Reynolds-averaged Navier-Stokes equations closed with the isotropic k-ω turbulence model of Wilcox, which resolves the near-wall flow and accounts for roughness effects in a straightforward manner. Calculations are carried out for flows through rectangular closed-duct and open-channel T-junctions. Comparisons of the predicted mean velocity field with laboratory measurements indicate that the model captures most experimental trends with reasonable accuracy. For the parametric study, flows are predicted for a range of discharge ratios, aspect ratios, and main channel-bed-roughness distributions. The numerical solutions are examined to elucidate the complex 3D flow patterns of lateral-intake flows, including zones of flow division, separation and reversal, vortices, and singular points within the bed-shear stress vector field. The model reproduces known 3D flow patterns and provides novel insights about the complex hydraulics and sediment transport processes encountered in lateral intakes at a level of detail that is not attainable by laboratory studies alone.     

Modeling 3D Supercritical Flow with Extended Shallow-Water Approach

Despite the three-dimensional (3D) nature of the flow, the classical shallow-water equations are often used to simulate supercritical flow in channel transitions. A closer comparison with experimental data, however, often shows large discrepancies in the height and pattern of the shock waves that increase with the Froude number. An extension to the classical shallow-water approach is derived considering higher-order distribution functions for pressure and horizontal and vertical velocities, therefore taking nonhydrostatic pressure distribution and vertical momentum into account. The approach is applied to highly supercritical flow in a channel contraction (F0 = 4.0), a channel junction (F0 ≈ 4.5 for both branches), and a channel expansion (F0 = 8.0). Specific problems of such flows—wetting and drying of computational cells and wave breaking due to steep free-surface gradients—are discussed and solved numerically. The solutions with the extended approach are compared both with experimental data and classical shallow-water computations, and the influence of the additional terms considering the 3D nature of such flows is illustrated.     

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.     

Hydraulic Jet Control for River Junction Design of Yuen Long Bypass Floodway, Hong Kong

The Yuen Long Bypass Floodway (YLBF) was designed to collect flows from the Sham Chung River (SCR) and the San Hui Nullah (SHN) and to serve as a diversion channel of the Yuen Long Main Nullah (YLMN). Under a 200-year return period design condition, the floodway was designed (1) to divert a flow of approximately 38?m3/s from the supercritical YLMN flow and (2) to convey a total combined flow of 278?m3/s to downstream within acceptable flood levels. The success of the design depends critically on complicated junction flow interactions that cannot be resolved by 1D unsteady flow models. These features include the supercritical-subcritical flow transition at the San Hui-Floodway (SHN-YLBF) junction and the diversion of part of the supercritical flow from the Main Nullah (YLMN). A laboratory Froude scale physical model was constructed to study water stages and flow characteristics in the floodway and to investigate optimal design arrangements at channel junctions and transitions. This paper summarizes the main features of the unique river junction network, in particular the use of the hydraulic jet principle at the SHN-YLBF junction to lower flood levels. In addition, a numerical flow model is employed to study flow details at the river junctions. The model is based on the general 2D shallow water equations in strong conservation form. The equations are discretized using the total variation diminishing finite-volume method which captures the discontinuity in hydraulic jumps. The numerical model predictions are well supported by the laboratory data, and the theoretical and experimental results offer useful insights for the design of urban flood control schemes under tight space constraints.     

Momentum Transport in Sharp Open-Channel Bends

Flow in open-channel bends is characterized by cross-stream circulation, which redistributes the velocity and the boundary shear stress and thereby shapes the characteristic bed topography. Besides a center-region cell, classical helical motion, a weaker counterrotating outer-bank cell often exists. In spite of its engineering importance, the mechanisms underlying distributions of the velocity and the boundary shear stress in open-channel bends, and especially the role of both circulation cells, are not yet fully understood. In order to investigate these mechanisms, an evaluation is made of the various terms in the momentum equations based on the data measured, which gave the following results. The outer-bank cell forms a buffer layer that protects the outer bank from any influence of the center-region cell and keeps the core of maximum velocity a distance from the bank. Advective momentum transport by the center-region cell is a dominant mechanism; it significantly contributes to the observed outward shift of the downstream velocity and the bed shear stress and to flattening of the vertical profiles of the velocity. This important advective momentum redistribution has to be included in the depth-integrated flow models often used in engineering practice. Commonly used linear models overpredict the effects of the center-region cell. Based on results of the analysis of experimental data, these models are extended by accounting for the feedback between the center-region cell and the downstream velocity. The nonlinear model obtained clearly reveals the mechanisms of the center-region cell and its advective momentum transport. An analysis of nonlinear model results confirms and complements the analysis of experimental data. A true quasithree-dimensional flow model is obtained by coupling this nonlinear model to depth-integrated flow models, thus providing an engineering tool for morphodynamical investigations.     

Role of Bed Discordance at Asymmetrical River Confluences

This paper studies laboratory open-channel confluences using a 3D, elliptic solution of the Reynolds-averaged Navier-Stokes equations, including a method for approximating the effects of water surface elevation patterns and a renormalization group modified form of the k-ε turbulence model. The model was tested by comparison with laboratory measurements of an asymmetric tributary junction. This suggests that although the model is unable to reproduce the quantitative detail (notably upwelling velocity magnitudes) of the flow structures as measured in laboratory experiments, statistically significant aspects of the experimental observations are reproduced. The model is used to (1) describe and explain the characteristic flow structures that form in a confluence with one of the tributaries angled at 45°, both with and without an elevation difference (bed discordance) in the angled tributary; and (2) investigate the relative importance of junction angles (30°, 45°, and 60°), bed discordance, and ratio of mean velocities in the tributary channels upon flow structures. This shows that bed discordance significantly enhances secondary circulation because of the effects of flow separation in the lee of the bed step, which significantly increases lateral pressure gradients at the bed and reduces water surface superelevation in the center of the tributary and water surface depression at the downstream junction corner. Extension to consideration of a number of junction angles, levels of bed discordance, and velocity ratios suggests that a small (10%) reduction in tributary depth can significantly increase the intensity of secondary circulation, albeit in a relatively localized manner. Simulations involving a numerical tracer illustrate the importance of bed discordance for mixing between the two flows and question the use of simple 2D parameterizations of mixing processes that do not consider bed discordance when the latter is present.     

Hydraulic Models of the Flow Distribution in a Four Branch Open Channel Junction with Supercritical Flow

Intense rainfall on urban areas can generate severe flooding in the city, and if the conditions are right, the flow in the streets can be supercritical. The redistribution of the flow in street intersections determines the flow rates and water levels in the street network. We have investigated the flow that occurs when two supercritical flows collide in a 90° junction formed by streets of identical cross section. Several flow configurations within the intersection are possible, depending on the position of the hydraulic jumps that form in and upstream of the intersection. Previous work has identified three flow types, with Type II flows being further classified into three subregimes. Hydraulic models have been developed, based on the principles of the conservation of flow and momentum flux in the intersection, which predict the angles at which the jumps will form. These models can be used to determine the flow type that will occur. Moreover, additional models have been developed for computing the outflow discharge distribution. For Type I flows, it has not been possible to develop such a hydraulic model for the discharge distribution, but some data are provided for one configuration to indicate the influence of different parameters. For Type II and Type III flows, such models are developed, and their predictions agree with data obtained from the channel intersection facility at the Laboratory of Fluid Mechanics and Acoustics in Lyon.     

Vertical Dispersion of Fine and Coarse Sediments in Turbulent Open-Channel Flows

Through using a kinetic model for particles in turbulent solid–liquid flows, underlying mechanisms of sediment vertical dispersion as well as sediment diffusion coefficient are investigated. Four hydrodynamic mechanisms, namely gravitational settling, turbulent diffusion, effect of lift force, and that of sediment stress gradient, coexist in two-dimensional (2D) uniform and steady open-channel flows. The sediment diffusion coefficient consists of two independent components: one accounts for the advective transport of sediment probability density distribution function due to sediment velocity fluctuations, and the other results from sediment–eddy interactions. Predictions of the kinetic model are in good agreement with experimental data of 2D open-channel flows. In such flows, it is shown that: (1) the parameter γ (i.e., the inverse of the turbulent Schmidt number) may be greater than unity and increases toward the bed, being close to unity for fine sediments and considerably large for coarse ones; (2) effects of lift force and sediment stress gradient become significant and need to be considered below the 0.1 flow depth; and (3) large errors may arise from the traditional advection–diffusion equation when it is applied to flows with coarse sediments and/or high concentrations.     

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.