## Assignment 10: River Discharge Function

River discharge function

The volume of water that passes by a point in a river over a certain time increment is the discharge (for example, the Tempe Town Lake can apparently withstand about 40,000 cubic feet per second with the dams up). The discharge can be simply expressed as the product of the velocity times the cross-sectional area of the flow:
Q = v * A, where Q is discharge [L^3/T], v is velocity [L/T], and A is area [L^2]
Many channels are effectively rectangular, so we can express the cross-sectional area as the product of width (w) times depth (d):
A = w * d So how to calculate the velocity?
We use a simple formulation called Manning's equation that relates the flow velocity to the local channel slope (s [dimensionless]), a channel shape parameter called the hydraulic radius [L]--see below, and a roughness parameter n [L^1/6]. NOTE THAT THIS IS FOR ENGLISH UNITS (feet, seconds): The hydraulic radius R is equal to the ratio of the channel area (A) to its wetted perimeter (P):
P = w + 2d
R = A/P
The Manning roughness coefficient (n) is estimated based upon inspection of the channel and its shape. Here is a table of roughnesses depending on the channel boundary type.

Here is a link with similar explanation and a few pictures and hydrographs from 2005 flow along the Salt River: pdf. A basic video explanation of the MATLAB implementation of these equations is here:

For this homework, I would like for you all to build a couple of functions that make some interesting geomorphic calculations. Build a web page for this assignment and present your results and answers to the following set of questions:

### Tasks

1) Develop a well documented Matlab function that will calculate the velocity and the discharge at a point given the channel shape, roughness, and slope parameters. Make sure that you use the fprintf command to write the variables back to the command line and let the user know what is going on.
2) Make a plot of the relationship of velocity versus channel slope for values of slope varying from 0 to 0.2 (about 12 degrees), all other variables constant. How does velocity vary with channel slope?
3) Make a plot of the relationship of discharge versus channel slope for values of slope varying from 0 to 0.2 (about 12 degrees), all other variables constant. How does discharge vary with channel slope?
4) Make a plot of the relationship between velocity and hydraulic radius, all other variables constant. Choose a range of radius values that are for shallow and deep channels. All things being equal, is the velocity faster in a deeper or shallower channel?
5) Make a plot of the relationship between velocity and channel roughness, all other variables constant. Vary the roughness from that of a smooth glass chute to a mountain stream. How does the velocity vary with roughness?
6) Make a plot of the relationship between discharge and channel roughness, all other variables constant. Vary the roughness from that of a smooth glass chute to a mountain stream. How does the discharge vary with roughness?
7) If you were to be watching the Salt River at the downstream end of the Tempe Town Lake and the water were flowing across the location of the downstream dam that is 19' high and the channel width is about 850' and the channel slope there is 0.0013 (Arrowsmith actually calculated based on surveyed elevations in the channel), what would the flow velocity and discharge be (you estimate the roughness)?
8) If you were watching the Salt River at the downstream end of Tempe Town Lake and the water were bank full (to the top of the levees=40'), the channel width were 900' and the slope was the same as above, what would the discharge and velocity be (you estimate the roughness)? In 1993, the river was supposedly flowing at 100,000 cubic feet per second. Based upon your calculations, was it at bankfull?

Grading Rubric (50 points):
• M-files: (15 points)
• Function for calculating velocity and discharge
• Script for plotting and calculations
• Five Plots (5 points each)
• #7, include justification of roughness coefficient based on inspection (5 points)
• #8, include justification of roughness coefficient based on inspection (5 points)

Assignment is due Wednesday, October 21, 2015 before class .

### GLG410 Computers in Earth and Space Exploration

Last modified: October 13, 2015