Routing MODULE - ROUTING METHODS Routing is the process of find the distribution of flow rate and depth in space and time along a river or storm sewer. Routing is also called Flow routing or flood routing. Two different types of routing models Hydraulic routing Hydrologic routing Hydrologlic Routing Involves the simultaneous solution of continuous equation and a storage volume equation written as a function of inflow and outflow flowrates. The basic equation used in hydrologic routing is called the continuity equation. ds = I(t) O(t) dt This is often applied to storage reservoirs and storm-water detention basins. Two types of methods are used: Modified Puls method (storage is a function of outflow rate, for example: storage basin) Muskingum Method (Storage is a function of both inflow and outflow, for example: channel routing) Modified Puls method Over a finite time interval between t and t+dt, the continuity equation can be written as: S S 1 = I + I 1 Δt O + O 1 Rewriting the above equation for convenience ( ) + S 1 I 1 + I Δt O 1 = S Δt + O The subscripts 1 and refer to the values at times t and t+dt. What are the known and unknown parameters in the above equation? All inflows are known (usually inflow hydrographs are given) For any reservoir, we can figure out the storage vs stage (Z) relationship Also outflow for given stage can be determined by equations. Steps involved to route the inflow hydrograph through a reservoir 1
S 1 1. Substitute the known values of I 1, I and Δt O 1. This gives the values of S Δt + O. From the reservoir and outlet structure characteristics, determine the discharge, O S corresponding to the calculated value of Δt + O. S 3. Subtract O from Δt + O S 1 to yield Δt O 1 at the end of this time step. 4. Repeat the steps 1 to 3 until the entire outflow hydrograph is calculated. Problem 1 A stormwater-detention basin is estimated to the have the following storage characteristics Stage (m) Storage (m 3 ) 5.0 0 5.5 694 6.0 155 6.5 507 7.0 365 7.5 4973 8.0 6484 The discharge weir from the detention basin has a crest elevation of 5.5m and the weir discharge is given by Q=1.83H 5/ The catchment runoff hydrograph is given by Time (min) 0 30 60 90 10 150 180 10 40 70 300 330 360 390 Runoff (m 3 /s) 0.4 5.6 3.4.8.4. 1.8 1.5 1. 1.0 0.56 0.34 0
ELEVATION-STORAGE-DISCHARGE RELATIONSHIP Stage (m) Storage (m 3 ) 5.0 0 5.5 694 6.0 155 6.5 507 7.0 365 7.5 4973 8.0 6484 3
Time (min) Runoff (m 3 /s) 0 0 30.4 60 5.6 90 3.4 10.8 150.4 180. 10 1.8 40 1.5 70 1. 300 1 330 0.56 360 0.34 390 0 4
Muskingum method Developed for Muskingum river basin in eastern Ohio by McCarthy in 1938. It is primarily used for routing water in drainage channels including rivers and streams. The routing expression for this method is given as: where C 1, C and C 3 are given by O j+1 = C 1 I j+1 + C I j + C 3 O j Δt KX C 1 = K(1 X) + Δt Δt + KX C = K(1 X) + Δt K(1 X) Δt C 3 = K(1 X) + Δt It should be apparent that C 1 +C +C 3 =1 The time step should be taken such that it lies between KX and K(1-X) to ensure that you do not get negative flows. A convenient value for time step can be between K/3 and K. Problem The flow hydrograph at a channel section is given by Time (min) 0 30 60 90 10 150 180 10 40 70 300 330 360 390 Runoff (m 3 /s) 10 10 5 45 31.3 7.5 5 3.8 1.3 19.4 17.5 16.3 13.5 1.1 Time 40 450 480 (min) Runoff (m 3 /s) 10 10 10 Use the Muskingum routing method to estimate the hydrograph 100 m downstream from the channel section. Assume that X=0. and K=40 minutes. 5
Time (min) Runoff (m 3 /s) Outflow (m 3 /s) 0 10 30 10 60 5 90 45 10 31.3 150 7.5 180 5 10 3.8 40 1.3 70 19.4 300 17.5 330 16.3 360 13.5 390 1.1 40 10 450 10 480 10 6
Problem 3 The flow hydrograph at a channel section is given below. Use Muskingum routing to find the outflow. The Muskingum parameters are as follows: X=0.15, K =.3 minutes Time (hr) I (m 3 /s) Outflow (O) (m 3 /s) 1 93 137 3 08 4 30 5 44 6 546 7 630 8 678 9 691 10 675 11 634 1 571 13 477 14 390 15 39 16 47 17 184 18 134 19 108 0 90 7
Problem 4 The flow hydrograph at a channel section is given below. Use Muskingum routing to find the outflow. The Muskingum parameters are as follows: X=0.3, K = 35 minutes Time (hr) I (m 3 /s) 0 0 30 1.5 60.1 90 15.4 10 13.6 150 1.4 180 11.7 10 10.8 40 9.9 70 8.4 300 8.1 330 7.5 360 4. 390 0 8