Chapter 2 Hydrology and Hydraulics: City of Raleigh Stormwater Design Manual January 2002
Chapter 2 Hydrology and Hydraulics: City of Raleigh Stormwater Design Manual January 2002
Chapter 2 Hydrology and Hydraulics: City of Raleigh Stormwater Design Manual January 2002
Every model has certain limitations that will effect its behavior for
different size drainage areas. The designer should be familiar with the
limitations of the method he is using. In general, street drainage and
small drainage areas (less than a hundred acres) can be modeled
using the rational equation. Larger areas can be modeled using
methods developed by the Natural Resources Conservation Service.
Distributed element models generally use one of these hydrologic
methods but also employ algorithms to account for reduction in the
peak discharge due to storage in the watershed from reservoirs or road
crossings. The hydrologic model used in the distributed element model
should be appropriate for each of the sub-watershed areas.
2
1
Method Size Limitations Comments
2.2.1 Equation The rational formula estimates the peak rate of runoff at any location in
a watershed as a function of the drainage area, runoff coefficient, and
mean rainfall intensity for a duration equal to the time of concentration
(the time required for water to flow from the most remote point of the
basin to the location being analyzed). The rational formula is
expressed as follows:
Where:
Q = Peak flow from the drainage area (cfs)
C = Coefficient of runoff (dimensionless)
I = Rainfall intensity for a given time to peak (in/hr)
A = Drainage area (acres)
Playgrounds 0.35
Lawns:
Sandy soil, flat, 2% 0.10
Sandy soil, average, 2 - 7% 0.15
Sandy soil, steep, > 7% 0.20
Clay soil, flat, 2% 0.17
Clay soil, average, 2 - 7% 0.22
Clay soil, steep, > 7% 0.35
Residential:
Single-family (R - 4) 0.50
Single-family (R - 6) 0.55
Multi-family (R - 10) 0.60
Multi-family (R - 20) 0.70
Multi-family (R - 30) 0.75
Business:
O & I (I, II, III) 0.85
I1 & I2 0.85 – 0.95
Shopping Centers 0.85 – 0.95
Streets:
Gravel areas 0.50
Drives, walks, and roofs 0.95
Asphalt and Concrete 0.95 – 1.00
2.2.3 Rainfall Intensity The rainfall intensity (I) is the average rainfall rate in in./hr for a duration
equal to the time of concentration for a selected return period. Once a
particular return period has been selected for design and a time of
concentration calculated for the drainage area, the rainfall intensity can
be determined from the intensity-duration-frequency (IDF) data for the
City of Raleigh given in Table 2.3.
Frequency (Yrs)
Duratio 2 5 10 25 50 100
n
5 mins 5.76 6.58 7.22 8.19 8.96 9.72
10 4.76 5.54 6.13 7.01 7.71 8.40
15 4.04 4.74 5.25 6.03 6.64 7.24
20 3.47 4.12 4.64 5.42 5.93 6.47
30 2.70 3.28 3.71 4.32 4.80 5.28
40 2.28 2.77 3.15 3.70 4.08 4.48
50 1.94 2.38 2.71 3.19 3.53 3.88
60 1.70 2.12 2.41 2.84 3.17 3.50
90 1.22 1.52 1.74 2.06 2.29 2.53
2 hr 0.95 1.20 1.37 1.62 1.81 2.00
3 0.71 0.89 1.02 1.21 1.35 1.50
6 0.44 0.56 0.65 0.77 0.86 0.96
12 0.26 0.33 0.39 0.46 0.52 0.57
24 0.15 0.19 0.22 0.27 0.30 0.33
Simple Nomograph
To obtain the total time of concentration, the pipe or open channel flow
time must be calculated and added to the inlet time. After first -
determining the average flow velocity in the pipe or channel, the travel
time is obtained by dividing velocity into the pipe or channel length.
Velocity can be estimated by using the nomograph shown on Figure
2.2. Note: time of concentration cannot be less than 5 minutes.
Kinematic Wave
2.2.5 Kirpich Equation The Kirpich equation is based on empirical data and observation.
Although it has no analytical basis, it has proven an effective method in
many years of use. It is therefore widely considered an acceptable
method for estimating time of concentration for small drainage areas.
The basic form of the equation is:
3 0.385
Equation 2.3 Kirpich Equation Tc = (L / H)
128
Where :
Tc =Time of Concentration (min)
H =Height of the most remote point on the watershed above
the outlet (ft)
L =Length of flow from the most remote point on the
watershed to the outlet (ft)
Example:
Height above outlet, H = 100’
Length of travel, L = 3000’
Time of concentration, Tc = 14 mins
In some cases runoff from a portion of the drainage area which is highly
impervi ous may result in a greater peak discharge than would occur if
the entire area were considered. In these cases, adjustments can be
made to the drainage area by disregarding those areas where flow time
is too slow to add to the peak discharge.
a. For most applications using the rational method the designer may
assume the time of concentration (Tc ) is equal to the time to peak
(Tp). Adjustments are commonly made to Kirpich Equation to
compensate for channelization.
b. For well-defined natural channels, use Tc .
Solution:
The initial pipe selection can then be made on the basis of a design
discharge of 156 cfs.
Soil types are divided into four major hydrologic soil groups denoted by
the letters A through D. A soils are those which have high infiltration
capacity and subsequently low runoff rates. D soils are those with very
low infiltration capacity and very high runoff rates. A list of soils
common in North Carolina can be found at the local NRCS office.
Those soils given dual notation represent hydrologic classifications for
drained and undrained conditions. The Wake County Soil Survey or
maps from the Central Engineering or Inspections Departments are
good sources of information on soil types for specific locations.
For the purposes of the NRCS method, antecedent moisture content
(AMC) is divided into dry, normal, and wet conditions based on the
rainfall in the prior five days. If the 5-day antecedent rainfall is greater
than 2.1 inches in the growing season or 1.1 inches in the dormant
season the moisture content is presumed wet (AMC III). If the 5 day
antecedent rainfall is less than 1.4 inches or 0.5 inches respectively, it
is presumed dry (AMC I). The runoff curve numbers presented here are
based on normal conditions (AMC II). For design of proposed facilities,
normal conditions are generally used.
The time of concentration can be broken into three types of flow, sheet
flow, shallow concentrated flow, and channel flow (or pipe flow). Sheet
flow is assumed to be no longer than a few hundred feet and can be
described by Manning's kinematic solution:
0.8
0.007 ( n L )
Equation 2.7 Time of Tc = --------------------
0.5 0.4
Concentration P2 S
For pipe systems the flow velocity can be estimated by the Manning
equation as well. Assuming the pipe is circular and is flowing just full,
the equation simplifies to the form:
0.667 0.5
0.59 D S
Equation 2.10 Flow in a Pipe V = ------------------
n
2.3.3 Peak Discharge The peak discharge equation used by the Natural Resources
Conservation Servi ces has the form:
Example:
Qp = Qu A Q Fp
2.5.1 Discharge When an impoundment is required, the peak discharge rate after
Limitations development must not exceed that specified by the particular
requirement for every point of discharge from the site. Please note that
the exemption for sites having less than 15% impervious area applies
to the entire site, not each individual outfall, and must include existing
as well as new impervious surfaces. For example, a wooded site has a
2-year discharge of 100 cfs. After development, the site has a peak
discharge rate of 200 cfs. If the site is required to control to the pre-
developed rate, then the impoundment facility must therefore be
designed to limit the discharge to 100 cfs instead of the 200 cfs that
would result after development. The 10-year design peak may also be
similarly impounded
2.5.2 Required Storage The quantity of water that must be detained in order to adequately
Volume reduce the peak discharge is referred to as the Required Storage
Volume. This is the volume that must be available in the facility without
exceeding the maximum permissible release rate. Although the
required volume can only be found by routing the design storm through
the proposed facility, for smaller volumes (less than 20,000 cf), it is
adequate to estimate the volume by subtracting the permissible outflow
from the basin from the peak inflow for the critical storm duration.
Equation 2.12 Storage Volume
S= (Qp - MPRR)*Tp
This may also be a good initial estimate for larger basins. Note that it is
only an estimate and should be verified by routing the design storm
through the proposed facility.
2.5.3 Types of Storage In general, the type of storage device selected depends on the quantity
of water to be stored and the associated cost of storage. Guidelines for
each are included herein. The selection of the type of storage used is
up the individual owner or engineer. Although all of the following types
of facilities will work, some will present more of a maintenance problem.
2.5.3.1 Dry Basins Detention basins should be designed such that the primary outlet
devices restrict the flow and allow water to pond in a safe contained
fashion. A properly designed emergency spillway should be provided
capable of passing the 100-year storm if the drainage area is greater
than 25 acres. Side slopes should be no steeper than 2:1 but if
vegetative groundcover is to be used, flatter slopes are highly
recommended. The basin should be constructed to insure positive
drainage. This will reduce the risk of mosquito problems and reduce
maintenance costs. In larger basins, a concrete low flow swale is
recommended since vegetation may be difficult to maintain with
frequent flow through the basin.
2.5.3.2 Wet Basins The same basic standards apply to wet basins as to dry impoundment
facilities. Outlet devices should be appropriately sized and an
emergency spillway provided. Because of their added benefit to water
quality, wet basins are highly encouraged by the City staff. Although
they may be impractical for smaller areas, their use in larger drainage
areas can provide improved water quality and an attractive, aesthetic
component to the development. Hydrologic modeling of the wet basin is
similar to that of other basins with only some minor changes to the
stage-storage curve.
2.5.3.4 Pipe Storage When space is severely limited on the site pipe storage may be an
option. Oversized pipes with a restricted outlet can provide storage but
usually only at a very high cost. Access to the pipe and outlet device
must be provided for adequate maintenance. Debris control should be
a prime consideration in designing pipe storage since the restricting
outlet device is generally much smaller than the storage or inflow pipes.
2.5.4 Outlet Devices The following sections are a general description of some common
outlet devices used in impoundment facilities. Other devices are
available. Because controlling multiple design storms may be required,
some rather imaginative outlet devices may result. To the extent
possible, outlet devices should be kept simple. This may require and
optimal design for one storm frequency and an over design for the other
storm event.
Example:
At HW /D = 1, h = 0.5 ft.
0.5
Q = 0.65 * 0.785 * [(2)(32.2)(0.5)] = 2.9 cfs
When the headwater depth is below the top of the orifice, h < 0.5 the
discharge can be described by a power curve, the constant for which
can be solved for by substituting in the value at h = 0.5.
The resulting power curve for 3/2 power would then be:
Qp = K(HW)
Therefore:
1.5
Qp = 2.9 (HW)
2.5.4.2 Weir Equation Most weirs used in impoundments will fall into one of two categories;
sharp-crested weirs, such as flow over a standpipe, or broad-crested
weirs such as emergency overflows in basins. Although considerable
research has been conducted in the modeling of weirs, a simple
expression can be applied to most weirs used in stormwater
impoundments. The equation is usually expressed as:
Or
1.5 1.5
Qtot = Cw La Ha + Cw Lb Hb
2.5.4.3 Riser Barrel Riser-barrel outlets act as a combination of several types of outlet
Outlets devices. At different stages the outlet may behave differently. At
shallow depths the riser may act as a weir. As the depth increases the
riser may begin to act as an orifice or the barrel may begin to control.
The controlling factor will be that with the smallest discharge at a given
depth.
Riser as an orifice:
Barrel as an orifice:
Equation 2.18 Barrel as Orifice 0.5
Q = Cd A (2gH)
Note that vortex action can and should be eliminated by using an anti-
vortex plate or designing a riser-barrel combination incapable of
producing a vortex.
Where:
dS
---- = Change in storage in the basin with respect to time
dT
I = Inflow
O = Outflow
2.5.5.1 Inflow The nature of impoundment basin routing is such that the inflow to the
Hydrograph basin must be described in small time increments. In order to
Formulation accomplish this, an inflow hydrograph must be formulated for each of
the post-development design storms. If one of the more complex
hydrologic methods is used to determine discharge from the drainage
area, the inflow hydrograph is already available and can be used for the
routing. If only the peak discharge has been determined, a hydrograph
must be formulated based on that peak.
1. Calculate the peak flow for the given design storm (rational method,
SCS methods, etc.)
2. Estimate the volume of runoff (Type II 6 hr. storm runoff for
example)
3. Given that
Qp Π×t
Equation 2.21 Step Function Q= 1 − cos in radians
2 tp
The resulting plot of stage vs. storage may be used for routing by hand
or a "best fit" equation of the points may be used. The best fit is usually
of the form
Storage = K × Stage b
Equation 2.23 Stage-Storage
Function
Where
3
Storage is in ft
Stage is in feet
K = 200
b=3
2Si 2Si
Ii + Ij + − Oi = + Oj
Tj − Ti Tj − Ti
In order to route the hydrograph, it is necessary to plot the relationship
between outflow and the right half of the equation above. Since the
stage-storage and stage-discharge relationships are known, this can be
easily accomplished. The routing is best explained by example.
Example:
Solution:
1. Select a time increment for the routing. Since the routing is actually
a numerical integration of the differential equation, the time
increment should be small enough to accurately reflect the area
under the hydrograph. Therefore, for any routing technique, the
time step should be no longer than 10 percent of the time to peak.
The time step for this routing will be 2 minutes.
2. Develop the inflow hydrograph. Using the shape of the step
function in Section 2.5.5.1, the inflow hydrograph is listed in column
2 of the routing worksheet.
3. Develop the stage-discharge function. The outlet device will be a
single 42 inch RCP with no tailwater constraints. The discharge
function can then be developed as discussed in Section 2.5.4.1.
The Cd is taken to be 0.65.
Therefore :
3/2
Q = 10.14 * stage
1/2
Q = 0.65 * 9.62 * (2 * 32.2 * (stage - 1.75))
Time Ii Ij 2S - O 2S + O O
Tj-Ti Tj-Ti
(min) (cfs) (cfs) (cfs) (cfs) (cfs)
0 0 10 0 0
2 10 38 10 10 0
4 38 82 58 58 0
6 82 138 46 178 66
8 138 200 122 266 72
10 200 262 288 460 86
12 262 318 550 750 100
14 318 362 906 1130 112
16 362 390 1344 1586 121
18 390 400 1840 2096 128
20 400 390 2358 2630 136
22 390 362 2866 3148 141
24 362 320 3328 3618 145
26 320 281 3712 4010 149
28 281 247 4011 4313 151
30 247 217 4233 4539 153
32 217 190 4389 4697 154
34 190 167 4486 4796 155
36 167 147 4533 4843 155
38 147 129 4537 4847 155
40 129 113 4503 4813 155
42 113 99 4437 4745 154
44 99 87 4341 4649 154
46 87 77 4221 4527 153
48 77 67 4081 4385 152
50 67 59 3925 4225 150
52 59 52 3753 4051 149
54 52 46 3570 3864 147
56 46 40 3376 3668 146
58 40 35 3174 3462 144
60 35 3249 142
The peak discharge will be about 155 cfs occurring about 38 minutes
into the storm. Note that there is some minor instability in the first few
routing steps. As long as these errors do not account for much change
in storage volume they can be ignored.
2.5.5.4 Short Cut An alternative to the storage indication method of routing is a method
Routing Method that may lend itself to computer spreadsheet application. (Developed
by Dr. H.R. Malcom, NCSU, Dept. of Civil Engineering) The
components required for the short cut method are similar to those of
storage-indication method. The short cut is an incremental tabular
application of the same differential equation but simplified to the form:
Si = (Ii-Oi)*(Ti-Tj)
Example:
Repeat the same example used in Section 2.5.5.3 for the storage
indication method using the short-cut method.
Solution:
1. The time step and inflow hydrograph are the same as those used in
the previous example and appear in columns 1 and 2 respectively.
2. The stage-storage curve is the same as the one developed
previously:
1
Storage 3
Stage =
2
3. The stage-discharge function is the same as the one developed
previously:
3/2
Q = 10.14 * stage
1/2
Q = 0.65 * 9.62 * (2 * 32.2 * (stage-1.75))
2.5.6 Ten Percent Rule The “ten percent rule” may be used to determine the downstream
extent of design considerations for new detention. This rule recognizes
that in addition to controlling the peak discharge from the outlet works,
storage facilities change the timing of the entire outflow hydrograph.
Where required, channel routing calculations must proceed
downstream to a confluence point where the drainage area being
analyzed represents ten percent or less of the total drainage area. At
this point, the effect of the hydrograph routed through the proposed
storage facility on the downstream hydrograph is assessed and shown
not to have detrimental effects on downstream hydrographs. If
detrimental impacts are suspected, then backwater calculations and
determination of flood elevations for the areas impacted by increased
flows, if any, must be prepared.