Learn R Programming

iemisc (version 0.5.0)

Manningcirc: Circular cross-section using the Gauckler-Manning-Strickler equation

Description

Manningcirc and Manningcircy solve for a missing variable for a circular cross-section.

Usage

Manningcirc(Q = NULL, n = NULL, Sf = NULL, y = NULL, d = NULL,
  units = c("SI", "Eng"))

Manningcircy(y = NULL, d = NULL, y_d = NULL, theta = NULL, Sf = NULL, Q = NULL, units = c("SI", "Eng"))

Arguments

Q
numeric vector that contains the discharge value [m^3/s or ft^3/s], if known.
n
numeric vector that contains the Manning's roughness coefficient n, if known.
Sf
numeric vector that contains the the bed slope (m/m or ft/ft), if known.
y
numeric vector that contains the flow depth (m or ft), if known.
d
numeric vector that contains the diameter value (m or ft), if known.
units
character vector that contains the system of units [options are SI for International System of Units and Eng for English units (United States Customary System in the United States and Imperial Units in the United Kingdom)]
y_d
numeric vector that contains the filling ration (y/d), if known.
theta
numeric vector that contains the angle theta (radians), if known.

Value

  • the missing parameter (Q, n, or Sf) & theta, area (A), wetted perimeter (P), top width (B), and R (hydraulic radius) as a list for the Manningcirc function.

    the missing parameter (d or y) & theta, area (A), wetted perimeter (P), top width (B), and R (hydraulic radius) as a list for the Manningcircy function.

encoding

UTF-8

source

r - Better error message for stopifnot? - Stack Overflow answered by Andrie on Dec 1 2011. See http://stackoverflow.com/questions/8343509/better-error-message-for-stopifnot.

Details

The Manningcirc function solves for one missing variable in the Gauckler- Manning equation for a circular cross-section and uniform flow. The possible inputs are Q, n, Sf, y, and d. If y or d are not initially known, then Manningcircy can solve for y or d to use as input in the Manningcirc function.

The Manningcircy function solves for one missing variable in the Gauckler- Manning equation for a circular cross-section and uniform flow. The possible inputs are y, d, y_d (ratio of y/d), and theta.

Gauckler-Manning-Strickler equation is expressed as

$$V = \frac{K_n}{n}R^\frac{2}{3}\sqrt{S}$$

[object Object],[object Object],[object Object],[object Object],[object Object]

This equation is also expressed as

$$Q = \frac{K_n}{n}\frac{A^\frac{5}{3}}{P^\frac{2}{3}}\sqrt{S}$$

[object Object],[object Object],[object Object],[object Object],[object Object],[object Object]

Other important equations regarding the circular cross-section follow: $$R = \frac{A}{P}$$

[object Object],[object Object],[object Object]

$$A = \left(\theta - \sin \theta\right) \frac{d^2}{8}$$

[object Object],[object Object],[object Object]

$$\theta = 2 \arcsin\left[1 - 2\left(\frac{y}{d}\right)\right]$$ [object Object],[object Object],[object Object]

$$d = 1.56 \left[\frac{nQ}{K_n\sqrt{S}}\right]^\frac{3}{8}$$ [object Object],[object Object],[object Object],[object Object],[object Object]

Note: This will only provide the initial conduit diameter, check the design considerations to determine your next steps.

$$P = \frac{\theta d}{2}$$

[object Object],[object Object],[object Object]

$$B = d \sin\left(\frac{\theta}{2}\right)$$

[object Object],[object Object],[object Object]

Assumptions: uniform flow and prismatic channel

Note: Units must be consistent

References

  1. Terry W. Sturm,Open Channel Hydraulics, 2nd Edition, New York City, New York: The McGraw-Hill Companies, Inc., 2010, page 8, 36, 102, 120, 123-125, 153-154.
  2. Dan Moore, P.E., NRCS Water Quality and Quantity Technology Development Team, Portland Oregon, "Using Mannings Equation with Natural Streams", August 2011,http://www.wcc.nrcs.usda.gov/ftpref/wntsc/H&H/xsec/manningsNaturally.pdf.
  3. Gilberto E. Urroz, Utah State University Civil and Environmental Engineering, CEE6510 - Numerical Methods in Civil Engineering, Spring 2006, "Solving selected equations and systems of equations in hydraulics using Matlab", August/September 2004,http://ocw.usu.edu/Civil_and_Environmental_Engineering/Numerical_Methods_in_Civil_Engineering/.
  4. Tyler G. Hicks, P.E.,Civil Engineering Formulas: Pocket Guide, 2nd Edition, New York City, New York: The McGraw-Hill Companies, Inc., 2002, page 423, 425.
  5. Wikimedia Foundation, Inc. Wikipedia, 26 November 2015, “Manning formula”,https://en.wikipedia.org/wiki/Manning_formula.

See Also

Manningtrap for a trapezoidal cross-section, Manningrect for a rectangular cross-section, Manningtri for a triangular cross-section, and Manningpara for a parabolic cross-section.

Examples

Run this code
library(iemisc)
library(iemiscdata)
# The equations used in this function were solved analytically after
# following these steps:
library(rSymPy) # review the package to determine its system dependencies
Q <- Var("Q")
n <- Var("n")
k <- Var("k")
Sf <- Var("Sf")
y <- Var("y")
d <- Var("d")
theta <- Var("theta")

# Simplify with rSymPy
eqsimp <- sympy("expr = n*Q*((theta*d)/2)**(2/3)-k*(theta-sin(theta))*(d**2/8)**(5/3)*sqrt(Sf)")
# eqsimp is "Q*n - k*Sf**(1/2)*d**2*(theta - sin(theta))/8"
# This is the equation that was used to solve for the missing variables




# Example 4.1 from Sturm (page 124-125)
Manningcircy(y_d = 0.8, d = 2, units = "Eng")
y <- Manningcircy(y_d = 0.8, d = 2, units = "Eng")
# defines all list values within the object named y
y$y # gives the value of y

Run the code above in your browser using DataLab