Effect of Sediment Size Scaling on Physical Modeling of Bridge Pier Scour

Local pier scour experiments were performed in the laboratory to investigate the effect of relative sediment size on pier scour depth using three uniform sediment sizes and three bridge pier designs at different geometric model scales. When the data from a large number of experimental and field investigations are filtered according to a Froude number criterion, the effect of relative sediment size on dimensionless pier scour depth is brought into focus. The choice of sediment size in the laboratory model distorts the value of the ratio of pier width to sediment size in comparison with the prototype which in turn causes larger values of scour depth in the laboratory than in the field. This model distortion due to sediment size is shown to be related to the scaling of the large-scale unsteadiness of the horseshoe vortex by studying the relevant time scales of its coherent structure upstream of a bridge pier using acoustic Doppler velocimeter measurements. Observations of sediment movement, probability distributions of velocity components, and phase-averaging of velocity measured upstream of a bridge pier reveal properties of coherent motions that are discussed in terms of their contribution to the relationship between dimensionless pier scour depth and the ratio of pier width to sediment size over a large range of physical scales.     

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.     

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.     

Detached Eddy Simulation Investigation of Turbulence at a Circular Pier with Scour Hole

This paper uses results from detached eddy simulation to reveal the dynamics of large-scale coherent eddies in the flow around a circular pier with an equilibrium scour hole. This is important for the sediment transport because the local scour process is controlled to a large extent by the large-scale coherent structures present in the near-bed region. The present paper investigates the dynamics of these coherent structures, their interactions and their role in entraining sediment in the later stages of the scour process when the horseshoe vortex system is stabilized by the presence of a large scour hole. The pier Reynolds number was 2.06×105, outside the range of well-resolved large-eddy simulation (LES). Additionally, scale effects are investigated based on comparison with LES results obtained at a much lower Reynolds number of 16,000 in a previous investigation. The paper provides a detailed study of the dynamics of the main necklace vortices of the horseshoe vortex system, including an investigation of the bimodal oscillations, their effect on the amplification of the turbulence within the scour hole and the interactions of the necklace vortices with the downflow. Several mechanisms for the growth of the downstream part of the scour hole in the later stages of the scour process are discussed. Similar to the low-Reynolds-number simulation, and consistent with experimental observations, the presence of strong upwelling motions near the symmetry plane resulted in the suppression of the large-scale vortex shedding in the wake. The fact that the nondimensional values of the turbulent kinetic energy and pressure RMS fluctuations in the higher Reynolds number simulation were consistently lower inside the regions of high turbulence amplification associated with the main necklace vortex, the separated shear layers and the near-wake shows that changes in the flow and turbulence due to the Reynolds number and scour hole geometry can be quantitatively significant over Reynolds numbers between 104 and 105.     

Particle Image Velocimetry Study of Flow near Trashrack Models

This paper provides results of an experimental study of turbulent flow near trashrack models that are comprised of an array of three rectangular bars. The bar thickness, bar depth, and center-to-center spacing were maintained constant. The flow characteristics were studied by aligning the bars with the approach flow and conducting measurements at three different approach freestream velocities. Subsequently, the freestream velocity was kept constant and detailed measurements were conducted for four different bar inclinations relative to the approach flow. For each test condition, a high-resolution particle image velocimetry (PIV) technique was used to conduct detailed velocity measurements in streamwise-spanwise planes at middepth of flow. From these measurements, isocontours and profiles of the mean velocities, turbulence intensities, Reynolds shear stress, and production term in the transport equation for the turbulent kinetic energy were obtained to study the flow characteristics around and downstream of the aligned and inclined bars. Flow characteristics near hydroelectric station trashracks are important for efficient turbine operation and reduction of fish entrainment.     

Three-Dimensional Numerical Model for Flow and Bed Deformation around River Hydraulic Structures

This paper describes a numerical model developed to simulate flow and bed deformation around river hydraulic structures. The model solved the fully three-dimensional, Reynolds-averaged Navier–Stokes equation expressed in a moving boundary-fitted coordinate system to calculate the flow field with water and bed surfaces varying in time. A nonlinear k-ε turbulence model was employed in order to predict flow near the structure where three-dimensional flow is dominant. The temporal change in bed topography was calculated by coupling a stochastic model for sediment pickup and deposition using a momentum equation of sediment particles in order to account for the effect of nonequilibrium sediment transport. In validating the numerical model, a spur dike and a bridge pier, which are considered to be typical river-engineering structures, were selected. By comparing the numerical results with observed laboratory experimental data, the model was found to reproduce flow and scour geometry around these structures with sufficient accuracy.     

Bridge Pier Scour under Flood Waves

The effect of a single-peaked flood wave on pier scour is investigated both theoretically and experimentally. The conditions considered involve clear-water scour of a cohesionless material of given median sediment size and sediment nonuniformity, an approach flow characterized by a flow depth and velocity, a circular-shaped cylindrical bridge pier, and a flood hydrograph defined by its time to peak and peak discharge. A previously proposed formula for scour advance under a constant discharge was applied to the unsteady approach flow. The generalized temporal scour development along with the end scour depth are presented in terms of mainly the densimetric particle Froude number based on the maximum approach flow velocity and the median sediment size. The effect of the remaining parameters on the end scour depth is discussed and predictions are demonstrated to be essentially in agreement with model observations.     

