Comparison of Three Commonly Used Equations for Calculating Local Scour Depth around Bridge Pier under Ice Covered Flow Condition

A precise prediction of maximum scour depth around bridge foundations under ice covered condition is crucial for their safe design because underestimation may result in bridge failure and over-estimation will lead to unnecessary construction costs. Compared to pier scour depth predictions within an open channel, few studies have attempted to predict the extent of pier scour depth under ice-covered condition. The present work examines scour under ice by using a series of clear-water flume experiments employing two adjacent circular bridge piers in a uniform bed were exposed to open channel and both rough and smooth ice covered channels. The measured scour depths were compared to three commonly used bridge scour equations including Gao’s simplified equation, the HEC-18/Jones equation, and the Froehlich Design Equation. The present study has several advantages as it adds to the understanding of the physics of bridge pier scour under ice cover flow condition, it checks the validity and reliability of commonly used bridge pier equations, and it reveals whether they are valid for the case of scour under ice-covered flow conditions. In addition, it explains how accurately an equation developed for scour under open channel flow can predict scour around bridge piers under ice-covered flow condition.


Introduction
Scouring can be defined as the process by which the particles of soil or rock around an abutment or pier of a bridge get eroded and removed to a certain depth (called the scour depth) [1].The U.S. Federal Highway Administration (FHWA) has estimated 60% of bridge failure cases in the USA are due to scour and on average, approximately 50 to 60 bridges collapse annually in the USA [2].
Wardhana & Hadipriono [3] studied 500 failures of bridge structures in the United States between 1989 and 2000 and reported that the most frequent causes of bridge failures were due to floods, scour, and their cumulative impact.The average age of the 500 failed bridges was 52.5 years but ranged from 1 to 157 years old [4].Brice & Blodgett [5] reported that damages to bridges and highways from major regional floods in 1964 and 1972 were equivalent to approximately 100 million US dollars per event.Bridge foundations should be designed to withstand the effects of scour.Bridge damage and failure have huge negative social and economic impacts in terms of reconstruction costs, maintenance and monitoring of existing structures, the disruptions of traffic flow, and in some life-threatening cases, the cost of human lives [4].Moreover, a precise prediction of scour depth will not only help to prevent those bridge failures which are the consequence of under-estimation of scour depth but also will efficiently reduce unnecessary construction cost of those bridge piers in which scour depths are over-estimated.To safely design bridges located on waterways under severe flooding conditions, many researchers have developed a number of laboratory-derived equations for predicting bridge pier scour depth [6] [11] [12].These equations are mostly empirical formulae which are usually based on regressional analysis of laboratory and/or field scour data.However, they differ from each other in terms of the factors considered in constructing the scour model, parameters used in the equation, laboratory and/or site conditions.
Since the number of these equations is relatively large, selection of the best performing equations for a special case is a difficult task.Comparison studies of scour formulae especially for different flow conditions might be helpful to select the one with the best performance.Additionally, many rivers become ice-covered during the winter months.However, the winter season is often overlooked even though most rivers in Canada and northern parts of the United States, Europe, and Asia are annually affected by ice.The relatively smaller number of studies on the scour around bridge pier under ice-covered flow condition is due to the inherent difficulty in collecting field data while ice is present and complications in lab-based measurements as a result of different scales and of temperature effects [13].Ice cover can significantly change the flow field and impact sediment transport in natural rivers.The formation of a stable ice cover effectively doubles the wetted perimeter compared to open channel conditions.This alters the hydraulics of the channel by imposing an extra boundary to the flow, causing the velocity profile to be shifted towards a smoother boundary (channel bed) and adding to the flow resistance [14].Furthermore, ice cover can lead to issues such as ice jamming, flooding, restricting the generation of hydro-power, blocking river navigation, and affecting the overall ecosystem bal-

Materials and Methods
Most of the equations for the prediction of bridge pier scouring express the final scour depth as a function of the flow characteristics (mean flow velocity at the approach section, water depth), flow properties (density and viscosity of the fluid), stream bed material properties (mean particle diameter, density) and bridge geometry (shape and dimension of the pier, angle of attack of the flow).In this paper, scour around circular bridge piers will be experimentally examined and subsequently the validity and reliability of three of the more commonly used and cited scour equations developed specifically for open channel flow condition will be investigated to see how accurately they predict scouring around bridge piers under ice-covered flow conditions.

