Skip to contentUnited States Department of Transportation - Federal Highway AdministrationFHWA HomeFeedback

Hydraulics Engineering

 

CHAPTER 5: LOCAL SCOUR AT PIERS, continued

Figure 7. Charts. Scatterplots of completed versus observed scour, in meters (M), for selected pier scour equations. These 26 scatterplots each show the observed scour from 0 to 9 meters on the X axis and the computed scour from 0 to 9 meters on the Y axis. A line of equality bisects the plot area diagonally; data points on this line indicate that the observed scour was exactly what was predicted. Data points above this line indicate overprediction of scour, whereas data points below the line indicate underprediction. Each scatterplot contains between 254 and 266 data points. On most of the scatterplots, the data points are clustered in the lower left corner, indicating that the majority of the observed and computed values are between 0 and 2 meters.  Plot 7A, the Ahmad equation, shows most data points in a vertical line between observed values of 0 and 2 meters. Of the remaining data points, approximately 10 are overpredicted by several meters and fewer than 10 are underpredicted by several meters.

(A) AHMAD EQUATION

Plot 7B, the Arkansas equation, shows most data points clustered between 0 and 2 meters observed and computed values. Of the remaining data points, approximately 20 are underpredicted by as much as 5 meters.

(B) ARKANSAS EQUATION

Plot 7C, the Blench-Inglis 1 equation, shows most data points clustered between 0 and 2 meters observed and computed values. Of the remaining data points, approximately 15 are overpredicted and 15 are underpredicted.

(C) BLENCH-INGLIS I EQUATION

Plot 7D, the Blench-Inglis 2 equation, shows most data points clustered between 0 and 2 meters observed and computed values. Of the remaining data points, approximately 10 are overpredicted by several meters and 10 are underpredicted by smaller amounts.

(D) BLENCH-INGLIS II EQUATION

Plot 7E, the Breusers equation, shows most data points clustered above the line of equality between 0 and 2 meters observed and computed values. Of the remaining data points, approximately 30 are overpredicted by several meters and 5 are underpredicted by smaller amounts.

(E) BREUSERS EQUATION

Plot 7F, the Breusers-Hancu equation, shows most data points clustered between 0 and 1.5 meters observed and 0 and 3 computed values. Of the remaining data points, approximately 20 are overpredicted by several meters and 10 represent observed values that were not predicted at all.

(F) BREUSERS-HANC EQUATION

Figure 7. Scatterplots of computed versus observed scour, in meters (m), for selected pier scour equations.

Plot 7G, the Chitale equation, shows most data points clustered between 0 and 1.5 meters observed and 0 and 4 computed values. Of the remaining data points, approximately 25 are overpredicted by several meters and 10 are underpredicted by smaller amounts.

(G) CHITALE EQUATION

Plot 7H, the Froehlich equation, shows most data points clustered between 0 and 2 meters observed and computed values and lying along the line of equality. Of the remaining data points, approximately 10 are overpredicted by several meters and 20 are underpredicted by several meters.

(H) FROEHLICH EQUATION

Plot 7I, the Froehlich design equation, shows most data points clustered above the line of equality between 0 and 2 meters observed and 1 and 3 meters computed values. Of the remaining data points, approximately 25 are overpredicted by several meters and 3 are underpredicted, all by less than 1 meter.

(I) FROEHLICH DESIGN EQUATION

Plot 7J, the Hydraulic Engineering Circular 18 equation, shows most data points clustered above the line of equality between 0 and 2 meters observed and 0.5 and 3 meters computed values. Of the remaining data points, approximately 25 are overpredicted by several meters and fewer than 10 are underpredicted, all by less than 1 meter.

(J) HEC-18 EQUATION

Plot 7K, the Hydraulic Engineering Circular 18 K4 equation, shows most data points clustered above the line of equality between 0 and 2 meters observed and 0.5 and 3 meters computed values. Of the remaining data points, approximately 25 are overpredicted by several meters and fewer than 10 are underpredicted, all by less than 1 meter.

