Numerical Investigation into Secondary Currents and Wall Shear in Trapezoidal Channels

This paper presents the use of computational fluid dynamics (CFD) to determine the distribution of the bed and sidewall shear stresses in trapezoidal channels. The impact of the variation of the slant angle of the side walls, aspect ratio, and composite roughness on the shear stress distribution is analyzed. The shear stress data can be directly output from the CFD models at the boundaries, but they can also be derived using the Guo and Julien equations for the average bed and side wall shear stresses. These equations compute the shear stress as a function of three components; gravitational, secondary flows, and interfacial shear stress, and are hence used to gauge the respective merits of the different components of wall shear. The results show a significant contribution from the secondary currents and internal shear stresses on the overall shear stress at the boundaries. This work also extends previous work of the authors on rectangular channels.     

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.     

Comparisons of Sidewall Correction of Bed Shear Stress in Open-Channel Flows

When investigating sediment transport in laboratory open-channel flows, it is often necessary to remove sidewall effects for computing effective bed shear stress. Previous sidewall correction methods are subject to some assumptions that have not been completely verified, and different values of the bed shear stress may be obtained depending on the approach used in making sidewall corrections. This study provides a quantitative assessment of the existing correction procedures by comparing them to a new sidewall correction model proposed in this study. The latter was derived based on the shear stress function and equivalent roughness size for both rigid and mobile bed conditions, which were obtained directly from experimental measurements. The comparisons show that the Einstein correction formula and the Vanoni and Brooks method generally predict relatively lower and higher bed shear stresses, respectively, while the Williams’ empirical function leads to more scatter. This study also demonstrates that the widely used Vanoni and Brooks approach can be well approximated by a simple formula derived based on the Blasius resistance function. The sidewall effects, when removed in the different ways, would consequently affect the presentation of the bedload function. Experimental results of bedload transport, when plotted as the dimensionless transport rate against the dimensionless shear stress with the latter being corrected using the present model, exhibit less scatter than those associated with the previous procedures.     

Flow Resistance over Mobile Bed in an Open-Channel Flow

This paper deals with the underlying mechanism of flow resistance in an alluvial channel: The effects of sidewall and bed form on flow resistance, Einstein’s divided hydraulic radius approach and Engelund’s energy slope division approach are reexamined. These two approaches assume that the shear stress on a mobile bed is the summation of shear stresses caused by skin friction and bed form. Using a different approach, this paper presents a theoretical relationship between the total bed shear stress with grain and bed-form shear stresses. The contribution of sidewall on the total bed shear stress is also discussed. The writers found that the size of bed form plays a significant role for the flow resistance, and developed relevant expressions for the length of the separation zone behind the bed forms. In addition, a systematical approach has been developed to compute the flow velocity in an alluvial channel. This approach is tested and verified against 5,989 flume and field measurements. The computed and measured discharge/velocity are in good agreement and 83.0% of all data sets fall within the ±20% error band.     

Turbulence Characteristics in Gradual Channel Transition

A field study was conducted to determine the effects of a channel transition on turbulence characteristics. Detailed three-dimensional (3D) flow measurements were collected at a cross section that is located downstream of a gradual channel expansion. These measurements were obtained via an acoustic doppler velocimeter and include the 3D velocity field, the mean local velocities, the turbulent intensities, the frictional characteristics of the flow, the secondary velocity along the transverse plane, and the instantaneous shear stress components in the streamwise and transverse directions. Analysis of the 3D flow data indicates that the turbulent flow on the outer bank of the channel is anisotropic. Such anisotropy of turbulence, which is attributed to the gradual expansion in the channel and bed roughness, yields the development of a secondary flow of Prandtl’s second kind as reported in 1952. In particular, it was found that turbulent intensities in the vertical and transverse directions on the outer bank section are different in magnitude creating turbulence anisotropy in the cross-sectional plane and secondary flows of the second kind. Turbulent intensities increase toward the free surface indicating the transfer of a higher-momentum flux from the channel bed to the free surface, which contradicts common wisdom. Results for the normalized stress components in the streamwise and transverse direction show similar behavior to the intensities. Moreover, the nonlinear distribution of stresses is indicative of the oscillatory nature of the flow induced by the secondary flows of Prandtl’s second kind. A similar behavior was found for flows in straight rectangular channels over different roughness. Finally, a comparison between the secondary current velocity with the mainstream velocity indicates that secondary flow of Prandtl’s second kind is present within the right half of the measured cross section.     

Distribution of Bed Shear Stress in Rotating Circular Flume

