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.     

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.     

Turbulence Characteristics in Flows Subjected to Boundary Injection and Suction

The influence of seepage (lateral flow) on the turbulence characteristics in free-surface flows over an immobile rough boundary is investigated. Steady flows having zero-pressure gradient over an immobile rough boundary created by uniform gravels of 4.1 mm in size were simulated experimentally with injection (upward seepage) and suction (downward seepage) applied through the boundary. A Vectrino (acoustic Doppler velocimeter) was used to measure the instantaneous velocities, which are analyzed to explore second- and third-order correlations, turbulent kinetic energy, turbulent energy budget, and conditional Reynolds shear stresses. It is observed that the second-order correlations decrease in presence of injection and increase in suction. The turbulent diffusivity and mixing length increase in presence of injection and decrease in suction. The third-order correlations suggest that the ejections are prevalent over the entire flow depth. The near-boundary flow is significantly influenced by the existence of upward seepage, which is manifested by a reduction in streamwise flux and the vertical advection of streamwise Reynolds normal stress. In addition, the upward flux and the streamwise advection of vertical Reynolds normal stress are also affected. The streamwise flux of turbulent kinetic is found to migrate upstream, while the vertical flux of turbulent kinetic energy is transported upward. The fluxes increase in presence of injection and decrease in suction. Energy budget evidences a lag between the turbulent dissipation and production and an opposing trend in the turbulent and pressure energy diffusions. A quadrant analysis for the conditional Reynolds shear stresses reveals that the ejection and sweep events are the primary contributions toward the total Reynolds shear stress production, with ejections dominating over the entire flow depth. The effect of seepage is shown to affect the magnitude of such events. However, in case of sweeps, this phenomenon is the opposite. The mean-time of occurrence of ejections and that of sweeps in suction are more persistent than those in no-seepage and injection.     

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.     

Turbulence Intensities of Shallow Rain-Impacted Flow over Rough Bed

A thorough understanding of the effect of raindrop impacts on the characteristics of a shallow flow is very helpful to properly predict and control water erosion. The turbulence properties of a shallow rain-impacted flow over a rough bed were studied using a two-dimensional fiber-optic laser Doppler velocimeter and an oscillating type raindrop generator. It was found that the value of momentum correction factor β for a rough boundary was much smaller than that for a smooth boundary at a given value of Reynolds number because of the vertical mixing effect of the boundary roughness. For the condition without rainfall, the vertical distribution of the longitudinal turbulence intensity became more uniform as the boundary roughness increased. For the “fully rough” condition with rainfall and constant roughness size, the effect of roughness on the distributions of turbulence intensities was small. Similar to the condition with smooth boundary, the rainfall-generated turbulence was greatest near the water surface and decreased downward at an exponential rate for the flow over rough bed. Two types of statistical correlations were developed to predict the distributions of both the longitudinal and the vertical turbulence intensities depending on whether the flow is completely penetrated by the raindrops.     

Velocity Distribution of Turbulent Open-Channel Flow with Bed Suction

This study investigates theoretically and experimentally the velocity distributions of turbulent open channel flow with bed suction. A velocity profile with a slip velocity at the bed surface and an origin displacement under the bed surface is proposed and discussed. Based on this assumption, a modified logarithmic law is derived. The measured experimental velocity distribution verifies the accuracy of the theoretically derived profile. The data show a significant increase in the near bed velocity and a velocity reduction near the water surface, resulting in the formation of a more uniform velocity distribution. The values of the origin displacement, slip velocity and shear velocity are found to increase with increasing relative suction. The measured data show the occurrence of two flow regions in the suction zone: a transitional region in which the velocity readjusts rapidly; and an “equilibrium” region.     

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.     

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.     

Experimental Study of Three-Dimensional Flow Field around a Complex Bridge Pier

In this paper, three-dimensional turbulent flow field around a complex bridge pier placed on a rough fixed bed is experimentally investigated. The complex pier foundation consists of a column, a pile cap, and a 2×4 pile group. All of the elements are exposed to the approaching flow. An acoustic-Doppler velocimeter was used to measure instantaneously the three components of the velocities at different horizontal and vertical planes. Profiles and contours of time-averaged velocity components, turbulent intensity components, turbulent kinetic energy, and Reynolds stresses, as well as velocity vectors are presented and discussed at different vertical and horizontal planes. The approaching boundary layer at the upstream of the pile cap separated in two vertical directions and induced an upward flow toward the column and a contracted downward flow below the pile cap and toward the piles. The contracted upward flow on the pile cap interacts with downflow in the front of the column and deflects toward the side of the pier, which in return produces a strong downflow along the side of the pile cap. The flow at the rear of the pile cap is very complex. The strong downward flow at the downstream and near the top of the pile cap in interaction with the reverse flow behind the column and upward flow near the bed produce two vortices close to the upper and lower corners of the pile cap with opposite direction of rotation. On the other hand, the back-flow from the wake of the pile cap is forced into the top region resulting in a secondary recirculation at the wake of the column. The contracted flow below the pile cap and toward the piles, a strong downflow along the sides of the pile cap at the upstream region, and a vortex flow behind the pile cap and an amplification of turbulence intensity along the sides of the pile cap at the downstream region are the main features of the flow responsible for the entrainment of bed sediments.     