(K) HEC-18-K4-EQUATION

Plot 7L, the Hydraulic Engineering Circular 18 K4 Molinas equation, shows most data points clustered between 0 and 2 meters observed and 0 and 2.5 meters computed values. Of the remaining data points, approximately 25 are overpredicted by several meters and fewer than 10 are underpredicted, all by less than 1 meter.

(L) HEC-18-K4Mo EQUATION

Figure 7. Scatterplots of computed versus observed scour, in meters (m), for selected pier scour equations, continued.

Plot 7M, the Hydraulic Engineering Circular 18 K4 Molinas greater than 2 millimeters equation, shows most data points clustered above the line of equality between 0 and 2 meters observed and 0 and 3 meters computed values. Of the remaining data points, approximately 25 are overpredicted by several meters and fewer than 10 are underpredicted, all by less than 1 meter.

(M) HEC-18-K4Mo (>2MM) EQUATION

Plot 7N, the Hydraulic Engineering Circular 18 K4 Mueller equation, shows most data points clustered above the line of equality between 0 and 2 meters observed and 0 and 3 meters computed values. Of the remaining data points, approximately 25 are overpredicted by several meters and fewer than 10 are underpredicted, all by less than 1 meter.

(N) HEC-18-K4Mu EQUATION

Plot 7O, the Inglis-Poona 1 equation, shows most data points in a vertical line between observed values of 0 and 2 meters, up to 8 meters computed values. Of the remaining data points, approximately 10 are overpredicted by several meters and 10 are underpredicted by several meters.

(O) INGLIS-POONA I EQUATION

Plot 7P, the Inglis-Poona 2 equation, shows most data points clustered between 0 and 2 meters observed and computed values. Of the remaining data points, approximately 25 are overpredicted by several meters and 15 are underpredicted by smaller amounts.

(P) INGLIS-POONA II EQUATION

Plot 7Q, the Larras equation, shows most data points clustered between 0 and 2 meters observed and 0.5 and 2 meters computed values. Of the remaining data points, approximately 15 are overpredicted by several meters and 15 are underpredicted by several meters.

(Q) LARRAS EQUATION

Plot 7R, the Laursen 1 equation, shows most data points clustered above the line of equality between 0 and 2 meters observed and 0 and 3 meters computed values. Of the remaining data points, approximately 25 are overpredicted by several meters and 3 are underpredicted by smaller amounts.

(R) LAURSEN I EQUATION

Figure 7. Satterplots of computed versus observed scour, in meters (m), for selected pier scour equations, continued.

Plot 7S, the Laursen 2 equation, shows most data points clustered above the line of equality between 0 and 2 meters observed and 0 and 3 meters computed values. Of the remaining data points, approximately 30 are overpredicted by several meters and 3 are underpredicted by smaller amounts.

(S) LAURSEN II EQUATION

Plot 7T, the Melville and Sutherland equation, shows most data points clustered above the line of equality between 0 and 2 meters observed and 0 and 3 meters computed values. Of the remaining data points, approximately 30 are overpredicted by several meters and none are underpredicted, except for those in the cluster.

(T) MELVILLE AND SUTHERLAND EQUATION

Plot 7U, the Mississippi equation, shows most data points clustered above the line of equality between 0 and 2 meters observed and 0 and 3 meters computed values. Of the remaining data points, approximately 25 are overpredicted by several meters and fewer than 10 are underpredicted by smaller amounts

(U) MISSISSIPPI EQUATION

Plot 7V, the Molinas equation, shows most data points clustered between 0 and 2 meters observed and computed values. Of the remaining data points, approximately 20 are overpredicted by several meters and 15 are underpredicted by smaller amounts.

(V) MOLINAS EQUATION

Plot 7W, the Shen equation, shows most data points clustered between 0 and 2 meters observed and 0 and 3 meters computed values. Of the remaining data points, approximately 15 are overpredicted by several meters and 15 are underpredicted by smaller amounts.