Distributions of bed shear stress across the width of a rotating circular flume with smooth and rough bed surfaces were obtained by measurement and model prediction. Results with flows over smooth beds showed that the flow in the central part may be considered to be two-dimensional and that effects of flow depth over the operating range of the flume are minor for flow depths not exceeding 0.14 m. For rough beds, the bed shear stress distributions were found to be skewed toward the inner wall. This can be corrected if a compensating roughness is added to the bottom of the ring. Such measures are also effective for flumes with smooth beds. Measured bed shear stress distributions agreed well with the predicted distributions for smooth beds and reasonably well for rough beds. The modified Preston tube, for measurement of bed shear stress in flows over rough beds, was found to give promising results. Further tests are required to completely define the uncertainty in bed shear stress measurements made with this instrument.     

Modeling Depth-Averaged Velocity and Boundary Shear in Trapezoidal Channels with Secondary Flows

The Shiono and Knight method (SKM) offers a new approach to calculating the lateral distributions of depth-averaged velocity and boundary shear stress for flows in straight prismatic channels. It accounts for bed shear, lateral shear, and secondary flow effects via 3 coefficients—f,λ, and Γ—thus incorporating some key 3D flow feature into a lateral distribution model for streamwise motion. The SKM incorporates the effects of secondary flows by specifying an appropriate value for the Γ parameter depending on the sense of direction of the secondary flows, commensurate with the derivative of the term Hρ(UV)d. The values of the transverse velocities, V, have been shown to be consistent with observation. A wide range of boundary shear stress data for trapezoidal channels from different sources has been used to validate the model. The accuracy of the predictions is good, despite the simplicity of the model, although some calibration problems remain. The SKM thus offers an alternative methodology to the more traditional computational fluid dynamics (CFD) approach, giving velocities and boundary shear stress for practical problems, but at much less computational effort than CFD.     

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.     

Boundary Shear Stress in Spatially Varied Flow with Increasing Discharge

The distribution of the wall shear stress on the bed and sidewalls of an open channel receiving lateral inflow was obtained from experimental measurements of the distribution of the velocity in the viscous sublayer using a laser doppler velocimeter. The experiments were conducted in a 0.4 m wide by 7.5 m long flume. Lateral inflow was provided into the channel from above via sets of nozzles positioned toward the downstream end of the flume. Lateral inflow was provided over a length of 1.9 m. The results indicate that the local boundary shear stresses are significantly influenced by lateral inflow. The significant variation occurs near and around the region where the lateral inflow enters the channel. At various measurement positions along the lateral inflow zone, mean boundary, mean wall, and mean bed shear stresses were obtained and compared. The results indicate that the mean boundary shear stresses increase from the upstream to the downstream ends of the lateral inflow zone. The results also indicate that the mean bed shear stress is always greater than the mean wall shear stress, which are approximately 30–60% of the mean bed shear stress. The friction factor in the Darcy–Weisbach equation was obtained from both the mean boundary shear stress and from the equation describing the water surface elevation in an open channel receiving lateral inflow (equation for spatially varied flow with increasing discharge). The results indicate that the estimated friction factors from the latter approach are significantly larger. Also, the estimated friction factors from both approaches are higher than the values predicted from the Blasius equation which describes the friction factor for wide uniform open channel flows. They were also higher than values predicted from the Keulegan equation, which is an empirically derived equation for flow in roof drainage gutters. The study highlights the deficiencies in the existing equations used to predict friction factors for spatially varied flow and that further research is required to explore the distribution of boundary shear stress in an open channel receiving lateral inflow.     

Semianalytical Model for Shear Stress Distribution in Simple and Compound Open Channels

Semianalytical equations were derived for distribution of shear stress in straight open channels with rectangular, trapezoidal, and compound cross sections. These equations are based on a simplified streamwise vorticity equation that includes secondary Reynolds stresses. Reynolds stresses were then modeled and their different terms were evaluated based on the work of previous researchers and experimental data. Substitution of these terms into the simplified vorticity equation yielded the relative shear stress distribution equation along the width of different channel cross sections. In compound channels the effect of additional secondary flows due to the shear layer between the main channel and the flood plain were also considered. Comparisons between predictions of the model and experimental data, predictions of other analytical or three dimensional numerical models with advanced turbulent closures, were made with good agreement.     

Lateral Depth-Averaged Velocity Distributions and Bed Shear in Rectangular Compound Channels

