manningc: Solves the Manning Equation for gravity flow in a circular pipe

Description

This function solves the Manning equation for water flow in a circular pipe at less than full. Uniform flow conditions are assumed, so that the pipe slope is equal to the slope of the water surface and the energy grade line. This is a modification of the code prepared by Irucka Embry in his iemisc package. The iemisc::manningcirc function was adapted here for more limited cases commonly used in classroom exercises, additional checks were included to ensure the pipe is flowing less than full, and a cross-section figure is also available. The iemisc::manningcirc and iemisc::manningcircy functions were combined into a single function. Manning n supplied is assumed to be that for full pipe flow; an optional argument may be supplied to account for n varying with depth.

Usage

manningc(
  Q = NULL,
  n = NULL,
  Sf = NULL,
  y = NULL,
  d = NULL,
  y_d = NULL,
  n_var = FALSE,
  units = c("SI", "Eng"),
  ret_units = FALSE
)

Value

Returns a list including the missing parameter:

Q - flow rate
V - flow velocity
A - cross-sectional area of flow
P - wetted perimeter
R - hydraulic radius (A/P)
y - flow depth
d - pipe diameter
Sf - slope
n - Manning's roughness (for full flow, or as adjusted if n_var is TRUE)
yc - critical depth
Fr - Froude number
Re - Reynolds number
Qf - Full pipe flow rate

Arguments

Q: numeric vector that contains the flow rate [ $m^{3} s^{- 1}$ or $f t^{3} s^{- 1}$ ]
n: numeric vector that contains the Manning roughness coefficient (for full flow or fixed).
Sf: numeric vector that contains the slope of the pipe [unitless]
y: numeric vector that contains the water depth [ $m$ or $f t$ ]
d: numeric vector that contains the pipe diameter [ $m$ or $f t$ ]
y_d: numeric vector that contains the ratio of depth to diameter [unitless]
n_var: If set to TRUE the value of n for full flow is adjusted with depth. [Default is FALSE]
units: character vector that contains the system of units [options are SI for International System of Units and Eng for English (US customary) units. This is used for compatibility with iemisc package]
ret_units: If set to TRUE the value(s) returned are of class units with units attached to the value. [Default is FALSE]

Author

Ed Maurer, Irucka Embry

Details

The possible applications of this function for solving the Manning equation for circular pipes are:

Given	Solve for
y_d, Q, Sf, n	d
d, Sf, Q, n	y
y, d, Q, n	Sf
y, d, Sf, n	Q
d, Q, Sf, y	n

The Manning equation (also known as the Strickler equation) describes flow conditions in an open channel under uniform flow conditions. It is often expressed as: $Q = A \frac{C}{n} R^{\frac{2}{3}} {S_{f}}^{\frac{1}{2}}$ where $C$ is 1.0 for SI units and 1.49 for Eng (U.S. Customary) units. Critical depth is defined by the relation (at critical conditions): $\frac{Q^{2} B}{g A^{3}} = 1$ where $B$ is the top width of the water surface. Since B equals zero for a full pipe, critical depth is set to the pipe diameter $d$ if the flow $Q$ exceeds a value that would produce a critical flow at $\frac{y}{d} = 0.99$ .

Examples

Run this code


#Solving for flow rate, Q: SI Units
manningc(d = 0.6, n = 0.013, Sf = 1./400., y = 0.24, units = "SI")
#returns 0.1 m3/s

#Solving for Sf, if d=600 mm and pipe is to flow half full
manningc(d = 0.6, Q = 0.17, n = 0.013, y = 0.3, units = "SI")
#returns required slope of 0.003

#Solving for diameter, d when given y_d): Eng (US) units
manningc(Q = 83.5, n = 0.015, Sf = 0.0002, y_d = 0.9, units = "Eng")
#returns 7.0 ft required diameter

#Solving for depth, y when given Q: SI units
manningc(Q=0.01, n=0.013, Sf=0.001, d = 0.2, units="SI")
#returns depth  y = 0.158 m, critical depth, yc = 0.085 m

#Solving for depth, y when given Q: SI units, and n variable with depth
manningc(Q=0.01, n=0.013, Sf=0.001, d = 0.2, n_var = TRUE, units="SI")
#returns depth  y = 0.174 m, critical depth, yc = 0.085 m

Run the code above in your browser using DataLab

Last chance! 50% off unlimited learning