Smada PDF
Smada PDF
The methods used in SMADA are outlined in Hydrology: Water Quantity and Water
Quality Control 2nd Edition, published by John Wiley and Sons.
PROGRAM OVERVIEW/OUTLINE
NOTE: The following information is a partial listing of the SMADA HELP files. Some
of the text, figures and equations have been omitted. Please use the SMADA ON-LINE
HELP for complete information.
I. WATERSHEDS
II. RAINFALL
III. HYDROGRAPHS
IV. PONDS
V. UTILITY PROGRAMS
A. MATRIX CALCULATOR
B. TIME OF CONCENTRATION CALCULATOR
C. POLLUTANT ANALYSIS
D. STATISTICAL REGRESSION ANALYSIS
E. STATISTICAL DISTRIBUTION ANALYSIS
F. STORM SEWER DESIGN
VI. ADDITIONAL PROGRAMS INSTALLED WITH SMADA (NOT PART OF SMADA MENU)
A. LOOPS (Simplified Hardy Cross Calculator)
B. CIRCULAR PIPE CALCULATOR
CREATING A WATERSHED
The two infiltration methods available to estimate the potential infiltration are:
1. Horton Method
More information about the use of the watershed parameters is contained in the
Water Budget section of the help system.
Watershed Files can be saved and opened. The files should be assigned a legal DOS 8
character name. A default extension of (. SHD) will automatically be placed on the
watershed. Watersheds can only be opened and saved by SMADA. Opening a watershed
file will automatically replace the existing watershed file in memory. If you wish
to delete old watershed files you may do so using File Manager.
Maximum Infiltration
Horton Limiting Infiltration Rate
Horton Initial Infiltration Rate
Horton Depletion Coefficient
Maximum Infiltration
Curve Number
Initial Abstraction Factor
Other data on soil types, hydrologic classification, and curve numbers are
available in appendix F(USDA-SCS, 1986). Judgment on the type of land use has to
be exercised and a conservative estimate (higher CN values) is prudent when doubt
exists.
Using the SCS-CN procedure, rainfall excess calculations are a function of rainfall
volume and curve number. Assuming that storage at any time is proportional to
maximum storage and rainfall excess is proportional to precipitation volume, (USDA-
SCS, 1986).
When using the SCS-CN method with composite curve numbers, the user does not have
to specify the three land categories. The percentage of impervious area should be
specified as 0.0. In place of the Pervious Curve Number the user should specify the
Composite Curve Number.
A rainfall file consists of a rainfall volume (in inches) for a series of time
increments. The maximum number of time increments allowed is 960. To determine the
number of time increments a rainfall will have divide the duration of the rainfall
by the time increment. For example a 1 hour long storm with 10 minute time
increments will have 6 steps (60 / 10).
If you choose to enter a incremental rainfall by hand, you must choose a duration
and a time step for the rainfall. If you choose to create a USER rainfall, you do
not need to enter a total rainfall volume. When you hit the OK button you will be
placed in the Rainfall Edit Dialog Box. There may be more edit spaces available
than the number of steps in your rainfall. Any extra values you enter outside of
the duration of the rainfall will be ignored. The total rainfall will be
calculated after the rainfall values have been entered.
SMADA has currently many types of dimensionless curves which can be entered; some
are SCS type II, SCS type II Florida modified, and SCS type III. (There may be
others available at the time this package arrives to you.) To use these curves you
must enter a duration, time step, and a total rainfall. Click the button
pertaining to the type of curve which you would like to use. A rainfall will be
created and you will automatically enter the rainfall View dialog box. You may
enter the rainfall edit routine to edit this rainfall. The edit routine
automatically truncates the rainfall to three decimal places and it may cause you
to lose the accuracy of the original rainfall create (8 places).
There should be no spaces and there should be a number on every line of the file.
From File/Import select rainfall. You will be asked to prompt the time step for
the rainfall data. The file will then be imported and the rainfall will be
assigned to the time steps specified by the rainfall interval.
The rainfall files saved to disk will automatically be assigned the default
extension (.RNF). Opening a rainfall file will automatically replace the rainfall
which is currently in memory. The program can then later retrieve these files.
SCS Method
Santa Barbara Urban Hydrograph Method
Unit Hydrograph Method
Selecting the desired method from the menu will generate the hydrograph using the
current rainfall and watershed.
Using these relationships and the summing of flows for each discrete time step,
SMADA creates a unit response function, similar to a unit hydrograph. An
individual response function is then generated for each of the steps of the
Instantaneous Hydrograph. These response functions are summed to produce an output
hydrograph.
All response functions and routings are performed using the default rainfall time
step. It is advisable that this time step be at least half the time of
concentration of the watershed.
Q2 = Q1*K[R1 + R2 - 2Q1]
where:
All routing involving this equation is done with the default time step used in the
rainfall entry. The instantaneous hydrograph is developed using the Infiltration
and abstraction information.
POND GEOMETRY
Edit Basic Pond Information - The basic information for pond geometry is entered
using this routine. The following diagram shows the relationship between the
different entries:
At the heart of all pond routings is the stage - storage discharge relationship.
This relationship relates the pond stage to the storage in the pond. The
relationship between these two parameters is dependent upon the pond geometry. One
method for the determination of this stage-storage relationship is using the stage-
area relationship. The area referred to is the surface area of the pond at any
given stage.
The Flow in (Qin) is received from a hydrograph or a pond. The Flow Out (Qout) is
calculated knowing the discharge which is found from the stage which is found from
the storage. The discharge relationship can be calculated from weir information or
input directly by hand.
Simple Hydrograph to Pond routing can be performed using the Quick Route option in
the Routing Menu. This option will route the current hydrograph to the current
pond. The results can then be displayed with either the Print Pond Routing or Plot
Pond Routing.
A. MATRIX CALCULATOR (EZMAT Version 1.2) Information Not Available See Program
The following formulas and equations are available with the TCCALC:
Izzard's Formula
Kerby's Equation
Kirpich's Equation
Kinematic Equation
Bransby Williams Equation
Federal Aviation Agency Equation
TR-55 Worksheet
Using Regress is a simple matter of getting the data into the program, selecting
the type of curve you wish to fit and then viewing the output. You can input the X-
Y data directly into the spreadsheet section on the regression form. You can
import the data from an ASCII file which has numbers in an X Y format. You can
paste two columns from a spreadsheet directly onto the spreadsheet section on the
regression form. All analyses will fit the data to a curve. The only difference
between the different analyses is the type of output which is given.
Iterative Search - In equations with a third parameter "n" allows you to try
numerous values for "n" and view the results.
Plot Regression - Performs the regression analysis and gives a plot of the
fitted data. This is a real handy option to get a general idea of the goodness
of fit.
WHAT IS A DISTRIBUTION?
We run DISTRIB and we find our data fits a 2 parameter log normal distribution very
nicely. Odds are any storms in the future will also fit nicely. We can predict
storms which occur every 200 storms, every 1000 storms, etc.. We have fit our data
to a distribution and we can make predictions based on that distribution. You may
be familiar with the most common distribution - the normal distribution. You know,
the one that is shaped like a typical bell curve. Well, surprise! there are many
distributions which fit many different types of data. This program will help you
fit your data to the distributions available in the program.
TYPES OF DISTRIBUTIONS
Normal
Log Normal
3 Parameter Log Normal
Pearson
Log Pearson
Gumbel
METHODS OF ANALYSIS
Maximum Likelihood
Method of Moments
IMPORTING DATA
Data can be imported in both DISTRIB 1.0 and DISTRIB 2.0. The data must be in ASCII
format with one number per line. There should be no extra spaces. For example if
you wish to import a set of data with 9 pieces of data then the file would look
like;
Data can be pasted from a spreadsheet directly into DISTRIB. To do this you must
have the data in columnar format in the spreadsheet. In your spreadsheet highlight
your data in columnar format and copy it to the clipboard. Switch to DISTRIB [Alt
Tab or Ctrl-Esc]. Select Edit, Paste from the DISTRIB menu. The data will be
sorted and displayed in the spreadsheet control in the distribution window. You
may edit the data within this display. Do not worry about blanks in the data, as
DISTRIB will remove these automatically. Click Ok from the Paste data dialog. You
may now analyze the data.
EDITING DATA
Data can be edited directly by changing the data on the spreadsheet in the
distribution window. When the distribution is selected the data will be sorted.
This data can also be saved to disk.
The distribution Analysis can be plotted. The plot will contain the actual data
and the prediction. The plot can be sent to the Windows clipboard using the copy
button.
DISTRIBUTION ANALYSIS
FID CURVES
The major design concern for a sewer system is the water flow rate and surcharge
condition which must be handled. Surcharging occurs when the pipe design cannot
pass the flow rate and water builds up before entering the conveyance system. All
water resulting from rainfall must be transported from one place to another either
above ground or below ground in pipes. A key concept in determining this volume
and flow rate of water is event average intensity. This intensity is the average
amount of rainfall which occurs per time period (hour) over the course of the
rainfall event. For example if 2 inches of rainfall occurs during a hour storm
then the intensity will be 4 in/hr.
This brings up another important consideration- the event duration. This duration
is the length of the storm being considered. A typical Frequency-Intensity-
Duration (FID) curve will display statistics for storm events which last from about
8 minute to 24 hours. When reading an FID curve you may determine that (for a given
In the last example both of these rainfalls are for a given return period. The
return period should be thought of as an inverse of an exceedence probability. A
five year return period means that on the average this rainfall will occur once
every five years. We should expect that a 100 year storm has more rainfall than a
2 year storm. Upon observing the curves we will find that this is the case. (Note:
FID curves do not tell us anything about the time variable rainfall distribution.
The rainfall distribution is a tabulation of how much rain falls for each time
period during the rainfall event.)
We would not expect that every part of the country or a state would have the same
probability of rainfall as every other part of the country or a state. This is why
any geographic area is divided into zones which have approximately the same
probability of the same rainfall events. The FID curves for the state of Florida
are included in the appendix. OPSEW supports Florida zones 1-11; Orlando is
located in zone 7, Miami in zone 10, Tallahassee in zone 2, etc.. Curve fits are
available for each of these zones, and these curves fits give reasonably good fit
for design purposes. (Note: The FID curves were developed using daily data with
assumed inter-event dry periods of approximately one day. With this assumption the
drainage system must drain in a 24 hour period or less. This includes time for the
downstream stage to return to the initial design stage.) For the state of Florida
these curve fits are published by the Florida Department of Transportation. The
form of the curve fit is:
where-
When designing a sewer system, the time of concentration is used in place of the
duration for the determination of any given average intensity. The time of
concentration is the maximum amount of time it takes water to travel from any point
in the watershed to a sewer inlet. By using the time of concentration and the FID
curves we can determine the average intensity of rainfall and thus the flow rate
(runoff) reaching the sewer inlet at any of the sewer nodes. The data used for the
calculation of intensity are shown below. These data are the contents of the file
FLORIDA.FID. The format of this file is:
ZONE_NUMBER
RETURN_PERIOD_2 A B C D
RETURN_PERIOD_50 A B C D
Where A, B, C, D are the constants used in the above equation for the given zone
and return period.
In the ROUTINES menu there is a routine which will allow the user to read the FID
curves automatically without performing any design analysis.
FID data for Florida zones 1-10 for 3rd Order Polynomial Fit
Another FID curve method used by OPSEW is the equation: I= A/(DR+ B)C
The file which contains the coefficients for this method is ZONE-DAT.ZN2 The format
for this file is,
NUMBER OF ZONES
NUMBER OF RETURN PERIODS
ZONE
RETURN PERIOD, A, B, C
<- Must reflect the number of return periods
ZONE
RETURN PERIOD, A, B, C
ZONE <- Must reflect the number of zones
etc...
DETERMINATION OF RUNOFF
After the rainfall reaches the ground there are four possible results; infiltrate
into the ground, flow along the ground, pool on the ground or evaporate back into
the air. Of these processes-, infiltration, runoff, surface abstraction, and
evaporation, we will be primarily concerned with the amount that is runoff. The
rate of rainfall which becomes runoff can be expressed as a fraction of the total
volume of rainfall. This fraction is called the Rational Coefficient and is
usually denoted with a capital "C". OPSEW uses the rate of rainfall which becomes
runoff to design sewer systems. The design is based upon the peak flow which can
be expected for a pipe. This peak flow can be determined from the equation:
Q = CiA
where-
C = Rational Coefficient
i = intensity (in/hr)
Business Pavement
Downtown 0.70-0.95 Asphalt/Concrete 0.70-0.95
Neighborhood 0.50-0.70 Brick 0.70-0.85
Residential Lawns, Sandy soils
Single Family 0.30-0.50 Flat, 2% 0.05-0.10
Multiunit detached 0.40-0.60 Average, 2-7%
0.10-0.15
Multiunit attached 0.60-0.75 Steep, > 7%
0.15-0.20
Suburban resident 0.25-0.40
apartment 0.50-0.70
Industrial Lawns, Heavy Soils
Light 0.50-0.80 Flat 2% 0.13-0.17
Heavy 0.60-0.90 Average, 2-7% 0.18-0.22
Steep >7% 0.25-0.35
Parks and Cemeteries 0.10-0.25 Railroad Yard
0.20-0.35
Unimproved 0.10-0.30
OPSEW uses one rational coefficient for each inlet to the system. The area in
concern may, however, contain multiple land uses. To handle this situation we can
calculate an effective rational coefficient. This coefficient is a weighted
average based on area of the rational coefficients of the different areas.
To determine the peak flow for the design we must first determine the time of
concentration for the watershed in question. The time of concentration is the time
it takes water to flow from the farthest point in the watershed to the outlet.
This can be accomplished by using any number of methods outlined in most hydrology
texts (see Wanielista, 1990). OPSEW contains a routine to assist you in
determining a correct choice for the time of concentration. In the Routines menu
the Calculate Time of Concentration option will allow you to calculate a time of
concentration for a watershed. If you do not know all the input parameters for
this routine then fill in only those parameters that you do know. OPSEW will
calculate the time of concentration based upon various equations from the
literature depending upon the input information.
The time of concentration is used to determine the average intensity used in the
rational formula for the sewer design. The duration is set equal to the time of
concentration and this number is used to find the intensity used in the rational
formula. When determining design flows for downstream segments of pipe we use an
effective time of travel to determine the intensity. This effective time of travel
is equal to the time of concentration of the most upstream node + the cumulative
time of travel in each of the upstream pipe sections. A flow at the inlet of the
invert in question is then calculated for this segment and then added to the flow
which is coming to the manhole through the pipe.
Where the concern of Hydrology was to determine the amount and rate of water our
conveyance system would have to convey; we are now interested in the study of
Hydraulics which will allow us to look at the details of the water flowing through
our system. One of the main objectives of our hydraulic design is to select invert
elevations and pipe diameters so to prevent excessive surcharging of our sewer
system. Surcharging can cause flooding of upstream areas if the design is not
proper.
FLOW CONDITIONS
OPSEW has a design approach that deals with the hydraulics of circular pipes.
Circular pipes can three possible flow conditions:
1. Partial Flow
2. Full Flow
3. Surcharged Flow
OPSEW is designed to perform calculations for all three of these flow conditions.
The automatic design routines for OPSEW will size pipes and place inverts so as to
prevent the last of these cases; surcharged flow. To determine velocities and flows
OPSEW uses Manning's formula.
OPSEW will calculate a velocity based on any of the three flow conditions.
OPSEW contains individual routines which will allow for the analysis of all three
of these conditions of flow. The automatic and manual design routines also make
the determination of the flow condition and base all calculations on the flow
condition.
Using the sciences of hydraulics and hydrology we can see that the design of a
sewer system is an iterative process. We will outline the steps of that process.
5. Starting at the most upstream pipe and using the time of concentration for
that watershed determine the flow through that pipe segment.
a. Use time of concentration as storm duration.
b. Determine storm intensity from FID curves.
c. Using "C" and intensity calculate peak flow.
7. Calculate the flow in this pipe segment in the same manner as step 5, but
using the effective time of concentration.
OPSEW considers the headloss in two places; headloss encountered in the length of
pipe and headloss in the manholes, joints, or junctions. The headlosses are used
in the calculation of the hydraulic grade.
The pipe costs are split into two parts- pipe cost and excavation cost.
The calculation of the excavation cost is slightly more complex; first the depth is
calculate. OPSEW calculates the average depth by the following equation:
Once the average depth is calculated this depth is interpolated against the cost
data depths. For example if the average depth is 5.6 ft and the available cost
data has excavation cost per foot for 4 ft. and 6 ft. the cost per foot is found
by interpolating the cost per foot from 4 and 6 ft. Once this cost per foot is
found the total excavation cost is found by multiplying the cost per foot by the
pipe length. The two costs are then added to give total cost.
Loops is a program which allows users to perform Hardy Cross method hydraulic loop
analysis (50 loops, 100 pipes, 80 nodes).
SMADA END