This paper presents a method to predict depth-averaged velocity and bed shear stress for overbank flows in straight rectangular two-stage channels. An analytical solution to the depth-integrated Navier–Stokes equation is presented that includes the effects of bed friction, lateral turbulence, and secondary flows. The Shiono and Knight method accounts for bed shear, lateral shear, and secondary flow effects via three coefficients, f, λ, and Γ, respectively. A novel boundary condition at the internal wall between the main channel and the adjoined floodplain is proposed and discussed along with other conventional boundary conditions. The analytical solution using the novel boundary condition gives good prediction of both lateral velocity distribution and bed shear stress when compared with experimental data for different aspect ratios.     

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.     

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.     

Unsteady Bed Shear Stresses Induced by Navigation: Laboratory Observations

Time-dependent bed shear stresses induced by the passage of a barge tow have been measured with hot film shear stress sensors in a 1:25 scale model. Conditions typical of those observed for Upper Mississippi River navigation traffic were simulated in the experimental facility. Two sets of experiments were carried out: the first set consisted of simultaneous shear stress measurements at different locations for a variety of flow depths and boat operating conditions, providing space-time distributions of ensemble averaged wall shear stresses. The second set included a large number of realizations gathered for one particular flow condition at a single position, allowing analysis of the time evolution of the turbulence characteristics (i.e., standard deviation) of the bed shear stresses. The results of the first set of experiments show that for all the experimental conditions the basic patterns of the shear stress are similar, with two regions of high shear stress associated with the passage of the bow and the stern of the barge tow, respectively. Analysis of the second set of experiments showed that, as a result of the passage of the barge tow, the bed-shear stress standard deviation departs from the values commonly observed under steady, uniform, open-channel flow conditions. This behavior has important implications for sediment transport.     

Boundary Shear in Circular Pipes Running Partially Full

The distribution of boundary shear stress in circular conduits flowing partially full, with and without a smooth flat bed simulating deposited sediments, has been examined experimentally ranging from 0.375 < F < 1.96 and 6.5 × 104 < R < 3.42 × 105, using the Preston tube technique The invert level of the flat bed and the water depth have been varied to simulate a wide range of possible flow conditions that may occur in culverts, sewers, and hydropower tunnels. The distribution of boundary shear stress around the wetted perimeter is shown to be highly sensitive to changes in cross-sectional shape. The results have been analyzed in terms of the variation of local∕global shear stress versus perimetric distance, and the percentage of the total shear force acting on the wall or bed of the conduit. The %SFw results have been shown to agree well with Knight's empirical formula for prismatic channels. The influence of secondary flows on the distribution of boundary shear stress and the implications of this for sediment transport have also been examined.     

Simplified Model of Tractive-Force Distribution in Closed Conduits

A simplified model for the computation of boundary shear stress distributions acting on the flow perimeter of closed ducts is presented. The model assumes that the surplus energy within any control volume in a three-dimensional flow will be transferred towards the nearest boundary to be dissipated. Based on this model, the flow cross sectional area in the closed duct is divided into subflow regions corresponding to the side walls and bed, and the shear distributions over the wetted perimeter within each subflow area are assessed. Analytical equations, valid for all channel aspect ratios, for the prediction of local and mean shear stresses along the bed and side walls in smooth rectangular duct flow are derived. The formulae give good predictions of the shear stress distributions when compared with existing experimental data in the literature. The possible applications of the model to nonrectangular duct flows are also discussed.     

Flow Field at a Vertical-Wall Abutment

Experiments were conducted to measure the three-dimensional turbulent flow field, using the acoustic Doppler velocimeter, at a short vertical-wall abutment (ratio of abutment length to approach flow depth less than unity) before and after the development of a scour hole under a clear water scour condition. In the upstream, the presentation of flow field through vectors at vertical sections shows a primary vortex associated with the downflow. In the downstream, the upward flow is comprised of with irregularities owing to the vortex shedding. The flow separation near the bed and within the scour hole is evident from the turbulent kinetic energy distribution. Using Reynolds stresses, the bed shear stresses are calculated.     

The state of residual stress in the near surface region of homogeneous and heterogeneous materials after grinding

The residual stresses near the surface were measured after grinding under different conditions with X-rays of various wavelengths using various peaks from homogeneous and heterogeneous materials with fcc-and bcc-structures. The strain distributions were subject toψ-splitting. A method was elaborated to measure the variation of stress as a function of depth with X-rays of various wavelengths. For two ground steels it was shown that this method gave stress distributions with depth which were in good agreement with those evaluated from stress measurements with one wavelength and layer removal. Normal-and shear stresses decreased absolutely with increasing depth. The existence of shear residual stresses was discussed and it was concluded that for heterogeneous materials these stresses can be compensated by those of opposing sign in the other phases.     

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.     

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.