Reynolds Stress Anisotropy in Open-Channel Flow

This paper reports the results of an experimental study characterizing turbulence and turbulence anisotropy in smooth and rough shallow open-channel flows. The rough bed consists of a train of two-dimensional transverse square ribs with a ratio of the roughness height (k) to the total depth of flow (d) equal to 0.10. Three rib separations (p/k) of 4.5, 9, and 18 were examined. Here, p is the pitch between consecutive roughness elements and was varied to reproduce the classical condition of d- and k-type roughness. For each case, two-component velocity measurements were obtained using a laser Doppler velocimetry system at two locations for p/k = 4.5 and 9: on the top of the rib and above the cavity, and an additional location for p/k = 18. The measurements allow examination of the local variations of the higher-order turbulent moments, stress ratios as well as turbulence anisotropy. Large variations of the turbulence intensities, Reynolds shear stress, turbulent kinetic energy and turbulence production are found for y1<3k. In this region, the flow is more directly influenced by the shear layers from the preceding ribs. The higher-order moments appear to be similar for all rough surfaces beyond y1 ≈ 7k. In the outer layer (y1>3k), all higher-order turbulent moments for the k-type roughness show a substantial increase due to the complex interactions between the roughness and the remnants overlying shear layers shed from succeeding ribs. Analysis of the components of the Reynolds stress anisotropy tensor shows that at p/k = 18, the flow at y1<5k tends to be more isotropic which implies that for this particular case, the effect of the roughness density could also be important. On the smooth bed, at the shallower depths, the correlation coefficient near the free surface increases and turbulence tends to become less anisotropic.     

Flow Turbulence over Fixed and Weakly Mobile Gravel Beds

Characteristics of turbulence structure in quasi-2D flows with static and weakly mobile gravel beds are presented. Three sets of measurements with acoustic Doppler velocimeters in an irrigation canal were used: two with subcritical bed shear stress (static beds) and one with the bed shear stress τo close to critical τoc (weakly mobile bed). The analyses included vertical distributions of local mean velocities, turbulence intensities, turbulent shear stresses, velocity auto- and cross-spectra, the quadrant method, and high-order velocity moments. A number of properties of turbulence intensities, high-order moments, streamwise bursting parameters, and velocity spectra appeared to be similar for all three flows, but some properties were different. The most important one was an observed reduction in the von Kármán constant for the flow with weakly mobile bed. Comparison of these results with other studies and analogies with drag-reducing flows suggest that at τo∕τoc ≈ 1 the drag on the bed for a given granular material should be minimized.     

Numerical Simulation of Flow over a Rough Bed

This paper presents results of a direct numerical simulation (DNS) of turbulent flow over the rough bed of an open channel. We consider a hexagonal arrangement of spheres on the channel bed. The depth of flow has been taken as four times the diameter of the spheres and the Reynolds number has been chosen so that the roughness Reynolds number is greater than 70, thus ensuring a fully rough flow. A parallel code based on finite difference, domain decomposition, and multigrid methods has been used for the DNS. Computed results are compared with available experimental data. We report the first- and second-order statistics, variation of lift/drag and exchange coefficients. Good agreement with experimental results is seen for the mean velocity, turbulence intensities, and Reynolds stress. Further, the DNS results provide accurate quantitative statistics for rough bed flow. Detailed analysis of the DNS data confirms the streaky nature of the flow near the effective bed and the existence of a hierarchy of vortices aligned with the streamwise direction, and supports the wall similarity hypothesis. The computed exchange coefficients indicate a large degree of mixing between the fluid trapped below the midplane of the roughness elements and that above it.     

Double-Averaged Open-Channel Flows with Small Relative Submergence

We investigate the turbulent structure of shallow open channel flows where the flow depth is too small (compared with the roughness height) to form a logarithmic layer but large enough to develop an outer layer where the flow is not directly influenced by the roughness elements. Since the log layer is not present, the displacement height d, which defines the position of the zero plane, and the shear velocity u* cannot be found by fitting the velocity data to the log law. However, these parameters are still very important because they are used for scaling flow statistics for the outer and roughness layers. In this paper we propose an alternative procedure for evaluating d in laboratory conditions, where d is found from additional experiments with the fully developed log layer. We also point out the appropriate procedure for evaluating the shear velocity u* for flows with low submergence. These procedures are applied to our own laboratory flume experiments with uniform sphere roughness, where velocities were measured using Particle Image Velocimetry. Results were interpreted within the framework of the double-averaged Navier–Stokes equations and include mean velocities, turbulence intensities, Reynolds stresses, and form-induced normal and shear stresses. The data collapse well and show that in flows without a developed log layer the structure of turbulence in the outer layer remains similar to that of flows with a log layer. This means that even though the roughness layer in the experiments reported herein was sufficiently high to prevent the development of the log layer, influence of the bed roughness did not spread further up into the outer layer. Furthermore, the results show that flow statistics do not depend on relative submergence except for the form-induced stresses which increase when relative submergence decreases.     