Near-Bed Turbulence Characteristics at the Entrainment Threshold of Sediment Beds

This experimental study is devoted to quantification of the near-bed turbulence characteristics at an entrainment threshold of noncohesive sediments. Near the bed, the departure in the distributions of the observed time-averaged streamwise velocity from the logarithmic law is more for immobile beds than for entrainment-threshold beds. In the Reynolds shear stress distributions, a damping that occurs near the bed for sediment entrainment is higher than that for immobile beds. Quadrant analysis reveals that in the near-bed flow zone, ejections and sweeps on immobile beds cancel each other, giving rise to the outward interactions, whereas sweeps are the dominant mechanism toward sediment entrainment. The bursting duration for entrainment-threshold beds is smaller than that for immobile beds. On the other hand, the bursting frequency for entrainment-threshold beds is larger than that for immobile beds. The third-order correlations indicate that during sediment entrainment, a streamwise acceleration associated with a downward flux and advection of streamwise Reynolds normal stress is prevalent. The streamwise and the downward vertical fluxes of turbulent kinetic energy (TKE) increase with sediment entrainment. The TKE budget proves that for sediment entrainment, the pressure energy diffusion changes drastically to a negative magnitude, indicating a gain in turbulence production.     

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.     

Velocity Pulse Model for Turbulent Diffusion from Flowing Water into a Sediment Bed

The “velocity pulse model” simulates the transfer of turbulence from flowing water into a sediment bed, and its effect on the diffusional mass transfer of a solute (e.g., oxygen, sulfate, or nitrate) in the sediment bed. In the “pulse model,” turbulence above the sediment surface is described by sinusoidal variations of vertical velocity in time. It is shown that vertical velocity components dampen quickly inside the sediment when the frequency of velocity fluctuations is high and viscous dissipation is strong. Viscous dissipation (ν) inside the sediment is related to the apparent viscosity depending on the structure of the sediment pore space, i.e., the porosity and grain diameter, as well as inertial effects when the flow is turbulent. A value ν/ν0 between 1 and 20 (ν0 is kinematic viscosity of water) has been considered. Turbulence penetration into the sediment is parametrized by the Reynolds number Re = UL/ν and the relative penetration velocity W/U, where U=amplitude of the velocity pulse; and W=penetration velocity; L = WT=wave length of the velocity pulse; and T is its period. Amplitudes of vertical velocity components inside the sediment and their autocorrelation functions are computed, and the results are used to estimate eddy viscosity inside the sediment pore system as a function of depth. Diffusivity in the sediment pore system is inferred by using turbulent or molecular Schmidt numbers. Turbulence penetration from flowing water can enhance the vertical diffusion coefficient in a sediment bed by an order of magnitude or more. Penetration depth of turbulence is higher for low frequency velocity pulses. Vertical diffusivity inside the pore system is shown to decrease more or less exponentially with depth below the sediment/water interface. Vertical diffusivities in a sediment bed estimated by the “velocity pulse model” can be used in pore water quality models to describe vertical transport from or into flowing surface water. The analysis has been conducted for a conservative material, but source and sink terms can be added to the vertical transport equation.     

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.     

Mean Flow and Turbulence Structure in Vertical Slot Fishways

This paper presents the results of an experimental study on the mean and turbulence structures of flow in a vertical slot fishway with slopes of 5.06 and 10.52%. Two flow patterns existed in the fishway and for each one, two flow regions were formed in the pools: a jet flow region and a recirculating flow region. The mean kinetic energy decays rapidly in the jet region and the dissipation rate in most of the areas in the pool is less than 200?W/m3. For the jet flow, the nondimensional mean velocity profile across the jet agrees very well with that of a plane turbulent jet in the central part of the jet with some scatter near its boundaries. Its maximum velocity decays faster compared to a plane turbulent jet in a large stagnant ambient. The jet presents different turbulence structure for the two flow patterns and for each pattern, the turbulence characteristics appear different between the left and right halves of the jet. However, the turbulence characteristics show some similarity for each case. The normalized energy dissipation rate shows some similarity and has a maximum value on the center of the jet. The results are believed to provide useful insight on the turbulence characteristics of flow in vertical slot fishways and can be used to verify numerical models and also for guidance in the design of fishways in the future.     

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.     

Turbulence Structure of Hydraulic Jumps of Low Froude Numbers