(W) SHEN EQUATION

Plot 7X, the Shen-Maza equation, shows most data points in a vertical line between observed values of 0 and 2 meters, up to 6 meters computed values. Of the remaining data points, approximately 20 are overpredicted by several meters and fewer than 10 are underpredicted by smaller amounts.

(X) SHEN-MAZA EQUATION

Figure 7. Scatterplots of computed versus observed scour, in meters (m), for selected pier scour equations, continued.

Plot 7Y, the Sheppard equation, shows most data points clustered above the line of equality between 0 and 2 meters observed and 0.5 and 4 meters computed values. Of the remaining data points, approximately 20 are overpredicted by several meters and 5 are underpredicted by smaller amounts.

(Y) SHEPPARD EQUATION

Plot 7Z, the Simplified Chinese equation, shows most data points clustered between 0 and 2 meters observed and computed values. Of the remaining data points, approximately 25 are overpredicted by several meters and 15 are underpredicted by smaller amounts.

(Z) SIMPLIFIED CHINESE EQUATION

Figure 7. Scatterplots of computed versus observed scour, in meters (m), for selected pier scour equations, continued.

Ranking the performance of scour-prediction equations is difficult because of the tradeoff between accuracy and underpredictions. If only accuracy is considered, the sum of squared errors can be used to evaluate the equations' performance (table 5). This statistic shows the Froehlich equation to be the most accurate; however, the Froehlich is a regression expression and underpredicted the depth of scour for 129 of 266 field observations. If the smallest number of underpredictions is used to evaluate the equations, the Froehlich Design equation is the best equation because it underestimated only four observations. The Froehlich Design equation, however, ranked 19th based on the sum of squared errors criteria. The magnitude of the underpredictions is just as important, if not more so, than the number of underpredictions; thus, the sum of squared errors for those observations that were underpredicted is another factor that should be considered. The Melville and Sutherland equation had the lowest sum of squared errors for the underpredicted observations, but this equation ranked 26th in overall sum of squared errors. The Melville and Sutherland equation did not underestimate scour by much, but grossly overestimated scour for many cases (figure 7T). The Froehlich Design, HEC-18-K4, HEC-18, HEC-18-K4Mu, and HEC-18-K4Mo (>2 mm) equations all had low sum of squared errors for the underpredicted observations. If the all the ranks are totaled, the Froehlich Design equation appears to be the top equation, followed by the HEC-18-K4Mu, HEC-18-K4, HEC-18, Mississippi, and HEC-18-K4Mo (>2 mm) equations; however, the Froehlich Design equation had the largest sum of squared errors for this group. If only the ranks based on the two sum of squared error categories are used, the HEC-18-K4Mu equation is favored and the Froehlich Design equation drops to a rank of 8.5.No single equation is conclusively better than the rest, but the top six equations generally appear to be the Froehlich Design, HEC-18-K4,

HEC-18-K4Mu, HEC-18-K4Mo (>2 mm), Mississippi, and HEC-18 equations.

Table 5. Summary of the performance of the selected pier scour equations.

Number of Underpredictions

Summation of Ranks

 

   

 

 

 

 

All Ranks

SSE Ranks

Equation

Number of Observations

SSE Magnitude

Rank

Count Number

Rank

SSE Magnitude

Rank

Total

Rank

Total

Rank

Ahmad

266

7536.86

27

61

14

159.48

22

63

23

49

25.5

Arkansas

266

239.52

4

74

20.5

165.61

23

47.5

20

27

16

Blench-Inglis I

266

265.83

5

74

20.5

52.14

17

42.5

18

22

11

Blench-Inglis II

266

954.55

17

174

27

824.60

27

71

25

44

23

Breusers

266

670.40

13

18

9.5

7.14

9

31.5

7.5

22

11

Breusers-Hancu

266

1205.60

21

77

22

201.18

25

68

