Scour Technology | Bridge Hydraulics | Culvert Hydraulics | Highway Drainage | Hydrology | Environmental Hydraulics |
FHWA > Engineering > Hydraulics > HEC 15 |
Design of Roadside Channels with Flexible Linings
|
(C.1) |
where,
V | = mean channel velocity, m/s (ft/s) |
V* | = shear velocity which is , m/s (ft/s) |
a, b | = empirical coefficients |
R | = hydraulic radius, m (ft) |
kS | = roughness element height, m (ft) |
g | = acceleration due to gravity, m/s2 (ft/s2) |
Manning's equation and Equation C.1 can be combined to give Manning's roughness coefficient n in terms of the relative roughness. The resulting equation is:
(C.2) |
where,
α | = unit conversion constant, 1.0 (SI) and 1.49 (CU) |
General form of the relative roughness equation (Kouwen and Unny, 1969; Kouwen and Li, 1981) is as follows.
(C.3) |
where,
f | = Darcy-Weisbach friction factor |
a, b | = parameters for relative roughness formula |
d | = depth of flow, m (ft) |
k | = roughness height (deflected height of grass stem), m (ft) |
Coefficients "a" and "b" are a function of shear velocity, V*, relative to critical shear velocity, V*crit, where the critical shear velocity is a function of the density-stiffness property of the grass stem. Table C.1 provides upper and lower limits of the relative roughness coefficients. Within these limits, values of the coefficients can be estimated by linear interpolation on V*/V*crit.
V*/V*crit | a | b |
---|---|---|
1.0 | 0.15 | 1.85 |
2.5 | 0.29 | 3.50 |
Critical shear velocity is estimated as the minimum value computed from the following two equations:
(C.4a) |
(C.4b) |
where,
V*crit | = critical shear velocity, m/s (ft/s) |
MEI | = density-stiffness parameter, N·m2 (lb·ft2) |
α1 | = unit conversion constant, 0.028 (SI) and 0.092 (CU) |
α2 | = unit conversion constant, 6.33 (SI) and 3.55 (CU) |
α3 | = unit conversion constant, 0.23 (SI) and 0.69 (CU |
Note: For MEI < 0.16 N·m2 (0.39 lb·ft2) the second equation controls, which is in the D to E retardance range.
The roughness height, k, is a function of density, M, and stiffness (EI) parameter (MEI). Stiffness is the product of modulus of elasticity (E) and second moment of stem cross sectional area (I).
(C.5) |
where,
h | = grass stem height, m (ft) |
τo | = mean boundary shear stress, N/m2, lb/ft2 |
Values of h and MEI for various classifications of vegetative roughness, known as retardance classifications, are given in Table C.2.
Average Height, h | Density-stiffness, MEI | |||
---|---|---|---|---|
Retardance Class | m | ft | N·m2 | lb·ft2 |
A | 0.91 | 3.0 | 300 | 725 |
B | 0.61 | 2.0 | 20 | 50 |
C | 0.20 | 0.66 | 0.5 | 1.2 |
D | 0.10 | 0.33 | 0.05 | 0.12 |
E | 0.04 | 0.13 | 0.005 | 0.012 |
Eastgate (1966) showed a relationship between a fall-board test (Appendix E) and the MEI property as follows:
(C.6) |
where,
hb | = deflected grass stem height resulting from the fall-board test, m (ft) |
α | = unit conversion constant, 3120 (SI) and 265 (CU) |
Kouwen collected additional data using the fall-board test (Kouwen, 1988). These data have been interpreted to have the following relationship.
(C.7) |
where,
Cs | = grass density-stiffness coefficient |
Combining Equations C.5 and C.7 gives:
(C.8) |
Over a range of shallow depths (y < 0.9 m (3 ft)), the n value is a function of roughness height as shown in Figure C.1. The linear relationships shown between roughness height, k, and Manning's n differ with vegetation condition.
(C.9) |
where,
Cl | = k-n coefficient |
α | = unit conversion constant, 1.0 (SI) and 0.213 (CU) |
As is apparent in Figure C.2, a relationship exists between Cl and Cs and is quantified in the following equation:
(C.10) |
Substituting Equations C.10 and C.8 in Equation C.9 yields the following:
(C.11) |
Defining a grass roughness coefficient, Cn as,
(C.12) |
yields the following relationship for Manning's n:
(C.13) |
Figure C.1a. Relative Roughness Relationships for Excellent Vegetated Conditions
Figure C.1b. Relative Roughness Relationships for Good Vegetated Conditions
Figure C.1c. Relative Roughness Relationships for Poor Vegetated Conditions
Figure C.2. Relationship between Cl and Cs
Most of the flow resistance in channels with large-scale (relative to depth) roughness is derived from the form drag of the roughness elements and the distortion of the flow as it passes around roughness elements. Consequently, a flow resistance equation for these conditions has to account for skin friction and form drag. Because of the shallow depths of flow and the large size of the roughness elements, the flow resistance will vary with relative roughness area, roughness geometry, Froude number (the ratio of inertial forces to gravitational forces), and Reynolds number (the ratio of inertial forces to viscous forces).
Bathurst's experimental work quantified these relationships in a semi-empirical fashion. The work shows that for Reynolds numbers in the range of 4 x 104 to 2 x 105, resistance is likely to fall significantly as Reynolds number increases. For Reynolds numbers in excess of 2 x 105, the Reynolds effect on resistance remains constant. When roughness elements protrude through the free surface, resistance increases significantly due to Froude number effects, i.e., standing waves, hydraulic jumps, and free-surface drag. For the channel as a whole, free-surface drag decreases as the Froude number and relative submergence increase. Once the elements are submerged, Froude number effects related to free-surface drag are small, but those related to standing waves are important.
The general dimensionless form of the Bathurst equation is:
(C.14) |
where,
V | = mean velocity, m/s (ft/s) |
V* | = shear velocity = (gdS)0.5 , m/s (ft/s) |
da | = mean flow depth, m (ft) |
g | = acceleration due to gravity, 9.81 m/s2 (32.2 ft/s2) |
n | = Manning's roughness coefficient |
Fr | = Froude number |
REG | = roughness element geometry |
CG | = dchannel geometry |
α | = unit conversion constant, 1.0 (SI) and 1.49 (CU) |
Equation C.14 can be rewritten in the following form to describe the relationship for n.
(C.15) |
The functions of Froude number, roughness element geometry, and channel geometry are given by the following equations:
(C.16) |
(C.17) |
(C.18) |
where,
T | = channel top width, m (ft) |
Y50 | = mean value of the distribution of the average of the long and median axes of a roughness element, m (ft) |
b | = parameter describing the effective roughness concentration |
The parameter b describes the relationship between effective roughness concentration and relative submergence of the roughness bed. This relationship is given by:
(C.19) |
where,
S50 | = mean of the short axis lengths of the distribution of roughness elements, m (ft) |
a, c | = constants varying with bed material properties |
The parameter, c, is a function of the roughness size distribution and varies
with respect to the bed-material gradation. σ, where:
(C.20) |
For standard riprap gradations the log standard deviation is assumed to be constant at a value of 0.182, giving a c value of 0.814.
The parameter, a, is a function of channel width and bed material size in the cross stream direction, and is defined as:
(C.21) |
In solving Equation C.15 for use with this manual, it is assumed that the axes of a riprap element are approximately equal for standard riprap gradations. The mean diameter, D50, is therefore substituted for Y50 and S50 parameters.
Dan Ghere
Resource Center (Olympia Fields)
708-283-3557
dan.ghere@fhwa.dot.gov