Drainage Manual Application Guide2013
Drainage Manual Application Guide2013
hed by The South Afriican Nation
nal Roads A
Agency SOC
C Ltd, © 20113
48 Tambbotie Avenuue, Val de GGrace, Pretooria, 0184, S
South Africa
a
P.O. Box
x 415, Preto
oria, 0001, South Africaa
[Link]
First edition in Afrrikaans calleed “Handleeiding vir Pa addreinerin ng” Pretoriaa, 1981
Second edition tran nslated into
o English caalled “Manu ual on Road Drainage” Pretoria a, 1983
Third partly reviseed edition Pretoriaa, 1986
Fourth p partly revissed edition Pretoriaa, 1993
Fifth Ed
dition – Fu ully revised
d and called the “Drain nage Manual” Pretoriaa, 2006
2 print (inclluding correections up tto 15 Octobe
n d er 2007) Pretoriaa, 2007
Sixth Eddition – Fuully revised
d and extend ded Pretoriaa, 2013
ISBN 9778‐0‐620‐55428‐2 (Drainage Manu ual)
ISBN 9778‐0‐620‐55429‐9 (Drainage Manu ual Applicattion Guide – – only availlable electro
onically)
Copyrig ght:
Apart ffrom any faair dealing for the pu udy, or criticism or
urpose of reesearch or private stu
review, as permittted under the
t Copyrig ght Act, thiis publicatio
on may not
ot be reprod
duced in
whole, o or in part, ffor commerrcial purposses. Additio onal copies for private use may, h however,
be printted off the PPDF version n, obtained from the w website, www [Link] .za
The coppyright doees not exten
nd to the m
methods or intellectual property rrights of anny of the
authorss.
In addittion, photoccopies of th he whole, orr parts of th he publicatiion may be made by a academic
institutiions for studdents.
Editor Edwin Kru ger
Sub‐ediitor Nuno Gom
mes
Cover d
design Nitrogen A
Advertising & Design, JJohannesbu
urg
Disclaimmer:
Althoug gh every effort has been mad de as to the
t accuraccy and appplicability
y of the
informaation contained in thiss publicatio on (includinng supporting flash d drive and sooftware),
the pub
blisher, The South Affrican Natiional Roads Agency SOC
S SANRAL), and the
Ltd (S
authorss, cannot acccept any le
egal responnsibility or liability fo
or any errorrs, omissionns or for
any othher reason w
whatsoever..
F
FOREWOR
RD
Foundedd by the Egyyptian ruler, during the PPtolemies dyynasty (323 BC to AD 330), the Mou useion in
Alexanddria hosted scholars such
h as Euclid aand Archimedes (287 to o 212 BC) w who made siignificant
advancess in mathemmatics of connes and cyliinders as weell as differeential equatioons leading to major
advancess in hydraullic engineeriing. These AAlexandrian scholars laiid the founddations of th heoretical
hydrologgy in conneection with practical
p appplications. Around
A the same time tthe Persians too had
already made an inngenious con ntribution too hydraulic engineering
e by developiing a water delivery
system kknown as qaanats – a suubterranean ssystem of tu unnels conneecting wells.. However, it is the
Romans who were innstrumental in expandingg the sciencee of hydraulic engineerinng to variouss parts of
their emppire.
Throughh the ages, civil engineers have alw ways had to o cope with h unforeseenn natural forrces. The
external forces creatted by climatic change, aand further exacerbated
e by human innduced variaables, can
unexpectedly and siggnificantly influence thee hydrologicaal cycle with h serious soccio-economic effects.
Althoughh mathematiical analysiss and modellling cannot cater for eveery eventuallity, we can certainly
attempt to scientificaally predict the behaviouur of these natural
n forcees and minimmise their im
mpacts on
our envirronment.
South AAfrica, for insstance, is knoown for its loow average annual rainfaalls and largee seasonal variations.
Despite the latter, abnormal
a raiinfalls have historically had disastro
ous consequuences. Altho ough our
problem
ms are not witthout precedeent, societiess are always inter-linked and local cattastrophes co ould have
serious rregional and national repeercussions.
Nazir Allli
Chief Exxecutive Offiicer
The Souuth African National
N Road
ds Agency L
Limited
i Drainage manual – Ap
pplication Gu
uide
ACK
KNOWLED
DGEMENT
TS AND ST
TRUCTUR
RE OF THE
E DRAINA
AGE MANU
UAL
Feedback, comments and suggeestions from users of thee previous editions of thhe manual have h been
incorporrated where possible.
p Thiis edition of tthe manual still
s covers alll the previouus backgroun nd theory
but has been extendded to includ de additionaal flood calcu ulation methhods, the annalyses and design
d of
stormwaater systems,, the hydraullic assessmeent of existin a the modeelling of free surface
ng culverts and
flows annd flood line calculations. With the fuurther expanssion of the manual
m it wass deemed neccessary to
separate the manual into two disttinct documeents; the firstt being the Drainage
D Maanual and th
he second
being thee Drainage Manual
M Appplication Gu uide.
Thank yyou to my feellow editors, Professor F Fanie van Vu uuren, Marco van Dijk aand Nuno Gomes
G for
their commmitment annd enthusiassm in updatiing the manu ual. The com
mpilation and
nd editing waas not an
easy taskk but has beeen completed
d with passioon and dedicaation.
Edwin KKruger
Editor
The Souuth African National
N Ro
oads Agencyy SOC Limitted
Feedbacck:
ii Drainage manual – Ap
pplication Gu
uide
DRAINA
AGE MAN UAL APPL
LICATION
N GUIDE
Tablle of Contents
Foreworrd i
Acknow wledgementss ii
Table off Contents iii
List of ssymbols v
AP
PPENDICE
ES
Appendiix 3A - STATIS
STICAL ANA ALYSIS
Appendiix 3B - STANDARD DESIG GN FLOOD METHOD
Appendiix 3C - STANDARD FLOO OD CALCULLATION FOR
RMS
Appendiix 3D - QT/QRMFF RATIOS FO
OR DIFFER
RENT CATC
CHMENT AR
REAS
Appendiix 3E - SCS-SA
A ADDITION NAL INFORRMATION
iv Drainage manual – Ap
pplication Gu
uide
LIST
T OF SYMB
BOLS
Chapter 2
F = fuuture value
i = annnual discou unt rate as a ddecimal fracttion
IRR = innternal rate of
o return techhnique
n = d
discount perio
od in years
NPV = n present vaalue
net
r = raate at which the left-handd and right-h hand sides off the equatioon are equal, resulting
inn a NPV of zero
z
Chapter 3
a = constant
A = arrea of catchm ment (km²)
ARF = arrea reduction n factor (%)
ARFiT = arrea reduction n factor (%)
b = constant
C = ruun-off coeffiicient (dimennsionless)
C = catchment parrameter withh regard to reeaction time
C1 = ruun-off coeffiicient for ruraal area with a value betw ween zero andd one
C1D = ruural run-off coefficient
c inncorporating the effect off dolomites
C1T = ruural run-off coefficient inncorporating g the effect ofo dolomites and initial saturation
s
faactor
C100 = calibration co oefficient (SD DF method)
C2 = calibration co oefficient (SD DF method)
C2 = ruun-off coeffiicient for urbban area with h a value betw ween zero annd one
C3 = ruun-off coeffiicient for lakkes with a vallue between zero and onee
CN = C
Curve Numbeer
CNf = F
Final Curve Number
N
CNw = C
Curve Numbeer for wet coonditions
CN-II = reetardance facctor approxim mated by thee initial Curve Number unnadjusted forr
anntecedent soil moisture
CP = ruun-off coeffiicient accordding to averag ge soil permeeability
CS = ruun-off coeffiicient accordding to averag ge catchment slope
CT = combined run n-off coefficiient for T-year return period (dimensiionless)
CV = ruun-off coeffiicient accordding to averag ge vegetal grrowth
D = sttorm duration (hours)
F = laag coefficiennt
FT = adjustment faactor for initi al saturation n for return peeriod T
fiT = fllood run-off factor (%)
H = h
height (m)
H = h
height of mosst remote poiint above outtlet of catchm ment (m)
H0,10L = ellevation heig ght at 10% off the length of o the waterccourse (m)
H0,85L = ellevation heig ght at 85% off the length of o the waterccourse (m)
heiT = effective rainffall (mm)
I = raainfall intenssity (mm/houur)
Ia = innitial losses (abstractions
( s) prior to thee commencem ment of storm
mflow, compprising of
d
depression sto
orage, interceeption and in nitial infiltrattion (mm)
IT = avverage rainfaall intensity ffor return peeriod T (mm//h)
K = reegional constant
KRP = constant for T-year
T returnn period
KT = constant for T-year
T returnn period
Ku = d
dimensionless s factor
l = h
hydraulic leng
gth of catchm ment along th he main chan nnel (m)
v Drainage manual – Ap
pplication Gu
uide
L = h
hydraulic leng
gth of catchm ment (waterccourse length
h) (km)
L = catchment lag g time (h)
LC = d
distance fromm outlet to cenntroid of catcchment area (km)
M = 2-year return period dailyy rainfall from m TR102
MAP = m
mean annual precipitation
p n (mm/a)
n = leength of reco ord (years)
P = m
mean annual rainfall
r (mmm/annum)
P = probability (% %)
P = d
daily rainfall depth (mm)), usually in nput as a onee-day designn rainfall forr a given
reeturn period
P1 = probability off at least one exceedence during the design
d life
PAvgT = avverage rainfaall over the ccatchment fo or the T-year return periodd (mm)
PAvgiT = avverage rainfaall for T-yearar storm duration (mm)
PiT = point intensity y for the retuurn period T (mm/h)
PT = point rainfall for the returnn period T (m mm)
Pt,T = thhe precipitatiion depth forr a duration of
o t minutes and
a a return period of T years
y
qp = peak discharg ge (m3/s)
Q = peak discharg ge (m³/s)
Q = sttormflow dep pth (mm)
Qe = peak discharg ge of unit hyddrograph (m³³/s)
QiT = peak discharg ge for T-yearr return perio
od (m³/s)
Qp = u hydrograaph peak disccharge (m³/s))
unit
QRMF = reegional maxiimum flood ppeak flow raate (m³/s)
QT = peak discharg ge for T-yearr return perio
od (m³/s)
r = rooughness coeefficient
R = avverage numb ber of days pper year on which
w thunderr was heard ((days/year)
S = potential max ximum soil w water retentio
on (mm),
Sav = avverage slopee (m/m)
T = tiime (hours)
T = reeturn period (years)
TC = tiime of conceentration (houurs)
TL = laag time L (ho ours)
Tp = tiime to peak (hours)
(
TSD = sttorm duration (hours)
t = d
duration (minnutes)
α = arrea distributiion factor
β = arrea distributiion factor
γ = arrea distributiion factor
qp = peak discharg ge of incremeental unit hyd drograph (m3/s)
Q = inncremental stormflow dep epth (mm)
D = u duration of time, usedd with the distribution off daily rainfalll to accountt for
unit
raainfall intenssity variationns (hours)
30 = 30-minute raiinfall intensitty for the 2-yyear return period (mm/hh)
Chapter 4
A = seectional areaa (m²)
A1 = uupstream secttional area (m m²)
A2 = d
downstream sectional
s areaa (m²)
B = frree surface width
w of crosss section (m))
C = C
Chézy constant
E = sppecific energgy (m)
Fr = F
Froude number
vi Drainage manual – Ap
pplication Gu
uide
g = ggravitational acceleration
a (m/s²)
hl = trransition losss (m)
hf = frriction lossess (m)
ks = m
measure of ab
bsolute roughhness (m)
L = d
distance (m)
P = w
wetted perimeeter (m)
q = d
discharge per unit width ((m³/s/m)
Q = d
discharge (m³³/s)
rc = centre line rad dius (m)
R = h
hydraulic radiius i.e. area ddivided by wetted
w perimeeter (m)
Re = R
Reynolds nummber
So = B slope (m//m)
Bed
S = ennergy slope, which is equual to bed slo ope only wheen flow is unniform (m/m)
v = u
uniform channel velocity (m/s)
v = avverage veloccity (m/s)
vc = crritical flow velocity
v (m/ss)
y = d
depth of flow measured peerpendicularr to the stream mbed (m)
y = d
distance betw
ween water suurface and ceentre of graviity of sectionn (m)
yc = crritical flow depth
d (m)
yn = n
normal/unifor rm flow deptth (m)
z = bed level at point where ddepth of flow w = y (m)
γ = sppecific weigh ht (value forr water 9,8 x 103 N/m3)
Δx = d
distance (m)
ρ = m density = 1 000 kg/m
mass m3 for water
υ = k
kinematic viscosity (≈ 1,114 x 10-6 m²/ss for water)
Chapter 5
A = effective cross-sectional pplan area of thet opening (m²)
(
A = crross sectionaal area (m²)
B = tootal flow wid dth (m)
C = innlet coefficieent (0,6 for shharp edges or
o 0,8 for rou
unded edges)
C = C
Chézy constant
CD = d
discharge coeefficient
d = fllow depth off water (mm))
D = d
depth of flow (m)
d1 = particle size (m)
d2 = siide slope parrticle size (m m)
E = sppecific energ gy (m)
F = blockage factor (say, 0,5)
Fr = F
Froude number
H = tootal energy head
h above ggrid (m)
H = ennergy head ≈ flow depthh for upstream m conditionss (m)
H = h
head (m)
I = raainfall intenssity (mm/h)
KL = d
discharge coeefficient
Lf = leength of floww path (m)
n = M
Manning rougghness valuee (s/m1/3)
n1 = rooad crossfalll (%)
n2 = rooad gradient (%)
P = w
wetted parammeter (m)
Q = d
discharge (m³³/s)
R = h
hydraulic radiius i.e. area ddivided by wetted
w perimeeter (m)
s = ennergy gradieent (m/m)
Chapter 6
Aover = arrea of flow over
o structuree at the flow depth selectted (m²)
Aeff = thhe effective inlet
i area thrrough the struucture (m²)
B = thhe width of the channel ( or the length h of the structure) (m)
d = d
depth of flow over the struucture (m)
D = thhe height of the
t soffit of tthe deck abo ove the river invert level ((m)
fi = a dimensionleess factor rel ated to the design
d level
Fr = F
Froude number
g = g
gravitational acceleration
a (9,81 m/s²)
LB = thhe total widthh of the deckk of the struccture (m)
n = M
Manning n-vaalue (s/m1/3)
nconcrete = M
Manning rougghness coeffi ncrete (s/m1/3)
ficient of con
nriver = M
Manning rougghness coeffi ficient of the river bed (s//m1/3)
Pcell = thhe total wetteed perimeterr of each cell (m)
Pconcrete = thhe part of thee wetted periimeter that haas a concretee surface perr cell (m)
Peff = Σ Pcell (effectiv t flow passsing through the structuree) (m)
ve wetted peerimeter for the
Pover = w
wetted perimeeter at the floow depth selected (m)
Priver = thhe part of thee wetted periimeter that iss made up by y the riverbedd per cell (m)
Q2 = d
discharge with a 1:2 year return period d (m³/s)
Qdesign = d
design dischaarge (m³/s)
Qover = d
discharge oveer the structuure within thee selected flo
ow depth (m³³/s)
Qunder = d
discharge cap
pacity of the oopenings thrrough the stru ucture (m³/s))
R = h
hydraulic radiius (m)
S0 = sllope in direction of flow (m/m)
v under = avverage veloccity of flow tthrough the structure
s (m/ss)
x = thhickness of th he deck (deppending on th he structural design outcoome) (m)
Chapter 7
A = crross sectionaal area (m²)
B = width (inside of culvert) ((m)
w
CB = innlet coefficieent for culverrts
Ch = innlet coefficieent for culverrts
D = innside diameter (m)
D = h
height (insidee of culvert) ((m)
Fr = F
Froude number
h f 1 2 = frriction lossess between crooss-section 1 and 2 (m)
h 1 1 2 = trransition lossses between cross-section n 1 and 2 (m)
H1 = u
upstream enerrgy level, rellative to the invert
i level (m)
(
H2 = d
downstream energy
e level,, relative to th
he invert lev
vel (m)
Kin = innlet secondarry loss coeffi ficients
Kout = outlet secondaary loss coeffficients
R = h
hydraulic radiius i.e. area ddivided by wetted
w perimeeter (m)
S0 = n
natural slope (m/m)
Sc = crritical slope (m/m), wherre Fr = 1
v = avverage veloccity (m/s)
yn = n
normal flow depth
d (m)
yc = crritical flow depth
d (m)
ix Drainage manual – Ap
pplication Gu
uide
v 2c = avverage critical velocity inn constriction
n (m/s)
v n2 = avverage flow velocity at ssection 2 baseed on An2
va = avverage veloccity in the maain channel
V* = shhear velocity y (m/s)
V*c = crritical shear velocity (m//s)
Vi = appproach velo ocity when pparticles at piier begin to move
m (m/s)
V = chharacteristic average veloocity in the contracted
c seection (m/s)
Vss = particle settlin ng velocity (m m/s)
VR = v
velocity ratio
Vc50 = crritical velociity for D50 beed material size (m/s)
Vc90 = crritical velociity for D90 beed material size (m/s)
y = m
mean depth of flow (m)
y = d
depth of flow in the contraacted bridge opening (m))
y = projected norm mal flow deppth in the con nstriction (m
m)
y0 = d
depth upstream of pier (m m)
y1 = fllow depth dirrectly upstreeam of pier (m m)
y2 = fllow depth un nder bridge (m m)
y2c = crritical depth in constricti on (m)
ys = sccour depth (m m)
Yt = tootal maximum m scour deptth (m)
Y0 = m
maximum genneral scour ddepth (m)
Ys = loocal scour deepth (m)
α1 = v
velocity coeffficient
α2 = v
velocity head coefficient ffor the consttriction
θ = bank angle wiith the horizoontal (°)
ρ = d
density of waater (kg/m3)
ρd = d bulk denssity (kg/m3)
dry
ρs = saaturated bulk k density (kgg/m3)
φ = riiprap angle of o repose (°)
τc = crritical tractiv
ve stress for sscour to occuur (N/m²)
ν = k
kinematic fluiid viscosity ((m²/s)
Chapter 9
A = fuull-flow areaa (m²)
A1, A2 = fuull-flow areaa for the inlett pipe and ou
utflow pipe (m
m²)
D = pipe inner diaameter (m)
g = g
gravitational acceleration
a (9,81 m2/s)
h = d
difference in elevation beetween the hiighest incom
ming pipe invvert and the centreline
c
of the outlet pipe
p (m)
hf1-2 = frriction lossess between crooss-section 1 and 2 (m)
hL = m
minor loss (m
m)
hl1-2 = seecondary lossses betweenn cross-sectio on 1 and 2 (m
m)
x Drainage manual – Ap
pplication Gu
uide
z1, z2 = innvert elevatio
ons of the infflow pipes reelative to thee outlet pipe iinvert (m)
γ = sppecific weighht (value forr water 9,8 x 103 N/m3)
ν = k
kinematic viscosity (m²/s))
Chapter 10
D = v
vertical dimen
nsion of the eexisting culv
vert.
dS
= thhe change in storage overr the time steep of dt (m3)
dt
I = avverage inflow w (m³/s)
N = auuxiliary funcction (m³/s)
O = outflow throu ugh culvert (m m³/s)
O = avverage outflo ow (m³/s)
QT0 = d
design flood for the desiggn return peeriod which was w obtainedd from the review
r of
thhe road classsification, RCC0 and the ind
dex flood, Q20
2 (m³/s).
QC1 = m
maximum caalculated exiisting inlet capacity
c of the culvert bby limiting the total
ennergy head tot 1,2D (m³/ss)
QC2 = m
maximum callculated currrent existing hydraulic caapacity of thhe culvert byy limiting
thhe total energgy head to bee equal to thee shoulder brrake point lev
evel (SBP) (mm³/s)
QT1 = d
design flood for the desiggn return peeriod which was w obtainedd from the review
r of
thhe road classsification, RCC-1and the ind
dex flood, Q20
2 (m³/s)
Q20 = inndex flood foor the contribbuting catchm ment with a return
r periodd of 20 years (m³/s)
Q2T0 = fllow rate relaated to a retuurn period tw wice that whhich was obbtained for th he design
fllood, QT0 (mm³/s)
Q2T1 = fllood rate related to a retturn period twice
t that which was obbtained for th he design
fllood, QT1 (mm³/s)
RC0 = original road classificationn
RC-1 = reeflects the selection
s of a road classsification which
w is onee class less than that
d
determined fo
or the road
S = teemporal storaage or pondiing volume (m³)
S = suum of the stoorage volumee of the prism m and the weedge (m3)
T = d
design return period
Tc = tiime of conceentration (h)
Ts = tootal time durring the routting of the fllood when th he upstream energy head d is more
thhan 1,2D (h)
VT1 = m
maximum stoorage volum me upstream of the culv vert assuminng level-poo ol routing
conditions an nd an inflow hydrograph with a peak flow rate off QT1 and a triangular
t
d
distribution with
w a base w a the peak discharge occcurring at Tc (m³)
width of 3Tc and
Vstorm = calculated sto orm volume based on thee assumption n of an infloow hydrograp ph with a
peak flow rate of QT1 andd a triangularr distribution n with a basee width of 3T Tc and the
peak discharg ge occurring at Tc (m³)
x = a dimensionleess weightinng factor ind dicating the relative
r impoortance of th
he inflow
(II) and the ou
utflow (O) in determining g the storage (S) in the reaach
ΔS = chhange in storrage volume (m³)
Δt = tiime step thatt is used (s)
I1 I 2
Δ
Δt = avverage volum
metric inflow
w (m³)
2
O1 O 2
Δt = avverage volum
metric outfloow (m³)
2
Chapter 11
A = crross sectionaal flow area ((m)
B = toop width (m))
dy = chhange in watter depth (m))
xi Drainage manual – Ap
pplication Gu
uide
dx = ddistance over which changge occurs (m m)
Esc = sppecific minimmum energy (m)
Fr = F
Froude number
g = g
gravitational acceleration
a (m/s²)
h = sttage height (m)
Q = fllow rate (m³//s)
S = bed slope (m//m)
Sc = crritical slope (m/m)
Sf = reepresents thee slope of thee total energy
y line
S0 = bed slope (m//m)
V = m
mean cross-seectional veloocity
y = fllow depth (mm)
yc = crritical flow depth
d (m)
yn = n
normal flow depth
d (m)
α = v
velocity coeffficient
Es = sppecific energgy (m)
Chapter 12
A = suurface area (m2)
Ag = g
geotextile areea available ffor flow (m²))
At = tootal geotextille area (m²)
AOS = appparent open ning size (mm m)
B = a coefficient (dimensionle
( ess)
B = w
width of colleector drain (m m)
Cu = thhe uniformity y coefficientt
d = d
diameter of piipe (m)
D85 = sooil particle size for whichh 85% of opeenings are sm maller (mm)
Dx = thhe sieve size through whiich x% of thee material paasses (mm)
g = loongitudinal slope
s of the rroad (m/m)
i = h
hydraulic graddient (m/m)
I = d
design infiltraation rate (mm
mm/h)
k = D
Darcy coefficcient of perm meability (m/ss) and
ks = permeability of o material (m m/day)
kb = permeability of o an open-ggraded layer (m/day)
kt = permeability of o the channnel backfill (m
m/day)
L = leength of the pipe
p (m)
L = leength of paviing (1 m widde) subject to o infiltration (m)
1/3
n = M
Manning’s n (s/m )
nb = porosity of an n open-gradeed layer
O95 = opening size in i geotextile for which 95% of openin ngs are smalller (mm)
P = 1h duration/1 year return pperiod rainfaall intensity (mm/h)
(
q = d
drainage rate (mm/day)
q = d
discharge per meter widthh (m3/s.m)
S = sppacing (m)
S = crross-slope off a drainage llayer (m/m)
So = sllope of the pipe
p (m/m)
t = d
depth of flow in material ((mm)
T = d
drainage perio od for layer ((h)
tb = thhickness of drainage
d layeer (mm)
tb = effective thick kness of draiin layer (mm
m)
W = w
width of the drainage
d layeer (m)
ψ = g
geotextile perrmittivity
The SAN NRAL Drainage Manuaal had a maj ajor update ini 2006 whiich includedd additional chapters,
worked examples annd links to ap pplicable draainage softw
ware. SANRA AL decided tto revise andd upgrade
the manuual to be useed as a user tool for all ppersons invo
olved in the design of dra
rainage structures and
systems.. The purposse of the Draainage Manuual (cover sh hown in Figuure 1.1) is too provide a reference
documennt with regarrd to drainagge and to dem nd reference some softwaare for the modelling
monstrate an m
and anallysis of drainnage problem
ms.
Figure 1.1: SA
ANRAL Dra
ainage Manu
ual
Reference too applicable software uti lities and useer manuals for
fo the prograams.
1-1 Introducttion
1.1 Layout of the Drainage Manual Application Guide
The Drainage Manual Application Guide contains twelve chapters. The focuses in the different
chapters are:
Chapter 3: Illustrate the various flood calculation methods for different recurrence intervals.
Chapter 7: Contains analysis and design details for lesser culverts and storm water pipes.
Chapter 8: Focuses on bridges and major culverts and scour at these structures.
Chapter 10: Assessment of hydraulic capacity of existing drainage structures and the application of
flood routing.
It is trusted that the document will provide valuable assistance in the design of drainage systems.
References to all the figures and literature can be found in the Drainage Manual
1-2 Introduction
2 ECONOMIC EVALUATION OF DRAINAGE SYSTEMS
The following three simple examples have been included to illustrate the use of the economic
evaluation procedures. The supporting software is capable of determining the NPV, IRR and LCA for
more complex income and expenditure streams. It is suggested that the supporting software be used to
conduct sensitivity analyses.
Evaluate which of the future income streams S1 or S2 is more favourable if the cost of capital is 10%
on a yearly basis and the amounts realize at the beginning of the year.
Year S1 S2
1 250 100
2 350 400
3 600 350
4 100 600
5 400 250
If you assume year 1 to be the base year then the NPV’s of the two income streams are:
NPVS1 = R1 412.39
NPVS2 = R1 374.43
These calculations reflect that the income stream S1 is more favourable when comparing the Net
Present Values (NPV).
Each future value was brought back to present values and accumulated to obtain the total NPV for
each income stream.
Determine the current investment that should be made for the replacement of a R1,5 million
installation (current cost) after 15 years, if the expected CPIX is 15 % and the return on a fixed
investment is 8% p.a.
Firstly the future value (F) of the investment should be determined. The current installation (P) is
worth R1 500 000 and the escalation will be 15 % for a 15-year period.
F P 1 i
n
…(2.2)
Now the current investment (P) should be calculated by discounting the future required value (F) by
8% per annum for the 15-year period.
F
P …(2.3)
1 i n
12 205 592
P R3 847 712
1 0,0815
2.3 Example 2.3 – Internal Rate of Return (IRR)
Determine the Internal Rate of Return (IRR) for the following cash flow.
The internal rate of return is the rate where the NPVincome = NPVexpenditure
1300
NPVexpenditure …(2.5)
1 i 0
Equation (2.11) = Equation (2.12)
IRR = 9,525%
In the tw
wo paragraphhs below, flo ood peaks w will be calculated for a smmall as well as a large catchment
c
using thee various releevant determ
ministic, statisstical and em
mpirical meth
hods.
Data reqquirements
Calculattion proceduure
H 0,855L H 0,10L
Sav … (3.1)
1 00000,75L
where:
Sav = averagee slope (m/m m)
H0,100L = elevatioon height at 10% of the length of the watercoursee (m)
H0,855L = elevatioon height at 85% of the length of the watercoursee (m)
L = length of watercourrse (km)
The elevatiion at 10% ofo the lengthh of the long ourse is H0,100L = 1 412,1 m and at
gest waterco
85% of thee length the elevation
e is H 0,85L = 1 528
8,8 m.
Step 4: Calculate the time off concentratiion from caatchment ch haracteristics . The recom
mmended
empirical formula
f for calculating the time off concentration in naturaal channels has been
developed by the US Sooil Conservaation Servicees.
0,385
00,87L2
Τ C … (3.2)
1 000 S av
where:
TC = time off concentratioon (hours)
L = length of watercourrse (km)
Sav = averagee slope (m/m m)
In most caases the long gest water ppath includess both overland and chaannel flows. In large
catchmentss the channel flow is uusually dom minant, but in small caatchments it may be
necessary to
t determinee TC as the ssum of the flow
f times fo
or both the ooverland and d channel
flow stretcches. To obttain a broadd indication, it may usu ually be acccepted that a defined
watercoursse exists wheen the averagge slope of th
he catchmentt is greater thhan 5 per cen
nt and the
catchment itself is largeer than 5 km²
m².
Figure 3.6:
3 Rainfalll stations ussed in determ
mining the representati
r ive MAP
The mean annual precipitation (MAP) for this catchment is 746,6 mm (see Table 3.1) and
the catchment is located in the inland region.
Historically there have been a number of ways in which the rainfall intensity could be
determined. These alternative methods have been retained in this document although the latest
method using the Design Rainfall Estimation Software is recommended (Alternative 3).
Alternative 2 – The TR102 representative rainfall data and the modified Hershfield equation
is used (in the past this was referred to as the Alternative Rational Method).
Alternative 3 – The design rainfall from the Design Rainfall Estimation software is used to
determine the point rainfall of the catchment.
Alternative 1
Step 6a: Determine the point rainfall values (PT) (mm) for the required return periods. Based on the
mean annual precipitation (MAP), the rainfall region, the time of concentration (TC) and the
required return period, Figure 3.7 can be used to determine the point rainfall. As shown in
Figure 3.7 the point rainfall for the 1:20 year and 1:50 year return periods are determined
using the co-axial Depth-Duration-Frequency diagram.
The point rainfall for the 1:20 and 1:50 year return periods is P20 = 78 mm and P50 = 104 mm
respectively.
The point intensity (PiT) is the point rainfall divided by the time of concentration (if TC >
0,25 hours). If TC ≤ 0,25 hours divide by 0,25 hours.
PT
PiT … (3.3)
TC
where:
PiT = point intensity for the different return periods (mm/h)
PT = point rainfall (mm)
TC = time of concentration (hours)
The point intensities for the 1:20 and 1:50 year return periods are Pi20 = 58,3 mm/h and
Pi50 = 77,7 mm/h respectively.
In this exam
mple the catcchment area is small and
d thus Figuree 3.8 is usedd. The resultiing ARFs
are ARF20 = 94% and ARF
A 50 = 91% % for the 1:20 and 1:50 years
y return pperiods respeectively.
Step 9a:: Determine the average rainfall intennsity or effecctive catchmeent precipitaation.
ARFT
I T PiT … (3.4)
100
w
where:
IT = rainfall inntensity averraged over th he catchmennt in millimeetres/hour
for the retturn period T.
T
ARFT = area reducction factor asa a percentaage for returnn period T (sshould be
smaller thhan 100%)
PiT = point intennsities for th
he different reeturn periodss (mm/h)
Step 6b: Determine the representative rainfall from the available TR102 South African Weather
Service stations in and around the catchment (see Table 3.2).
Step 7b: Based on the calculated time of concentration and representative rainfall, determine the
precipitation depth. In this example the time of concentration is 80 minutes, in other words
less than 6 hours, and thus the modified Hershfield relationship will be used.
Pt,T 1,130,41 0,64lnT 0,11 0,27lnt 0,79M 0,69 R 0,20 … (3.5)
where:
Pt,T = precipitation depth for a duration of t minutes and a return period of T
years (mm)
t = duration (minutes)
T = return period
M = 2-year return period daily rainfall from TR102
R = average number of days per year on which thunder was heard (days/year)
(Figure 3.9)
The average number of days on which thunder was heard (R) is equal to 61 and M is
60 (from Table 3.2). The calculated precipitation depths are:
Pt20 = 85,55 mm and Pt50 = 107,10 mm
The point intensity (PiT) is the point rainfall divided by the time of concentration (if TC >
0,25 hours). If TC ≤ 0,25 hours, divide by 0,25 hours.
PtT
PiT … (3.6)
TC
where:
PiT = point intensity for the different return periods (mm/h)
PtT = precipitation depth for a duration of t minutes and a return period of T
years (mm)
TC = time of concentration (hours)
The point intensity for the 1:20 and 1:50 year return periods is Pi20 = 63,93 mm/h and
Pi50 = 80,04 mm/h respectively.
where:
ARF = area redduction factorr as a percentage (shouldd be less than
n 100%)
A = catchmeent area (km²²)
TC = time of cconcentration (hours)
Step 10b
b: Determinee the averagee rainfall inteensity or effeective catchm
ment precipittation
ARFT … (3.8)
I T PiT
100
w
where:
IT = rainfall intensity av veraged overr the catchm ment (mm/h) for the
return peeriod T.
ARFT = area redduction facto or as a percenntage for retturn period T (should
be smalller than 100 %)
PiT = point inttensities for the different return perioods (mm/h)
The averrage rainfall intensities arre I20 = 61,4 mm/h and I50
5 = 76,9 mm
m/h.
This is the
t preferred method off obtaining point
p rainfa
all data for uuse in the
Rationall Method prrocedure.
Step 6c:: The designn rainfall frrom the Dessign Rainfalll Estimation n software iis the recom mmended
method to determine th he point desiggn rainfall of a catchmen nt (see Figurre 3.11). Utilising the
Design Raiinfall Estimaation in Soutth Africa software application the reepresentativee weather
station or coordinates
c as
a shown in F Figure 3.12 can be enterred. A summ mary of all th
he closest
rainfall stattions as welll as the n-dayy rainfall vallues as shown in Figure 3.13 is obtaiined. The
software ennables the esstimation of design rainffall for durations rangingg from 5 min nutes to 7
f 2 to 200 year return periods at any
days and for a 1' latitu ude x 1' longgitude point in South
Africa.
Step 7c:: Based on thhe calculated a representtative rainfallll, determine the point
d time of conncentration and
rainfall values (see Figu
ure 3.14) forr the catchment area.
The run-offf coefficientt in the ratioonal method is an integraated value reepresenting the
t many
factors influencing the rainfall run-ooff relationsh
hip.
There is noo objective theoretical m method for determining
d C and as a result the subjective
elements ofo experiencce and engiineering jud dgement playy a very im mportant rolle in the
successful application of
o this methood.
Based on the available data from m the catchmment, the folllowing tablees were com
mpiled to
characterise the catchm
ment (Table 33.4 and Tablle 3.5).
Utilising Table 3.4 the run-off coefficient for the rural area is calculated using the following formula:
C1 C S C P C V … (3.10)
where:
C1 = run-off coefficient with a value between zero and one
CS = run-off coefficient according to average catchment slope
CP = run-off coefficient according to average soil permeability
CV = run-off coefficient according to average vegetal growth
The average rainfall falls between 600 and 900 mm, and thus:
If it is estimated that up to 10% (D%) of the area could be dolomitic, then the run-off factor should be
reduced as described earlier in this chapter. Based on the defined slopes, the following factors (Dfactor)
are used to adjust the run-off coefficient.
Vleis and pans (slopes <3%) - 0,10
Flat areas (3 to 10%) - 0,20
Hilly (10 to 30%) - 0,35
Steep areas (slopes >30%) - 0,50
C1D C1 1 D % C1 D % D factor x C S%
C1D 0,36651 0,1 0,36650,10,10 x 0,20 0,20 x 0,70 0,35 x 0,1
C1D 0,337
C1D is the rural run-off coefficient that incorporates the effect of the dolomitic area.
The influence of initial saturation is incorporated by means of an adjustment factor. Using these
adjustment factors (FT) for rural areas, the run-off coefficients (C1D) for the 1:20 and 1:50 year return
periods are adjusted as follows with F20 = 0,67 and F50 = 0,83.
C1T C1D x FT … (3.11)
The run-off coefficient for the urban area (C2) is calculated using the defined catchment
characteristics (Table 3.5) and the most conservative coefficients of Table 3.5 (for this
example).
Step 11: Determine the peak flow for each of the required return periods utilising the simple linear
relationship and for this example the various methods used to calculate the average rainfall
intensity (Alternatives 1 to 3):
C I A
QT T T … (3.13)
3,6
where:
QT = peak flow rate for T-year return period (m³/s)
CT = combined run-off coefficient for T-year return period
IT = average rainfall intensity over catchment for a specific return
period (mm/hour)
A = effective area of catchment (km²)
3,6 = conversion factor
The peak flow rates based on the rational method for the 1:20 year and 1:50 year return periods
(calculated by means of equation 3.13) are:
Step 1: The first three steps of the Rational method described above are also applicable to the Unit
Hydrograph method, thus A = 28,5 km2, L = 7,25 km and S = 0,02146 m/m.
Step 2: Determine the veld-type zone in which the catchment is located from Figure 3.15. The
catchment of the Moretele Spruit falls in Zone 8.
Step 7: Obtain the mean annu ual precipitattion (MAP), as described for the Raational meth
hod. The
his catchmentt is 746,6 mm
determinedd MAP for th m/a.
For this example the point rainfalls for the 0,25 hour, 0,5 hour, 1 hour and 2 hour
storms have been determined for the different return periods 1:20 and 1:50 year (see
Table 3.7).
Step 8.2: Calculate the point rainfall intensity (mm/hour). The point intensity (PiT) is the point
rainfall divided by the storm duration (TSD).
P
PiT T … (3.16)
TSD
where:
PiT = point intensities for the different return periods (mm/h)
PT = point rainfall (mm)
TSD = storm duration (hours). If duration < 0,25 hours use 0,25 hours.
See solution in Table 3.7.
Step 8.3: Determine the area reduction factors (ARFiT) for the different return periods based on
the catchment area and different storm durations from Figure 3.8. The determined
ARFiT values are shown in Table 3.7.
Step 8.4: Calculate the average rainfall (PAvgiT) for the different return periods and storm
durations. This is the area reduction factor (ARFiT) multiplied by the point rainfall
(PT). The average rainfall values are shown in Table 3.7.
Step 8.5: Determine the flood run-off factor from Figure 3.17. This factor is based on the
average rainfall, veld-type zone and catchment area. The flood run-off factors (fiT) are
given in Table 3.7.
Step 8.6: Calculate the effective rainfall (heiT) for each return period and selected storm duration
by multiplying the flood run-off factors (fiT) with the average rainfall values (PAvgiT).
Step 9: The maximmum flood peeak is obtainned by multip plying the efffective rainffall for speciific storm
durations with
w the unnit hydrograpph peak flo ow. The durration of sto torms that cause
c the
maximum peak
p discharrge is obtaineed by trial an
nd error. Sin
nce the standdard duration n of a unit
hydrographh is one houur (from thee one-hour raainfall), the duration sho hould be incrreased or
decreased to
t make prov vision for othher rainfall durations.
d
The dimennsionless onee-hour unit hhydrograph forfo veld-type zone 8 is shhown in Figu ure 3.19,
obtained frrom Table 3.8, and the cconstructed S-curve
S in Fiigure 3.20. TThe rising annd falling
limbs of thhe unit hydro
ograph used tto construct the S-curve are not equaal. If the ord dinates of
the staggerred unit hydrographs aare summed,, the constructed S-curv rve is not constantly
c
increasing until it reacches the maxximum valuee thereof. Ass illustrated in this exam mple, this
leads to ann uneven S-curve as show wn in Figuree 3.20. It is suggested
s thhat this be rectified as
shown by preventing
p th
he S-curve vvalues for ex
xample (Q/Q Qp)t being lesss than (Q/Qp)t-1. This
approach iss conservativve, which coould probably y lead to an over-estimat
o tion in calcullating the
b should prrovide a concceptually corrrect answer in terms of the flood
volume of discharge, but
peak value.
Figure 3.19:
3 Dimenssionless one--hour unit hydrograph
h for veld typpe Zone 8
Step 11:: The last step is to caalculate the peak flows for the diffferent returnn periods an nd storm
durations. The maxim mum peak floow for each return
r period
d is then useed as the dessign peak
flood for thhat specific return
r periodd as shown in
n Table 3.9. The peak vvalues are addjusted as
indicated for
fo QPiT/QP < 1; in this casse 0,9923.
Figu
ure 3.20: S-ccurve for veeld type Zon
ne 8 and lagg
ged by 0,25 hhour
The calcculation sequuence to deterrmine the floood peaks baased on the SDF method iis as followss:
Step 1: Identify thee drainage baasin in whichh the site is located from Figure 3.211. The catchm
ment falls
in drainagee basin number 1.
Step 4: Determine the average slope of thhe catchmentt as describeed for the R
Rational method. The
calculated average
a slope (S) for thiss example is 0,02146 m/m
m.
Pt,T 1,130,41 0,64lnT 0,11 0,27lnt 0,79M 0,69 R 0,20 … (3.18)
where:
M = mean of the annual daily maxima from Table 3B.1 (in Appendix B)
equals 56 mm.
R = average number of days per year on which thunder was heard from
Table 3B.1 equals 30.
Step 7: Multiply the point precipitation depth PtT (mm) by the area reduction factor ARF (%) to
determine the average rainfall over the catchment for the required return period (PAvgT). The
corresponding rainfall intensity IT (mm/h) is obtained by dividing this value by the time of
concentration.
ARF 90000 12800lnA 9830lnt
0,4
… (3.19)
Step 8: The above steps constitute the standard procedures used in the conventional rational method.
The SDF uses calibrated run-off coefficients C2 (2-year return period) and C100 (100-year
return period) from Table 3B.1 instead of determining these from catchment characteristics.
The run-off coefficients for the range of return periods are derived by applying the return
period factors YT in Table 3.12, using the equation below:
C 2 YT C100 C 2
CT … (3.20)
100 2,33 100 100
With the calibrated coefficients being C2 = 10 % and C100 = 40 %, and the return period factors Y20 =
1,64 and Y50 = 2,05 the run-off coefficients are calculated for the 1:20 and 1:50 year return periods:
Step 9: Finally, the flood peak QT (m³/s) for the required return period is calculated as follows:
C I A
QT T T … (3.21)
3,6
Step 10: The SDF hydrograph is triangular in shape with the duration of the rising limb equal to the
time of concentration TC (hours), and that of the falling limb equal to twice the time of
concentration. Use linear interpolation between these two values.
The SCS method requires a minimum input of catchment area, the catchment response time (Time of
Concentration), design rainfalls and soils and land cover classification. The catchment can be divided
into sub-areas or Hydrological Response Units (HRUs) which are areas with similar soils and land
cover and hence have a relatively similar hydrological response to rainfall.
From the ISCW images of soil Land Types, the following Land Types are found in the
catchment: Ib7, Ba3 and Ba9. For each of the Land Types, the dominant terrain unit was
assumed to represent the entire HRU. Similarly, the dominant land class and soil series
within a terrain unit was assumed to respresent the entire unit and, where more than one soil
series is grouped to form the land class, the median SCS soil classification was used. For
each soil series representing each HRU, the appropriate SCS soil classification was
determined using tables from the Visual SCS-SA User Manual. The classification determined
as described above is summarised in Table 3.15, with the catchment divided into four HRUs.
For each HRU in Table 3.15 the soil and land cover classification is used to determine the
initial CN using information contained in Table 3E.3 (in Appendix 3E). From the initial
CN, values of S can be determined using Equation 3.22.
25400
S 254 …(3.22)
CN
The catchment lag is computed from the time of concentration for the catchment using
Equation 3.23.
Given that the SCS method is primarily intended to estimate peak discharge from small
catchments, over which a uniform daily rainfall may be assumed, the reduction of point to
catchment rainfall by means of an ARF is generally not applied. However, to be consistent
with the example, an ARF20 = 94% and ARF50 = 91% for the 20 and 50 year return periods
respectively is used.
The SCS method requires the input of one-day design rainfalls which can be estimated using
a number of different methods. The Visual SCS-SA software includes the following options:
(i) User estimated design rainfalls computed directly from raingauge data.
(ii) Design rainfall extracted for a selected station from Adamson’s TR102 report.
(iii) Using at-site design rainfall calculated from the representative rainfall stations used to
represent the rainfall in each of the 712 hydrological zones.
(iv) Design rainfall estimated using the regional method developed by Smithers and
Schulze (2003).
Given the regional approach used and the longer periods of record from all available stations
used in the analysis compared to all previous studies, it is recommended that the regional,
scale invariance approach to design rainfall estimation developed for South Africa by
Smithers and Schulze (2003) be used to estimate the one-day design rainfall.
The results from the application of this approach at the location of the three rain gauges
located in and close to the catchment are summarised in Table 3.16. Given the relatively
small differences in design rainfall between the values computed from the at-site data and
the RLMA&SI methodology, as well as between the different locations, the RLMA&SI
values for Garsfontein were used in the estimation of the design floods. These were then
reduced by the areal reduction factors of 94 % and 91 % for the 20 and 50 year return period
events respectively in order to determine the catchment design rainfalls.
The ARF is used to calculate the catchment design rainfall, as shown in Table 3.17.
For each sub-area or HRU which the catchment has been divided into, the stormflow is
calculated using Equation 3.24, as shown in Table 3.18. For example, the computation for
the 20 year return period stormflow from HRU 1 is:
(P Ia) 2
Q for P > Ia …(3.24)
P Ia S
where
Q = stormflow depth (mm),
P = daily rainfall depth (mm), usually input as a one-day design rainfall for a
given return period,
S = potential maximum soil water retention (mm),
index of the wetness of the catchment's soil prior to a rainfall event,
Ia = initial losses (abstractions) prior to the commencement of stormflow,
comprising of depression storage, interception and initial infiltration (mm)
= 0,1S (recommended for use in South Africa)
Q 20
111,5 (0,1 75,9)2 60,0 mm
111,5 (0,1 75,9) 75,9
The stormflow from each sub-area is area weighted to calculate an equivalent stormflow
depth from the entire catchment using Equation 3.25 with the results shown in Table 3.19.
N
Q Q A
i 1
i i … (3.25)
where
Q = average stormflow depth from entire catchment (mm)
Qi = stormflow depth (mm) from i-th sub-area calculated using Equation 3.24
Ai = fraction of sub-area of total catchment area for i-th sub-area
N = number of sub-areas
The stormflow volume from the catchment is calculated using Equation 3.26 and is given in
Table 3.20.
QA … (3.26)
V
1000
where
V = stormflow volume (m3x106)
Q = stormflow depth (mm)
A = catchment area (km2)
The peak discharge from the catchment is calculated using Equation 3.27.
0,2083AQ
q … (3.27)
p 1,83L
The computed peak discharge for the 20 and 50 year return periods are summarised in Table 3.21 for
both the manual calculation method, as summarised above, and computed by the Visual SCS-SA
software. The differences in the calculated peak discharge by the two methods are a result of:
(i) No adjustment of the initial CNs into final CNs is included in the manual method. In the
Visual SCS-SA method, the median condition method was used to adjust the initial CNs into
final CNs, thus accounting for typical soil moisture conditions.
(ii) The manual method does not use the incremental unit hydrograph approach and does not take
into account the regionalised typical rainfall hyetographs, both of which are used in the Visual
SCS-SA software.
A standard calculation sheet is included in Appendix 3C for the manual method. The output from the
Visual SCS-SA software is shown in Figure 3.22 and Figure 3.23. The manual method outlined in this
guide is a simplified version of the adaptations to the SCS method for South Africa. The full
adaptations can be performed manually using the SCS-based design runoff user manual developed by
Schmidt and Schulze (1987b). The full adaptations are incorporated into the Visual SCS-SA software,
including the estimation of design rainfall using the RLMA&SI developed by Smithers and Schulze
(2003).
300
250
Discharge (m3/s)
200
150
100
50
0
600 700 800 900 1000 1100 1200
Time (minutes)
20 Year 50 Year
Peak discharges for return periods less than or equal to 100 years can be determined by means of an
empirical deterministic method developed by Midgley and Pitman. The formula reads:
where:
QT = peak flow for T return period (m³/s)
KT = coefficient based on veld-type region (see Figure 3.24 and Table 3.22).
P = mean annual precipitation over catchment (mm/a) (see Figure 3.5 or utilise Design
Rainfall software for the catchment (see Figure 3.11).
Α S
and C (Catchment parameter with regard to reaction time) … (3.29)
L LC
With the catchment being in veld-type zone 8, the KT values are 0,57 and 0,79 for the 1:20 and 1:50
year return periods respectively. The calculated catchment parameter C with regard to reaction time is
0,1238.
The flood peaks as determined with the empirical method for the 1:20 and 1:50 year return periods are
78,85 m³/s and 109,28 m³/s respectively.
Table 3..23 summarises the resultts for the fivve methods used in this ex
xample.
Table
T 3.23: C n of solution
Comparison ns
Return period
p
Method
1:20
1 1:50
Raational (Alterrnative 1) 164 224
Raational (Alterrnative 2 184 243
Raational (Alterrnative 3) 128 166
Unnit Hydrograph 129 172
SD
DF 125 183
SC
CS-SA 205
2 269
Em
mpirical 79 109
The secoond worked example is a design floood calculation n for a new low-level
l briidge across the
t Tsitsa
River, w
which runs through
t the Eastern Cappe in a soutth-easterly direction
d (seee Figure 3..25). The
position of the propoosed Tsitsa lo
ow-level riveer bridge is shown on Figgure 3.25.
F
Figure 3.25: Proposed T
Tsitsa low-leevel river brridge positioon
Table
T 3.24: C
Catchment characteristi
c ics
Descriptiion of charaacteristic Determ mined valuee Com mment
Catchmennt area (see Figure
F 3.25) 43318 km² Catchmment area mayy be clearly defined
d
Starts at Anteloppe Spruit, joins j the
Length off longest wattercourse Tsitsanaa River andd further do ownstream
1779,5 km
(see Figu
ure 3.27) joins th
he Ixnu Riveer to form the Tsitsa
River
Height diifference Total height
h differe
rence equals 1 814 m,
(1085-meethod) 5500 m very steeep slopes allong the upper reaches
(See Figu
ure 3.28) of the water
w course
Average ccatchment sllope 00,37% See detaailed descripption below
Distance to catchmennt centroid 885 km
SDF Draiinage basin number
n 23
Based on calculatted averagee from a
Average rrainfall 8 60 mm numberr of weatherr stations in n the T35
drainage basin
Catchmment area fallls within reegions K5
RMF K-ffactor 5,,0 – 5,2 and K6 (assume higghest value).
(See Figgure 3.26)
Figure 3.26:
3 Maxim
mum flood peeak regions in southern Africa from
m Kovács
Figure 3.28:
3 Waterccourse profiile and avera
age slope
Another peak flow event, during which the m measuring weeir was overtopped, occurrred in Octob
ber 1976.
It is estim
mated that thhis peak flow
w was 1347 m
m³/s.
Based on the historical flow records the results from statistical analyses are shown in Table 3.26.
Equations and statistical tables are included in Appendix 3A.
The calcculation sequuence to determine the floood peaks haas been descrribed in the ffirst worked example,
and will therefore noot be repeated d here. The main resultss are presenteed in Table 33.27. In thiss example
the pointt precipitatioon is obtained
d from the w
weather service station selected for thhis basin from
m TR102;
i.e. Station nr. 1804339 @ INSIZW WA.
Peak discharges for return periods less than or equal to 100 years could be determined by means of an
empirical deterministic method developed by Midgley and Pitman. The formula reads:
In this example the catchment falls within Zone 2 (All year) with KT values of K10 = 0,83,
K20 = 1,04, K50 = 1,36 and K100 = 1,60
The flood peaks as determined using the empirical method (Equation 3.30) for the different return
periods are shown in Table 3.28.
Step 2: Identify the region in which the site is located (Figure 3.26). In this example, as shown in
Table 3.24, the region is K6 (higher K-value). Note that the regions on the map refer to the
location of the site and not to the catchment. Only if the site is located near a boundary
between regions would it be necessary to consider adjusting the K-factor.
Step 3: Utilise the equation provided in Table 3.29 to calculate the RMF as 8059 m³/s
Comparing the calculated flood peaks below (Table 3.30) provides an overview of the range of
expected floods. Based on the flood calculations above, a structurally sound low-level bridge structure
was designed to withstand the high floods and provide a safe crossing during the lower floods (Figure
3.31 and Figure 3.32). Further analysis was also performed on the historical flow data to ensure that
the bridge will not be inundated for prolonged periods, cutting off communities from each other.
m description
Problem n Example 4.1
4
Figure
F 4.1: C
Cross-sectio
on of channeel
Determiine:
(i) T
The energy gradient
g (Sf)..
(ii) W
Whether the flow is sub-- or supercrittical.
(iii) T
The average velocity throough sectionn 3.
(iv) Whether the flow is lamiinar or turbullent.
W
(v) T
The normal flow
f depth.
n Example 4.1
Solution 4
Table 4.1:
4 Channeel characteriistics (Exammple 4.1)
Section
n
Paraameter
1 2 3
Areea (A) 28,,5 m² 9,0 m² 499,25 m²
Wetted perimeter (P) 12,,24 m 3,0 m 1 5,42 m
Hydraulic raadius (R = A//P) 2,333 m 3,0 m 33,19 m
Absolute roughness
r (ks) 0,,3 m 0,7 m 0,7 m
Chézy
12R
l
C 18 log 35,44 m½/s 30,8 m½/s 311,3 m½/s
ks
T
The energy gradient (Sf )
BBy assumingg uniform flo ow conditionns the local slope of the channel, S0, may be sett equal to
tthe energy slope, Sf. Con ntinuity of m
mass and eneergy (Chézy equations) ccombined provide the
ffollowing rellationships:
4-48 Hydra
aulic calculatiions
Q total Q
Q total A1 C1 R 1Sf A 2 C 2 R 2Sf A 3 C 3 R 3Sf
Q total 28,535,4 2,33 9,030,8 3,0 49,2531,3 3,19 Sf
Q total (1541 480 2754) Sf
Q 2total
Sf 0,01 m/m (Energy gradient)
4776 2
Q 2 B 477 25,5
2
(Froude)² = Fr 0,952 1,0 thus subcritical
gA 3 9,81(86,75) 3
Q 3 2754 477
v3 5,59 m/s
A 3 4776 49,25
v
Q 477
5,50 m/s
A 28,5 9 49,25
R
A 86,75
2,83 m
P 12,24 3,0 15,42
υ = 1,14 x 10-6 m²/s (kinematic viscosity of water)
vR 5,502,83
Re 13,65 10 6
υ 1,14 x 10 6
Re >> 2000 Highly turbulent
The normal (uniform) flow depth for a given discharge is calculated by the same procedure,
except that the flow depth is the unknown quantity and the energy gradient, Sf, is equal to the
average (near constant) bed slope, S0.
As shown above, the area, wetted perimeter and hydraulic radius can be written in terms
of the unknown depth, Y. If Y is the depth in section 3, then the variables may be written
as shown in Table 4.2.
With a known slope, S0, and flow rate, Q, Y can be solved. In this case Y = 5 m.
4.2 Example 4.2 - Gradually varying river flow (backwater calculation – simple sectional
details)
Determine the flood level at section 3 for a river of trapezoidal section with side slopes 1:2 and
varying bed width. The characteristics of the cross-sections are reflected in Table 4.3.
Using the principle of conservation of mass and energy would solve this problem. It is assumed that
the flow rate is constant at 43,3 m³/s. It is necessary to determine the type of flow to establish the
control, and then to work away from the control.
Although the assumption that uniform flow will be present at the cross-sections is incorrect,
calculation of the “normal flow depth” at each section will give an indication of the type of flow. In
Table 4.4 the “normal flow depths” have been calculated. This is not the solution to the problem but
merely a way to establish the type of flow!
From Table 4.4 it could be concluded (Fr < 1) that the flow will be subcritical and hence that the
control will be downstream.
Now start with the assumption of a flow depth at section 1 (downstream) and work upstream by
applying the continuity of energy as shown in Table 4.5. Assume that the secondary losses will be
negligible.
The flow depth at section 1 is assumed to be 2,258 m. The total energy level at this section is then
equal to 2,427 m.
Problem
m description
n Example 4.4
4
Solution
n Example 4.3
4
2g 29,81
4-52 Hydra
aulic calculatiions
v 22
Consequently H 3,0 0,212 y 2 can be assumed.
2g
q v 2 y 2 0,167 m 3 /sm
v2
v 2 = 7,9 m/s; y2 = 0,021 m and Fr2 17
gy 2
The actual velocity will be slightly lower. (If H>> y2, then v2 2gH )
The normal (uniform) flow depth in a long 2 m wide, rectangular canal is 2 m and the normal flow
velocity 2 m/s. The Manning n-value is 0,02 s/m1/3. There is a 90° bend with a centre-line radius of
7 m. Calculate the Froude number for uniform flow conditions.
Calculate the flow depths just upstream and just downstream of the 90° bend.
vn 2,0
Frn 0,452 (based on the normal flow depth)
gy n 9,812,0
Downstream control and hence the depth just downstream of the bend will be 2 m.
H1
2,02 2,0 0,1165 2,321 m
29,81
v12
y1 2,321 m
2g
Write the velocity in terms of the upstream flow depth (y1) using continuity:
Problem
m description
n Example 4.5
4
Water flows across a 16 m wide road. The rooad has cross-falls of 2%%. Calculate tthe dischargee per unit
width thhat would floow across th
he road whenn the adjaceent level risees 0,5 m aboove the shouulder, see
Figure 44.3.
Solution
n Example 4.5
4
Assume that the conntrol (point of release) occurs at th he crown (B). If there are no energgy losses
between A and B, thhe specific energy
e at thee top of the crown will be 0,5 – 0, 16 = 0,34 m and the
critical ddepth yc = 0,2227 m (two thirds
t of the specific enerrgy).
q v c y c g y c2 0,338
0 m 3 /s.m
m Assuume that the actual disch
harge is 0,3 m
m³/s.m
q2 and 3
yc 3 0,209 m Ec y c 0,314 m
g 2
5 1
y 3S 2
q
n
Manning n-vaalue = 0,013 s/m1/3
With a M
5 1
y 3 0,02 2
0,3 y n 0,116 m y c 0,209 m
0,013
vy 0,3
Re 6
2,6 10
0 5 2000 , indicating the
t flow is tu
urbulent.
υ 1,14 100
4-54 Hydra
aulic calculatiions
v A2 v2
0,5 y A 0,35 A (Energy equuation with prrovision for transition
t lossses)
2g 2g
Thus yA = 0,472 m
dE
So Sf (Energy equation for prismaticc channels)
Δx
dE v2 v2
x dE y 1 1 y 2 2
So Sf 2g 2g
v2n 2
So = -0,002 m/m (uphill) Sf 4
y 3
ΣΔx = 8,64 > 8,0 thuus the dischaarge per unit width should be less ([Link]. q should bee smaller).
Choose ssmaller q andd repeat until ΣΔx ≈ 8,0 m for a moree accurate annswer.
4-55 Hydra
aulic calculatiions
5 SURFAC
CE DRAINA
AGE
m description
Problem n Example 5.1
5
The aim
m should alw ways be to liimit the flow
w depth on the maximum of 6 mm to
t road surrface to a m
prevent hhydroplaningg.
n Example 5.1
Solution 5
F
Figure 5.1: Depth
D of sheeet flow on rroad surfacee (Laminar flow
f conditioons assumed
d)
5-56 Surface
S drain
nage
Sf 2 2 6 2
Sf 6,322%
d 4,6xx102 L f I
0,55
S
0.2
f …(5.3)
…
d 4,6xx102 31,62100 0,5
0,066320.2 4,499 mm
5.2 Worked Example
E 5.2 – Capacity of side chan
nnel
m description
Problem n Example 5.2
5
n Example 5.2
Solution 5
Y – Y1 )
X = ZA(Y …(5.4)
…
X = (20))(100 – 40)
X = 12000
ZBY1 = ((40)(40)
ZBY1 = 11600
Top widdth, T = 12000 + 1600 = 28
800 mm
T
Z
Y
5-57 Surface
S drain
nage
Q = QA+ QB see calculations
c in Table 5.11.
S = 0,05 m/m
s 1/3
Manningg n = 0,015 s/m
Ta
able 5.1: Ch
hannel flow characterist
c tics
Parameterr Section A SSection B
A A 0,,50,1 0,041,2 A B 00,50,040 1,6
Crosss-sectional arrea (m²)
A A 0,0084 A B 00,032
m description
Problem n Example 5.3
5
Determinne the flow capacity of a drop grid inlet, dimennsions of 0,9 by 0,6 m an
and a submerrgence of
0,2 m iff the approacching flow is subcriticall. For compaarison with Figure
F 5.6 aassume that the inlet
coefficieent = 0,8 andd the blockag
ge factor, F = 0,5.
n Example 5.3
Solution 5
Figure 55.3 provides the relationsship of flow rrate for an orrifice controll or a broad-ccrested weir..
Q CFA 2gH
H …(5.5)
…
w
where:
C = inlet coefficient (0,6 forr sharp edges or 0,8 for rounded
r edgees)
F = blockage faactor (say 0,55)
A = effective crross-sectionaal plan area of
o the openinng (m²)
H = total energyy head abovee grid (m)
5-58 Surface
S drain
nage
Fig
gure 5.4: Exxample of ty
ype of grid in
nlet
Q 1,777A H
A 0,90,6 0,54 m
m²
Q = 0,4227 m3/s
Problem
m description
n Example 5.4
5
Solution
n Example 5.4
5
5-59 Surface
S drain
nage
P Y Y 2 T 2 0,1 0,12 5,02 5,101 m
4
S 0,04 m/m
100
5
Q
1 0,25 3 0,04 0,446 m³/s .
0,015 5,101 23
Strictly speaking, the very wide section with variable velocities should be subdivided into narrower
sections.
Determine the maximum flow depth and velocity in a wide channel with a slope, S is 2%. The channel
is lined with stones (relative density 2,65) and representative size (more than 50% by mass) of
250 mm.
R=D=y
where:
R = hydraulic radius (m)
D = y = flow depth (m)
From Chezy:
12R
V 18log RS
k
s
…(5.6)
d1 > 11 DS …(5.7)
Dy
0,25 1,14 m
110,02
Now the velocity can be calculated using Equation 5.6:
121,14
v 18 log 1,140,02
0,25
v = 4,71 m/s
Problem
m description
n Example 5.6
5
Solution
n Example 5.6
5
For stablle side slopes the particlee size (d2) shhould at least be:
8,3Ds
d2
tan 2 θ
cossθ 1
tan 2 φ
…(5.8)
…
where:
F
Figure 5.5: Required
R sizzes of the stoone for erossion protectiion of loose bbed channells
(Thhe side slopee, θ, should always
a be sm
maller than th
he angle of reepose, Ø, to eensure stability.)
d2
8,3 1,8 0,001
0 0,028
0 m
tan 2 25
cos 25 1
tan 2 30
5-61 Surface
S drain
nage
6 LOW LE
EVEL CRO
OSSINGS
Problem
m description
n Example 6.1
6
The apprroach gradieents of the ro oad are modeerate and therre is no horizzontal curvat
ature. The preeliminary
design oof the verticall alignment of
o the road aacross the stru
ucture to be provided hass also been done.
d The
straight section in thhe middle (L2) has a lenggth of 20 m, and K1 and K3 are both 4 m (refer to t Figure
6.1). Thee slope of thee road on thee southern baank is - 5,6%
%, and on the northern bannk 7,0%.
Solution
n Example 6.1
6
The design level is determined as per Secttion 6.3.3. Design D level 1 is taken as the initiaal choice.
Because of the expected traffic volumev of 3000 vehicles per
p day exceeeding the suuggested 250 0 vehicles
per day, the design level
l is increeased to leveel 2. This is supported
s by
y the availabbility of an allternative
route of length less than
t the sugg gested 50 kmm. As the critteria suggestted for designn level 3 aree not met,
design leevel 2 is seleected.
6-62 Low
w level crossiings
The design flow rate is determined from Equation 6.1:
Q design f i Q 2 …(6.1)
From Table 6.1 follows that f2 is 0,50 and Q2 is the discharge with a 1:2 year return period (120 m³/s).
Cross-section
With a design period of 20 years and 2% growth in traffic per year, the anticipated 300 vehicles per
day is expected to increase to 446 vehicles per day after 20 years. Because of this being less than the
suggested 500 vehicles per day and visibility being good, a single-lane structure is opted for. The
cross-fall in the direction of flow is taken as 2%.
Selection of structure
Because of good, but uneven founding conditions a low-level bridge is opted for. Six spans of 6 m
each fit the river cross-section well. Piers are 300 mm thick.
Hydraulic calculations
The capacity of the structure is determined as the sum of the flow that can be accommodated over the
structure and through the structure.
Assume supercritical flow and decide on a maximum flow depth of 0,1 m (d). The flow that can be
accommodated over the structure is determined from Equation 6.2.
A 5/3 1/2
over S 0
Q over 2/3
…(6.2)
n Pover
A over A 1 A 2 A 3 , or
1 1
A over d 800K1d dL 2 d 800K 3 d …(6.3)
3 3
and
1 1
Pover 800K 1d L 2 800K 3 d …(6.4)
2 2
where:
Qover = the discharge that could be accommodated over the structure within the
selected flow depth (m3/s)
Aover = area of flow over structure at the flow depth selected (m²)
S0 = slope in direction of flow, for example 0,02 or 0,03 m/m
n = Manning n-value. For a concrete deck nconcrete
Pover = wetted perimeter at the flow depth selected (m)
A1, A2, A3 = the areas defined in Figure 6.1 (m2)
d = depth of flow over the structure (m)
K1 = the geometric K value for vertical curve 1
K2 = the geometric K value for vertical curve 3
With K being a vertical road alignment parameter, defined as the horizontal length of
road required for a 1% change in the gradient of the road.
S0 is 0,02 (2% as above) and Manning n roughness parameter for concrete is 0,016 s/m1/3. The cross-
section area of flow is determined using Equation 6.2.
1 1
A over d 800K1d dL 2 d 800K 3 d where K1 and K3 are both 4 and L2 is 20 m (Equation 6.3)
3 3
Aover = 3,19 m2
1 1
Pover 800K 1d L 2 800K 3 d = 37,89 m (equation 6.4)
2 2
Establish whether flow is indeed supercritical by calculating the Froude number (Equation 6.5):
2
Q over B …(6.5)
Fr
g A 3over
where:
B = L1 + L2 + L3 (m), the width of the channel (or the length of the structure)
g = 9,81 m/s², the gravity constant
1 1
and L1 800K 1d and L 3 800K 3 d
2 2
Fr = 1,87, which is > 1,0 m, confirming supercritical flow over the deck of the structure.
H1 H 2 …(6.7)
…
v under 2
C n eff LB
4/3
2gg R
w
where:
Aeff = the effective inllet area throu ugh the struucture (m²) = ΣAcell (the effective
inlett area througgh the structu
ure)
LB = the total
t width o f the deck off the structure (m)
v undeer = the velocity
v of fllow through the structuree (m/s)
C = facto
or that refleccts the transittion losses (E
Equation 6.8))
6-65 Low
w level crossiings
By applying the conservation of energy principle, determine the depth upstream of the structure, h as
shown in Figure 6.2, that is required to pass the flow rate, Qover:
v 22
h d …(6.9)
2g
H1 = h + x + D …(6.10)
where:
x = the thickness of the deck (depending on the structural design outcome) (m)
D = the height of the soffit of the deck above the river invert level (m)
H 2 D L BS0 …(6.11)
where:
LB = the total width of the deck of the structure (m)
S0 = slope of the conduit underneath the structure (m/m)
Determine the depth upstream of the structure, h, that is required to pass the flow rate, Qover using
Equation 6.9.
v2
h 2 d , where v 2 = Qover/Aover = 1,70 m/s
2g
h = 0,247m
Assume Kinl = 0,5 and Kout = 1,0 (both sudden transitions), then
C = 6 x (0,5 + 1,0), see Equation 6.8.
C=9
n eff
n cell Pcell
Peff
n eff
60,02214,2
85,2
neff = 0,022 s/m1/3
R = 0,562 m
0,562 4/3
4/3
2g R 2(9,81)
v under = 1,344 m/s
Design discharge
The capacity of the structure at the design level Qover + Qunder = 69,8 m³/s
As Qover + Qunder is larger than Qdesign (60 m³/s), the design is complete as the structure is adequate. If
this was not the case, the level of the deck would have to be adjusted, and the calculation be redone.
Problem
m description
n Example 7.1
7
Figu
ure 7.1: Desiign flood freequency estimate
No signiificant debris is anticipatted since thee catchment area consists mainly of grassland. The T client
favours tthe use of cirrcular (pipe)) culverts. Thhe final levell of the roadw
works acrosss the river will be at a
level of 22,5 m above the riverbed d.
Solution
n Example 7.1
7
First calcculate the noormal flow depth, Yn andd flow condittions in the channel upstrream and dow
wnstream
of the planned culverrt.
D
Determine the
t upstream
m normal floow depth
QC RS A
…(7.1)
…
where:
Q = flow raate (m3/s)
C = Chèzyy constant
R = hydrauulic radius (m
m)
A = m2)
area (m
Yn upstream = Y = uppstream norm mal flow dep
pth (notation used here) (m
m)
12R
44,5188 log RS A
ks
12A AS
44,5 188 log A
k s P P
Solve th
he upstream normal flow
w depth Y in
n the above equation.
In a similar way as above, the downstream normal flow depth can be determined. In this case
the downstream normal flow depth, Yds = 1,541 m, Frds = 0,841 and the flow is subcritical.
2
Total energy head upstream of the culvert, H 1 v Y 2,175 m
2g
Determine the size of the culverts to manage the flow
The height difference between the river bed and the final road level is 2,5 m. If the optimum H/D ratio
of 1,2 is used the maximum vertical dimension of the culvert (D) is 2,5/1,2 = 2,08. Based on Figure
7.4 multiples of 1,8 m diameter pipe culverts will be used.
For a culvert with a diameter of 1,8 m and the downstream flow depth of 1,541 m, the flow will
probably be inlet controlled (to be verified) and the flow can be evaluated based on the relationship for
inlet control (Table 7.1).
Table 7.1: Relationships for the flow rate under inlet control
0,05 1,9
Q S H1 2 2
0,48 0 D Q C B BH1 gH1
D 2
gD 0,4 3 3
Where: CB = 1,0 for rounded inlets (r > 0,1B)
CB = 0,9 for square inlets
And for: 0,8 < H1/D ≤ 1,2 * And for: H1/D > 1,2
Q C h BD 2gH 1 C h D
0,05
1,5
Q S0 H1
0,44
0,4 D
2
D gD
Where: Ch = 0,8 for rounded inlets
(S0 = slope of culvert bed with slight
Ch = 0,6 for square inlets
effect on capacity)
Note:
* For H1/D > 1,2, the orifice formulae
D
applies Q C D A 2 g H 1 with
2
CD ≈ 0,6
H 1 2,5
In this example, the maximum 1,39 and the flow rate through a culvert can be determined.
D 1,8
Assume the distancee between thee pipes is 1000 mm and th he wall thickness of the ppipes is abou
ut 78 mm,
then the total width of six culveerts will be = (6)[1,8+(2 2)(0,078)] +(5)(0,1) = 1 2,236 m. With
W some
groundw n in Figure 7.3
work it is posssible to place the culvertts as is shown 7 .
Figure 77.4 reflects the requiredd culvert sizee for a given n (design) flow
fl rate andd a H/D ratiio of 1,2.
Referencce to pipe annd portal (box
x) culverts arre reflected here.
h
By meanns of Figure 7.4 and by assuming thaat a portal (b box) culvert could be useed as an alterrnative to
the calcuulation abovee, the requireed culvert sizze for Inlet Control
C condiitions could bbe obtained.
Q C h BD 22gH1 C h D
…(7.3)
…
Where Ch = 0,8 for rounded inleets and Ch = 0,6 for squarre inlets.
It follow
ws from the new
n slope deetails (S0 upstrream and S0 doownstream is 0,0
0015 m/m) tthat the upstrream and
downstreeam normal flow depths will be 2,0 m, as was determined
d beefore. If the upstream water level
is limited to a maxim
mum of 2,5 m to prevennt the inundaation of the road, r the cullvert flow rate can be
determinned as follow
ws.
For inlett control connditions the length, rougghness, slopee and hydrau
ulic radius off the culvertt have no
F outlet coontrol these variables do influence thhe flow rate and have
influencee on the disccharge rate. For
to be connsidered.
7,4417
v 2,9155 m/s
π( 0,9) 2
h l1- 2 (K
K inlet K outleet )
v2
0,5 1
2,915 = 0,649 m
2
2g 29,81
1 3,7D
2log( ) …(7.7)
λ ks
λ 0,02015
h f1 2 = 0,121 m
For outlet control conditions, the upstream conditions can now be determined (Equation 7.4):
H1 H 2 Z1 h l1 2 h f1 2
H1 = 2,908 m
The value of H1 for outlet control is greater than the maximum allowable damming height of 2,5 m,
hence outlet control will be maintained through the culvert. The only way to reduce the upstream flow
depth is to consider more culverts of similar dimension in parallel or larger culverts.
Assume that seven culverts will be used. With the flow of 44,5/7 = 6,357 m³/s ( v = 2,498 m/s) and the
value of the losses, H L h l1 2 h f1 2 = 0,477 + 0,089 = 0,566 m, the upstream conditions can now be
determined.
H1 = H2 - Z1 + HL
H1 = 2,704 m (which is still greater than the allowable upstream damming height)
To illustrated the use of Figure 7.4, start at Q = 6,357 m³/s (flow rate in each of the 7 culverts) and
HL = 0,566 m, then it can be seen that the 1 800 mm pipe culvert is still insufficient to transport the
flow (Figure 7.7). Figure 7.4 could thus be used to consider other culvert sizes.
Placing the seven culverts in parallel will also however result in a section width of about
(7)[1,8+(2)(0,078)] +(6)(0,1) = 14,3 m, which is much wider than the river base of 8,0 m.
The number of 1,8 m pipe culverts required to prevent overtopping of the road will be eight (although
still marginally undersize). Alternatively portal culverts with a larger vertical dimension (2,1 m) could
be used following the same procedure as above.
An economic/technical assessment of the alternatives i.e. re-alignment of the road against provision of
a practical culvert design have to be conducted to select the solution for implementation without
changing the risk of failure.
Problem
m description
n Example 7.2
7
The flow
w releases intto a natural trapezoidal
t rriver section with a base width of 2,00 m and side slopes of
1V:2H. T The natural slope of the river is 0,0004 m/m and the roughnesss is 0,05 m.. Details of the
t cross-
section aare given in Figure
F 7.8.
Fig
gure 7.8: Up
pstream view
w of the culv
vert
Solution
n Example 7.2
7
vC R
RS
…(7.8)
…
122R
Q vA 18 log RSA
k s …(7.9)
…
where:
Q = flow w rate (m3/s)
v = averrage velocity y (m/s)
R = hydrraulic radius (m)
ks = absoolute roughneess (m)
S = slope of the riverr section (m//m)
Table 7.22 reflects thee cross-sectio
onal parametters for a circcular pipe.
vC R
RS
…(7.8)
…
12RR
Q vA 18 log RSA
k s …(7.9)
…
where:
Q = flow w rate (m3/s)
v = averrage velocity y (m/s)
R = hydrraulic radius (m)
ks = absoolute roughneess (m)
S = slope of the riverr section (m//m)
T
Table 7.2: Sectional parrameters forr a circular cross-sectionn
V
Variable Y < D/2 (θ < π radials)) Y > D/2 (θ > π radia
als)
Cross--sectional viiew
12R
8,9 18 log R 0,004 A
0,05
with:
A 2y n 20,5 y n 2y n
P 2 2 5 y n
A
R
P
where:
yn = the unknown
u norrmal flow deepth (m)
The critical slope in the culvert can be determined. Critical conditions will occur when Fr = 1, and the
critical slope Sc can be determined for full flow conditions as follows:
v2 Q2
Sc 2 2 2 …(7.10)
C R C A R
Sc = 0,00792 m/m
This indicates that the flow depth in the culvert will reduce downstream from the position where
critical flow occurs near the inlet of the culvert because S0 > Sc.
If it is assumed that the flow depth in the culvert Y is more than D/2, solve for a potential flow depth,
Y.
It is found that:
y = 1,330 m; A = 2,015 m² ; B = 1,582 m (top width of flow); v = 4,417 m/s and Fr = 1,25
(subcritical).
Figure 7.9 can now be used to select the appropriate erosion protection.
In this example the appropriate protection falls in the Type III. Figure 7.10 can now be used to obtain
dimensions for outlet erosion protection.
Problem
m description
n Example 8.1
8
Fig
gure 8.1: Up
pstream view
w of the brid
dge
Determiine:
(i) C
Characteristiics of the uncconstricted fflood state
(ii) T
The flow typpe
(iii) Bridge opening ratio
B
(iv) V
Velocity heaad coefficients
(v) C
Calculate baackwater
Solution
n Example 8.1
8
Utilising the Manning equation the normal flow depth can be calculated, yn = 3,157 m.
150 2 36,06
1 1
Q2B
2 2
Frn =
gA 3 3
n 9,8186,30
= 0,359 < 1
Q 150
vn =
A n 86,30
= 1,738 m/s
Esn = yn
vn
2
3,157
1,738
2
2g 29,81
= 3,311 m
150 2
1 1
Q2
3 3
y2c = 2 = 1,995 m
9,81177 2
gb
Q 150
v 2c = = 4,424 m/s
y 2c b 2,817 17
Esc = y 2c
v 2cc
2
1,995
44,4242 = 2,992 m < Esn indicaating Type I flow (see
2gg 29,81
Figure 8..2).
Because the values of Esn and Esc are fairlyy close, and other lossess are so far ignored, it would
w be
prudent tto check Typpe I and Type II flow.
17 Q b 1110,25
Qb = 120,69 M =
17 1,61 Q 1150
= 110,25 m3/s = 0,735
α1 = qv 2
2
Qv n
= 1,20
α2 = om Figure 8..3)
1,15 (fro
For Typ
pe I flow:
F
Figure 8.4: Chart
C to deteermine the backwater
b coefficient,
c K
K*
Ap = Wp y n An2 = (bcosø)(yn)
= (2)(3,157) = 17cos15 3,157
o
= 6,314 m² = 51,84 m²
Ap 6,314
J =
A n2 51,84
= 0,122 (use J = 0,1 from Figure 8.4)
Eccentricity
e = Qa 12,14
1 1
Qc 17
17,17 120,69 1
17 1,61
= 0,56
K* = 0,80
= Q 150
v n2
A n2 51,84
= 2,894 m/s
2
v
h1*1 = *
K α 2 n2 …(8.2)
2g
0,801,15 2,894
2
h1*1 =
29,81
= 0,393 m
bc = b W 17,0 2,0
p
= 15,0 m
Cb = om Figure 88.5
0,152 fro
Figure 8.5
5: Estimation
n of the Bacckwater Coeefficient, Cb
1 1
Q 2 3 150 2 3
y2c =
gb 2 9,811152
c
= 2,168 m
A n2 51,84
y =
b 17
= 3,050 m
=
1,154,6122 0,152 1 1,201,7382 2,168 3,050
29,81 29,81
= 1,436 - 0,185 + 2,168 - 3,050
= 0,371 m
150
v1 =
99,67
= 1,505 m/s
h1*1 = 1,436
1,201,5052 2,168 3,050
29,81
= 0,417 m
Although the difference in this case is negligible, to be conservative, the higher value should be used.
From the calculations h1*1 for Type II flow was 0,417 m which is less than the backwater calculated
for Type I flow, thus Type I flow prevails i.e. h1*1 = 0,441 m.
Note that this example was also modelled in HEC-RAS (provided on the supporting flash drive) and
that the highest backwater was obtained by the Standard Step Energy Method. The backwater is
300 mm, which is less than the value of 441 mm obtained above. However, in this model the
ineffective flow area option had been used. The model was then re-run, with this option switched off
and a higher backwater of 580 mm was obtained in the revised model.
Users of HEC-RAS should therefore carefully consider the option where the bridge approach
conditions are “smoothed”, thereby reducing the backwater.
Consideration is being given to construct a bridge across the Sand River, which is some 730 m wide at
the proposed bridge site. The potential scour at the bridge should be determined. This problem was
also evaluated with HEC-RAS, and the data files are included with the data files as [Link].
Figure 8.6 shows a plan view (obtained from the HEC-RAS problem evaluation) and the position of
the bridge relative to the other cross-sections. The cross-sectional information for all the sections is
available. The bridge will be positioned at cross-section 6.5 (downstream from cross-section 7 and
upstream of cross-section 6).
The briddge opening between the sloping abuutments is ap pproximately 126,61 m w wide and the bridge is
supporteed by five piiers, each with
w a width of 1,5 m (eequally spacced). The higgh (road surface) and
low cordd (bridge sofffit) values fo
or the bridge deck on the upstream sid
de are 6,7 annd 5,5 m resp
pectively.
The userr can open thet project (Example2.p
( prj) in HEC-RAS and by y selecting thhe appropriaate icons,
review thhe bridge datta which is not
n repeated hhere in detail.
D50 = 0,00200 m
D90 = 0,00455 m
Cross-seection detailss
Determiine
Solution
n Example 8.2
8
It is estimated that the bed roughness under flood conditions will be 0,002m, equal to D50 the
representative sediment material size.
Top flow width = 698,2 m for the calculated normal flow depth of 2,98 m.
It is assumed that the bed material consists of deep alluvial sand with no cohesion.
The regime equations are applied to establish equilibrium conditions at the design flow:
B 14Q 0,5 D 50
0,25
Fs0,5 …(8.3)
where:
B = mean channel width (m)
y = mean depth of flow (m)
Q = equivalent steady discharge which would generate the channel geometry (m3/s)
q = discharge per unit width (Q/B) (m3/s.m) (Note: To estimate channel geometry
conditions under flood conditions the design flood flow may be used.) (8.9)
D50 = median size of bed material (m)
Fs = side factor to describe bank resistance to scour (Table 8.4)
With Fs = 0,1 from Table 8.4 for sandy loam, the width B can be calculated.
Use Equation 8.4 to determine the mean flow depth at the equilibrium width:
0,17
y 0,38q 0,67
D 50 …(8.4)
where:
y = mean depth of flow (m)
D50 = median size of bed material (m)
The maximum live-bed depth, Ymax, is slightly less than the fixed bed depth, yn of 2,98 m, which
reflects that no general short term scour will occur.
It has been indicated that contraction scour can be determined by applying either the regime
equations (Equations 8.4 for the case of an alluvial channel) or the contraction equations (Equations
8.5 with 8.6 or 8.7).
v 22 v12
d s y 2 y 1 1 K
…(8.5)
2g
Sediment-laden flow
6/7 2/3 1/3
y2 Qt B1 n2
…(8.6)
y1 Q c B2 n1
3/7
Q2
y2 ...(8.7)
2/3
40D m B 22
First apply the regime equation on the reduced width of 126,61 m. In this case q = 850/126,61 =
6,714 m³/s.m leading to a mean flow depth, y of 3,915 m. From Table 8.5, ymax can be determined as
follows: ymax = 1,25 x 3,915 = 4,893 m. This reflects a scour depth, ds = 4,893 – 2,98 = 1,913 m
Table 8.5: Factors to convert mean flow depth (y) to maximum channel depth
Description Multiplying factor
Straight reach of channel 1,25(*)
Moderate bend 1,50
Severe bend 1,75
Right-angled abrupt turn 2,00
Note: * Neill recommends that this factor be increased to 1,50 in cases where dune movement
takes place on the riverbed.
Secondly the contraction equations are used to determine the scour depth after it has been established
if the flow will be sediment laden or not.
and the term, V*D50/ν = 483 >> 13, thus in turbulent flow region (see Figure 8.8).
12R
Vc 5,75V*cc log …(8.10)
…
ks
12R 0,937
From Eqquation 8.10:: Vc 5,75V
V*c log 5,75 0,029log12 0,6255 m/s
ks 0,00
02
Figure 8.10:
8 Long constriction
c in sedimentt-laden flow
w: definition of terms
6/7
y 2 850
1,469 (widths and n-values are equal for these sections)
y1 542,6
y2 = (1,469)(2,98) = 4,378 m
Assuming a level bed with a total depth of 4,378 m, the velocity in the contraction can be determined:
850
v2 1,63 m/s
4,378126,61 51,5
Note that in this case the downstream area is 521,5 m², calculated as follows (4,378 x (126,61 –
5(1,5))). This is larger than the upstream main channel area of 258,73 m² (Table 8.3), and thus the
flow is expanding. Equation 8.12 is used to determine the contraction scour depth.
v 2 v12
d s y 2 y1 1 K 2 and with K = 1 for a sudden transition
…(8.12)
2g
1,63 2 1,299 2
d s 4,378 2,98 1 1,0
29,81
ds = 1,468 m
This scour depth (1,468 m) is less than that obtained with the regime theory (1,913 m).
Use Equations 8.13 and 8.14 to compute local scour in two different ways for alluvial channels
(cohesionless material). Obtain the factors needed for Equation 8.13 from Table 8.6 and Table 8.7.
Obtain the factors needed for Equation 8.14 from Table 8.6 to Table 8.8 and K4 in cases where
armouring is applicable. Compare answers obtained from Equations 8.13 and 8.14 and select a
conservative answer using good engineering judgement.
Table 8.7: Correction factor K2, for angle of attack of the flow
Angle (skew angle
L/b = 4 L/b = 8 L/b = 12
of flow)
0 1,0 1,0 1,0
15 1,5 2,0 2,5
30 2,0 2,75 3,5
45 2,3 3,3 4,3
90 2,5 3,9 5,0
Note: In the case of L/b larger than 12, the ratio’s for L/b = 12 should be used.
From Equation 8.13, with depth y0 in the bridge section as determined from regime:
d s 1,8y 0,75
0 b
0,25
y0 …(8.13)
where:
ds = local scour depth at pier (m)
y0 = depth upstream of pier (m)
b = pier width (m)
Note that the scour level is (3,915 + 1,629) = 5,544 m below the design flood level.
0,35
ys y
2,0 K 1 K 2 K 3 K 4 1 Fr10,43 …(8.14)
b b
with
b = 1,5 m
y1 = 2,98 m (normal flow depth upstream of the bridge, Table 8.3)
Fr1 = 0,429 based on main channel data directly upstream of the pier
K1 = 1,0 for round nose
K2 = 1,0 for zero skew angle
K3 = 1,1 for small dunes
K4 = 1,0 for uniform sediment (no armouring), see Drainage Manual for details of
calculating K4 correction factor, then
0,35
ys 2,98
2,0 1,0 1,0 1,1 1,0 0,429 0,43 1,943
b 1,5
ys = 2,915 m
Note that the scour level is (2,98 + 2,915) = 5,895 m below the design flood level associated with
the normal flow depth and a fixed bed level.
Apply factors in Table 8.9 to the general short-term average scour depth obtained from Equation 8.13.
From Table 8.9 the factor for flow that impinges at right angles on bank = 2,25; hence the scour at the
abutments can be determined as shown below:
Total scour is the sum of the long and short-term general scour, contraction scour and local scour.
Table 8.9 reflects a summary of all the calculated scour depths.
Table 8.9: Summary of the calculated scour depths for Worked Example 8.2
Scour type Calculated scour depth, ds (m)
Short term general scour No scour
Regime equation 1,913
Contraction scour
Contraction equation 1,468
Piers 2,915
Local scour
Abutments 2,102
Piers 4,828
Total expected scour
Abutments 4,013
The potential general scour at the bridge has been determined in (i) using Equations 8.5. A more
correct approach is to estimate contraction scour separately for the main channel and over banks, as is
done in HEC-RAS, where the over bank flows may then reflect clear water scour. The scour depth in
the channel calculated in the approach used in HEC-RAS is less than 3,1 m. With the regime theory
reflecting a scour depth of 1,91 m and the HEC-RAS result of 3,1 m, the contraction scour calculation
of 1,468 m using Equation 8.12 is too conservative and hence discarded.
Based on the summary in Table 8.9 the total scour can be determined as follows.
Total scour level at piers, below design flood level (not accounting for backwater)
= 1,913 + 2,915 = 4,828 m
With the right abutment on the edge of the main channel, the scour would be the sum of main channel
contraction scour plus abutment scour, thus:
Total scour at abutments = 1,913+ 2,1 = 4,013 m below design flood level
Verify the scour depth with the method based on the principle of applied stream power
For total scour at piers in alluvial rivers, check the answer against values obtained by means of
Equation 8.15 that is based on the principle of applied stream power.
With F = 0,8; ks = 0,002 m ; vss = 0,24 m/s; q = 542,6/126,61 = 4,29 m³/s.m and C calculated from
the Chezy relationship for total section, C = 67,5; it follows that:
Yt = 2,03 m below design flood level, which is substantially less than obtained above.
The designer will experience these conflicting results, which reflect amongst other the complexities
involved in the mathematical description of scour estimates and the shortcoming in the assumption
that the material is cohesionless.
Considering the risk of failure of the structure due to scour and the potential consequences, these cases
require further evaluation by experienced persons.
The problem is also evaluated using HEC-RAS, and the data files are contained as Example 2 in the
supporting data files.
Water flows from a submerged catch pit at a kerb inlet through a 20 m concrete drain pipe where it
releases into a river stream. The inside diameter of the pipe is 0,3 m. The pipe has a slope of 0,5 % and
the depth of the water in the catch pit above the top of the pipe is 0,2 m (i.e. available head).
To solve this problem the conservation of energy principle will be used i.e. Equation 9.1. The entrance
(point 1) and exit (point 2), are open to atmosphere and thus p1/γ = p2/γ = 0 if the stream line is
selected along the top of the water surface.
v12 p 1 v 22 p 2
z1 z 2 h f1 2 h ι 1 2
2g γ 2g γ ...(9.1)
The velocity in the catch pit ( v1 ) is negligible (large catch pit area in comparison with pipe area). If
the datum line is selected through the invert of the pipe outlet (point 2) z2 will be zero and z1 can be
calculated as follows:
The first term is the height above the pipe inlet and the second term the length of the pipe multiplied
by the slope of the pipe to obtain the vertical height difference between inlet and outlet.
The last two terms in the energy equation are the loss terms. The friction loss can be calculated using
the Chezy equation (Equation 9.2) with the C-value determined by Equation 9.3. The absolute
roughness value of the pipe was estimated as ks = 0,0004 m.
v2L
h f 1- 2 …(9.2)
C2R
D
12
12R 4
C 5,75 g log or C 5,75 g log for circular pipes …(9.3)
ks
ks
0,3
12
4
C 5,75 9,81 log 60,37
0,0004
v 22 20
h f1- 2 0,07317v 22
60,37 2 0,3
4
Well-rounded 0,2
θ = 15° 0,2
Converging sections
kv 2 Sudden 0,5
h1
2g 0,25
Cone
Bends
kv 22 θ = 90° 0,4
h1
2g θ = 45° 0,3
Outlets
2
kv 2 A Sudden 1,0
h1 1 1 1
2g A2
The diameter does not change and the velocity where the water enters the pipe can be assumed to be
equal to v 2 .
kv 22 0,5 v 22
h1 0,02548 v 22
2g 29,81
v 2 = 1,416 m/s
0,3 2
Q v 2 A 2 1,416 π 0,1 m /s
3
2
This pipe can discharge 100 l/s if it is allowed to dam at the entrance up to 500 mm above the pipe
inlet level.
The EPA Storm Water Management Model (SWMM) is a dynamic rainfall-runoff simulation model
that computes runoff quantity and quality from primarily urban areas. The runoff component of
SWMM operates on a collection of subcatchment areas that receive precipitation and generate runoff
and pollutant loads. The routing portion of SWMM transports this runoff through a system of pipes,
channels, storage/treatment devices, pumps and regulators. SWMM tracks the quantity and quality of
runoff generated within each subcatchment and the flow rate, flow depth, and quality of water in each
pipe and channel during a simulation period comprised of multiple time steps.
SWMM was first developed in 1971 and since then has undergone several major upgrades. It
continues to be widely used throughout the world for planning, analysis, and design related to storm
water runoff, combined sewers, sanitary sewers, and other drainage systems in urban areas and has
also been used for modelling non-urban areas. The most current implementation of the model is
version 5.0 which was released in 2005. It has modernized both the model’s structure and its user
interface, making SWMM easier to use and more accessible to a new generation of hydrologists,
engineers, and water resources management specialists.
All the practical exercises after this first tutorial are developed for the same catchment area and each
one builds in some degree on the results of a previous example. Therefore, it is recommended to begin
with Exercise 1 and work sequentially through up to Exercise 6, while hopefully building the required
input data files and running them with SWMM for each exercise. These files, as well as the backdrop
image files that are needed to complete the exercises, are available on the flash disk.
Goal: This first tutorial provides an introduction to using EPA SWMM, Version 5,
for modeling the quantity of storm water runoff produced from an urban area. The
topics to be covered include:
Project Setup
Constructing a SWMM Model
Saving and Opening Projects
Setting the Properties of SWMM Objects
Running a Single Event Analysis
Viewing Simulation Results
In the exercises the decimal separator will be the decimal comma (,) whilst in the EPA SWMM
program the decimal point (.) is used. The EPA SWMM program is constantly upgraded and the
screen layouts on newer versions may differ from what is shown in this example.
PROJECT SETUP
To begin this tutorial, start the EPA SWMM program by double clicking the EPA SWMM icon on the
desktop. The main window should appear as shown in Figure 9.1.
In this ttutorial a draainage system serving a 0,48 ha residential areea will be m modelled. Th he system
layout iss shown beloow in Figuree 9.2 and coonsists of sub bcatchment areas
a S1 throough S3, storrm sewer
conduitss C1 through C4, and con nduit junctionns J1 through
h J4. The system discharrges to a stream at the
point labbelled Out1. The first steeps are to creeate the objects shown inn this diagramm on SWMM M's Study
Area Maap and setting the variouss properties of these objeects. The water quantity rresponse to a 75 mm,
6-hour raainfall eventt, as well as a continuous,, multi-year record will th
hen be simullated.
Go to File menu
1. G m on the main w and select New
m window N .
Figure 9.3: P
Project Defau
ults: ID labeels
T
This will make
m EPA SWMM
S bjects with cconsecutive numbers
autoomatically laabel new ob
ffollowing thee designated prefix.
On the Sub
4. O bcatchments page of thhe dialog sett the following default values (as shown
s in
F
Figure 9.4):
A
Area: 1,619
W
Width: 120
% Slope: 0,5
% Imperv: 50
N
N-Imperv: 0,01
N
N-Perv: 0,10
D
Dstore-Impeerv: 1,3
D
Dstore-Perv:: 1,3
%
%Zero-Impeerv: 25,0
6. Select Save as
a defaults for
fo all new prrojects and click
c OK to accept
a these choices and close the
ddialog.
Setting M
Map Option
ns
Next wee will set somme map displlay options sso that ID lab
bels and sym
mbols will bee displayed as
a we add
objects tto the study area
a map, and links will hhave directioon arrows.
5. Finally, seleect the Flow Arrows pagge; select thee Filled Arro
ow style, andd set the Arro
ow size
to 7 (see Figgure 9.11).
Fig
gure 9.11: M
Map Optionss: Flow Arro
ows
2. Y
You can leavve the dimen nsions at theeir default vaalues for thiss example. SSet the Map Units to
M
Meters by seelection this option
o and cclick on the OK
O button.
Finally, look in the status bar at a the bottomm of the main window and a check thhat the Auto o-Length
feature is Off. If it is on, then clicck the down arrow button
n and select "Auto-Lengtth: Off" from
m the pop-
up menuu that appearss.
We are nnow ready too begin addinng componeents to the Sttudy Area Map
M to createe the model layout as
shown inn Figure 9.12. We will sttart with the subcatchments.
3. DDo the samee for the nextt three corneers and then right-click
r th
he mouse (orr hit the Enteer key) to
cclose up the rectangle th
hat representts subcatchmment S1. You u can press th
the Esc key if instead
yyou wanted to cancel yo our partially drawn subcaatchment and d start over again. Don'tt worry if
tthe shape or position of the
t object isnn't quite righ
ht. We will gog back laterr and show how
h to fix
tthis.
4. R
Repeat this process
p ubcatchmentts S2 and S3. Note.
for su
1. T
To begin addding junction
ns, click the o the Objecct Toolbar.
button on
Now wee will add thhe storm sew wer conduits that connecct our drainaage system nnodes to onee another.
(You muust have creaated a link's end
e nodes ass described inn the previouus topic befoore you can create
c the
We will begin with conduitt C1 which cconnects juncction J1 to J2
link.) W 2.
1. C
Click the t Object T
button on the Toolbar. Thee mouse curssor changes sshape to a crosshair.
2. C
Click the moouse on juncttion J1. Notee how the mo
ouse cursor changes
c shappe to a pencil.
1. C
Click the Raain Gage buttton on thhe Object Toolbar.
Note in F Figure 9.15 that the subcatchments aare not yet liinked to a sp
pecific junctiion. This will be done
later in tthe exercise.
Re-Posittioning Objects
1. C
Clicking the button to map in Objeect Selection
t place the m n mode;
2. cclicking on the
t object to be moved;
3. ddragging thee object with the left mouuse button held down to itts new positiion.
2. T
Then click thhe button on the Maap Toolbar to put the map into Verteex Selection mode.
3. SSelect a verttex point on n the subcatcchment outliine by clicking on it (nnote how thee selected
vvertex is indiicated by a filled
f solid sqquare).
4. D
Drag the verrtex to its new
w position w
with the left mouse
m button
n held down.
6. W
When finisheed, click the button tto return to Object
O Selecction mode.
4. In the Save
S As dialog that apppears (Figure 9.18), seleect a folder and file nam me under
which too save this project. We ssuggest namiing the file [Link]
E . (An eextension off .inp will
be addedd to the file name
n if one iis not supplieed.)
5. Click Sa
ave to save th
he project to file.
Setting P
Properties
D
Double-clickk the object on
o the map.
Select the obbject from the Data pagge of the Bro ow and thenn click the Browser's
owser windo B
button.
Whenever the Prop perty Editorr has the foccus you can press the F1
F key to obbtain a moree detailed
descriptiion of the prooperties listeed.
Setting S
Subcatchmeent Propertiies
3. SSelect Subcaatchments asa the class oof object to edit, Rain Gage
G as the property to edit, and
ttype in Gagee1 as the new
w value (as shhown in Figu
ure 9.20).
4. CClick OK too change the rain gage off all subcatchments to Gage1
G . A connfirmation diialog will
aappear notinng that 3 sub
bcatchments have chang ged (Figure 9.21). Selecct No when asked to
ccontinue editting.
To set tthe outlet noode of our subcatchmentts we have to proceed one by one,, since thesee vary by
subcatchhment:
1. D hment S1 orr select it from the Data Browser annd click the Browser's
Double clickk on subcatch B
he Property Editor.
button too bring up th
2. TType J1 in the
t Outlet field
f and preess Enter. Note
N how a dotted
d line iss drawn between the
ssubcatchmennt and the node (Figure 99.22).
Setting N
Node/Link Properties
P
Only one of the condduits in our example sysstem has a no on-default prroperty valuee. This is con
nduit C4,
the outleet pipe, whosse diameter should
s be 0,445 m. instead
d of 0,3 m. To
T change itss diameter (m maximum
depth), sselect conduiit C4 into thee Property E
Editor and seet the Max. Depth
D value to 0,45.
Setting R
Rain Gage Properties
P
In order to provide a source of rainfall inpuut to our prooject we need to set the rain gage properties.
Select G
Gage1 into thhe Property Editor
E wing propertiies (see Figuure 9.24):
and sset the follow
As menttioned earlierr, a simulatioon of the respponse of the study area to o a 75 mm, 66-hour design
n storm is
requiredd. A time seriies named TS1 will conttain the hourly rainfall in ntensities that
at make up thhis storm.
Thus a tiime series obbject is createed and popullate with datta. To do thiss:
From the Da
1. F ata Browser select the Tiime Series category
c of objects.
2. C
Click the he Browser w
button on th which will bring
b up a Time Series E
Editor form.
Click the OK
6. C K button to accept
a the new
w time series.
RUN
NNING A SIINGLE EVE
ENT ANALYSIS
From the Da
1. F ata Browser, select the O
Options categ
gory and clicck the buutton.
2. OOn the Geneeral page off the Simulaation Optionns dialog that appears, seelect Kinema
atic Wave
aas the flow routing
r method. The Infiiltration model should already
a be sett to Green-A
Ampt. The
AAllow Ponding option sh hould be uncchecked (see Figure 9.27
7).
On the Dates page of thee dialog, set the End Analysis time to [Link] ( Figure 9.28
3. O 8) and the
Dates and Tiimes as indiccated on Figu
D ure 9.28.
Figure 9.28: S
Simulation options
o – Da
ates
mulation. Too do so, selecct Project | Run Simulaation (or clicck the
We are nnow ready too run the sim
button). After the runn was completed the Run
n Status willl be indicatedd as shown inn Figure 9.3
30.
o V
Viewing the Status Repo ort
o V
Viewing resuults on the map
m
o V
Viewing a time series ploot
o V
Viewing a prrofile plot
The Staatus Report contains useeful summarry information about thee results of a simulation
n run. To
view thee report, selecct Report | Status
S wn in Figure 9.31)
(show
2. OOf the 75 mmm of rain th hat fell on tthe study areea, 43,75 mm
m infiltratedd into the gro
ound and
eessentially thhe remainderr became runnoff.
The Conduiit Surchargee Summary table shows that Conduit C2, just dow
4. T ownstream off node J2,
w
was at full caapacity and therefore
t apppears to be sllightly underrsized.
2. SSelect the vaariables to view for Subccatchments, Nodes, and Links from the dropdow wn combo
bboxes in the Themes pan nel. Try for iinstance Run
noff, Depth and Flow resppectively as shown in
FFigure 9.32.
3. T
The colour coding
c used for a particuular variable is displayedd with a legeend on the sttudy area
map. To togggle the displaay of a legennd, select Vieew | Legendss (see Figuree 9.33).
m
Figu
ure 9.33: Leg
gends
4. T
To move a leegend to another locationn, drag it with
h the left mo
ouse button hheld down.
1. Select Report
R | Grap
ph | Time Seeries or simp
ply click on the Stanndard Toolb
bar.
2. A Time Series Plot dialog will appear. It is used to seleect the objectts and variab
bles to be
plotted.
For this example, thee Time Seriees Plot dialogg can be used
d to graph th
he flows in coonduits C1 and
a C2 as
follows ((shown in Fiigure 9.35):
1. Select Links
L as the Object
O Cate gory.
3. Click onn conduit C1 (either on thhe map or in n the Data Browser) andd then click in the
dialog too add it to thee list of linkss plotted. Do
o the same for conduit C22.
4. Press OK
K to create th
he plot as shhown in Figu
ure 9.36.
SWMM can generatee profile plots showing hhow water su urface depth varies acrosss a path of connected
c
nodes annd links. Creeate such a plot necting junction J1 to thhe outfall Ou
p for the cconduits conn ut1 of our
examplee drainage syystem. To do this:
Either enter J1
2. E J in the Sta
art Node fielld of the Pro
ofile Plot dialog that appeears or selectt it on the
m
map or fromm the Data Brrowser and cclick the button next to
t the field.
4. CClick the Finnd Path buttton. An ordeered list of thhe links whicch form a connnected path
h between
tthe specifiedd Start and End nodes w will be displlayed in the Links in Prrofile box seee Figure
99.38. You caan edit the en
ntries in this bbox if need be.
b
5. CClick the OK K button to create the pllot (Figure 9.39), showiing the waterr surface pro
ofile as it
eexists at the simulation tiime currentlyy selected in the Map Brrowser.
Figure
F 9.39: Profile – No
ode J1 – Out1
In the annalysis just run we chosse to use thee Kinematic Wave meth hod of routinng flows thrrough our
drainagee system. Thiis is an efficient but simpplified appro
oach that can
nnot deal with
th such phenoomena as
backwatter effects, prressurized flow, flow revversal, and non-dendritic
n c layouts. SW
WMM also in ncludes a
Dynamicc Wave rouuting proced dure that cann represent these conditions. This procedure, however,
requires more compuutation time, due to the nneed for smalller time stepps to maintainn numerical stability.
From the Da
1. F ata Browser, select the O
Options categ
gory and clicck the buttton.
2. OOn the Geneeral page of the Simulattion Optionss dialog that appears, seleect Dynamicc Wave as
tthe flow routting method (see Figure 9.41).
If you loook at the Status Reportt for this runn you will seee that there is
i no longer aany flooding
g and that
the peakk flow carriedd by conduitt C2 has beeen increased from 95,17 l/s to 112,566 l/s (the con nduit now
flows pressurized).
1. W
What is the maximum
m deepth at junctiion J2 and when
w does this occur?
2. W
What is the determined
d system
s runofff coefficientt for this exerrcise? Hint: ssee Status Reeport
3. W
What is the maximum
m flo
ow at the sysstem outlet?
4. W
What is the peak
p runoff from
f catchm
ment S2 and when
w does it occur?
5. W
What is the head
h nce between junctions J1 and J2 at 02
differen 2:45?
7. W
What capaciity is still avaailable in linkk C1?
8. D
Draw the floow versus tim
me graph for link C1.
Problem
m description
n Example 10.1
1
Outflow
w stage relatiionship of th
he spillway oof the dam is
i given by:
Q C d LH 1,55
w
where:
Q = discharge (m3/s)
Cd = discharge coefficient
L = length of the spillway ((m)
H = total energ
gy head (meaasured abovee the spillway
y level) (m)
In this ccase the outfflow can be determined by the follo onship: Q 110H 1,5
owing relatio
Figure 10
0.1: Section through thee spillway of the dam
i Figure 100.2.
The infloow hydrograaph is given in
Solution
n Example 10.1
1
H H
S AddH 10 6 (7,,5 1,5H)dH
O O
S 10 6 (7,5 H 0,755 H 2 k)
S O
N :
Δt 2
I1 I 2
N 2 N1 O1 and by su ubstituting thee known valu ues, it follow
ws:
2
A distannce away froom the spill section
s wherre the velociity approachees zero in thhe dam the difference
d
between the water level and the spillway
s leveel reflects thee total energy
y, i.e. h = H
100 6
N (7,5 H 00,75 H 2 ) 55
5 H1,5
7 200
N 10441,7H 104 ,17 H 2 55 H 1, 5
N 1044,17H(10 H 0,53 H )
F
Figure 10.3: Graphical p
presentation
n of the auxiliary functioon
If the infflow and outtflow hydrographs are plootted (Figuree 10.4) it willl be observeed that:
F
Figure 10.4: The inflow and calcula
ated outflow hydrographhs
Problem
m description
n Example 10.2
1
L
Longitudina
al profile forr longest wa
ater course in
i catchmennt
0,385
0,87L2
ΤC
1 0000 S av …(10.1)
…
W
Where: TC = timee of concentrration (hours))
L = hydrraulic length
h of the catchment, measuured along th he
flow
w path from the
t catchmen nt boundary tto the point
wheere the flood needs to be determined
d ((km)
Sav = averrage slope (m
m/m)
The userr may calculaate the averaage slope as ddefined in Fiigure 10.6.
H 0,85LL H 0,10L
S av
1 00000,75L …(10.2)
…
W
Where: Sav = averrage slope off the catchmeent (m/m) (seee Figure 100.6)
H0,10L = elevvation height at 10% of thhe length of tthe watercou
urse (m)
H0,85L = elevvation height at 85% of the
t length off the waterco ourse (m)
L = lenggth of the wattercourse (km
m)
16 450000
14 400000
350000
12
300000
10
Volume (m3)
Area vs. Height Relationship
Height (m)
4
100000
2 50000
0 0
0 20000 40000 60000 80000 100000 120000 140000
Area (m2)
Figure 10.7: Area-Height and Area-Volume relationship for the upstream side of the culvert
The discharge capacity of culverts operating as inlet controlled systems and varying upstream water
levels, are reflected in Table 10.5. The outlet characteristics can be defined for a rectangular culvert
by the equations as given below in column 2 for the different ratios of H1/D
For : For:
0 < H1/D < 0,8 0 < H1/D ≤ 1,2
0,05 1,9
Q S H1 2 2
0,48 0 D Q C B BH1 gH1
D 2
gD 0,4 3 3
And for: 0,8 < H1/D ≤ 1,2 * And for: H1/D > 1,2
Q C h BD 2gH 1 C h D
0,05 1,5
Q S H1
0,44 0 D
0,4
2
D gD
Where: Ch = 0,8 for rounded inlets
Ch = 0,6 for square inlets
(S0 = slope of culvert bed with slight
effect on capacity)
Note:
The estimated peak flow rate for the catchment was anticipated to be 76 m3/s for a return period of
Q2T. By assuming a triangular distribution of inflow, the inflow hydrograph is given in Figure 10.8.
70 Inflow (m³/s)
60
50
Flow rate (m3/s)
40
30
20
10
0
0 50 100 150 200 250
Time (min)
Step 3: Enter the site specific catchment and culvert parameters as given in the problem statement
above in the allocate blue cells in the Routing Utility:
The auxiliary function used to estimate the outflow from the culvert can be viewed under the
Auxiliary Function tab (Figure 10.9). Notice the vertical increase for outflows exceeding
63 m3/s. This was due to the Area-Height relationship that was not completed for all 16
available increments. However, this does not have any influence on the outcome of the results
since the maximum damming depth is only 4,64 m and that corresponds to an outflow of
45,18 m3/s which is lower than the allowable damming depth of 5 m to the SBP.
100
80
Outflow(m3/s)
60
40 Auxiliary Function
20
0
0.0 2000.0 4000.0 6000.0 8000.0 10000.0 12000.0 14000.0 16000.0
Auxiliary Function (m3/s)
Figure 10.9: Graphical presentation of the auxiliary function
Depicted on the inflow (blue line) and outflow (red line) hydrographs (Figure 10.10) are also
the change in volume with time (dV/dt) as well as the duration for which the upstream energy
head (H1) exceeds a ratio of 1,2D. These values are reflected by the dashed purple line and
green solid line respectively.
If the inflow and outflow hydrographs are plotted (Figure 10.10) it will be observed that:
the intersect of the hydrographs coincides with the maximum storage; and
the maximum outflow rate will be associated with the time of the intersect.
50
Flow rate (m3/s)
0.20
-0.20
20
10 -0.40
0 -0.60
0 50 100 150 200 250 300 350 400
Time (min)
Step 6: The results are summarized under the Summary of Results tab and can be printed by clicking
on the Print Summary Sheet button in the upper right corner of the screen (Figure 10.11).
Included on the summary sheet is a table that depicts the outcome of the assessment for
existing culvert structures. In this case the culvert has not met the criteria as described in
Chapter 10 and thus Fails the hydraulic criteria for existing culverts.
HEC-RAAS is constantly being upgraded and the scrreen layoutss on newer versions ma
ay differ
from wh
hat is shown
n in this chap
pter.
STARTIN
NG A NEW PROJECT
P
To beginn this exerciise, start thee HEC-RAS program by y double cliccking the HE
EC-RAS ico
on on the
desktop. The main window
w shouuld appear ass shown in Fiigure 11.1.
Figure
F 11.1: HEC-RAS main windo
ow
The firstt step in deveeloping a HEEC-RAS appplication is to w project. Goo to File men
o start a new nu on the
main winndow and seelect New Project. Thee New Projecct window sh hould appearr as shown in i Figure
11.2. Set the drive and
a directory y you wouldd like to work k in. Enter the project ttitle and filee name as
shown inn Figure 11.2. Once yo ou have enterred the inforrmation, presss the OK bbutton to accepted the
title and file name annd create the new projectt.
Once baack at the HEEC-RAS Maain window select from the menu baar Options, aand set the units u that
you wouuld like to work in to be metric unitss as well as be
b the defaullt setting for all new projjects (see
Figure 111.3 and Figu
ure 11.4).
ENTERING
G GEOMET
TRIC DATA
A
Figure
F 11.5: Geometric Data window
Figure 11.6: P
Plan view of river sectio
on
HHEC-RAS will
w indicate that the pictuure extents past
p the curreent schematiic boundaries and ask
yyou if you would
w o increase thee schematic extents (Figu
want to ure 11.8). SSelect Yes.
PPress the lefft mouse bu utton once too start drawiing the reachh. Move thhe mouse pointer and
ccontinue to press
p the lefft mouse butt
tton to add additional
a points to the liine segment.. To end
tthe drawing of the reach,, double clickk the left mo
ouse button and
a the last ppoint on the reach
r will
bbe placed at the current mouse pointter location (right
( click will
w remove tthe last poin nt drawn).
AAll reaches must be draw wn from thee upstream too downstream m (in the poositive flow direction)
d
ii.e. start at cross
c sectionn 12 down too cross sectiion 1 (Decreaasing numeriic values).
OOnce the reaach is drawn face will prompt you to enter an ideentifier for the River
n, the interfa
nname and thhe Reach nam me. The Rivver identifieer can be up to 32 characcters, while the reach
w be called,, Tsitsa and the reach
nname is limiited to 12 characters. In this exercisee the river will
LLower reachh (see Figuree 10).
Figure
F 11.100: River and
d reach namees
GGo to the Options menu u and select Add a new Cross Sectiion. An inpput box will appear to
pprompt you to
t enter a riv he new crosss section (seee Figure 11.14).
ver station ideentifier for th
Figure
F 11.14 : Add a new
w river statio
on
E
Enter the:
Description: Upstream boundary
D b of tthis river secction
D
Downstream m reach lengthhs: LOB = 1 00, Channel = 95 and RO OB = 90
Manning n-vvalues: LOB = 0,05, Chan
M annel = 0,035 5 and ROB = 0,05
Station and elevation
e dettails:
Nr S
Station Elevatioon
1 100 100,3
2 103 97,3
3 125 96,1
4 133 94,4
5 135 94,3
6 138 95,6
7 148 96,0
8 154 99,5
M L bank = 1 25 and Righ
Main channeel stations: Left ht bank = 138
8
C
Cont\Exp cooefficients: Contraction
C = 0,1 and Exp pansion = 0,3
3
P g the Plot
Plot the crosss section to visually insspect the dataa. This is acccomplishedd by pressing
C
Cross Sectioon option un nder the Plott menu on the Cross Secttion Data Edditor or by cliicking on
tthe Expand XS editor to t include a XS plot . The cro
oss section sshould look similar
s to
tthat shown inn Figure 11..16 below.
NNote the maanning rough hness valuess at the top of the crosss section annd the speciffied bank
sstations (red dots) (Figurre 11.16).
Two op
ptions of coompleting cross
c sectioon data entrry:
Method
d 2: # Expeert users: Skip
S the nexxt two pagess.
F
Figure 11.199: Cross secttion 9 data Figure 11
1.20: Cross ssection 8 datta
F
Figure 11.211: Cross secttion 7 data Figure 11
1.22: Cross ssection 6 datta
F
Figure 11.255: Cross secttion 3 data Figure 11
1.26: Cross ssection 2 datta
F
Figure 11.277: Cross secttion 1 data
M
METHOD 2
Press the OK
P K button and
d the new crooss section will
w appear in the editor.
C
Change the cross
c n description to “Tsitsa River 11”.
section
A
Adjust the crross section stationing
s too reduce the overbanks
o by
y 10%.
TThis is accommplished by selecting thhe Adjust Stations featurre from the O Options men nu on the
C or, then selecct Multiply by
Cross Sectioon Data Edito b a Factor. When the input box ap ppears for
tthis option, three data entry fieldss will be av vailable to adjust
a the sstationing off the left
ooverbank, chhannel, and the right overrbank separaately. Enter values
v of 0,99 for the righ
ht and left
ooverbanks, but
b leave thee main channnel field blaank. This willw reduce thhe stationing g of both
ooverbanks byy 10%, but leeave the maiin channel unnchanged.
PPress the Appply Data bu hat shown inn Figure 11.17). Plot
utton (The daata should bee similar to th
tthe cross secction to visuaally inspect itt (see Figuree 11.29).
These seeven steps abbove should be repeated to enter all the data for Tsitsa Riverr (Lower Reaach). The
d in Table 111.1. Perform
necessarry adjustmennts are listed m the cross section
s dupliications in order
o that
they are listed in thee table. Mak
ke sure to chhange the desscription of each
e cross seection, and also
a press
the Applly Data button after mak king the adjuustments for each
e cross seection.
This com mpletes all thhe cross secttion data forr the Tsitsa River
R (Loweer reach) savve the data fiile before
continuinng. Saving thet data to a file is achieeved by exiting the Crosss Section Dat ata editor win
ndow and
selectingg the Save Geometry
G Data
D As optiion from thee File menu on the Geom metric Data window.
After sellecting this option
o you will
w be prom mpted to enteer a Title forr the geomettric data (Figgure 30).
Enter “N Natural Tsittsa River” forf this exeercise, and then t t OK butttton. A file name is
press the
automatiically assignned to the geeometry dataa based on whatw you enteered for the project file name i.e.
Exercisee1.g01.
Figure
F 11.30:: Save Geom
metry Data As
A
RReturn to thee Cross Secttion Data ediitor window and from the Plot menuu select Plot Profile...
tto view the longitudinal
l profile
p of thee entered cro ( Figure 111.31).
oss sections (see
We havee completedd the required geometricc data and caan now conttinue and ennter the Steaady Flow
Data.
ENTERING
E STEADY FLOW
F DAT
TA
E
Enter the 1:220, 1:50 and 1:100 year fflood peak values of 35, 80 and 150 m
m³/s respectiively (see
F
Figure 34).
TThe next steep is to enteer any requirred boundaryy conditions. To enter bboundary coonditions,
ppress the Reeach Bounda ary Conditiions button at
a the top of the Steady FFlow Data ed
ditor (see
FFigure 11.344). The boun w appear as shown in Fiigure 11.35.
ndary conditiions editor will
A pop up boox will appeaar (see Figu re 11.36) reequesting you u to enter ann average slo
ope at the
ddownstream end of the riiver reach. E m/m) then clicck on the OK
Enter a valuee of 0,001 (m K button.
T
This completes all of thee necessary bboundary con ( Figure 111.37).
ndition data (see
Click the OK
C K button on the Boundarry Condition
ns window to
o accept the ddata and retu
urn to the
Steady Flow
w Data screen
n.
O
Once the datta has been saved, you caan close the Steady
S Flow Data Editor..
PERFORM
MING THE HYDRAUL
LIC CALCU
ULATIONS
Now thaat all of the data has beeen entered, we can callculate the steady
s waterr surface pro
ofiles. To
perform the simulations, go to th
he HEC-RAS S main window and selecct Steady Floow Analysis from the
Run meenu or click on the Stea ady Flow A Analysis buttton t menu baar. The Stea
on the ady Flow
Analysiss window shhould appear as shown in Figure 11.3
39, except yo
ours will not have any plaan title or
Short ID
D yet.
You will be prompted to enter a shhort identifierr. Enter a tittle of Run1 in the Short ID box
Y
(Figure 11.4
41) and click
k on the OK bbutton.
TThe next sttep is to seelect the deesired flow regime for which the model will perform
ming Subcritical flow caalculations only
ccalculations.. For this exaample we wiill be perform o since
oonly a downnstream boun ndary condittion was speecified. Makee sure Subcrritical is thee selected
fflow regime..
Now thaat everythingg has been saved and sset, the steady flow commputations ccan be perfoormed by
pressing the Compuute button att the bottom of the Stead
dy Flow Simmulation winndow (Figurre 11.39).
Once thee computatioons have beeen completedd, the compu
utation wind
dow can be cclosed, as well as the
Steady F
Flow Simulattion windoww.
VIEW
WING RESU
ULTS
Profile plotts
Rating curvves
Profile p
plot
ot as shown in Figure 11
This shoould give youu a profile plo 1.44.
The thirdd plot optionn, which is of value, is thhe General Prrofile Plot. Select
S Generral Profile Plot
P from
the View
w menu. Thiss should givee you a profille plot as sho own in Figurre 11.45.
From thee Standard Plots menu various otheer useful plo b made succh as Froude numbers
ots can also be
(on left bbank, main channel
c and right
r bank), ssee Figure 11.46.
Figure 11
1.46: Generral profile pllot (Froude numbers)
n
Figu
ure 11.48: Geeneral profiile plot (Wriiting data to
o a file (csv ttype))
This cop
pying of picttures/graphs/data can b
be done thro
oughout HEC-RAS.
Rating ccurve
Select R
Rating Curvees from the View
V menu pplot a compu uted rating cu
urve. A ratinng curve based on the
computeed water surfface profiles will appear as shown in Figure 11.4 49. You can look at the computed
c
rating cuurve for any location
l (cro
oss section) bby selecting the
t appropriaate river, reac
ach and river station.
X-Y-Z P
Perspective Plot
P
Select X
X-Y-Z Persp pective Plotss from the V
View menu to o plot a 3D view
v n (Figure
of the river section
11.50). S
Set the Rotattion Angle to –70 and thhe Asimuth Angle
A to 19 to obtain thee same view as shown
in Figurre 11.50. Thhis type of view
v gives a clear view of
o the wideniing of the rivver at cross section 5
and 6.
Tabularr output
Detailed
d tabular ou
utput
Now loook at some taabular outpu ut. Go to thee View menu u on the HEC C-RAS mainn window. There
T are
two typees of tables available, a detailed outtput table an nd a profile summary tab able. Select Detailed
Output Tables to geet the first tab i Figure
ble to appearr. The table should be siimilar to the one shown in
11.51 (nnotice the warning at the bottom off the table fo or river statiion number 3). This tab
ble shows
detailed hydraulic innformation at a single cross sectio on. Other cross sectionns can be viiewed by
selectingg the approprriate reach an
nd river from
m the table.
The HE EC-RAS softtware has a system of Errors, Warrnings, and Notes that are passed from the
computaation program ms to the useer interface. During the computation
c s, the compuutation progrrams will
set flagss for at a paarticular node (nodes aree cross sections, bridgess, culverts, oor multiple openings)
o
wheneveer it is necesssary. These message flaags are writteen to the staandard outpuut file, along with the
computeed results forr that node. When the uuser interfacee reads the computed
c ressults from thhe output
file, if anny errors, waarnings, or no
otes exist, thhey are interp
preted and diisplayed in vvarious locatiions from
the interface.
The useer can requeest a summaary of all thhe errors, waarnings, and
d notes that occurred du
uring the
i accomplished by seleccting Summ
computaations. This is mary Errors,, Warnings,, and Notes from the
View meenu on the main
m HEC-RA
AS window oor clicking th
he short cut button
b .
ERROR RS: Error meessages are only sent wheen there are problems
p thaat prevent thee program fro
om being
able to ccomplete the run.
Cross seections spaceed to far apart. This cann cause severral warning messages
m to bbe set.
Cross seections startting and end ding station ns not high enough
e . If a computed w
water surface is higher
than eithher end pointt of the cross section, a w
warning messsage will appear.
Bad Staarting Waterr Surface Ellevation. If tthe user speccifies a bound dary conditioon that is nott possible
for the sppecified flow
w regime, thee program wiill take action and set an appropriate warning messsage.
Bad Crooss Section Data. This can cause seeveral problems, but mo ost often the program wiill not be
able to bbalance the ennergy equation and will ddefault to criitical depth.
Profile S
Summary Table
T
There arre several typpes of profile tables listeed under the Std. Tabless menu (see Figure 11.5 53) of the
profile taable window
w. Each one of o these tablees shows typpical detail reelevant to thee specific structure.
User defined tables are limited to 15 variables. Once you have selected all of the variables that you
want, press the OK button at the bottom of the window. The profile table will automatically be
updated to display the new table
Once you have the table displayed in the profile table window, you can save the table headings for
future use. To save a table heading, select Save Table from the Options menu on the profile table
window. When this option is selected, a pop up window will appear, prompting you to enter a name
for the table. Once you enter the name, press the OK button at the bottom of the pop up window. The
table name will then be added to a list of tables included under the User Tables menu on the profile
table window.
Create a user defined table as detailed above with the following columns:
Q Total (Total flow rate)
W.S. Elev (Water surface elevation)
Crit W.S. (Critical water surface elevation)
Vel Chnl (Velocity in the main channel)
Froude # Chl (Froude number in the main channel)
Save the table (Exercise 1 – Table) and view the summary table thereof with 1:50 and 1:100 year
profiles (selected from the Options menu) (Figure 11.55).
At the ennd of this exeercise the folllowing objeectives shoulld have been met:
B
Be able to seet-up a HEC--RAS projectt
K
Know how to
t enter geom
metric data
U
Understand the
t setting off boundaries and controlss
K
Know how to
t enter stead
dy flow data
K
Know how to
t analyse a river
r system
K
Know how to
t extract info
formation
3. What is the kinetic energy component at cross section 11 for the 1:50 year flood?
4. What is the kinetic energy component at cross section 5 (the proposed site) for the 1:50 year
flood?
5. Why is there a difference between the kinetic energy component at cross section 5 and 11?
6. What are the flood levels (1:20, 1:50 and 1:100) at the proposed site (cross section 5)?
7. What will the water level be at cross section 5 if the flood peak is 100 m³/s?
8. What is the energy weighting coefficient (alpha) for cross section 5 for the 1:50 year flood?
9. How wide is the river flowing at cross section 5 during the 1:100 year flood?
10. What is the hydraulic depth in the main channel during the 1:100 year flood at cross section 5?
11. How would you know if cross section 5 or 6 was functioning as a control point in the river
section?
12. Are you in a position to draw in the 1:20, 1:50 and 1:100 year flood lines for the proposed
site?
STARTIN
NG A NEW PROJECT
P
To beginn this exerciise, start thee HEC-RAS program by y double cliccking the HE
EC-RAS ico
on on the
desktop. The main window
w shouuld appear ass shown in Fiigure 11.56.
The firstt step in deveeloping a HE EC-RAS appplication is to w project. Goo to File men
o start a new nu on the
main winndow and seelect New Project. Thee New Projecct window sh i Figure
hould appearr as shown in
11.57. Set the drive and directorry you wouldd like to worrk in. Enter the project title and filee name as
typicallyy shown in Figure
F 11.57
7. Once yoou have entered the information, preess the OK button to
acceptedd the title andd file name and
a create thee new projecct.
Once baack at the HEEC-RAS Maain window select from the menu baar Options, aand set the units u that
you wouuld like to work
w in to be
b metric uniits as well asa be the default settingg for all new
w projects
(assignm
ments). In thee right hand corner of thee main screen it will now
w indicate SI units.
PPress the lefft mouse bu utton once too start drawiing the reach h. Move thhe mouse pointer and
ccontinue to press
p the lefft mouse butt
tton to add additional
a points to the liine segment.. To end
tthe drawing of the reach,, double clickk the left mo
ouse button and
a the last ppoint on the reach
r will
bbe placed at the current mouse
m pointeer location (rright click wiill remove thhe last point drawn).
d
AAll reaches must be draw wn from thee upstream too downstream m (in the poositive flow direction)
d
ii.e. start at cross
c section 100 down too cross sectiion 70, 65 do
own to crosss section 40 and from
ccross sectionn 8 to 3 (see Figure
F 11.600).
OOnce a reachh is drawn, the interface will prompt you to enterr an identifieer for the Rivver name
aand the Reacch name. Th he River ideentifier can be up to 32 chharacters, whhile the reach
h name is
vers will be called, Riet and Blesbokk and the
llimited to 122 characters.. In this exeercise the riv
rreaches Uppper and Loweer reach (see Figure 11.6 60).
W
When you first
fi draw thee schematic tthere will bee no tic mark ks representinng the cross sections.
T
The tic markks only show
w up after youu have entereed cross sectiions.
Figure 11.60:
1 Geom
metric Data window (sch
hematic)
GGo to the Options menu u and select Add a new Cross Sectiion. An inpput box will appear to
pprompt you to
t enter a riv he new crosss section (seee Figure 11.62).
ver station ideentifier for th
Figure
F 11.62 : Add a new
w river statio
on
E
Enter the:
Description: Upstream River
D R Station
D
Downstream m reach lengthhs: LOB = 225, Channel = 28 and ROB = 30,5
Manning n-vvalues: LOB = 0,035, Chhannel = 0,02
M 25 and ROB = 0,035
Station and elevation
e dettails:
Nr Statioon Elevation
1 0,00 26,82
2 3,05 25,60
3 4,57 24,69
4 8,53 23,77
5 9,75 22,25
6 10,67 21,03
7 12,19 20,42
8 15,24 19,51
9 18,29 19,45
10 19,81 20,73
11 21,34 22,25
12 22,56 24,23
13 24,08 24,38
14 25,30 25,91
15 26,21 28,04
M L bank = 8,
Main channeel stations: Left 8,53 and Righ
ht bank = 22,,56
C
Cont\Exp cooefficients: Contraction
C = 0,1 and Exp
pansion = 0,3
3
Figure
F 11.64:: Save Geom
metry Data As
A
o Stepp 5: Return too HEC-RAS,, and select same s numberr (or more) nnumber of ro
ows in the
Cross Section Data
D editor wwindow (see Figure 11.6 67). From thhe Edit menuu click on
Pastte, to insert the data in thhe table.
o Stepp 6: Enter thee rest of the oother charactteristic data:
Desccription: “”
Dowwnstream reaach lengths: L LOB = 22,9, Channel = 29,0 2 and ROB B = 33,5
Mainn channel staations: Left bbank = 8,32 anda Right baank = 22,56
Cont\Exp coefficcients: Contrraction = 0,1 and Expansion = 0,3
All the
t other chaaracteristic daata is providded in Table 11.2.
The spreadsheet priintout of alll the cross section dataa is attachedd in Appendix A for
reference.
Upper
Riet
85 7,92 22,56 23,8 28,0 31,4
80 7,71 22,56 25,6 27,4 30,5
75 7,53 22,56 24,4 29,0 33,5
70 7,32 22,56 0,0 0,0 0,0
65 8,78 31,64 25,0 29,0 33,5
60 8,53 31,64 25,9 30,5 36,6
Lower
Lower
Figure 11.69: A
Automaticallly Backup Data
D
T
There is onlyy one junctio
on in the riveer system. For this junctiion add the ffollowing description:
D
Division of Upper Reacch and Loweer Reach at Confluence with Blesbook River (see Figure
11.70)
RReach lengthhs across thee junction arre entered in n the junctioon editor, ratther than in the cross
ssection data.. This allows for the lenggths across very
v compliccated conflueences (i.e. floow splits)
tto be accomm modated. Th his is the reaason why the reach length hs in the crosss section daata for the
llast cross secction of each
h reach was lleft blank or set to zero. Enter the junnction length hs of 24,4
m and 21,3 m for the Rieet River – U Upper reach and
a Blesbok River – Low wer reach resspectively
((see Figure 11.70).
O
Once you haave entered all n the Apply button and close the
a the data ffor the juncttion, click on
w
window by pressing
p the OK
O button.
R
Remember to king on Savee Geometryy Data from
t save the Geometry ddata by click m the File
m
menu.
Fig
gure 11.71: U
Unsteady Flow Data Ed
ditor
Bounndary condiitions are en ntered by ffirst selecting the Boun ndary Condditions tab from the
Unstteady Flow Data editorr. River, Reaach, and Riveer Station locations of thhe external bounds b of
the ssystem will automaticallly be shownn in the table. Boundarry conditionss are entered d by first
seleccting a cell inn the table fo
or a particulaar location, th
hen selecting
g the boundarary condition type that
is deesired at thatt location. Not
N all bounddary conditio on types are available foor use at all locations.
l
The program will automaticaally gray-outt the boundaary condition n types that aare not relev
vant when
the uuser highlighhts a particulaar location inn the table.
Bou
undary Cond
ditions
TThere are several differeent types of bboundary conditions available to the user. The following
f
iis a short discussion of eaach type:
FFlow Hydroograph (this is the type tthat you can n enter for this exercise))
A flow hydrrograph can be used ass either an upstreamu bouundary or do downstream boundary
b
ccondition, buut is most commonly ussed as an up pstream boun ndary condittion. When the flow
hhydrograph button
b is preessed, the winndow shown n in Figure 11.72
1 will apppear. As sh
hown, the
uuser can eithher read the data from a HEC-DSS (HEC ( Data Storage
S Systeem) file, or can enter
tthe hydrograaph ordinatess into a tablee (Since Riv ver station 8 was selectedd as shown in i Figure
111.71 the floow hydrograp ph detail for this stations is first requiired).
Figu
ure 11.72: Flow
F Hydroggraph (River station 8 – Blesbok Riiver)
T
The user also has the opption of enterring a flow hydrograph directly intoo a table, as shown in
F
Figure 11.772. The firsst step is to select a Da ata time interval from m the drop down
d list.
C
Currently thee program on
nly supports regular interrval time seriies data.
T
To enter dataa into the tab o select eitheer Use Simullation Time or Fixed
ble, the user iis required to
S
Start Time. For this exeercise select Use Simulattion Time.
Figu
ure 11.73: Flow
F Hydroggraph (River station 8 – Blesbok Riiver)
AAn option liisted at the bottom of thhe flow hyd drograph bou undary conddition is to make
m this
bboundary a Critical Boundary Con ndition. When you seleect this optiion, the prog gram will
mmonitor the inflow
i hydro
ograph to seee if a changee in flow ratee from one tim
ime step to th
he next is
eexceeded.
IIf the changge in flow rate does exceed the user entereed maximum m, the program will
aautomaticallyy cut the tim
me step in haalf until the change in flow rate doess not exceed
d the user
sspecified maax.
T
The flow hydrograph forr this exercisse can be plo king on the P
otted by click Plot Data bu
utton (see
F
Figure 11.744).
T
The maximuum peak flow
w value is 31,,0 m³/s (timee 2:30)
T
The completted Flow Hy ydrograph edditor screen and plot screeen is shownn in Figure 11.75
1 and
F
Figure 11.766 respectively. The maxiimum peak flow
f value is 84,0 m³/s (ttime 2:00).
O
Other bounddary types not used in thiss example:
SStage Hydroograph
A stage hyddrograph can n be used ass either an upstream
u or downstream boundary condition.
c
TThe editor foor a stage hyydrograph iss similar to the flow hydrograph edittor. The useer has the
cchoice of eitther attaching a HEC-DS SS file and path
p name orr entering thhe data directly into a
ttable.
Inittial Conditioons
IIn addition tot the bounddary conditioons, the user must estab blish the inittial conditionns of the
ssystem at thee beginning of
o the unsteaady flow sim mulation. Inittial conditionns consist off flow and
sstage inform
mation at eacch of the crooss sections,, as well as elevations ffor any storaage areas
ddefined in thhe system. In
nitial conditiions are estab m within the U
blished from Unsteady Fllow Data
eeditor by selecting the In
nitial Condittions tab. Affter the Initial Conditionns tab is seleected, the
UUnsteady Fllow editor will
w appear as shown in Fiigure 11.78.
AAs shown inn Figure 11..78, the user has two opttions for estaablishing thee initial cond
ditions of
tthe system. The
T first optiion is to enteer flow data for each reacch and have the programm perform
a steady flow
w backwater run to compuute the correesponding staages at each cross section n.
T
The second option
o is to read
r in a file of stages and d flows that were writtenn from a prev
vious run,
w
which is callled a “Restarrt file”. The ffirst option iss used in thiss exercise.
EEnter the iniitial flows ass indicated inn Figure 11.78, 0,5 m³/ss, 1,0 m³/s aand 1,5 m³/s for cross
ssections 8, 100 and 65 reespectively.
TThe last stepp in developin ng the unsteaady flow datta is to save the
t informatition to a file. To save
tthe data, seleect the Save Unsteady F Flow Data AsA from the FileF menu onn the Unstea ady Flow
DData editor. A pop-up window wiill appear prrompting you to enter a title for thee data as
sshown in Fig gure 11.79. Enter “Base Unsteady Flow F ” for thiss exercise, annd then press the OK
bbutton. A fiile name is automatically
a y assigned to
o the Unstead dy Flow Datta based on what
w you
eentered for thhe project fille name i.e. E
Exercise1.u001.
Fig
gure 11.79: S
Saving Unsteeady Flow Data
D
Oth
her unsteadyy Flow Data
a Options
OObserved Data
D In DSS
TThis option allows the user to attaach observed
d data pathn
names from a HEC-DS
SS file to
sspecific riverr stations witthin the moddel.
WWhen an observed data pathp name iss attached too a specific river station llocation, thee user can
gget a plot of the observed
d flow or stagge hydrograp ph on the samme plot as thhe computed flow and
sstage hydroggraphs. Add ditionally, thee observed data
d will show w up on proofile and crosss section
n, the user seelects Obserrved Data In
pplots. To usse this option n DSS from the Optionss menu of
tthe Unsteady y Flow Data a editor.
MMinimum Flow
F and Floow Ratio Taable
TThis option brings up a global e ditor that willw show all a the locattions in wh hich flow
hhydrographs have been attached
a as bboundary connditions. The editor allow
ws the user to
t enter a
mminimum floow or a flow factor for eaach flow hyd
drograph bou
undary condittion.
UNSTEA
ADY FLOW
W ANALYSIIS
P
Performing Unsteady Flow
F Calculaations
TThe first stepp is to put together a Plan n. The Plan n defines whiich geometryy and flow data
d are to
bbe used, as well
w as provid ding a title an
and short iden
ntifier for thee run.
Y
You will be prompted to o enter a shoort identifierr. Enter a tiitle of Base in the Short ID box
(Figure 11.8
82) and click on the OK bbutton.
o Post-Prrocessor
The posst-processor is used to ccompute detailed hydrau ulic informattion for a seet of user
specifieed time lines during the uunsteady flow
w simulation periods. In general, the unsteady
flow computations only
o computte stage and flow at all ofo the compuutation nodess, as well
as stagee and flow hyydrographs aat user speciffied locationss. If the postt-processor is not run,
then thee user will only
o be ablee to view thhe stage and flow hydroographs and no other
output from
f HEC-R RAS. By ruunning the post-processo
p or, the user will have alla of the
availablle plots and tables
t for unnsteady flow that HEC-R RAS normallyy produces for f steady
flow.
SSimulation time
t window w
TThe user is required
r to enter
e a time window thaat defines the start and eend of the siimulation
pperiod. The time window w requires a starting datee and time an
nd an ending date and tim
me.
CComputatioonal settings
T ngs area conntains the Co
The computaational Settin omputationa al Interval, H
Hydrograph h Output
IInterval, Deetailed Outpput Interval, the name an nd path of the output DSSS file, and whether
w or
nnot the proggram is run in
n a mixed fl
flow regime. The computtation intervval is probably one of
tthe most impportant param
meters entereed into the model.
m It shou
uld be small enough to accurately
a
ddescribe the rise and fall of the hyddrographs being routed but not smaall to take forever
f to
ccompute.
FFor this exerrcise set the Computatioonal Interva al at 5 minutees (from the drop down list). Set
tthe Hydrogrraph Outpu ut Interval tto 15 minutees and the Detailed
D Outtput Interva
al also at
115 minutes.
Figure 11.83:
1 Unsteeady Flow Analysis
A (com
mpleted)
Clicck on the Co
ompute butto
on to run the Unsteady Flow
F Analysiis.
Oncce the modeel has finisheed all of thee computatio ons successffully, you caan begin vieewing the
resuults. Severaal output opttions are avaailable from the View menum bar on the HEC-RRAS main
winndow. Thesee options incllude:
a) CCross sectionn plots
b) PProfile plots
c) GGeneral proffile plot
d) RRating curvees
e) XX-Y-Z persppective plots
f) DDetailed tabuular output at
a a specific ccross section (cross sectio
on table)
g) LLimited tabuular output att many cross sections (profile table)
BBegin by plootting a crosss section. Sellect Cross Sections fromm the View m menu bar on the
t HEC-
RRAS main window.
w Anyy cross sectioon can be plo otted by seleecting the apppropriate riv
ver, reach
aand river staation (See Figure 11.84) . Several plo otting featurees are availabble from the Options
mmenu bar onn the cross section ploot window. These T optionns include: zoom in; zo oom out;
sselecting whhich plans, profiles,
p varriables to pllot; and con
ntrol over linnes, labels, symbols,
sscaling etc.
SSet the Zeroo Delay horizontal scrolll bar as indiccated (in Fig
gure 11.85) aand click on
n the play
bbutton, , to view the changing w water surfacee levels and energy gradde lines at this
t cross
ssection. Clicck the Stop button,
b , too stop the an
nimation.
N
Next plot a water surfacce profile. S elect Waterr Surface Prrofiles from the View menu
m bar.
T g you a prrofile plot ass shown in Figure 11.88.
This should give
Figu
ure 11.88: Water
W Surfacce Profile (fo
or Riet Riveer: Upper reeach)
Figure
F 11.899: X-Y-Z Perrspective plo
ot
BRIDGE
ES
In the neext section a bridge struucture will bee added dowwnstream of the confluennce of the twwo rivers.
HEC-RA AS computess energy losses caused bby structuress such as briidges and cuulverts in thrree parts.
One partt consists off losses that occur
o in the reach immeediately down nstream fromm the structu
ure where
an expannsion of flow w takes placce. The seconnd part is th
he losses at the
t structuree itself, whicch can be
modeledd with severral different methods. T The third parrt consists of
o losses thaat occur in the
t reach
immediaately upstream m of the stru
ucture where the flow is contracting
c to
o get throughh the opening
g.
Cross seection 2 is located imm mediately dow wnstream frrom the brid dge (i.e. withhin a short distance).
d
This cross section shhould represeent the naturaal ground jusst outside thee bridge.
Cross seection 3 shoould be locatted just upstrream of the bridge. The distance bettween cross section 3
and the bbridge shoulld be relativeely short. Thhis distance should
s only reflect the leength requireed for the
abrupt accceleration and
a contraction of the flow w that occurrs in the imm
mediate area oof the openin ng.
Cross seection 4 is ana upstream cross sectionn where the flow lines arre approxim
mately paralleel and the
cross secction is fully effective.
T
To enter briidge data thhe user pressses the Brid
dge/Culvert button on tthe Geomettric data
window (Fig
w gure 11.60). Once the B ridge/culverrt button is pressed,
p the B
Bridge/Culv
vert Data
E
Editor will appear
a own in Figurre 11.90.
as sho
F
From the Op ptions menuu, select Addd a Bridge and/or Culv vert from thhe list. An input
i box
w
will appear prompting
p yo
ou to enter a river station
n identifier fo
or the new brridge.
Enter 52 as shown
E s in Fig
gure 11.91.
F
Figure 11.91: Bridge rivver station (R
Riet River: Lower
L reachh)
E
Enter the Deescription off the bridge: SStephnie brid
dge (see Figure 11.92)
F
From the Bridge/Culv
B ert Data E Editor selectt the Deck//Roadway iicon to actiivate the
D
Deck/Roadw
way Data Ed wn in Figure 11.93.
ditor as show
TThe bridge deck will be 0,9 m higgh and will have a slop
pe across it of 0,05 m (for road
ddrainage).
TThe US and DS Embank kment SS (uupstream and d downstreamm embankmeent side slop pe) values
sshould be enntered as 2 (horizontal to 1 verticaal). These values
v are ussed for the graphical
rrepresentatioon on the pro
ofile plot.
AAt the bottoom of the Deeck/Roadwaay Data Ediitor, there are a three addditional fieldds of data
fi is the Max
eentry. The first M Allowab ble Submerg gence. This input ratio oof downstreaam water
ddepth to upstream energy y, as measurred above thee minimum weir
w elevatioon. When th he ratio is
eexceeded, thhe program will
w no longerr consider th he bridge decck to act as a weir and will switch
tthe computattion mode to o energy (staandard step) method. For this exercisee the defaultt value of
00,95 (95%) should
s be sellected.
TThe second field at the bottom of tthe editor is the Min Weir W Flow E Elevation. Th his is the
eelevation thaat determiness when weir flow will staart to occur over
o the briddge. If this field is left
bblank (as in this exercisee), the progrram will deffault to use th
he lowest hiigh cord valu ue on the
uupstream sidde of the bridge. The lasst field at th
he bottom of the editor iss the selectio on of the
WWeir Crest Shape. Thiis selection w will determinne the reducction of the w weir flow co oefficient
ddue to submeergence. Forr this exercisse, a Broad Crested
C weir shape shouuld be selecteed.
Click the OK
C K button
TThis bridge has three piiers that shoould be entered. From the Bridge/C Culvert Datta Editor
sselect the Pieer icon to activate the Pi er Data Edittor as shownn in Figure 111.95. The thhree piers
aare entered by specifyiing the Cen nterline staation at thee upstream side as welll as the
ddownstream side. The fiirst pier is poositioned at 151 m. It has a width of 00,5 m and it starts
s at a
llevel below the ground profile
p and eends at a lev ined cords [Link]. 20 and
vel inside thee bridge defin
223 m.
C
Click on thee Add button d pier is at ceenterline 20 m, has a
n to add a piier (Pier #2).. The second
width of 0,5 m and startss at elevationn 19 m and en
w nds at elevatiion 23 m.
C
Click on thee Add button n to add a ppier (Pier #3). The third pier is at ceenterline 25 m, has a
width of 0,5 m and startss at elevationn 19 m and en
w nds at elevatiion 23 m.
Figure 11
1.96: Bridgee data (View
w of bridge with
w piers)
TThis will poosition the neew cross secction downsttream of the newly addedd bridge (wh hich is at
ccross sectionn 52). The reach
r lengthss should alsoo be corrected. This new on is 8 m
w cross sectio
ddownstream of the bridg ge, thus the distance froom River station 55 to thhe bridge is 5 m, the
wwidth of the bridge is 7,5 m and thiss cross sectioon is a furtheer 8 m downnstream of thhe bridge.
TThe total disstance from River
R stationn 55 to Riverr station 51 is
i 20,5 m. G Go to River station
s 55
a the reach lengths to 200,5 m. Click
aand change all k on the App ply button to accept the ch
hanges.
TThe average slope betweeen River stattions 55 and 50 is 0,0006656 m/m. Thhe newly add ded River
sstation 51 cuurrently has the same eleevations as River
R station 55 (since itt is a copy of RS 55).
TTo adjust thhe elevation of this Riverr station (RS
S 51) with a value of thhe distance multiplied
m
wwith the aveerage slope i.e. 20,5 m x 0,000656 = ± 13 mm click on Adjuust Elevatio ons under
tthe Options menu. Enteer an adjustm ment of -0,01 13 m (as showwn in Figurre 11.98) and d click on
tthe OK buttoon.
BBefore we continue
c T is done by clicking on Save
we need to savee the Geomeetric data. This
GGeometry Data
D As undeer the File m
menu on the Geometric
G Data
D window w. After selecting this
ooption you will
w be prom mpted to enteer a Title forr the geometric data (Fiigure 11.100
0). Enter
“Base Geom metry + bridg ge” for this exercise, an
nd then press the OK buutton. A filee name is
aautomatically assigned to
t the geomeetry data baased on whatt you entereed for the prroject file
nname i.e. Exxercise2.g02.
TThe next stepp is to enter the ineffectivve flow areas. Any inefffective flow aareas that ex
xist due to
tthe bridge shhould be enteered. At a bbridge ineffecctive flow arreas normallyy occur just upstream
aand downstrream of the road embannkment, awaay from the bridge openning. It wass for this
rreason that River
R station 51 was addeed just down nstream of thee bridge.
NNo ineffectivve flow areaas were speccified for thee newly creaated downstrream river sttation 51.
TThe ineffective flow areas could alsoo be set by clicking
c t Boundinng XS’s 55 button
on the b on
tthe Bridge Culvert
C Data a Editor screeen (see Figu
ure 11.96).
The entered bridge shoulld now look similar to that shown in Figure
T F 11.1002.
B
Bridge Mod deling Appro oach
T
The bridge routines
r allow
w the modeller to analysee the bridge flows by usiing differentt methods
wwith the sam
me geometry y. The diffeerent methodds are: low flow,
f high fl
flow and com mbination
fflow.
F
For this exerrcise select th
he energy, m
momentum an
nd Yarnell eq
quations.
TThe Momenntum Balancee method perrforms a mo omentum balance throughh the bridge area and
rrequires the selection off drag coeffi g coefficientt to 2,0 (Squ
ficient, Cd. Set the drag uare nose
ppiers).
YYarnell Class A flows exists wheen the waterr surface th hrough the bbridge is coompletely
ssubcritical (ii.e. above the critical deppth). The flow regime without
w the bbridge was su
ubcritical
aand thus Yarrnell should also be sele cted. Enter the Yarnell pier
p coefficieent K as 1,255 (Square
nnose and taill).
F ws Energy Only
For high flow O (Standaard Step) sh
hould be selected.
Go to the Hy
G ydraulic Tab
ble Parametters Editor by o the HTabb Param icon
b clicking on n .
B
Before we continue
c we need to savve the Geommetric data. This
T is done by clicking on Save
G
Geometry Data
D he File menuu on the Geometric Data
under th a window.
P
Performing Unsteady Flow
F Calculaations with the
t Bridge
TThe first stepp is to put to w Plan. Thee Plan defines which geoometry and flow
ogether a new f data
aare to be useed, as well as providing a title and sho
ort identifier for the run.
TTo establishh a plan, seleect New Plaan from the File menu on o the Unstteady Flow Analysis
wwindow. Ennter the plan title as Basee analysis + bridge and then
t he OK button
press the n (Figure
111.107).
Y
You will be prompted too enter a shoort identifier.. Enter a titlle of Bridge in the Shorrt ID box
(Figure 11.1
108) and click on the OKK button.
S
Select the coorrect Geomeetry file and U
Unsteady flo
ow file from the drop dow
wn list.
SSimulation time
t window w
TThe user is required
r to enter
e a time window thaat defines the start and eend of the siimulation
pperiod. The time window w requires a starting datee and time an
nd an ending date and tim
me.
CComputatioonal settings
TThe computaational Settinngs area conntains the Co
omputationa al Interval, HHydrograph h Output
IInterval, Deetailed Outpput Interval, the name an nd path of the output DSSS file, and whether
w or
nnot the progrram is run in
n a mixed fllow regime. The compu utation intervval is probably one of
tthe most im mportant paraameters enteered into the model. This T should be small en nough to
aaccurately describe
d the rise and falll of the hyd
drographs being routed bbut not smalll to take
fforever to coompute.
Figure 11.109:
1 Unstteady Flow Analysis
A (completed)
Oncce the modeel has finisheed all of thee computatio ons successffully, you caan begin vieewing the
resuults. Severaal output optiions area avvailable fromm the View menu
m bar onn the HEC-R
RAS main
winndow. Thesee options incllude:
a) CCross sectionn plots
b) PProfile plots
c) GGeneral proffile plot
d) RRating curvees
e) XX-Y-Z persppective plots
f) DDetailed tabuular output at
a a specific ccross section (cross sectio
on table)
g) LLimited tabuular output att many cross sections (profile table)
BBegin by plotting a cross section. Select Crosss Sections from the Viiew menu baar on the
HHEC-RAS main m window w. Any crosss section can n be plotted by
b selecting the appropriiate river,
rreach and rivver station (See Figure 11.110). Seeveral plottin ng features aare available from the
OOptions mennu bar on th he cross sectiion plot winndow. Thesee options incclude: zoom in; zoom
oout; selectingg which planns, profiles, variables to plot; and co
ontrol over liines, labels, symbols,
sscaling etc.
N
Next plot a water surfacce profile. S elect Waterr Surface Prrofiles from the View menu
m bar.
C
Click on thee Play button
n, , to vieew the Anim
mation Contrrol window (see Figuree 11.111).
Click on the Expand buttton, , to ssee the entiree control.
C
SSet the Zeroo Delay horizzontal scroll bar as indiccated (in Figu
ure 11.111) and click on n the play
bbutton, , to t view the changing
c watter surface leevels and eneergy grade liines of the river reach
a shown in Figure 11.1112 for the maximum
pprofile plot. This should give you a pprofile plot as m
wwater surfacee.
F
Figure 11.112: Water Su
urface Profille (for Riet River:
R Loweer reach- M
Maximum WS)
AD
DDING AN IINLINE STR
RUCTURE (WEIR)
In the neext section ann inline weir will be addeed in the Blesbok River – Lower reacch.
T
To enter inline weir dataa click on thee Inline stru
ucture button
n on the Geoometric data a window
(Figure 11.558). Once thhe Inline stru ucture button is pressed, the Inline sttructure Da
ata Editor
w a shown in Figure
will appear as F 11.1114.
FFrom the Op ptions menuu, select Add d an Inline structure from f the listt. An input box will
aappear prommpting you to enter a riverr station iden
ntifier for thee inline structture.
EEnter 6.5 as shown in Fig
gure 11.115 .
Figurre 11.115: In
nline weir rivver station (Blesbok
( Riv
ver: Lower rreach)
E
Enter the Deescription off the inline sttructure: Mayya weir (see Figure 11.1 16).
TThe first inpput at the topp of the editoor is the distance from th
he upstream cross section
n and the
ddeck (i.e. rivver station 7). This distannce is 10 m.
T
The weir flow nt selected forr this analysiis is 1,44.
w coefficien
Click the OK
C K button.
Figurre 11.118: In
nline weir (V
Viewing the weir)
EEnter the Pilot flow off 50 l/s i.e. 0,05 m³/s (Pilot dischaarge for leak
akage or to keep the
ddownstream channel wett at low flow s) (see Figu
ure 11.118).
IIn addition to
t uncontrolled overflow weirs, the user u can add gated spillw ways (this is optional).
o
TTo add gatedd spillways to he Gate butto
t the structuure, press th on on the Innline Structuure Data
eeditor. Once this button n is pressed, the gated editor will apppear as show wn in Figurre 11.119
((except yourrs will be blaank until youu have enteredd some data)).
Figure
F 11.1199: Gated Sp
pillway Edito
or
TThe Gated Spillway ed ditor is simillar to the Cuulvert editorr in concept.. The user enters
e the
pphysical desscription of the
t gates, ass well as the required co oefficients, inn the Gated Spillway
eeditor. Thee functionality of the gaates is defin ned as part of o the Unsteeady Flow data.
d The
ffollowing is a list of the data
d containeed on this ed ditor:
GGate Group p - The Gate Group is auutomatically y assigned to "Gate #1" thhe first time you open
tthe editor. The
T user can enter up to 10 different Gate Groupss at each part rticular river crossing,
aand each gatte group can have up to 25 identical gate openings. If all of tthe gate openings are
eexactly the same,
s then only one gatee group need ds to be enterred. If the usser has gate openings
tthat are diffeerent in shap
pe, size, elevvation, or hav
ve different coefficients, then additio onal Gate
GGroups mustt be added fo or each Gate type.
W
Width - This field is used for enterinng the width of the gate in
n meters (Ennter 1,2 m).
Gate Type - This field is used for selecting the type of gate. A number of different gate types
are available. Select from the drop down list Sluice.
Sluice Discharge Coefficient - This field is used for entering the coefficient of discharge for
the gate opening. This coefficient ranges from 0,5 to 0,7 for sluice gates. For this sluice gate
enter 0,6.
Orifice Coefficient - This field is used to enter an orifice coefficient, which will be used for
the gate opening when the gate becomes more than 80 percent submerged. Between 67 percent
and 80 percent submerged, the program uses a transition between the fully submerged orifice
equation and the free flow equations. When the flow is less than 67 percent submerged, the
program uses the free flow gate equations. Enter an Orifice Coefficient of 0,8.
Head Reference – This field is used to select the reference point for which the upstream
energy head will be computed from. The default is the gate sill (invert), which is normally
used when the flow through the gate goes out into a channel. If the gate causes the flow to jet
out freely into the atmosphere, then the head reference should be selected as the centerline
elevation of the gate opening. If the gate crest is an ogee spillway crest, then the center of the
gate opening should be used. Ogee spillway crests are normally designed to follow the shape
of water jetting freely into the atmosphere. For this exercise select from the drop down list Sill
(invert).
Weir Shape - This parameter allows the user to select between a Broad Crested, Sharp
Crested shape weir and an Ogee shaped weir. Select Broad Crested.
Weir Coefficient - This field is used for entering a weir coefficient that will be used for the
gate opening. This coefficient will only be used when the gate is opened to an elevation
higher than the upstream water surface elevation. When this occurs, the flow through the gate
is calculated as weir flow. Enter the Weir Coefficient of 1,67.
Centerline Stations - This table is used for entering the centerline stationing of the identical
gate openings. The user should enter a different centerline stationing for each gate opening
that is part of the current gate group. All gate openings within the same gate group are exactly
identical in every way, except their centerline stationing. As a user adds new centerline
stationing values, the number of identical gates in the group is automatically incremented and
displayed in the field labelled "# Openings". Enter the Centerline Station as 16,25 m.
Once all of the data for the gates has been entered (as shown in Figure 11.119), the user needs
to press the OK button for the data to be accepted. If the user presses the OK button, this does
not mean that the data is saved to the hard disk; it is only stored in memory and accepted as
being good data. This data is part of the geometry data, and is stored in the geometric data file.
The data can be stored to the hard disk by selecting from the File menu of the Geometric
Data window Save Geometry Data As.
After selecting this option you will be prompted to enter a Title for the geometric data (Figure
11.120). Enter “Base Geometry + bridge + weir” for this exercise, and then press the OK
button. A file name is automatically assigned to the geometry data based on what you entered
for the project file name i.e. Exercise2.g03.
T
The newly added
a weir will wn on the Geometric
w be show G Data
D editor w
window as shown
s in
F
Figure 11.1222.
TThis option allows the user to controol the openinng and closin
ng of gates bbased on the elevation
oof the waterr surface upsstream of thee structure. A gate begin
ns to open w when a user specified
eelevation is exceeded inn this case 233,5 m (see Figure
F 11.12
25). The gatte will beginn to close
t water surrface elevatioon reaches 22,5 m.
aagain when the
T
The gate opeens at a rate specified
s by the user (0,1
1 m/min).
T
The closing of the gate iss at a user sppecified rate (also
( 0,1 m/m
min.).
T
The user must
m also entter a maxim mum and miinimum gatee opening. FFor this exeercise the
M
Maximum Gate
G Openinng is set at 1,,1 m and the Minimum Gate
G Openinng is set at 0 m.
Figure 11
1.125: Elevattion Contro
olled Gates data
d editor
Figu
ure 11.126: S
Saving Unstteady Flow Data
D
TThe completted Unsteady y Flow Dataa screen shou ude the addittional boundary at the
uld now inclu
iinline structuure (RS 6.5) as shown in Figure 11.1
127.
Figure
F 11.1227: Unstead
dy Flow Data
a
IInitial condditions. Theere is no neeed to add an ny initial flow conditionns at the new
w internal
ccontrol sincee the flow en
ntered at Rivver station 8 will
w also flow w over the w
weir at the sttart of the
aanalysis.
P
Performing Unsteady Flow
F Calculaations with the
t Inline Sttructure andd the Bridgee
TThe first stepp is to put to w Plan. Thee Plan defines which geoometry and flow
ogether a new f data
aare to be useed, as well as providing a title and sho
ort identifier for the run.
YYou will be prompted to o enter a shoort identifier.. Enter a title of Base+bb+w in the Short
S ID
bbox and clickk on the OK
K button.
S metry file annd Unsteady flow file from the dropp down list i.e. Base
Select the coorrect Geom
G
Geometry + bridge + weeir and Base Unsteady Fllow + gates respectively.
r .
SSimulation time
t window w
TThe user is required
r to enter
e a time window thaat defines the start and eend of the siimulation
pperiod. The time window w requires a starting datee and time an
nd an ending date and tim
me.
CComputatioonal settings
TThe computaational Settinngs area conntains the Co
omputationa al Interval, HHydrograph h Output
IInterval, Deetailed Outpput Interval, the name an nd path of the output DSSS file, and whether
w or
nnot the progrram is run in
n a mixed fllow regime. The compu utation intervval is probably one of
tthe most im mportant paraameters enteered into the model. This T should be small en nough to
aaccurately describe
d the rise and falll of the hyd
drographs being routed bbut not smalll to take
fforever to coompute.
A
Also select thhe Mixed Fllow Regimee option sincee mixed flow
w might occuur at the flow
w over the
weir. The coompleted Un
w nsteady Floww Analysis screen is showwn in Figuree 11.130.
Clicck on the Co
ompute butto
on to run the Unsteady Flow
F Analysiis.
VIEW
WING THE RESULTS
Oncce the modeel has finisheed all of thee computatioons successffully, you caan begin vieewing the
resuults. Several output optiions are avaiilable from the View menu
m bar on the HEC-R RAS main
winndow. Thesee options incllude:
a) C
Cross sectionn plots
b) P
Profile plots
c) G
General proffile plot
d) R
Rating curvees
e) X
X-Y-Z persppective plots
f) D
Detailed tabuular output at
a a specific ccross section (cross sectio
on table)
g) L
Limited tabuular output att many cross sections (profile table)
N
Next plot a water surfacce profile. S elect Waterr Surface Prrofiles from the View menu
m bar.
C
Click on thee Play button
n, , to vieew the Anim
mation Contrrol window (see Figuree 11.133).
Click on the Expand buttton, , to ssee the entiree control.
C
SSet the Zeroo Delay horizzontal scroll bar as indicated (in Figu ure 11.133) and click on n the play
bbutton, , to t view the changing
c watter surface leevels and eneergy grade liines of the river reach
a shown in Figure 11.1134 for the maximum
pprofile plot. This should give you a pprofile plot as m
wwater surfacee.
At the ennd of this exeercise the folllowing objeectives shoulld have been met:
B
Be able to a river system
m containing m
more than on
ne river reach
h
K
Know how too enter unsteeady flow daata and definiing and enterring the bounndary controls
B
Be able to annalyse the riv
ver system (uunsteady flow
w)
K
Know how too define a brridge structurre and a weirr
K
Know how too set the bou
undaries for uunsteady flow
w with openiing and closiing gates
K
Know how too obtain resu
ults from an uunsteady flow
w analysis
Question
ns
1. D
Define the fllow type in th
he river.
2. D
Does the flow t cross secttion flow ontto the banks of the river ssystem?
w at any of the
3. A
At what tim
me does the maximum
m floow at cross section 40 occur
o and whhat is the peeak flow?
W
What is the normal
n flow depth duringg this peak fllow?
5. Describe the reason for the differences in the Rating curves (Flow rate versus Water Surface
Elevation) for cross section 40 and cross section 100.
6. What effect does the bridge have on the maximum water surface level at cross section 70?
9. What will happen if the flow hydrographs of the Upper reach of the Riet river is increased
with 100%?
10. What increase in flow velocity is experienced due to the bridge compared to the before
scenario?
11. What is the maximum flow through the gate at the inline structure and at what time step does
this occur?
12. What effect does the added inline structure have on the flow conditions at the bridge (flow
depth, velocity etc.)?
13. What will be the effect if the Computational Interval, Hydrograph Output Interval and
Detailed Output Interval are also set to 15 minutes and the rate of opening and closing of the
gates are set to 1,2 m/min.
Problem
m description
n Example 12.1:
1
Calculatte the maxim mum infiltrattion rate (mm m/day), whiich may be discharged
d vvia the depiccted sub-
surface hherringbone system to a main drainaage pipe. Th he diameter of
o the centrall pipe is 150 mm and
its slope 1:500. The diameter off the laterals is 100 mm and a their slop
pes 1: 100. TThe Manning g n-value
for the ppipes is 0,0111s/m1/3. Figure 12.1 refleects the layou
ut.
Solution
n Example 12.1:
1
For each
h lateral, floowing 70% full:
f
q
26,922 x10 d S 0,7
6 8/3
0
1/2
nA
S
A S (L ) Refer
R to Figu
ure 12.2
2
L 152 152 21,21
1m
L
S 100,61m (V
Valid in thiss example since laterals are placed aat 45º and dim
mensions
2
are
a equal 15 m x 15 m)
10,61
A 10,61 (221,21)
2
281,3 m 2
q
26,92 x 10 0,1 0,01
6 8/3
0 0,7 1312 mm/daay
1/2
0,011281,3
12-242 Su
ub-surface drrainage
F
Figure 12.2: General vieew of a herringbone dra
ainage system
m
Q = (1,312)(281,3))(14)
Q = 51167 m³/day
2/3 1/2
π
4
0,15 2
0,155
4
1
Q 500
0,011
Q = 0,000805 m³/s
Q ≈ 6955 m³/day
695
q max
281,314
qmax = 1177 mm/day (<< 1312 mm
m/day)
12-243 Su
ub-surface drrainage
APPENDICES
A-244 Appendices
APPENDIX 3A
STATISTICAL ANALYSIS
A-245 Appendices
Single station direct statistical analysis
The following frequency distributions are discussed for untransformed and log10-transformed data:
Untransformed data
Normal, Extreme Value Type 1 and General Extreme Value
Log10-transformed data
Log-Normal, Log-Gumbel and Log-Pearson Type III
Step 1: Determine the mean, standard deviation and skewness coefficient of the raw data and the log10
transformed data as follows:
Mean x
x …(3A.1)
N
0,5
x x
2
Skewness coefficient
g
N
x x
3
…(3A.3)
N 1N 2 s3
s
Coefficient of variation cV …(3A.4)
x
where:
x = observed value (or of the logarithm of the observed value for the log
distributions)
x = mean of observed values (or of the logarithm of the observed value for the
log distributions)
N = the total number of observations
s = the standard deviation of the observed values (or of the logarithm of the
observed values)
g = skewness coefficient
cV = coefficient of variation
Step 2: The peak value for the desired return period and assumed distribution function can be derived
for each of the frequency distributions as follows:
Normal distribution
The normal distribution is applicable where the observed values represent the effects of a large
number of independent processes. The distribution is symmetrical about the mean and is
therefore only suitable for data where the skewness coefficient (g) is equal to, or close to zero.
The spread about the mean is a function of the coefficient of variation. For high coefficient of
variation values, the bottom tail may extend below zero and may result in negative flows
being generated when the distribution is applied to untransformed data.
y
x x …(3A.6)
s
A-246 Appendices
The value of y for a given value of G(y) cannot be solved directly from Equation 3A.4, and
hence published tables have to be used. Based on the return period, read from Table 3A.1b
the value of G(y) and obtain y from Table 3A.1a.
Q T x sy …(3A.7)
From Table 3A.4 (g = 1,14) read the value of WT for the required return period. Calculate QT
directly using:
For the known value of the skewness coefficient (g) read off the value WT from Table 3A.4
and the values of k, E(y) and var(y) from Table 3A.2 by using linear interpolation.
For this distribution, the logarithms of the data are assumed to be normally distributed. Based
on the skewness coefficient (g), obtain the value WT for the required return period from
Table 3A.3.
Q T antilog log(x) s log WT …(3A.11)
where:
slog = the standard deviation of the logarithms of the observed values as shown in
equation 3A.11
0,5
logx logx
2
s log …(3A.12)
N 1
and
log x = the logarithm of the mean of the observed values
Confidence bands
The confidence, with which the values of the magnitude-return period relationships are
estimated, depends on the number of observations contained in the data set. The greater the
number of observations, the greater the degree of assurance, and subsequently the narrower
the confidence band.
A-247 Appendices
The displacement of the two-sided 95% confidence band about the estimated value can be
read from Table 3A.3 where N is the number of observations. The 95% confidence limits are:
Q T(95%) antilog log(x) s log WT Wα …(3A.13)
where:
From Table 3A.4 (g = 1,14) read the value of WT for the required return period. Calculate QT
directly using:
Q T antilog log(x) s log 0,780WT 0,450 …(3A.14)
From Table 3A.3 determine the value of WT for the known skewness coefficient (g) of the
log-transformed data by linear interpolation.
Q T antilog log(x) s log WT …(3A.15)
Based on an example from Flood Risk Reduction Measures by WJR Alexander the incorporation of
the influence of historical information, missing data and outliers is required to determine the
confidence of the results. It is thus required to calculate the historically weighted mean ( x h ), standard
deviation ( s h ) and skewness coefficient ( g h ).
WT x b x a
xh …(3A.16)
YT WT LW
WT d 2b d a2
0,5
sh …(3A.17)
YT WT LW 1
YT WTLWWT d 3b d 3a
3
gh s …(3A.18)
YT LWWT 1YT WTLW 2
where:
YT = total time span (= NA + NB + NC)
WT = weight applied to data = (YT – NA) / NB
NA = floods equal to or above the high threshold
NB = floods between high and low thresholds
NC = missing data
LW = low outliers including zero flows
ZR = zero flows
and where: xa = is the value of a peak equal to or above the high threshold
xb = is the value of a peak below the high threshold
da and db = are deviations of xa + xb from x h
All values being the logarithms of the data.
A-248 Appendices
These historically weighted values of the mean, standard deviation and the skewness coefficient are
then used in the equations for the LN/MM, LP3/MM, EV1/MM and GEV/MM distribution in the
usual way (3.1).
Evaluating Example 3.2 for the Tsitsa River, utilizing Equations 3A.1 to 3A.15 (if the missing data is
not included in the statistical analysis) will provide the following results:
The next set of results is based on the historically weighted mean ( x h ), standard deviation ( s h ) and
skewness coefficient ( g h ) (Table 3A.6a and b) which incorporates the missing data.
A-249 Appendices
Return period N/MM EV1/MM GEV/MM LN/MM LEV1/MM LP3/MM
2 477 416 418 348 304 355
5 795 749 726 702 635 706
10 962 971 934 1012 1034 999
20 1099 1183 1137 1367 1651 1323
50 1254 1457 1407 1924 3022 1790
100 1357 1664 1615 2429 4764 2202
Arrange the observed data in descending order of magnitude and assign to each value a rank
number starting from one. Determine the plotting position (return period) for each value using
the Weibull formula. The general equation is given below and the values for the constants a
and b are provided in Table 3A.7.
nl a
Τ … (3A.19)
m-b
where:
T = return period in years
nl = length of record in years
m = number, in descending order, of the ranked annual peak floods
a = constant (see Table 3A.7)
b = constant (see Table 3A.7)
If the horizontal axis has a probability classification, the probability (P) is calculated as:
1
… (3A.20)
Some of the commonly used plotting positions recommended for use in hydrological analyses
are given in Table 3A.7. If several distributions are plotted on a single graph, then the general
purpose Cunane plotting position should be used.
Plot the values against their estimated return periods on log-probability paper; draw the best
fitting straight line through the plotted points and extrapolate to determine the estimated
maximum value for the required return period. Alternatively utilize software such as Utility
Programs for Drainage or HEC-SSP included on the supporting flash drive.
A-250 Appendices
Table 3A.1: Properties of the standardized normal distribution
Table 3A.1a Table 3A.1b
Standardized normal distribution Standardized normal distribution
y G(y)% y G(y)% T G(y)% WT
0,00 50,00 -0,00 50,00 1000 0,1 -3,09
0,05 51,99 -0,05 48,01 500 0,2 -2,88
0,10 53,98 -0,10 46,02
0,15 55,96 -0,15 44,04
200 0,5 -2,58
0,20 57,93 -0,20 42,07 100 1,0 -2,33
0,25 59,87 -0,25 40,13 50 2,0 -2,05
0,30 61,79 -0,30 38,21 20 5,0 -1,64
0,35 63,68 -0,35 36,32
0,40 65,54 -0,40 34,46 10 10,0 -1,28
0,45 67,36 -0,45 32,64 5 20,0 -0,84
0,50 69,14 -0,50 30,86 2 50,0 0,00
0,55 70,88 -0,55 29,12
0,60 72,57 -0,60 27,43
5 80,0 0,84
0,65 74,22 -0,65 25,78 10 90,0 1,28
0,70 75,81 -0,70 24,19 20 95,0 1,64
0,75 77,34 -0,75 22,66 50 98,0 2,05
0,80 78,81 -0,80 21,19
0,85 80,24 -0,85 19,76 100 99,0 2,33
0,90 81,59 -0,90 18,41 200 99,5 2,58
0,95 82,89 -0,95 17,11 500 99,8 2,88
1,00 84,13 -1,00 15,87
1,05 85,31 -1,05 14,69
1000 99,9 3,09
1,10 86,43 -1,10 13,57 5000 99,98 3,55
1,15 87,49 -1,15 12,51 10000 99,99 3,72
1,20 88,49 -1,20 11,51
1,25 89,44 -1,25 10,56
1,30 90,32 -1,30 9,68
1,35 91,15 -1,35 8,85
1,40 91,26 -1,40 8,08
1,45 92,65 -1,45 7,35
1,50 93,32 -1,50 6,68
1,55 93,94 -1,55 6,06
1,60 94,52 -1,60 5,48
1,65 95,05 -1,65 4,95
1,70 95,54 -1,70 4,46
1,75 95,99 -1,75 4,01
1,80 96,41 -1,80 3,59
1,85 96,78 -1,85 3,22
1,90 97,13 -1,90 2,87
1,95 97,44 -1,95 2,56
2,00 97,72 -2,00 2,28
2,05 97,98 -2,05 2,02
2,10 98,21 -2,10 1,79
2,15 98,43 -2,15 1,57
2,20 98,61 -2,20 1,39
2,25 98,78 -2,25 1,22
2,30 98,93 -2,30 1,07
2,35 99,06 -2,35 0,94
2,40 99,18 -2,40 0,82
2,45 99,29 -2,45 0,71
2,50 99,38 -2,50 0,62
2,55 99,46 -2,55 0,54
2,60 99,53 -2,60 0,47
2,65 99,60 -2,65 0,40
2,70 99,65 -2,70 0,35
2,75 99,70 -2,75 0,30
2,80 99,74 -2,80 0,26
2,85 99,78 -2,85 0,22
2,90 99,81 -2,90 0,19
2,95 99,84 -2,95 0,16
3,00 99,86 -3,00 0,14
3,05 99,88 -3,05 0,16
3,10 99,90 -3,10 0,10
3,15 99,92 -3,15 0,08
3,20 99,93 -3,20 0,07
3,25 99,94 -3,25 0,06
3,30 99,95 -3,30 0,05
3,35 99,96 -3,35 0,04
3,40 99,97 -3,40 0,03
3,45 99,97 --3,45 0,03
3,50 99,98 --3,50 0,02
3,55 99,98 --3,55 0,02
3,60 99,98 --3,60 0,02
3,65 99,99 --3,65 0,01
3,70 99,99 --3,70 0,01
3,75 99,99 --3,75 0,01
A-251 Appendices
Table 3A.2: Parameters of the standardized general extreme value distribution
Standardized general extreme value distribution
g k E(y) var(y)
-2,000 1,406 -1,247 3,204
-1,900 1,321 -1,182 2,505
-1,800 1,240 -1,127 1,984
-1,700 1,163 -1,080 1,590
-1,600 1,089 -1,041 1,287
-1,500 1,018 -1,008 1,052
-1,400 0,950 -0,980 0,868
-1,300 0,885 -0,957 0,721
-1,200 0,824 -0,938 0,602
-1,100 0,765 -0,922 0,507
-1,000 0,708 -0,910 0,428
-0,900 0,655 -0,901 0,362
-0,800 0,604 -0,894 0,307
-0,700 0,555 -0,889 0,261
-0,600 0,509 -0,887 0,222
-0,500 0,465 -0,886 0,188
-0,400 0,424 -0,886 0,159
-0,300 0,384 -0,888 0,134
-0,200 0,346 -0,892 0,112
-0,100 0,311 -0,896 0,094
0,000 0,277 -0,901 0,077
0,100 0,245 -0,907 0,063
0,200 0,215 -0,914 0,050
0,300 0,187 -0,922 0,039
0,400 0,160 -0,930 0,030
0,500 0,134 -0,938 0,022
0,600 0,110 -0,947 0,016
0,700 0,088 -0,956 0,010
0,800 0,067 -0,966 0,006
0,900 0,047 -0,975 0,003
1,000 0,028 -0,985 0,001
1,100 0,010 -0,994 0,000
1,200 -0,006 1,004 0,000
1,300 -0,022 1,013 0,001
1,400 -0,037 1,023 0,002
1,500 -0,050 1,032 0,005
1,600 -0,063 1,041 0,008
1,700 -0,075 1,049 0,011
1,800 -0,086 1,058 0,016
1,900 -0,097 1,066 0,021
2,000 -0,107 1,074 0,026
2,100 -0,116 1,082 0,032
2,200 -0,125 1,089 0,038
2,300 -0,133 1,097 0,044
2,400 -0,140 1,104 0,051
2,500 -0,148 1,110 0,058
2,600 -0,154 1,116 0,065
2,700 -0,160 1,123 0,072
2,800 -0,166 1,128 0,080
2,900 -0,172 1,134 0,087
3,000 -0,177 1,139 0,094
3,100 -0,182 1,145 0,102
3,200 -0,187 1,150 0,110
3,300 -0,191 1,154 0,117
3,400 -0,195 1,159 0,125
3,500 -0,199 1,163 0,132
3,600 -0,203 1,168 0,140
3,700 -0,207 1,172 0,148
3,800 -0,210 1,176 0,155
3,900 -0,213 1,180 0,163
4,000 -0,217 1,183 0,170
4,100 -0,220 1,187 0,178
4,200 -0,223 1,191 0,186
4,300 -0,225 1,194 0,193
4,400 -0,228 1,197 0,201
4,500 -0,231 1,201 0,208
4,600 -0,233 1,204 0,215
4,700 -0,236 1,207 0,223
4,800 -0,238 1,210 0,230
4,900 -0,240 1,213 0,237
5,000 -0,242 1,215 0,244
A-252 Appendices
Table 3A.3a: Values of the standardized variate WT for the normal and exponential
distributions
Return Non- Normal distribution
Exponential
period exceedance Confidence limits Wα
WT distribution
(years) probability 75% 95%
2 0,50 0,00 1,63 2N 2,77 2N 0,69
5 0,80 0,84 1,89 2N 3,23 2N 1,61
10 0,90 1,28 2,20 2N 3,74 2N 2,30
20 0,95 1,64 2,49 2N 4,25 2N 3,00
50 0,98 2,05 2,87 2N 4,89 2N 3,91
100 0,99 2,33 3,13 2N 5,34 2N 4,61
200 0,995 2,58 3,38 2N 5,76 2N 5,30
500 0,998 2,88 3,69 2N 6,27 2N 6,21
1000 0,999 3,09 3,91 2N 6,66 2N 6,91
10000 0,9999 3,72 4,58 2N 7,80 2N 9,21
Table 3A.3b: Values of the standardized variate WT for the Pearson Type III distribution
Return Pearson Type III distribution (Values of WT)
period g
(years) -1,0 -0,8 -0,6 -0,4 -0,2 0,0 0,2 0,4 0,6 0,8 1,0
2 0,16 0,13 0,10 0,07 0,03 0,00 -0,03 -0,70 -0,10 -0,13 -0,16
5 0,85 0,87 0,86 0,86 0,85 0,84 0,83 0,82 0,80 0,78 0,76
10 1,13 1,17 1,20 1,23 1,26 1,28 1,30 1,32 1,33 1,34 1,34
20 1,32 1,39 1,46 1,52 1,59 1,64 1,70 1,75 1,80 1,84 1,88
50 1,49 1,61 1,72 1,83 1,94 2,05 2,16 2,26 2,36 2,45 2,54
100 1,59 1,73 1,88 2,03 2,18 2,33 2,47 2,62 2,76 2,89 3,02
200 1,66 1,84 2,02 2,20 2,39 2,58 2,76 2,95 3,13 3,31 3,49
500 2,88
1000 1,79 2,02 2,27 2,53 2,81 3,09 3,38 3,67 3,96 4,24 4,53
10000 1,88 2,18 2,53 2,90 3,30 3,72 4,15 4,60 5,05 5,50 5,96
Table 3A.4a: Values of the standardized variate WT for the general extreme value distribution
(EV1 & EV2)
General extreme value (Values of WT)
Return
g
period
(years) 1,14 1,2 1,4 1,6 1,8 2,0 2,5 3,0 3,5 4,0 4,5 5,0 5,5 6,0
EV1 EV2
2 0,37 0,37 0,37 0,37 0,37 0,37 0,38 0,38 0,38 0,38 0,38 0,38 0,38 0,38
5 1,50 1,51 1,55 1,58 1,60 1,63 1,68 1,72 1,75 1,77 1,79 1,80 1,82 1,83
10 2,25 2,28 2,35 2,43 2,49 2,55 2,67 2,76 2,84 2,90 2,94 2,98 3,01 3,04
20 2,97 3,01 3,15 3,28 3,40 3,50 3,73 3,91 4,05 4,16 4,25 4,33 4,39 4,45
50 3,90 3,97 4,22 4,45 4,66 4,86 5,28 5,62 5,89 6,12 6,30 6,46 6,59 6,70
100 4,60 4,71 5,05 5,38 5,68 5,97 6,59 7,10 7,52 7,86 8,15 8,39 8,59 8,77
200 5,30 5,44 5,90 6,34 6,76 7,16 8,04 8,77 9,38 9,89 10,31 10,67 10,97 11,24
500 6,21 6,41 7,05 7,68 8,29 8,87 10,19 11,32 12,26 13,06 13,74 14,32 14,81 15,24
1000 6,91 7,15 7,95 8,75 9,53 10,29 12,02 13,53 14,82 15,92 16,86 17,66 18,36 18,96
10000 9,21 9,64 11,13 12,68 14,25 15,82 19,65 23,19 26,34 29,13 31,58 33,73 35,63 37,31
A-253 Appendices
Table 3A.4b: Values of the standardized variate WT for the general extreme value distribution
(EV3)
General extreme value (Values of WT)
Return
g
period
(years) -1,0 -0,8 -0,6 -0,4 -0,2 0,0 0,2 0,4 0,6 0,8 1,0
EV3
2 0,33 0,34 0,34 0,34 0,35 0,35 0,36 0,36 0,36 0,37
5 1,01 1,06 1,12 1,17 1,23 1,28 1,33 1,38 1,43 1,47
10 1,28 1,37 1,46 1,57 1,67 1,78 1,89 2,00 2,10 2,19
20 1,44 1,57 1,71 1,86 2,02 2,19 2,37 2,54 2,71 2,86
50 1,58 1,74 1,93 2,15 2,38 2,64 2,90 3,18 3,45 3,72
100 1,64 1,83 2,05 2,31 2,60 2,92 3,26 3,62 3,99 4,35
200 1,68 1,89 2,14 2,44 2,77 3,16 3,58 4,02 4,49 4,97
500 1,72 1,95 2,23 2,56 2,96 3,42 3,94 4,51 5,13 5,76
1000 1,74 1,98 2,27 2,64 3,07 3,59 4,19 4,86 5,58 6,35
10000 1,76 2,02 2,36 2,78 3,32 4,00 4,83 5,81 6,96 8,24
A-254 Appendices
APPENDIX 3B
STANDARD DESIGN FLOOD METHOD
A-255 Appendices
Table 3B.1: Information required for the calculation of the SDF
SAWS M R C100 MAP MAE
C2
Basin station SAWS site
(mm) (days) (%) (%) (mm) (mm)
number
1 546 204 Struan 56 30 10 40 550 1800
2 675 125 Autoriteit 62 44 5 30 450 1900
3 760 324 Siloam 64 28 5 40 470 1700
4 553 351 Waterval 58 20 10 50 630 1600
5 680 059 Leydsdorp 78 10 15 70 620 1700
6 369 030 Siloam 51 54 15 60 670 1500
7 328 726 Olivine 49 39 15 60 510 1700
8 322 071 Danielskuil 47 39 5 20 380 2100
9 258 452 Jacobsdal 43 47 15 60 380 1800
10 233 049 Wonderboom 54 55 10 50 560 1600
11 236 521 Mashai 39 66 40 80 430 1400
12 143 258 Scheurfontein 39 52 5 30 290 2100
13 284 361 Wilgenhoutsdrif 40 55 5 15 70 2600
14 110 385 Middelpos 25 13 10 30 140 2400
15 157 874 Garies 22 11 5 20 130 2100
16 160 807 Loeriesfontein 28 11 10 40 210 1900
17 84 558 Elandspoort 45 1 40 80 500 1500
18 22 113 La Motte 59 4 30 60 810 1400
19 69 483 Letjiesbos 34 16 10 35 160 2200
20 34 762 Uitenhage 53 12 15 60 480 1600
21 76 884 Albertvale 45 23 10 35 460 1700
22 80 569 Umzoniana 84 26 15 60 820 1200
23 180 439 Insizwa 60 45 10 80 890 1200
24 240 269 Newlands 76 15 15 80 910 1200
25 239 138 Whitson 55 9 10 80 830 1200
26 336 283 Nqutu 61 17 15 50 760 1500
27 339 415 Hill Farm 85 17 30 80 890 1400
28 483 193 Maliba Ranch 75 54 15 60 740 1400
29 556 088 Mayfern 66 11 15 50 740 1600
A-256 Appendices
Table 3B.2: D
Daily rainfalll from TR102
A-257 Appendices
Table 3B
B.2: Daily rrainfall from
m TR102 (continued)
A-258 Appendices
Table 3B
B.2: Daily rrainfall from
m TR102 (continued)
A-259 Appendices
APPENDIX 3C
STANDARD FLOOD CALCULATION FORMS
A-260 Appendices
RAT
TIONAL ME
ETHOD (AL
LTERNATIIVE 1)
A-261 Appendices
RAT
TIONAL ME
ETHOD (AL
LTERNATIIVE 1)
Table 3C.55
Return p period (years)) 2 5 10 20 50
0 100
Adjustmeent factor (Ft) for
steep andd impermeablee 0,7
75 0,80 0,85 0,90 0,9
95 1,00
catchmennts
Adjustmeent factor (Ft) for
flat and ppermeable 0,5
50 0,55 0,60 0,67 0,8
83 1,00
catchmennts
R1 - Page 2/2
A-262 Appendices
RAT
TIONAL ME
ETHOD (AL
LTERNATIIVE 2)
A-263 Appendices
RATIONAL M
METHOD (A
ALTERNAT
TIVE 2)
LEGGEND TABL LE
Ration
nal method (A Alt 2)
ID Referen nce
F
Figure 3.5 or SA
S
W
Weather Serviices
T
TR102 or otheer
F
Figure 3.9
T
Table 3C.1
T
Table 3C.2
T
Table 3C.3
T
Table 3C.4
T
Table 3C.5
T
Table 3C.6
F
Figure 3.10
T
TR102
Table 3C.6
Sellection criteriia Calculation
C method
m
Modiified Hershfield equation
TC < 6 hours Pt, T 11,130,41 0,6
64lnT 0,11 0,27lnt 0,779M 0,69 R 0,20
Linear interpolation between calcu ulated modifieed Hershfield
6 hours ≤ TC < 24 hours
equatioon point rainfaall and 1-day point
p rainfall from TR102
R2 - Page 2/2
A-264 Appendices
RAT
TIONAL ME
ETHOD (AL
LTERNATIIVE 3)
A-265 Appendices
RAT
TIONAL ME
ETHOD (AL
LTERNATIVE 3)
Table
T 3C.7
LEEGEND TABLE Rural
R (C1)
Ratioonal method (Alt
( 3) Componeent Classification Mean annua
al rainfall (mm
m)
ID Referen nce 600 600
0 - 900 9000
FFigure 3.5 or SA
S Wetlaands and pans (<3%) 0,01 0,03
0 0,05
WWeather Serviices
Surface sslope Flat arreas (3 to 10%
%) 0,06 0,08
0 0,11
TTable 3C.7
TTable 3C.8 (CS) Hilly (10 to 30%) 0,12 0,16
0 0,220
TTable 3C.9 Steep areas (>30%)) 0,22 0,26
0 0,30
TTable 3C.10 Very permeable
p 0,03 0,04
0 0,05
TTable 3C.11 Permeeable 0,06 0,08
0 0,10
Permeabiility
# FFigure 3.12 annd Figure (CP) Semi--permeable 0,12 0,16
0 0,220
33.13
FFigure 3.8 (oor Figure Imperrmeable 0,21 0,26
0 0,30
33.26 DM) Thickk bush and planntation 0,03 0,04
0 0,05
FFigure 3.10, Figure Vegetatioon Light bush and farm m-lands 0,07 0,11
0 0,15
33.13 or other (CV) Grasslands 0,17 0,21
0 0,225
No veegetation 0,26 0,28
0 0,30
Table 3C.111
Return p period (years)) 2 5 10 20 50
0 100
Adjustmeent factor (Ft) for
steep andd impermeablee 0,7
75 0,80 0,85 0,90 0,9
95 1,00
catchmennts
Adjustmeent factor (Ft) for
flat and ppermeable 0,5
50 0,55 0,60 0,67 0,8
83 1,00
catchmennts
Note: # CCalculate the point intensity making use off the provided Design Rainffall estimation
y rainfall by m n software.
The exacct point intenssity can be caalculated by mmeans of lineear interpolatiion between ttwo consecutiive values
considering the time off concentration.
R3 - Page 2/2
A-266 Appendices
UNIT HY
YDROGRAP
PH METHO
OD
Description
n of catchmen
nt
River detaiil
Calculated
d by Date
Physsical characteeristics
Size of catcchment (A) kkm² Veldd type
Longest waatercourse (L) kkm Lagg (TL)
Average sloope (Sav) mm/m Coeefficient (KU)
Length to catchment centtroid (LC) kkm Peakk discharge K A m³/s
Mean annuaal precipitatioon(MAP)# mmm of unit
u Qp U
TL
Catchment L LC hydrograph (QP)
IC
index S av
LEGEEND TABLE
Unit Hydrrograph method
ID Reference ID Referencee
F
Figure 3.5 or SA
S Weather Services
S Figure 3.8 (or Figure
F 3.26 D
DM for large arreas)
F
Figure 3.15 Figure 3.17
F
Figure 3.16 Table
T 3.9
T
Table 3.6 Paragraph [Link] in DM
D
Design rainfalll Database or Figure 3.7 QPiT = Qp x [(S1 – S2T)/TSD]mmax
U - Page 1/1
A-267 Appendices
ST
TANDARD DESIGN FLOOD MET
THOD
Description
n of catchmen
nt
River detaiil
Calculated
d by Date
Physiical characteristics
Size of catcchment (A) km
m² Time of 00,385
0,87L2
Longest waatercourse (L) km
m concentration TC hours
Average sloope (Sav) m//m (TC) 1000Sav
SDF basin# Time of
o concentratio
on, t (= 60TC) minutes
2-year returrn period rainffall (M) mmm Days of
o thunder per year (R) days/year
TR1022 n-day rainfa
all data
Weather Seervice station Mean annual precipitation (MA AP) mm
Weather Seervice station no.
n Cooordinates &
Return period
p (yearss)
Duration (d
days)
2 5 10 20 550 100 200
1 day
2 days
3 days
7 days
Rainfall
Return perriod (years), T 2 5 10 20 550 100 200
Point precippitation depth (mm), Pt,T
Area reducttion factor (%
%), ARF
90000 12800lnA 99830lnt
0,4
Average inttensity (mm/hhour), IT
(= Pt,T x AR
RF / TC)
Run n-off coefficients
Calibration factors C2 (2-year retu
urn period) (%
%) C100 (100-year return pperiod) (%)
Return perriod (years) 2 5 10 20 550 100 200
Return period factors (YT) 0 0,84 1,28 1,64 2, 05 2,33 2,58
Run-off coeefficient (CT),
C Y C C
C T 2 T 100 2
100 2,33 100 100
CT IT A
Peak flow ((m³/s), Q T
3,6
L
LEGEND TAABLE
Standaard Design Flo
ood method
ID Referencce ID
D R
Reference
Figure 3.211 Table 3C.12
Table 3B.11
Tablee 3C.12
Criteria Calculaation method
Modified Heershfield equaation
TC < 6 houurs
Pt, T 1,130,441 0,64lnT 0,11 0,27llnt 0,79M 0,69 R 0,20
Linear
L interpoolation betweeen calculated modified
m Hersshfield
6 hours ≤ TC < 24 hourss
equation poinnt rainfall and 1-day point raainfall from T
TR102
A-268 Appendices
SC S-SA METH
HOD
C - Page 1/2
A-269 Appendices
SC
CS-SA MET
THOD
LE
EGEND TABL LE
SC
CS-SA methood
ID Reference
Figure 3.122 or Figure
3E.1 to Figgure 3E.6
Figure 3.255 or 3.26
Table 3E.1 or 3E.2
Table 3E.33
Table 3E.33
Table 3C.88
Table 3C.99
Table 3C.110
Table 3C.111
Table 3C.112
Table
T 3C.8
Adjustmentt of Curve Nu
umbers
1100
CN f
C
Median Coondition Methhod 00
110 ΔS
CN II 25,4
C II
CN
Wet/saturaated Conditioons CN w
0,4036 0,0059CN II
Table 3C.9
254400
Potential m
maximum soill water retentiion S 254
CNN
T
Table 3C.10
(P Ia) 2
Stormflow
w depth Q for P > Ia
P Ia S
Table 3C.11
QA
Stormflow
w volume V=
1000
Tablee 3C.12
0,2083AQ
0
Peak disccharge estimaation qp =
1,83L
C - Page 2/2
A-270 Appendices
EMP
PIRICAL METHODS
M
Description
n of catchmen
nt
River detaiil
Calculated
d by Date
Physiical characteristics
Size of catcchment (A) km² Veld tyype
Longest waatercourse (L) km Catchmment parameteer (C) Α S
C
Length to catchment centtroid (LC) km with reegard to reactio
on time L LC
Average sloope (Sav) m/m
m Kovácss region
Mean annuaal precipitatioon (P) mm
Return perriod (years), T 10 20 50 100
Constant vaalue for KT
Peak flow ((m³/s), QT based
Q T 0,0377K
0 0,66 0,2
T PA C
on Midgleyy & Pitman
Peak flow ((m³/s), QRMF based
b on Kováács
Return perriod (years), T 50 100 200
QT/QRMF rattios
Peak flow ((m³/s), based on
o QT/QRMF raatios
LEGEN ND TABLE
Empiriccal methods
ID Refereence ID Referencee ID
D R
Reference
Figure 3.5 or SA
Fiigure 3.26 Tablee 3C.14
Weather Seervices
Figure 3.155 Taable 3C.13 Tablee 3D.1 or 3D.2
2
Table 3C.133
Con stant values of
o KT
Return Veld type (Figure
( 3.15)
period T 2 4& 6
1 All
A 3 5 All 7 8 9
in yearss Winter
yeear
5A Winter
year
10 0,17 0,42 0,83 0,29 0,59 0,5
59 0,33 0,67 0,67 0,42
2 0,50
20 0,23 0,52 1,04 0,40 0,68 0,8
80 0,45 0,91 0,91 0,57
7 0,68
50 0,32 0,68 1,36 0,55 0,95 1,11 0,63 1,26 1,26 0,79
9 0,95
100 0,40 0,80 1,60 0,70 1,20 1,4
40 0,80 1,60 1,60 1,00
0 1,20
Tabble 3C.14
RMF reegion classificcation in soutthern Africa
Kovács K* N
Number of Trannsition zone Flood zzone
region floods # Area rangee QRMMF Arrea range QRMF
(km²) (m³/s) (km²) (m³/s)
K1 2,8 6 1 – 500 30A0,262
0
500
0 – 500 000 1,74A0,72
0,265
0
K2 3,4 12 1 – 300 50A 300
0 – 500 000 5,25A0,66
0,34
K3 4,0 26 1 – 300 70AA 300
0 – 300 000 15,9A0,60
K4 4,6 55 1 – 100 100AA0,38 100
0 – 100 000 47,9A0,54
0,50
K5 5,0 155 1 – 100 100AA 100
0 – 100 000 100A0,50
0,56
K6 5,2 61 1 – 100 100AA 100 – 30 000 145A0,48
K7 5,4 34 1 – 100 100AA0,62 100 – 20 000 209A0,46
0,68
K8 5,6 25 1 – 100 100AA 100 – 10 000 302A0,44
Notes:
* RMF
MF K value as used
u in Equation 3.32
# Recoorded flood data
d are refleccted in the DW
WAF report TR105
T – Maximum flood peeak discharges
in SSouth Africa: An
A empirical approach
a
E -Page 1/1
A-271 Appendices
APPENDIX 3D
QT/QRMF RATIOS FOR DIFFERENT CATCHMENT AREAS
A-272 Appendices
Table 3D.1: QT/QRMF ratios for different catchment areas in South Africa, Lesotho and
Swaziland (3.13)
Return Effective catchment area - Ae (km2)
Region period KT
(years) 10* 30* 100 300 1 000 3 000 10 000 30 000 100 000 300 000
50 5,06 0,537 0,508 0,474 0,503 0,537 0,570 0,607
K8
100 5,25 0,668 0,645 0,617 0,640 0,668 0,695 0,724
(5,6)
200 5,41 0,803 0,788 0,769 0,784 0,803 0,821 0,838
50 4,70 0,447 0,416 0,380 0,411 0,447 0,482 0,523
K7
100 4,89 0,556 0,525 0,492 0,523 0,556 0,588 0,623
(5,4)
200 5,04 0,661 0,635 0,607 0,633 0,661 0,687 0,716
50 4,50 0,447 0,416 0,380 0,411 0,447 0,482 0,526 0,566
K6
100 4,69 0,556 0,528 0,494 0,524 0,556 0,588 0,626 0,660
(5,2)
200 4,86 0,676 0,650 0,624 0,650 0,676 0,701 0,733 0,758
K5 50 4,30 0,447 0,416 0,380 0,411 0,447 0,482 0,525 0,567 0,617
(5 - except in 100 4,48 0,550 0,521 0,488 0,517 0,550 0,582 0,619 0,657 0,699
SW Cape) 200 4,64 0,661 0,636 0,608 0,633 0,661 0,687 0,718 0,748 0,780
K5 50 4,45 0,531 0,502 0,468 0,497 0,531 0,564
(5 - G, H in 100 4,63 0,654 0,629 0,600 0,625 0,654 0,680
SW Cape) 200 4,78 0,777 0,758 0,738 0,757 0,777 0,795
50 3,84 0,416 0,385 0,350 0,381 0,416 0,453 0,496 0,541 0,591
K4
100 4,04 0,524 0,495 0,462 0,491 0,524 0,558 0,597 0,636 0,679
(4,6)
200 4,20 0,629 0,603 0,576 0,602 0,629 0,660 0,692 0,724 0,758
50 3,26 0,426 0,426 0,426 0,390 0,426 0,463 0,506 0,548 0,602 0,651
K3
100 3,50 0,562 0,562 0,562 0,529 0,562 0,595 0,631 0,666 0,710 0,749
(4)
200 3,68 0,692 0,692 0,692 0,665 0,692 0,718 0,745 0,771 0,804 0,831
50 2,40 0,317 0,317 0,317 0,281 0,317 0,353 0,398 0,444 0,500 0,560
K2
100 2,66 0,428 0,428 0,428 0,391 0,428 0,463 0,506 0,549 0,598 0,651
(3,4)**
200 2,91 0,570 0,570 0,570 0,536 0,570 0,600 0,638 0,672 0,710 0,753
Note: * Estimated ratios
** Ratios of this region may also be used in region K1 (2,8)
Table 3D.2: QT/QRMF ratios for different catchment areas in Namibia and Zimbabwe (3.13)
Return Effective catchment area - Ae (km2)
Region period KT
(years) 10* 30* 100 300 1 000 3 000 10 000 30 000 100 000 300 000
Namibia
K5 50 4,50 0,562 0,534 0,501 0,529 0,562 0,594 0,631
(5) 100 4,70 0,708 0,686 0,661 0,683 0,708 0,732 0,759
200 4,85 0,841 0,828 0,813 0,826 0,841 0,855 0,871
K4 50 4,14 0,589 0,561 0,530 0,558 0,589 0,620 0,654 0,690 0,727
(4,6) 100 4,34 0,741 0,721 0,699 0,719 0,741 0,763 0,787 0,811 0,835
200 4,48 0,871 0,860 0,848 0,860 0,871 0,883 0,895 0,909 0,920
K3 50 3,50 0,562 0,562 0,562 0,529 0,562 0,595 0,631 0,666 0,710
(4) 100 3,66 0,676 0,676 0,676 0,648 0,676 0,703 0,731 0,759 0,793
200 3,77 0,767 0,767 0,767 0,746 0,767 0,788 0,809 0,829 0,856
K2 50 2,88 0,550 0,550 0,550 0,517 0,550 0,585 0,619 0,656 0,696
(3,4) 100 3,01 0,639 0,639 0,639 0,610 0,639 0,669 0,698 0,729 0,762
200 3,13 0,733 0,733 0,733 0,711 0,733 0,758 0,779 0,803 0,828
Zimbabwe
K6 50 4,65 0,531 0,502 0,468 0,497 0,531 0,564 0,603 0,640
(5,2)** 100 4,86 0,676 0,652 0,625 0,649 0,676 0,702 0,731 0,759
200 5,03 0,822 0,807 0,791 0,806 0,822 0,838 0,855 0,871
Note: * Estimated ratios
** In region K5 use the same ratios as those applicable to South Africa
A-273 Appendices
APPENDIX 3E
SCS-SA ADDITIONAL INFORMATION
A-274 Appendices
Table 3E.1: Example of classification of soils in southern Africa into hydrological soil
groups by soil form, family and textural class (taxonomic classification)
Typical Typical
Soil SCS SCS
Code Soil Family Textural Soil Form Code Soil Family Textural
Form Grouping Grouping
Class Class
ADDO Ad 1111 Glenconnor LmSa A/B ARCADIA Ar 1100 Lonehill Cl C/D
B Ad 1111 Glenconnor SaLm B C/D Ar 1100 Lonehill Cl C/D
Ad 1111 Glenconnor SaClLm B Ar 1200 Rustenburg Cl C/D
Ad 1111 Glenconnor SaCl B/C Ar 1200 Rustenburg Cl C/D
Ad 1112 Dalby LmSa A/B Ar 2100 Minerva Cl C/D
Ad 1112 Dalby SaLm B Ar 2100 Minerva Cl C/D
Ad 1112 Dalby SaClLm B Ar 2200 Diepsloot Cl C/D
Ad 1112 Dalby SaCl B/C Ar 2200 Diepsloot Cl C/D
Ad 1121 Centlivres LmSa B Ar 3100 Bospoort Cl C/D
Ad 1121 Centlivres SaLm B/C Ar 3100 Bospoort Cl C/D
Ad 1121 Centlivres SaClLm B/C Ar 3200 Deercroft Cl C/D
Ad 1121 Centlivres SaCl C Ar 3200 Deercroft Cl C/D
Ad 1122 Kentvale LmSa B ASKHAM Ak 1000 Aroab LmSa A/B
Ad 1122 Kentvale SaLm B/C B Ak 1000 Aroab SaLm B
Ad 1122 Kentvale SaClLm B/C Ak 1000 Aroab SaClLm B
Ad 1122 Kentvale SaCl C Ak 1000 Aroab SaCl B/C
Ad 1211 Spekboom LmSa A/B Ak 2000 Noenieput LmSa B
Ad 1211 Spekboom SaLm B Ak 2000 Noenieput SaLm B/C
Ad 1211 Spekboom SaClLm B Ak 2000 Noenieput SaClLm B/C
Ad 1211 Spekboom SaCl B/C Ak 2000 Noenieput SaCl B/C
Ad 1212 Gorah LmSa A/B AUGRABIE Ag 1110 Hefnaar LmSa A/B
Ad 1212 Gorah SaLm B B Ag 1110 Hefnaar SaLm B
Ad 1212 Gorah SaClLm B Ag 1110 Hefnaar SaClLm B
Ad 1212 Gorah SaCl B/C Ag 1110 Hefnaar SaCl B/C
Ad 1221 Walkraal LmSa B Ag 1120 Giyani LmSa B
ADDO Ad 1221 Walkraal SaClLm B/C Ag 1120 Giyani SaLm B/C
B Ad 1221 Walkraal SaCl C Ag 1120 Giyani SaClLm B/C
Ad 1222 Sylvania LmSa B Ag 1120 Giyani SaCl C
Ad 1222 Sylvania SaLm B/C Ag 1210 Khubus LmSa A/B
Ad 1222 Sylvania SaClLm B/C Ag 1210 Khubus SaLm B
Ad 1222 Sylvania SaCl C Ag 1210 Khubus SaClLm B
Ad 2111 Maurmond LmSa A/B Ag 1210 Khubus SaCl B/C
Ad 2111 Maurmond SaLm B Ag 1220 Shilowa LmSa B
Ad 2111 Maurmond SaClLm B LEGEND
Ad 2111 Maurmond SaCl B/C A - low runoff potential
Ad 2112 Airedale LmSa A/B B - moderately low potential
Ad 2112 Airedale SaLm B C - moderately high potential
Ad 2112 Airedale SaClLm B D - high runoff potential
Ad 2112 Airedale SaCl B/C Sa - sand
Ad 2121 Felsenheim LmSa B Cl - clay
Ad 2121 Felsenheim SaLm B/C Lm - loam
Ad 2121 Felsenheim SaClLm B/C
Ad 2121 Felsenheim SaCl C
Ad 2122 Longhill LmSa B
Ad 2122 Longhill SaLm B/C
Ad 2122 Longhill SaClLm B/C
Ad 2122 Longhill SaCl C
Ad 2211 Mimosa LmSa A/B
Ad 2211 Mimosa SaLm B
Ad 2211 Mimosa SaClLm B
Ad 2211 Mimosa SaCl B/C
Ad 2212 Peperboom LmSa A/B
Ad 2212 Peperboom SaLm B
Ad 2212 Peperboom SaClLm B
Ad 2212 Peperboom SaCl B/C
Ad 2221 Suttondale LmSa B
Ad 2221 Suttondale SaLm B/C
Ad 2221 Suttondale SaClLm B/C
Ad 2221 Suttondale SaCl C
Ad 2222 Tregaron LmSa B
Ad 2222 Tregaron SaClLm B/C
Ad 2222 Tregaron SaCl C
A-275 Appendices
Table 3E.2: Example of classification of soils in southern Africa into hydrological
soil groups by soil form and series (binomial classification)
Typical Typical
SCS SCS
Soil Form Code Soil Series Textural Soil Form Code Soil Series Textural
Grouping Grouping
Class Class
ARCADIA Ar 40 Arcadia Cl C/D BAINSVLEI Bv 14 Wykeham SaLm A/B
C/D Ar 11 Bloukrans Cl C/D BONHEIM Bo 41 Bonheim LmSa C/D
Ar 21 Clerkness Cl C/D C Bo 20 Bushman SaClLm C
Ar 41 Eenzaam Cl C/D Bo 30 Dumasi SaClLm C
Ar 20 Gelykvlakte Cl C/D Bo 31 Glengazi SaCl C/D
Ar 10 Mngazi Cl C/D Bo 10 Kiora SaClLm C
Ar 32 Nagana Cl C/D Bo 21 Rasheni SaCl C/D
Ar 12 Noukloof Cl C/D Bo 11 Stanger SaCl C/D
Ar 31 Rooidraai Cl C/D Bo 40 Weenen SaClLm C
Ar 30 Rydalvale Cl C/D CARTREF Cf 10 Amabele LmSa B/C
Ar 42 Wanstead Cl C/D C Cf 12 Arrochar SaClLm C
Ar 22 Zwaarkrygen Cl C/D Cf 13 Byrne SaCl C/D
AVALON Av 13 Ashton SaLm A/B Cf 21 Cartref SaLm C
B Av 26 Avalon SaClLm B Cf 22 Cranbrook SaClLm C
Av 12 Banchory Sa A Cf 30 Grovedale Sa B/C
Av 27 Bergville SaCl B/C Cf 31 Kusasa SaLm B/C
Av 37 Bezuidenhout SaCl C Cf 32 Noodhulp SaClLm C
Av 33 Bleeksand SaLm B/C Cf 11 Rutherglen SaLm C
Av 34 Heidelberg SaLm B/C Cf 20 Waterridge LmSa B/C
Av 20 Hobeni LmSa A/B CHAMPAGNE Ch 11 Champagne SaLm D
Av 14 Kanhym SaLm A/B D Ch 21 Ivanhoe SaClLm D
Av 24 Leksand SaLm B Ch 10 Mposa SaLm D
Av 10 Mastaba LmSa A Ch 20 Stratford SaClLm D
AVALON Av 32 Middelpos Sa B CLOVELLY Cv 33 Annandale SaLm B
B Av 31 Mooiveld LmSa B A/B Cv 18 Balgowan Cl B
Av 25 Newcastle SaLm A/B Cv 40 Bleskop LmSa A
Av 17 Normandien SaCl B Cv 36 Blinkklip SaClLm B
Av 22 Rossdale Sa A/B Cv 17 Clovelly SaCl B
Av 16 Ruston SaClLm B Cv 28 Clydebank Cl B
Av 36 Soetmelk SaClLm B/C Cv 35 Denhere SaLm A/B
Av 21 Uithoek LmSa A/B Cv 46 Dudfield SaClLm A/B
Av 30 Viljoenskroon LmSa B Cv 11 Geelhout LmSa A
Av 23 Villiers SaLm B Cv 25 Gutu SaLm A
Av 11 Welverdien LmSa A Cv 47 Klippan SaCl B
Av 35 Windmeul SaLm B Cv 38 Klipputs Cl B/C
Av 15 Wolweberg SaLm A Cv 10 Lismore LmSa A
BAINSVLEI Bv 23 Ashkelon SaLm A/B LEGEND
A/B Bv 36 Bainsvlei SaClLm B A - low stormflow potential
Bv 12 Camelot Sa A B - moderately low potential
Bv 20 Chelsea LmSa A C - moderately high potential
Bv 30 Delwery LmSa A/B D - high stormflow potential
Bv 13 Dunkeld SaLm A/B Sa - sand
Bv 16 Elysium SaClLm A/B Cl - clay
Bv 10 Hlatini LmSa A Lm - loam
Bv 34 Kareekuil SaLm B
Bv 31 Kingston LmSa A/B
Bv 26 Lonetree SaClLm A/B
Bv 25 Maanhaar SaLm A
Bv 11 Makong LmSa A
Bv 27 Metz SaCl B
Bv 22 Oosterbeek Sa A
Bv 37 Ottosdal SaCl B/C
Bv 24 Redhill SaLm A/B
Bv 32 Trekboer Sa A/B
Bv 15 Tygerkloof SaLm A
Bv 33 Vermaas SaLm B
Bv 21 Vungama LmSa A
Bv 35 Wedgewood SaLm A/B
Bv 17 Wilgenhof SaCl B
A-276 Appendices
Table 3E.3 Initial Curve Numbers for selected land cover and treatment classes,
stormflow potentials and hydrological soil groups (various sources)
Land Cover Land Treatment/ Practice/Description Stormflow Hydrological Soil Group
Class Potential A A/B B B/C C C/D D
1 = Straight row 77 82 86 89 91 93 94
Fallow 2 = Straight row + conservation tillage High 75 80 84 87 89 91 92
3 = Straight row + conservation tillage Low 74 79 83 85 87 89 90
1 = Straight row High 72 77 81 85 88 90 91
2 = Straight row Low 67 73 78 82 85 87 89
3 = Straight row + conservation tillage High 71 75 79 83 86 88 89
4 = Straight row + conservation tillage Low 64 70 75 79 82 84 85
5 = Planted on contour High 70 75 79 82 84 86 88
Row Crops 6 = Planted on contour Low 65 69 75 79 82 84 86
7 = Planted on contour + conservation tillage High 69 74 78 81 83 85 87
8 = Planted on contour + conservation tillage Low 64 70 74 78 80 82 84
9 = Conservation structures High 66 70 74 77 80 82 82
10 = Conservation structures Low 62 67 71 75 78 80 81
11 = Conservation structures + conservation tillage High 65 70 73 76 79 80 81
12 = Conservation structures + conservation tillage Low 61 66 70 73 76 78 79
1 = Straight row Low 45 56 66 72 77 80 83
Garden Crops
2 = Straight row High 68 71 75 79 81 83 84
1 = Straight row High 65 71 76 80 84 86 88
2 = Straight row Low 63 69 75 79 83 85 87
3 = Straight row + conservation tillage High 64 70 74 78 82 84 86
4 = Straight row + conservation tillage Low 60 67 72 76 80 82 84
5 = Planted on contour High 63 69 74 79 82 84 85
6 = Planted on contour Low 61 67 73 78 81 83 84
Small Grain
7 = Planted on contour + conservation tillage High 62 68 73 77 81 83 84
8 = Planted on contour + conservation tillage Low 60 66 72 76 79 81 82
9 = Planted on contour - winter rainfall region Low 63 66 70 75 78 80 81
10 = Conservation structures High 61 67 72 76 79 81 82
11 = Conservation structures Low 59 65 70 75 78 80 81
12 = Conservation structures + conservation tillage High 60 67 71 75 78 80 81
13 = Conservation structures + conservation tillage Low 58 64 69 73 76 78 79
1 = Straight Row High 66 72 77 81 85 87 89
Close Seeded 2 = Straight Row Low 58 65 72 75 81 84 85
Legumes or 3 = Planted on contour High 64 70 75 80 83 84 85
Rotational 4 = Planted on contour Low 55 63 69 74 78 81 83
Meadow 5 = Conservation structures High 63 68 73 77 80 82 83
6 = Conservation structures Low 51 60 67 72 76 78 80
1 = Straight row: trash burnt 43 55 65 72 77 80 82
2 = Straight row: trash mulch 45 56 66 72 77 80 83
3 = Straight row: limited cover 67 73 78 82 85 87 89
Sugarcane 4 = Straight row: partial cover 49 60 69 73 79 82 84
5 = Straight row: complete cover 39 50 61 68 74 78 80
6 = Conservation structures: limited cover 65 70 75 79 82 84 86
7 = Conservation structures: partial cover 25 46 59 67 75 80 83
8 = Conservation structures: complete cover 6 14 35 59 70 75 79
A-277 Appendices
Table 3E.3 Initial Curve Numbers for selected land cover and treatment classes,
stormflow potentials and hydrological soil groups (various sources) (continued)
Land Cover Land Treatment/ Practice/Description Stormflow Hydrological Soil Group
Class Potential A A/B B B/C C C/D D
1 = Veld/pasture in poor condition High 68 74 79 83 86 88 89
2 = Veld/pasture in fair condition Moderate 49 61 69 75 79 82 84
Veld (range) 3 = Veld/pasture in good condition Low 39 51 61 68 74 78 80
and Pasture 4 = Pasture planted on contour High 47 57 67 75 81 85 88
5 = Pasture planted on contour Moderate 25 46 59 67 75 80 83
6 = Pasture planted on contour Low 6 14 35 59 70 75 79
A-278 Appendices
Figu
ure 3E.1: One-day desig
gn rainfall d
distribution over southeern Africa foor 2 year return
period
Figu
ure 3E.2: One-day desig
gn rainfall d
distribution over southeern Africa foor 5 year return
period
A-279 Appendices
Figuree 3E.3: One--day design rainfall disttribution over southernn Africa
for 10 year returnn period
A-280 Appendices
Figuree 3E.5: One--day design rainfall disttribution over southernn Africa
for 50 year returnn period
A-281 Appendices