24

46

24

Chitale

266

2299.40

25

90

23

169.37

24

72

26

49

25.5

Froehlich

266

160.67

1

129

26

98.24

21

48

21

22

11

Froehlich Design

266

1067.77

19

4

1

1.51

2

22

1

21

8.5

HEC-18

266

822.38

15

13

7

2.16

4

26

4.5

19

4.5

HEC-18-K4

262

791.54

14

15

8

1.93

3

25

3

17

2

HEC-18-KMO (All)

266

495.18

11

65

16

17.01

13

40

15.5

24

13

HEC-18-KMO (> 2 mm)

266

608.79

12

21

11

2.47

6

29

6

18

3

HEC-18-K4Mu

266

448.53

9

18

9.5

2.23

5

23.5

2

14

1

Inglis-Poona I

266

1758.81

24

119

25

597.74

26

75

27

50

27

Inglis-Poona II

266

229.68

3

72

19

45.67

16

38

12

19

4.5

Larras

266

311.13

7

48

13

72.09

20

40

15.5

27

16

Laursen I

266

1277.71

23

6

2

5.20

8

33

10

31

21

Laursen II

266

930.57

16

9

3.5

10.95

12

31.5

7.5

28

18

Laursen-Callander

266

960.55

18

9

3.5

10.39

11

32.5

9

29

19.5

Melville and Sutherland

262

3092.08

26

28

12

1.45

1

39

13.5

27

14

Mississippi

266

465.05

10

12

6

7.90

10

26

4.5

20

6

Molinas

262

199.79

2

103

24

55.96

18

44

19

20

7

Shen

266

300.77

6

69

18

37.00

15

39

13.5

21

8.5

Shen-Maza

266

1133.23

20

67

17

36.90

14

51

22

34

22

Sheppard

262

1276.04

22

11

5

3.89

7

34

11

29

19.5

Simplified Chinese

254

344.46

8

62

15

56.21

19

42

17

27

16

SSE-sum of squared errors

Sheppard's equations have been revised several times as reported by Alkhalidi, Sheppard et al., and Sheppard.(65,66,67) The 2001 version of the equation as listed in table 5 was used for comparisons in their analysis because it was the version provided by Sheppard to the author at the time the analysis was conducted. That version of the Sheppard equations ranked 19th based on the two sum of squared error categories.

Since no single equation was superior to the others and none of the equations accurately predicted the scour for all conditions, it is important to assess where the equations failed. Residuals of selected equations were plotted against Froude number (Vo/(gyo)0.5), relative velocity (Vc/Vo), median grain size (D50), pier width (b), relative bed material size (b/D50), and relative depth (yo/b) to assess where the equations may fail to properly account for the scour processes. Figure 8 shows that the Froehlich equation, which is a regression equation, has no significant patterns; it fits the data reasonably well. However, to convert the Froehlich equation from a regression equation to a design equation, Froehlich added the pier width as a factor of safety. The factor of safety increases the scatter in the data significantly. The plot of residuals versus pier width shows that factor of safety becomes too large as the pier width increases (figure 9). The HEC-18-K4 equation shows patterns of increasing overprediction as Froude number (0-0.4), median grain size, and pier width increase (figure 10). The K4, proposed by Mueller, reduces the effect of the Froude number and median grain size, but patterns are still evident in the pier width (figure 11).(5) Only pier width displays a pattern in the residuals of the Mississippi equation (figure 12). The revised HEC-18 equation, HEC-18-K4Mo, also shows patterns in the residuals with Froude number and median grain size, and the most dominant pattern is the bottom envelope on the pier width (figure 13).(1) Most underpredictions seem to occur for grain sizes less than 2 mm. Table 5 shows that two thirds of the underpredictions by HEC-18-K4Mo occur at grain sizes less than 2 mm. Thus, limiting the Ki and K4 corrections to grain sizes greater than 2 mm improves the performance of the Molinas correction.