Turbulence characteristics of hydraulic jumps with Froude numbers of 2.0, 2.5, and 3.32 are presented. A Micro Acoustic Doppler velocimeter was used to obtain measurements of the velocities, turbulence intensities, Reynolds stresses, and power spectra. The maximum turbulence intensities and Reynolds stress at any section were found to decrease rapidly from the toe of the jump towards downstream within the jump and then gradually level off in the transition region from the end of the jump to the friction dominated open channel flow downstream. The maximum turbulence kinetic energy at each section decreases linearly with the longitudinal distance within the jump and gradually levels off in the transition region. The Reynolds stress and turbulence intensities within the jump show some degree of similarity. The dissipative eddy size was estimated to vary from 0.04 mm within the jump to 0.15 mm at the end of the transition region. The dominant frequency is in the range from 0 to 4 Hz for both horizontal and vertical velocity components.     

Characteristics of Submerged Jets in Evolving Scour Hole Downstream of an Apron

This paper reports an experimental investigation on the velocity and turbulence characteristics in an evolving scour hole downstream of an apron due to submerged jets issuing from a sluice opening detected by an acoustic Doppler velocimeter. Experiments were carried out for the conditions of submerged jets, having submergence factors from 0.96 to 1.85 and jet Froude numbers from 2.58 to 4.87, over sediment beds downstream of a rigid apron. The distributions of time-averaged velocity vectors, turbulence intensities, and Reynolds stress at different streamwise distances are plotted for the conditions of initial flat bed, intermediate scour holes, and equilibrium scour hole downstream of an apron. Vector plots of the flow field show that the rate of decay of the submerged jet velocity increases with an increase in scour hole dimension. The bed-shear stresses are determined from the Reynolds stress distributions. The flow characteristics in evolving scour holes are analyzed in the context of self-preservation, growth of the length scale, and decay of the velocity and turbulence characteristics scales. The most significant observation is that the flow in the scour holes (intermediate and equilibrium) is found to be plausibly self-preserving.     

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.     

Application of Several Depth-Averaged Turbulence Models to Simulate Flow in Vertical Slot Fishways

Vertical slot fishways are hydraulic structures which allow the upstream migration of fish through obstructions in rivers. The velocity, water depth, and turbulence fields are of great importance in order to allow the fish swimming through the fishway, and therefore must be considered for design purposes. The aim of this paper is to assess the possibility of using a two-dimensional shallow water model coupled with a suitable turbulence model to compute the flow pattern and turbulence field in vertical slot fishways. Three depth-averaged turbulence models of different complexity are used in the numerical simulations: a mixing length model, a k?ε model, and an algebraic stress model. The numerical results for the velocity, water depth, turbulent kinetic energy, and Reynolds stresses are compared with comprehensive experimental data for three different discharges covering the usual working conditions of vertical slot fishways. The agreement between experimental and numerical data is very satisfactory. The results show the importance of the turbulence model in the numerical simulations, and can be considered as a useful complementary tool for practical design purposes.     

Characteristics of Velocity and Excess Density Profiles of Saline Underflows and Turbidity Currents Flowing over a Mobile Bed

Turbidity currents in the ocean and lakes are driven by suspended sediment. The vertical profiles of velocity and excess density are shaped by the interaction between the current and the bed as well as between the current and the ambient water above. We present results of a set of 74 experiments that focus on the characteristics of velocity and fractional excess density profiles of saline density and turbidity currents flowing over a mobile bed. The gravity flows include saline density flows, hybrid saline/turbidity currents and a pure turbidity current. The use of dissolved salt is a surrogate for suspended mud that is so fine that it does not settle out readily. Thus, all the currents can be considered to be model turbidity currents. The data cover both Froude-subcritical and Froude-supercritical regimes. Depending on flow conditions, the bed remains flat or bed forms develop over time, which in turn affect vertical profiles. For plane bed experiments, subcritical flow profiles have velocity peaks located higher up in the flow, and display a sharper interface at the top of the current, than their supercritical counterparts. The latter have excess density profiles that decline exponentially upward from the bed, whereas subcritical flows show profiles with a region near the bed where excess density varies little. Wherever bed forms are present, they have a significant effect on the profiles. Especially for Froude-supercritical flow, bed forms push the location of peak velocity upward, and render the near-bed fractional excess density more uniform. In the case of subcritical flow, bed forms do not significantly affect fractional excess density profiles; velocity profiles are pushed farther upward from the bed than in the case of a plane bed, but to a lesser extent than for supercritical bed forms. Overall, the relative position of the velocity peak above the bed shows a dependence upon flow regime, being lowered for increasing Froude number Fd. Gradient Richardson numbers Rig in the near-bed region increase with increasing Fd, but are lower than the critical value of 0.25, indicating that near-bed turbulent structures are not notably suppressed. At the top interface, values of Rig are above the critical value for subcritical and mildly supercritical Fd, effectively damping turbulence. However as Fd increases, Rig goes below the critical value. Shape factors calculated from the profiles for use in the depth-averaged equation of motion are evaluated for different flow and bed conditions. Normalized experimental profiles for supercritical currents scale up well with observations of field-scale turbidity currents in the Monterey Canyon, and the range of average bed slopes and Froude numbers also compare favorably with estimated field-scale flow conditions for the Amazon canyon and fan. This suggests that the experimental results can be used to interpret the kinds of flows that are responsible for the shaping of major submarine canyon-fan systems.