Culvert formula (Water Overlay)
Revision as of 12:38, 14 June 2019 by Bobby@tygron.nl (talk | contribs)
Flow through culverts is based on an open channel flow calculation.
The actual height of the culvert is at least the height of the terrain on either end of the culvert and the provided threshold height:
- B_{c} = max( T_{c} , B_{l} , B_{r} )
The radius of the culvert:
- r = D / 2.
The height of the water column at either end of the culvert, relative to the culvert, is calculated:
- w_{l} = max(0, w_{l} - B_{c})
- w_{r} = max(0, w_{r} - B_{c})
Flow depth d is:
- d = max ( 0 , min ( D , max ( w_{l} , w_{r} ) - B_{c} ) ) )
The loss coefficient for the culvert is calculated:
- U = sqrt( 1 / ( 1 + 2 * g * n^{2} * length / R_{h}^{4/3} ) )
The hydraulic radius R_{h} is calculated using the formula's in the image on the right.
The potential flow through the culvert is then calculated:
- C = U * K * sqrt( 2 * g * abs( w_{l} - w_{r} ) )
Finally the actual amount of water flow is calculated:
- Δf = Δt * C / Δx
Where:
- D = The CULVERT_DIAMETER attribute of the culvert.
- T_{c} = The CULVERT_THRESHOLD attribute of the culvert.
- B_{c} = The surface height of the base of the culvert.
- w_{left} = The water level on the left side of the culvert, relative to datum.
- w_{right} = The water level on the right side of the culvert, relative to datum.
- R_{h} = The hydraulic radius in the culvert^{[1]}.
- K = Flow area, based on the height of the water in the (circular) culvert.
- g = Acceleration factor of GRAVITY, defined for the Water Overlay.
- L = The length of the culvert, calculated as the distance between the culvert's endpoints.
- U = Loss coefficient for culverts.
- n = The CULVERT_N attribute of the culvert.
- C = The potential rate of water flow through the culvert.
- Δf = The water flow which takes place.
- Δt = Computational timestep.
- Δx = Cell size.
See also
References
- ↑ Hydraulic Radius Equations Formulas Calculator ∙ found at: https://www.ajdesigner.com/phphydraulicradius/hydraulic_radius_equation_pipe.php ∙ (last visited 2019-02-11)