Although many equations have been proposed for predicting the depth of local pier scour, no equation presented here accurately predicts the depth of scour for the wide variety of conditions represented in the field data set. Six equations are identified as being better than the others when assessed for their value as design equations; however, even these display patterns in their residuals that indicate they are not properly accounting for the scour processes. Additional research and analysis are needed to develop a better local pier scour prediction equation.

Figure 8. Charts. Evaluation of residuals for the Froehlich equation. These six scatterplots show the residuals from the observed and computed values for the Froehlich equation graphed against six potential explanatory variables. For each scatterplot, the Y axis is the residuals from negative 7 to positive 5 meters. Positive residual values indicate underprediction, and negative values indicate overprediction.

(A) FROUDE NUMBER

Plot 8B, relative velocity, shows the residuals plotted against the relative velocity from 0 to 5 on the X axis. The data points are clustered around the horizontal line at 0 between values of 0.5 and 3, with some scatter out to positive 4 and negative 3.

(B) RELATIVE VELOCITY

Plot 8C, 50 percent finer grain size, in millimeters, shows the residuals plotted against the median grain size on a logarithmic scale from 0.1 to 100.0 millimeters on the X axis. The data points are clustered all along the horizontal line at 0, with some scatter out to positive 4 and negative 3.

(C) 50% FINER GRAIN SIZE, IN MM

along the horizontal line at 0, with some scatter out to positive 4 and negative 3.

(D) PIER WIDTH, IN M

Plot 8E, relative bed material size, shows the residuals plotted against the relative bed material size on a logarithmic scale from 10 to 10,000 on the X axis. The data points are clustered all along the horizontal line at 0, with scatter out to positive 4 and negative 3 only after 5,000.

(E) RELATIVE BED MATERIAL SIZE

Plot 8F, relative depth, shows the residuals plotted against the relative depth on a logarithmic scale from 0.1 to 10.0 on the X axis. The data points are clustered all along the horizontal line at 0, with some scatter out to positive 4 and negative 3.

(F) RELATIVE DEPTH

Figure 8. Evaluation of residuals for the Froehlich equation.

Figure 9. Charts. Evaluation of residuals for the Froehlich Design equation. These six scatterplots show the residuals from the observed and computed values for the Froehlich equation graphed against six potential explanatory variables. For each scatterplot, the Y axis is the residuals from negative 7 to positive 5 meters. Positive residual values indicate underprediction, and negative values indicate overprediction.

(A) FROUDE NUMBER

Plot 9B, relative velocity, shows the residuals plotted against the relative velocity from 0 to 5 on the X axis. The data points are clustered between residual values of 0 and negative 3, and relative velocity values of 0.5 to 2. There is some scatter out to negative 7; only 2 data points are positive, and 2 more are exactly 0.

(B) RELATIVE VELOCITY

Plot 9C, 50 percent finer grain size, in millimeters, shows the residuals plotted against the median grain size on a logarithmic scale from 0.1 to 100.0 millimeters on the X axis. The data points are scattered between residual values of 0 and negative 3, with some scatter out to negative 7. Only 3 data points are positive, and 1 more is exactly 0.

(C) 50% FINER GRAIN SIZE, IN MM

Plot 9D, pier width, in meters, shows the residuals plotted against pier width from 0 to 6 meters on the X axis. The data points are clustered between residual values of 0 and negative 2.5 and pier width of 0.5 to 2 meters. There is some scatter after 2 meters between negative 1 and 7. Only 3 data points are positive, and 1 more is exactly 0.

(D) PIER WIDTH, IN M

Plot 9E, relative bed material size, shows the residuals plotted against the relative bed material size on a logarithmic scale from 10 to 10,000 on the X axis. The data points are scattered between residual values of 0 and negative 3, with scatter out to negative 7 mostly after 5,000. Only 3 data points are positive, and 1 more is exactly 0.

(E) RELATIVE BED MATERIAL SIZE