Transverse Dispersion Caused by Secondary Flow in Curved Channels

A new theoretical equation is proposed to describe the streamwise variations of the transverse velocity along a curved channel with a constant curvature. Furthermore, based on this theoretical equation for the transverse velocity, a new equation for the transverse dispersion coefficient is developed to incorporate the effect of the secondary flow on the transverse dispersion in curved channels. The new equations for the transverse velocity and dispersion coefficient are verified with experimental data sets that were obtained from laboratory experiments conducted in two different channels. The results show that the proposed velocity equation properly describes the streamwise variations of the secondary flow developed in the curved channels. The reach-averaged values of the transverse dispersion coefficient calculated by the new equation are in relatively good agreement with the observed values from the laboratory channels. Sensitivity analysis reveals that both the secondary flow and the transverse dispersion coefficient are proportional to the roughness factor, and in inverse proportion to the aspect ratio of the channel.     

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.     

Unobstructed and Obstructed Turbulent Flow in Gravel Bed Rivers

An acoustic Doppler velocimeter was used to characterize turbulence in two gravel bed rivers. Data were collected in unobstructed flow and compared to recent investigations. Additional data collected in the wake of emergent boulders indicate that mean flow velocity, turbulent kinetic energy, gradients in the streamwise velocity, and Reynolds stress downstream from large rocks deviate from unobstructed flow results, but similar turbulence patterns are found behind each boulder. Results of this study are discussed with regard to natural channel design and fish habitat.     

Numerical Study of the Transition Regime between the Skimming and Wake Interference Flows in a Water Flume by Using the Lattice-Model Approach

A numerical study to describe the transition regime between the skimming and wake interference flows due to the influence of an idealized bed roughness in a water flume was carried out here using the lattice model approach. The model reproduced the skimming, transition, and wake interference regimes for different aspect ratios that determine the bed roughness geometry. The simulated turbulent structures were visualized by drawing the trajectories of a large number of passive tracer particles released in the computational domain, and the results agreed with those reported by the research works. The dimensionless streamwise and vertical turbulent intensities were calculated at five test sections. The results obtained supported the visualized flow patterns permitting us to detect the presence of a shear layer developed at the top of the roughness element, whose strength varied according to the flow regime simulated.     

Response of Velocity and Turbulence to Sudden Change of Bed Roughness in Open-Channel Flow

The experimental study shows how an open-channel flow would respond to a sudden change (from smooth to rough) in bed roughness. Using a two-dimensional acoustic Doppler velocimeter and a laser Doppler velocimeter, the velocity, turbulent intensities, and Reynolds stress profiles at different locations along a laboratory flume were measured. Additionally, the water surface profile was also measured using a capacitance-type wave height meter. The experimental data show the formation of an internal boundary layer as a result of the step change in bed roughness. The data show that this boundary layer grows much more rapidly than that formed in close-conduit flows. The results also show that the equivalent bed roughness, bed-shear stress, turbulent intensities, and Reynolds stress change gradually over a transitional region, although the bed roughness changes abruptly. The behavior is different from that observed in close-conduit flows, where an overshooting property—which describes the ability of the bed-shear stress to attain a high-peak value over the section with the larger roughness, was reported. A possible reason for the difference is the variation of the water surface profile when an open-channel flow is subjected to a sudden change in bed roughness.     

Wall-Wake Flows Downstream of a Sphere Placed on a Plane Rough Wall