Characteristics of Horseshoe Vortex in Developing Scour Holes at Piers

The outcome of an experimental study on the turbulent horseshoe vortex flow within the developing (intermediate stages and equilibrium) scour holes at cylindrical piers measured by an acoustic Doppler velocimeter (ADV) are presented. Since the primary objective was to analyze the evolution of the turbulent flow characteristics of a horseshoe vortex within a developing scour hole, the flow zone downstream of the pier was beyond the scope of the investigation. Experiments were conducted for the approaching flow having undisturbed flow depth ( = 0.25?m) greater than twice the pier diameter and the depth-averaged approaching flow velocity ( = 0.357?m/s) about 95% of the critical velocity of the uniform bed sand that had a median diameter of 0.81?mm. The flow measurements by the ADV were taken within the intermediate (having depths of 0.25, 0.5, and 0.75 times the equilibrium scour depth) and equilibrium scour holes (frozen by spraying glue) at a circular pier of diameter 0.12?m. In order to have a comparative study, the ADV measurements within an equilibrium scour hole at a square pier (side facing the approaching flow) of sides equaling the diameter of the circular pier were also taken. The contours of the time-averaged velocities, turbulence intensities, and Reynolds stresses at different azimuthal planes (0, 45, and 90°) are presented. Vector plots of the flow field at azimuthal planes reveal the evolution of the characteristics of the horseshoe vortex flow associated with a downflow from intermediate stages to equilibrium condition of scour holes. The bed-shear stresses are determined from the Reynolds stress distributions. The flow characteristics of the horseshoe vortex are discussed from the point of view of the similarity with the velocity and turbulence characteristic scales. The imperative observation is that the flow and turbulence intensities in the horseshoe vortex flow in a developing scour hole are reasonably similar.     

Evaluation of Mean Velocity and Turbulence Measurements with ADCPs

To test the ability of acoustic Doppler current profilers (ADCPs) to measure turbulence, profiles measured with two pulse-to-pulse coherent ADCPs in a laboratory flume were compared to profiles measured with an acoustic Doppler velocimeter, and time series measured in the acoustic beam of the ADCPs were examined. A four-beam ADCP was used at a downstream station, while a three-beam ADCP was used at a downstream station and an upstream station. At the downstream station, where the turbulence intensity was low, both ADCPs reproduced the mean velocity profile well away from the flume boundaries; errors near the boundaries were due to transducer ringing, flow disturbance, and sidelobe interference. At the upstream station, where the turbulence intensity was higher, errors in the mean velocity were large. The four-beam ADCP measured the Reynolds stress profile accurately away from the bottom boundary, and these measurements can be used to estimate shear velocity. Estimates of Reynolds stress with a three-beam ADCP and turbulent kinetic energy with both ADCPs cannot be computed without further assumptions, and they are affected by flow inhomogeneity. Neither ADCP measured integral time scales to within 60%.     

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.     

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.     

Aspects of Inlet Boundary Prescription for a Turbulent Flow Field

In the literature on turbulent flow various combinations of velocity, turbulence kinetic energy, and eddy viscosity models have been proposed for the inlet boundary to a flow field. There appears to be no rational criterion for specifying inlet boundary conditions. The present study proposes a criterion to select the inlet boundary conditions by treating the inlet boundary as a part of the flow field. Using this criterion, any prescribed variation of velocity, turbulence kinetic energy k, and rate of dissipation of kinetic energy ε, must satisfy the governing flow field equations at the inlet boundary. Analysis of previously used profiles of velocity, kinetic energy, and dissipation of kinetic energy in the literature indicates that most of these do not satisfy the flow field equations. To substantiate the importance of inlet boundary conditions, several reported numerical simulations using k-ε turbulence models are reconsidered to determine if there is any linkage between the residual errors at the inlet boundary and the errors in the flow field simulation. Based on such an analysis, it is both logical and practical to hypothesize the inlet boundary as a part of flow field.     

Turbulent Flow near Vertical Angled Fish Screen

Mean and turbulent flow characteristics on the upstream and downstream sides of the screen in a flow diversion channel have important implications for operation and maintenance (e.g., sedimentation) and for assessing fish behavior related to flow turbulence. This technical note extends an earlier study on mean flow near screens to turbulence characteristics. Acoustic Doppler velocimeter was used to explore three-dimensional mean and turbulent flow characteristics on the upstream and downstream sides of vertical angled fish screens. The present study confirms the two-dimensional mean velocity observations of the previous experimental work and shows that the vertical mean velocities are less than 10% of the local magnitudes of longitudinal velocity and hence can be ignored. Horizontal components of the mean velocity on the downstream side of the screen were relatively small, but the turbulent velocity fluctuations were two to three times as intense as those measured on the upstream side.     

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.     