Plot 9F, relative depth, shows the residuals plotted against the relative depth on a logarithmic scale from 0.1 to 10.0 on the X axis. The data points are scattered between residual values of 0 and negative 3, with some scatter out to negative 7 around depth between 1 and 3. Only 3 data points are positive, and 1 more is exactly 0.

(F) RELATIVE DEPTH

Figure 9. Evaluation of residuals for the Froehlich design equation.

Figure 10. Charts. Evaluation of residuals for the Hydraulic Engineering Circular 18 K4 equation. These six scatterplots show the residuals from the observed and computed values for the Hydraulic Engineering Circular 18 K4 equation graphed against six potential explanatory variables. For each scatterplot, the Y axis is the residuals from negative 7 to positive 5 meters. Positive residual values indicate underprediction, and negative values indicate overprediction.

(A) FROUDE NUMBER

Plot 10B, relative velocity, shows the residuals plotted against the relative velocity from 0 to 5 on the X axis. The data points are clustered between residual values of 0 and negative 3, and relative velocity values of 0.5 to 2. There is some scatter out to negative 6; fewer than 10 data points are positive.

(B) RELATIVE VELOCITY

Plot 10C, 50 percent finer grain size, in millimeters, shows the residuals plotted against the median grain size on a logarithmic scale from 0.1 to 100.0 millimeters on the X axis. The data points are scattered between residual values of 0 and negative 3, with some scatter out to negative 6. Fewer than 10 data points are positive.

(C) 50% FINER GRAIN SIZE, IN MM

Plot 10D, pier width, in meters, shows the residuals plotted against pier width from 0 to 6 meters on the X axis. The data points are clustered in a downward-sloping shape between residual values of 0 and negative 2.5 and pier width of 0.5 to 2 meters. There is some scatter after 2 meters between 0 and negative 7. Fewer than 10 data points are positive.

(D) PIER WIDTH, IN M

Plot 10E, relative bed material size, shows the residuals plotted against the relative bed material size on a logarithmic scale from 10 to 10,000 on the X axis. The data points are scattered between residual values of 0 and negative 3, with scatter out to negative 6 around relative bed material size of 50 to 70 and after 5,000. Fewer than 10 data points are positive.

(E) RELATIVE BED MATERIAL SIZE

Plot 10F, relative depth, shows the residuals plotted against the relative depth on a logarithmic scale from 0.1 to 10.0 on the X axis. The data points are scattered between residual values of 0 and negative 3, with some scatter out to negative 6 around depth between 0.8 and 3. Fewer than 10 data points are positive.

(F) RELATIVE DEPTH

Figure 10. Evaluation of residuals for the HEC-18-K4 equation.

Figure 11. Charts. Evaluation of residuals for the Hydraulic Engineering Circular 18 K4 Mueller equation. These six scatterplots show the residuals from the observed and computed values for the Hydraulic Engineering Circular 18 K4 Mueller equation graphed against six potential explanatory variables. For each scatterplot, the Y axis is the residuals from negative 7 to positive 5 meters. Positive residual values indicate underprediction, and negative values indicate overprediction.

(A) FROUDE NUMBER

Plot 11B, relative velocity, shows the residuals plotted against the relative velocity from 0 to 5 on the X axis. The data points are clustered in a slightly downward-sloping shape between residual values of 0 and negative 2, and relative velocity values of 0.5 to 2. There is some scatter out to negative 5; fewer than 10 data points are positive.

(B) RELATIVE VELOCITY

Plot 11C, 50 percent finer grain size, in millimeters, shows the residuals plotted against the median grain size on a logarithmic scale from 0.1 to 100.0 millimeters on the X axis. The data points are scattered between residual values of 0 and negative 3, with some scatter out to negative 5. Fewer than 10 data points are positive.

(C) 50% FINER GRAIN SIZE, IN MM