The time-averaged characteristics of turbulent wall-wake flows downstream of a sphere placed on a rough wall are studied. The profiles of the defect of streamwise velocity, Reynolds shear stress, and turbulence intensities exhibit some degree of similarities when they are scaled by their respective peak defect values. For the velocity defect profiles, the vertical distances are scaled by the height of the location of the half-peak velocity defect. However, for the defect profiles of the Reynolds shear stress and the turbulence intensities, the vertical distances are scaled by the height of the location of the half-peak Reynolds shear stress defect. The magnitudes of the peak defect of all the quantities diminish with the distance downstream of the sphere characterizing the recovery of their undisturbed profiles. Additionally, the theoretical similarity solution for the velocity defect profiles is obtained. The third-order correlations imply that in the inner layer of wall wakes, a streamwise acceleration is prevalent and associated with a downward flux, suggesting sweeps. In contrast, in the outer layer, a streamwise deceleration exists and is associated with an upward flux, suggesting ejections. The profiles of the energy budget show that the turbulent and pressure energy diffusions oppose each other. The turbulent production has a positive peak, and the pressure energy diffusion has a negative peak, indicating a large gain in turbulence production in the wall-wake flows. The quadrant analysis confirms that in wall-wake flows, sweeps are the governing mechanism resulting from an inrush of fluid streaks. The bursting events have shorter duration, but they are more frequent than those in upstream.     

Validation of a Three-Dimensional Numerical Code in the Simulation of Pseudo-Natural Meandering Flows

Validation of a three-dimensional finite volume code solving the Navier–Stokes equations with the standard k-ε turbulence model is conducted using a high quality and high spatial resolution data set. The data set was collected from a large-scale meandering channel with a self-formed fixed bed, and comprises detailed bed profiling and laser Doppler anemometer velocity measurements. Comparisons of the computed primary and secondary velocities are made with those observed and it is found that the lateral momentum transfer is generally under predicted. At the apices this results in the predicted position of the primary velocity maximum having a bias towards the channel center, compared to the position where it has been measured. Using a simplified two zone roughness distribution whereby a separate roughness height was prescribed for the channel center and channel sides relative to a single distributed roughness height, generally led to a slightly improved longitudinal velocity distribution; the higher velocities were located nearer to the outside of the bend. Improving both the free surface calculation and scheme for discretization of the convection terms led to no appreciable difference in the computed velocity distributions. A more detailed study involving turbulence measurements and bed form height distribution should discriminate whether using distributed roughness height is a precursor to using an anisotropic turbulence representation for the accurate prediction of three-dimensional river flows.     

Turbulence Characteristics and Interaction between Particles and Fluid in Particle-Laden Open Channel Flows

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.     

Three-Dimensional Numerical Simulation of Compound Channel Flows

A three-dimensional numerical study is presented for the calculation of turbulent flow in compound channels. The flow simulations are performed by solving the three-dimensional Reynolds-averaged continuity and Navier–Stokes equations with the k?ε turbulence model for steady-state flow. The flow equations are solved numerically with a general-purpose finite-volume code. The results are compared with the experimental data obtained from the UK Flood Channel Facility. The simulated distributions of primary velocity, bed shear stress, turbulent kinetic energy, and Reynolds stresses are used to investigate the accuracy of the model prediction. The results show that, using an estimated roughness height, the primary velocity distributions and the bed shear stress are predicted reasonably well for inbank flows in channels of high aspect ratio (width/depth ≥ 10) and for high overbank flows with values of the relative flow depth greater than 0.25.     

Rough Wall Boundary Layer on Plates in Open Channels

Experiments were performed to measure the characteristics of a turbulent boundary layer developing on a rough surface placed in an open channel flow at close proximity to the free surface. Streamwise velocity measurements were made with a one-component laser Doppler velocimeter system at the top of the spherical roughness elements. Measurements at three stations downstream of the plate leading edge show the growth of the boundary layer on the rough wall and its interaction with the exterior open-channel flow and the free surface. Resorting to the turbulence profile provides an alternative definition of the boundary layer thickness. The near-wall flow follows the well-known logarithmic law with a shift due to roughness. In the outer layer, there are two opposing effects: the free surface tends to decrease the wake component while the roughness tends to increase it. The streamwise turbulence intensity is affected by the shear and turbulence in the exterior flow, the effect of the free surface being greater than that of wall roughness.     

Numerical Modeling of Three-Dimensional Flow Field Around Circular Piers