HEC-18/Jones Equation
The most commonly used pier scour equation in the United States is the Colorado State University (CSU) equation proposed by Richardson and Davis [18] and is recommended by the U.S. Department of Transportation's Hydraulic Engineering (HEC-I8) (1993).It was developed from laboratory data and is recommended for both live-bed and clear-water conditions.The HEC-18/Jones equation is based on the Colorado State University (CSU) equation: where y s = scour depth; y 0 = the approach flow depth; y s /y 0 is a dimensionless expression of the relative scour depth with respect to flow depth; K 1 = correction factor for pier nose shape which is unity for circular cylinder; K 2 = correction factor for angle of attack flow which is unity for 90 0 ; K 3 = correction factor for bed condition which is 1.1 for clear water scour; b = nominal pier width; and Fr = approach flow Froude number.K 4 is a correction factor to account for armoring World Journal of Engineering and Technology of the scour hole: where V R is the velocity ratio and is dimensionless: where V o is the approach velocity directly upstream from the pier and V i50 is the approach velocity, in feet per second, required to initiate scour at the pier for the particle size D 50 .V i50 is calculated as follows: where D 50 is the particle size for which 50 percent of the bed material is finer, in units of feet and V c50 is the critical velocity, in feet per second, for incipient motion of the particle size D 50 .V c50 is defined as follows: 1/6 1/3 50 0 50 11.21 while D 90 is the particle size for which 90 percent of the bed material is finer, in units of feet, and V c90 is the critical velocity, in feet per second, for the incipient motion of the particle size as given by:

Gao's Simplified Equation
Gao's simplified pier scour equation is based on laboratory and field data from China [19].This equation has different forms depending upon whether the scour condition is live-bed scour (bed material upstream from bridge is in motion) or clear-water scour (bed material upstream from bridge is not in motion) as discussed in Landers & Mueller [20].The Gao's simplified equation for clear-water pier scour is defined as [19]: where y s is the depth of pier scour below the ambient bed, in feet; K s is the simplified pier shape coefficient which is 1.0 for cylinders; b is the width of bridge pier, in feet; y 0 is the depth of flow directly upstream from the pier, in feet; D m is the mean particle size of the bed material, in feet (for this study D 50 was used as D m ); V o is the approach velocity directly upstream from the pier, in feet per second; and V c is the critical (incipient motion) velocity, in feet per second, for the D m -sized particle.V ic is the approach velocity, in feet per second, corresponding to critical velocity at the pier.V ic can be calculated using the following equation: If the density of water is assumed to be 62.

Froehlich Design Equation
The Froehlich design equation is included as a pier-scour calculation option within the computer model HEC-RAS, Version 3.1 [21].The Froehlich's [22] design equation is defined as: where ϕ is a dimensionless coefficient based on the shape of the pier nose, and is 1.0 for round-nosed piers; Fr 1 is the Froude Number directly upstream from the pier; b e is the width of the bridge pier projected normal to the approach flow, in feet; b is the width of the bridge pier, in feet; D 50 is the particle size for which 50 percent of the bed material is finer, in feet and y 0 is the depth of flow directly upstream from the pier, in feet.

Experiment Setup
Experiments Four different pairs of bridge piers with diameter of 6 cm, 9 cm, 11 cm and 17 cm were used (Figure 1(a)).Bridge piers were constructed from PVC plumbing pipe and were circular in shape.A pair of bridge piers were placed inside both sand boxes at a distance of 50 cm from each other and were fixed to the bottom of the flume.One pair of bridge piers were located in each sand box so two experiments were carried out simultaneously in each experimental run.Each pier was offset from the centre line by 25 cm, as illustrated in Figure 1(b).The water depth in the flume was adjusted by the position of the tailgates.In front of the first sand box, a SonTek Incorporated 2D flow meter was installed to measure the approaching flow velocity, water depth, and inflow discharge during the experiment.A staff gauge was also installed in the middle of each sand box smooth and rough ice conditions.Figure 2 shows the rough ice-covered flow around a bridge pier in experiment as well as ADV measurement around scour depth.The experimental runs were 24 hours long which allowed the scour hole to reach an equilibrium depth as noted in previous experiments conducted by Hirshfield [16].After 24 hours, the flume was gradually drained.The scour depth was manually measured along the outside lines of the circular bridge piers.a comparison of the predicted maximum scour depth from the three equations  1).Likewise, Pier-scour depths calculated using the HEC-18/Jones equation similarly exceeded measured pier scour for all 36 observations (Figure 3(b) and Table 1).

Results and Discussion
However, overestimations were larger for the Froehlich design than for the HEC-18/Jones equation.Overall, the most reliable and accurate equation which has predicted the pier scour depths under open channel and ice-covered flow to a very good extent was Gao's simplified equation.Pier-scour depths calculated using the Gao's simplified equation were smaller than measured scour depths for 6 of the 12 measurements for open channel flow and for 9 of the 12 measurements for smooth ice cover.However, it completely underestimated Pier-scour depths for the rough ice cover flow (Figure 3(c) and Table 1).Statistics for calculated and measured pier scour are summarized in Table 2.The averages of pier-scour depths for open channel flow calculated from Froehlich equation was 0.46 ft, from the HEC-18/Jones was 0.38 ft, and from the Gao's simplified equation was 0.15 ft.The average of the measured pier-scour depths was actually 0.15 ft.This result indicates Gao's simplified equation was entirely successful in predicting the scour depth under open channel conditions.Both the Froehlich and HEC-18/Jones equations resulted in significant over-estimation of the average scour depth.Furthermore, from  Mathematically, Gao's simplified equation is the most successful of the three equations at calculating scour depth but it does underestimate the extent of scour.In practical terms, this could lead to a false sense of security in practical situations as increased scour depth could lead to premature failure of a pier.