Plot 11D, pier width, in meters, shows the residuals plotted against pier width from 0 to 6 meters on the X axis. The data points are clustered between residual values of 0 and negative 2 and pier width of 0.5 to 2 meters. There is some scatter after 2 meters out to negative 5. Fewer than 10 data points are positive

(D) PIER WIDTH, IN M

Plot 11E, relative bed material size, shows the residuals plotted against the relative bed material size on a logarithmic scale from 10 to 10,000 on the X axis. The data points are scattered between residual values of 0 and negative 3, with scatter out to negative 5 after relative bed material size 5,000. Fewer than 10 data points are positive.

(E) RELATIVE BED MATERIAL SIZE

Plot 11F, relative depth, shows the residuals plotted against the relative depth on a logarithmic scale from 0.1 to 10.0 on the X axis. The data points are scattered between residual values of 0 and negative 3, with some scatter out to negative 6 around depth between 0.8 and 3. Fewer than 10 data points are positive.

(F) RELATIVE DEPTH

Figure 11. Evaluation of residuals for the HEC-18-K4Mu equation.

Figure 12. Charts. Evaluation of residuals for the Mississippi equation. These six scatterplots show the residuals from the observed and computed values for the Mississippi equation graphed against six potential explanatory variables. For each scatterplot, the Y axis is the residuals from negative 7 to positive 5 meters. Positive residual values indicate underprediction, and negative values indicate overprediction.

(A) FROUDE NUMBER

Plot 12B, relative velocity, shows the residuals plotted against the relative velocity from 0 to 5 on the X axis. The data points are clustered between residual values of 0 and negative 2, and relative velocity values of 0.5 to 2. There is some scatter out to negative 5; fewer than 10 data points are positive.

(B) RELATIVE VELOCITY

Plot 12C, 50 percent finer grain size, in millimeters, shows the residuals plotted against the median grain size on a logarithmic scale from 0 to 100.0 millimeters on the X axis. The data points are scattered between residual values of 0 and negative 2, with some scatter out to negative 5. Fewer than 10 data points are positive.

(C) 50% FINER GRAIN SIZE, IN MM

Plot 12D, pier width, in meters, shows the residuals plotted against pier width from 0 to 6 meters on the X axis. The data points are clustered between residual values of 0 and negative 2 and pier width of 0.5 to 2 meters. There is some scatter after 2 meters out to negative 5. Fewer than 10 data points are positive.

(D) PIER WIDTH, IN M

Plot 12E, relative bed material size, shows the residuals plotted against the relative bed material size on a logarithmic scale from 10 to 10,000 on the X axis. The data points are scattered between residual values of 0 and negative 2, with scatter out to negative 5 after relative bed material size 5,000. Fewer than 10 data points are positive.

(E) RELATIVE BED MATERIAL SIZE

Plot 12F, relative depth, shows the residuals plotted against the relative depth on a logarithmic scale from 0.1 to 10.0 on the X axis. The data points are scattered between residual values of 0 and negative 2, with some scatter out to negative 5 around depth between 1 and 3. Fewer than 10 data points are positive.

(F) RELATIVE DEPTH

Figure 12. Evaluation of residuals for the Mississippi equation.

Figure 13. Charts. Evaluation of residuals for the Hydraulic Engineering Circular 18 K4 Molinas equation. These six scatterplots show the residuals from the observed and computed values for the Hydraulic Engineering Circular 18 K4 Molinas equation graphed against six potential explanatory variables. For each scatterplot, the Y axis is the residuals from negative 7 to positive 5 meters. Positive residual values indicate underprediction, and negative values indicate overprediction.

(A) FROUDE NUMBER

Plot 13B, relative velocity, shows the residuals plotted against the relative velocity from 0 to 5 on the X axis. The data points are clustered between residual values of 1 and negative 3, and relative velocity values of 0.5 to 2. There is some scatter out to negative 5.

(B) RELATIVE VELOCITY

Plot 13C, 50 percent finer grain size, in millimeters, shows the residuals plotted against the median grain size on a logarithmic scale from 0.1 to 100.0 millimeters on the X axis. The data points are scattered between residual values of 1 and negative 3, with some scatter out to negative 5.

