Culvert formula (Water Overlay)

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={\frac {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 {\frac {1}{1+2g\cdot n^{2}\cdot length/R_{h}^{\frac {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:

Failed to parse (unknown function "\abs"): {\displaystyle C = U \cdot K \cdot \sqrt{ 2g \cdot \abs{ w_l - w_r } }}

Finally the actual amount of water flow is calculated:

Failed to parse (syntax error): {\displaystyle Δ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.