Summary and Conclusion
Three pier-scour equations, namely the Froehlich design equation, the HEC-18/ Jones equation and Gao's simplified equation, were evaluated against the data from 36 flume experimental runs obtained for open channel, rough and smooth ice-covered flow under uniform bed sediment type.The most important result that can be obtained from the comparison of these three equations is that under nearly the same flow depth and approach flow velocity but different flow cover, the average calculated values from all the three equations stayed nearly constant.On the other hand, in terms of ice cover, the rougher the ice surface, the more turbulent the flow and the deeper the scour depth generated.Although Gao's simplified equation was reasonably successful in prediction of pier scour depth for both the open-channel and smooth-ice conditions, it underestimated the pier-scour depth under rough ice-covered flow conditions.Therefore, it can be concluded that none of the equations adequately model scour depth under rough ice conditions and the equations are in need of another term to make them more suitable to be used for the ice-covered flow conditions.
How to cite this paper: Namaee, M.R., Li, Y.Q., Sui, J.Y. and Whitcombe, T. (2018) Comparison of Three Commonly Used Equations for Calculating Local Scour Depth around Bridge Pier under Ice Covered Flow World Journal of Engineering and Technology were carried out at the Quesnel River Research Centre, Likely, BC, Canada in a large-scale flume.The flume measures 40 m long, 2 m wide and 1.3m deep.The longitudinal slope of the flume bottom was 0.2 percent.A holding tank with a volume of nearly 90 m 3 was located at the upstream end of the flume and provided a constant head in the experimental zone.Two valves were connected to the holding tank to allow for control of the flow velocity.At the end of the holding tank and upstream of the main flume, water overflowed from a rectangular weir into the flume.Since the flow of water was turbulent while entering the flume, a flow diffuser was placed downstream of the rectangular weir to dissipate the turbulence in the flow of water.Two sand boxes with the depth of 0.30 m were filled with a uniform sediment having a median particle size (D 50 ) of 0.47 mm.The first sand box was 5.6 m in length and the second sand box was 5.8 m in length.The distance between the sand boxes was 10.2 m.

Figure 2 . 10 -
Figure 2. 10-Mhz Acoustic Droppler Velocimeter in use to measure the velocity field around bridge piers under rough ice-covered flow conditions.
to maximum measured scour depth under smooth and rough ice-covered flow condition.It should be mentioned the residual values are defined as predicted scour depth minus measured scour depth.A negative residual value indicates over-estimated scour depth.From the data provided, it is possible to determine which equation is most useful under various conditions.Evaluation of these equations are important to bridge design especially for the case of flow under ice-covered condition.The scatterplots in Figure 3(a)-(c) compare predicted pier scour depth for each of the three equations to the measured scour depths obtained experimentally under the different flow conditions.Pier-scour depths calculated using the Froehlich design equation exceeded measured scour depths for every measurement of both open channel and ice-covered flows (Figure 3(a), Table

Figure 3 .
Figure 3. (a) Comparison of calculated to measured pier scour for Froehlich equation under open, smooth and rough flow condition.(b) Comparison of calculated to measured pier scour for HEC-18/Jones equation under open, smooth and rough flow condition.(c) Comparison of calculated scour depths to the measured scour depth for Gao's simplified equation under open, smooth and rough flow condition.
4 pounds per cubic foot and the bed World Journal of Engineering and Technology material is assumed to have a specific gravity of 2.65, the equation for V c can be

Table 1
demonstrates the performance of the Gao's simplified equation (SCE), the HEC-18/Jones equation (HJE) and Froehlich Design Equation (FDE) under smooth and rough flow conditions.The measured scour in Table1stands for the maximum scour depth between left and right bridge pier.Table2 provides

Table 1 .
Comparison of calculated pier scour to measured pier scour from three equations.

Table 2 .
Error values of the three equations with respect to measured data under open, smooth and rough flow.

Table 2 ,
the highest index of agreement value (Ia) and the lowest RMSE and MAE values are for the Gao's simplified equation for the case of open channel flow.
In terms of smooth ice-covered flow, the averages of pier-scour depths for open channel flow calculated from Froehlich equation was 0.46 ft, from the M. R. Namaee et al.DOI: 10.4236/wjet.2018.62B00659 World Journal of Engineering and Technology