(C) 50% FINER GRAIN SIZE, IN MM

Plot 13D, pier width, in meters, shows the residuals plotted against pier width from 0 to 6 meters on the X axis. The X axis values are unevenly spaced; the amount of space given to each value decreases after 1, getting smaller for each value up to 6. The data points are clustered between residual values of 1 and negative 2 and pier width of 0.5 to 2 meters. There is some scatter after 2 meters out to negative 5.

(D) PIER WIDTH, IN M

Plot 13E, relative bed material size, shows the residuals plotted against the relative bed material size on a logarithmic scale from 10 to 10,000 on the X axis. The data points are scattered between residual values of 0 and negative 2, with scatter out to negative 5 between relative bed material size 50 and 80 and after 5,000.

(E) RELATIVE BED MATERIAL SIZE

Plot 13F, relative depth, shows the residuals plotted against the relative depth on a logarithmic scale from 0.1 to 10.0 on the X axis. The data points are scattered between residual values of 1 and negative 3, with some scatter out to negative 5 around depth between 0.8 and 4.

(F) RELATIVE DEPTH

Figure 13. Evaluation of residuals for the HEC-18-K4Mo equation.

EVALUATION OF LABORATORY RESEARCH

General

Laboratory research has been the primary tool for defining the relations among variables affecting the depth of pier scour. The validity of these relations has not been proven in the field. Landers and Mueller evaluated many relations developed in the laboratory by use of transformed data (to obtain a more normal distribution) and smoothing techniques to assess general trends in the data.(21) They found only minimal agreement between the field data and laboratory-based relations. The assessment presented here investigates the relations in the field data for variable combinations commonly reported by laboratory investigations. Unlike the data set used by Landers and Mueller, all data at skewed piers were removed to prevent bias by these data, as previously discussed.(21) No transformations were applied unless necessary for consistency with published relations. Using the basic data without transformation results in a less uniform distribution of the data, but it provides for a more direct comparison with laboratory work.

Pier Geometry

For piers aligned with the approach flow, laboratory research indicates that streamlining the pier nose reduces the depth of scour. Many shapes have been tested in the laboratory (table 6) to determine what effect the pier shape has on the depth of scour, relative to a circular pier. The pier shapes for the field data are classified as unknown, cylinder, round-nose, sharp, or square. Landers and Mueller did not account for variables known to affect the depth of scour other than pier shape and found no correlation between pier shape and scour depth.(21) It is necessary to remove the effects of pier width, velocity, depth, and bed material before comparing the depth of scour among different pier shapes. Figures 14 and 15 show the effect of pier nose shape on scour. First, only the effect of pier width is accounted for by dividing the depth of scour (ys) by the pier width (b). Figure 14 shows that the median of the relative depth of scour decreases as the pier is more streamlined; however, there is overlap in the interquartile range, and the differences are not significant. Second, the effects of pier width, velocity, depth, and bed material are removed by regressing these variables with the depth of scour and by plotting the residuals of a regression against pier shape, which is not included in the regression. The residuals from the regression analysis (figure 15) show the same trend as observed in figure 14. The pier shape does not affect the depth of scour in the field as much as in the laboratory. In the field, flow directions are variable, pier shapes vary with depth, and the effect of submerged debris is not easily accounted for; these combine to reduce the effect of pier shape on the depth of scour.

 

Back to Top

Previous | Table of Contents | Next

Related Links

Hydraulics Laboratory Overview

Hydraulics & Hydrology Research at TFHRC

Hydraulics of Iowa DOT Slope-Tapered Pipe Culverts

South Dakota Culvert Inlet Design Coefficients

For More Information:

Kornel Kerenyi
Turner Fairbank
202-493-3142
kornel.kerenyi@fhwa.dot.gov

 

FHWA
United States Department of Transportation - Federal Highway Administration