A three-dimensional numerical model FLUENT is used to simulate the separated turbulent flow around vertical circular piers in clear water. Computations are performed using different turbulence models and results are compared with several sets of experimental data available in the literature. Despite commonly perceived weakness of the k-ε model in resolving three-dimensional (3D) open channel and geophysical flows, several variants of this turbulence model are found to have performed satisfactorily in reproducing the measured velocity profiles. However, model results obtained using the k-ε models show some discrepancy with the measured bed shear stress. The Reynolds stress model performed quite well in simulating velocity distribution on flat bed and scour hole as well as shear stress distribution on flat bed around circular piers. The study demonstrates that a robust 3D hydrodynamic model can effectively supplement experimental studies in understanding the complex flow field and the scour initiation process around piers of various size, shape, and dimension.     

Outer Scaling for Open Channel Flow over a Gravel Bed

Similarity analysis is performed for hydraulically rough open channel flow over a gravel bed to provide mixed outer scaling of the mean-velocity profile. The analysis is based on equilibrium turbulent boundary-layer theory derived using the asymptotic invariance principle. Outer scaling based on the similarity theory is validated with velocity measurements from the laboratory and field, having a Reynolds number range that includes 1×104, 1×105, and 1×106 and a Froude number range from 0.26 to 0.83. The results show that the free-stream velocity is an appropriate outer scale for gravel-bed river flows at moderate and bankfull stage. The results agree well with the velocity measurements and collapse laboratory and field data, which allow an important connection between open channel research in the laboratory and the applications for which the research is performed in the field. The results show that the R/aD84 roughness parameter is consistent with the mixed scale used in boundary-layer velocity scaling. This is in agreement with the consistent turbulent structure of the flow for both flat plate boundary-layer and open channel flow scenarios. While R/D84 has been used empirically with depth-averaged velocity and roughness laws for many years, this roughness parameter is shown in a theoretical context due to its influence on the turbulent structure of the flow. The results are applicable to modeling the velocity distribution under fundamental gravel-bed flow cases that span to the bankfull flow regime, which provides a contribution to stream engineering.     

Characteristics of Turbulent Flow in Submerged Jumps on Rough Beds

The results of an experimental investigation on the flow field in submerged jumps on horizontal rough beds, detected by an acoustic Doppler velocimeter, are presented. Experiments were conducted for the conditions of submerged jumps, having submergence factors from 0.96 to 1.85 and jet Froude numbers from 2.58 to 4.87, over rough beds of Nikuradse’s equivalent sand roughness equaling 0.49, 0.8, 1.86, and 3?mm. The vertical distributions of time-averaged velocity components, turbulence intensity components, and Reynolds stress at different streamwise distances from the sluice opening and the horizontal distribution of bed-shear stress are plotted. Vector plots of the flow field show that the rate of decay of jet velocity in a submerged jump increases with increase in bed roughness. The flow characteristics on rough beds, being different from those on smooth bed, are discussed from the point of view of similarity, growth of the length scale, and decay of the velocity and turbulence characteristics scales. The most important observation is that the flow in the fully developed zone is found to be self-preserving.     

Nonlinear Mixing Length Model for Prediction of Secondary Currents in Uniform Channel Flows

A nonlinear turbulence model for numerical solution of uniform channel flow is presented. Turbulent stresses are evaluated from a nonlinear mixing length model that relates turbulent stresses to quadratic products of the mean rate of strain and the mean vorticity. The definition of the mixing length, based on a three-dimensional integral measure of boundary proximity, eliminates the need for solution of additional transport equations for the turbulence quantities. Experimental data from the literature for closed and open-channel flows are utilized to validate the model. The model produced the secondary flow vortices successfully. Velocity field and wall shear stresses affected by secondary flow vortices are accurately computed. Bulging of velocity contour lines toward the corners and dipping phenomena of maximum velocity are successfully simulated.     

Spatially Averaged Open-Channel Flow over Rough Bed

In this paper it is suggested that the double-averaged (in temporal and in spatial domains) momentum equations should be used as a natural basis for the hydraulics of rough-bed open-channel flows, especially with small relative submergence. The relationships for the vertical distribution of the total stress for the simplest case of 2D, steady, uniform, spatially averaged flow over a rough bed with flat free surface are derived. These relationships explicitly include the form-induced stresses and form drag as components of the total stress. Using this approach, we define three types of rough-bed flows: (1) Flow with high relative submergence; (2) flow with small relative submergence; and (3) flow over a partially inundated rough bed. The relationships for the double-averaged velocity distribution and hydraulic resistance for all three flow types are derived and compared with measurements where possible. The double-averaged turbulent and form-induced intensities and stresses for the case of regular spherical-segment-type roughness show the dominant role of the double-averaged turbulence stresses and form drag in momentum transfer in the near-bed region.