Temporal Evolution of Plunge Pool Scour

The temporal development of plunge pool scour was investigated using a novel experimental approach. Longitudinal profiles along the scour hole were recorded with an optical method to allow its definition at any time, from the initiation of scour to nearly the end-scour condition. The characteristics of the scour hole geometry were investigated, namely the maximum scour hole depth, the maximum ridge height, and their locations relative to the scour hole origin. It is demonstrated that the evolution is logarithmic, similar to that found for bridge pier and abutment scour. A distinction is further made between the developing and the developed scour hole phases. The main issue of the present research was to define the developed scour hole characteristics because the developing scour phase is influenced by turbulence features that may be difficult to assess. This work therefore allows for an appreciation of the temporal evolution of a scour process of engineering interest.     

Local Scour and Riprap Stability at Bridge Piers in a Degrading Channel

The experimental study examines local scouring and riprap stability at bridge piers in rivers subject to bed degradation. The data show that the equilibrium bed profile associated with that with or without a pier is essentially the same, except for the obvious section around the pier. Total scour depth is shown to be the sum of bed degradation and pier scour depth. The latter can be computed from the time-average live-bed scour depth associated with the undisturbed velocity ratio before bed degradation. The experimental data also show that pier-scour depth is invariant with time, for t ≥ 24?h. In a degrading channel, riprap around a pier will eventually develop into a stable mound when the bed shear stresses reduce with bed degradation. An auxiliary test shows that the mound is very vulnerable to another designed flood flow accompanied by large dunes. This type of riprap instability may be called bed-degradation induced failure.     

Numerical Simulation of Swirling Flow in Complex Hydroturbine Draft Tube Using Unsteady Statistical Turbulence Models

A numerical method is developed for carrying out unsteady Reynolds-averaged Navier-Stokes (URANS) simulations and detached-eddy simulations (DESs) in complex 3D geometries. The method is applied to simulate incompressible swirling flow in a typical hydroturbine draft tube, which consists of a strongly curved 90° elbow and two piers. The governing equations are solved with a second-order-accurate, finite-volume, dual-time-stepping artificial compressibility approach for a Reynolds number of 1.1 million on a mesh with 1.8 million nodes. The geometrical complexities of the draft tube are handled using domain decomposition with overset (chimera) grids. Numerical simulations show that unsteady statistical turbulence models can capture very complex 3D flow phenomena dominated by geometry-induced, large-scale instabilities and unsteady coherent structures such as the onset of vortex breakdown and the formation of the unsteady rope vortex downstream of the turbine runner. Both URANS and DES appear to yield the general shape and magnitude of mean velocity profiles in reasonable agreement with measurements. Significant discrepancies among the DES and URANS predictions of the turbulence statistics are also observed in the straight downstream diffuser.     

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.     

Flow Velocity and Pier Scour Prediction in a Compound Channel: Big Sioux River Bridge at Flandreau, South Dakota

The two-dimensional (2D) depth-averaged river model Finite-Element Surface-Water Modeling System (FESWMS) was used to predict flow distribution at the bend of a compound channel. The site studied was the Highway 13 bridge over the Big Sioux River in Flandreau, South Dakota. The Flandreau site has complex channel and floodplain geometry that produces unique flow conditions at the bridge crossing. The 2D model was calibrated using flow measurements obtained during two floods in 1993. The calibrated model was used to examine the hydraulic and geomorphic factors that affect the main channel and floodplain flows and the flow interactions between the two portions. A one-dimensional (1D) flow model of the bridge site was also created in Hydrologic Engineering Centers River Analysis System (HEC-RAS) for comparison. Soil samples were collected from the bridge site and tested in an erosion function apparatus (EFA) to determine the critical shear stress and erosion rate constant. The results of EFA testing and 2D flow modeling were used as inputs to the Scour Rate in Cohesive Soils (SRICOS) method to predict local scour at the northern and southernmost piers. The sensitivity of predicted scour depth to the hydraulic and soil parameters was examined. The predicted scour depth was very sensitive to the approach-flow velocity and critical shear stress. Overall, this study has provided a better understanding of 2D flow effects in compound channels and an overall assessment of the SRICOS method for prediction of bridge pier scour.     

Predicting Equilibrium Scour-Hole Geometry near Angled Stream Deflectors Using a Three-Dimensional Numerical Flow Model

Stream rehabilitation projects often involve the installation of instream structures such as flow deflectors. The objective of this study is to use the output of a three-dimensional numerically simulated flow field over a flat, predeformation bed to predict the planform extent of the equilibrium scour hole near stream deflectors of varying angles. It is shown that the upstream extent of the simulated flow separation zone determines the upstream limit of scour, whereas the lateral scour extent at the nose of the deflector is determined by the width of the separation zone. Further, scour depths are greatest in regions where strong downwelling (negative vertical velocity) exists, and the position of the local minimum dynamic pressure point in the simulated flow field defines the downstream limit of scouring. The results have implications for future design of habitat improvement structures where different angles and lengths of structures could potentially be tested prior to their implementation to determine the resultant scour geometry.