Breach flow formula (Water Overlay)

Flow through breaches is calculated based on the weir formula, including the consideration between free flow and submerged flow situations. The weir formula is supplied with the breach width (as stand-in for the weir width), which is always calculated with the breach growth formula before calculating the weir flow. The downstream water level can optionally be measured with an additionally configured Breach Level Area, or an automatically placed

For the entry area, the area where the water originates from, the following also applies:

• The water level used is defined by either an input area or the external water level.
• In case of an external area, the surface height used is defined by external surface level. This height limits how far the water level can be lowered.
• Optionally a second value can be supplied for the external surface level. In that case, the second value represents the area of the bottom.
• In case of an input area, the average water level is based on the water level on the individual grid cells of the input area.

Each timestep, the external water level is changed based on the amount of water flowing in or out.

$h_{e}=w_{e,t}-z_{b}$

$\Delta V=min(h_{e}\cdot A_{s},Q_{w,t}\cdot \Delta t)$

In case $A_{s}=A_{b}$:

$w_{e,t+1}=w_{e,t}-{\frac {\Delta V}{A_{b}}}$

Where:

$z_{b}$ Datum height of the bottom of the area water flows from.

$w_{e,t}$ = The water level of the entry area at time $t$; In case of an external area: EXTERNAL_WATER_LEVEL of the breach.

$Q_{w,t}$ = The calculated water flow which takes place at time $t$, based on the weir formula.

$\Delta t$ = The timestep in seconds.

$\Delta V$ = The limited volume of water transported, based on the remaining water.

$A_{s}$ = The size of in the entry area; In case of an external area: EXTERNAL_AREA of the breach.

$A_{b}$ = The (optional) size of the bottom of the entry area; In case of an external area: the second value provided for the EXTERNAL_AREA of the breach. In case only one value is provided, $A_{b}=A_{s}$

Trapezoidal prism

The external water area can also be represented as a trapezoidal prism. The new water level then has to be calculated differently.

The length and width of the surface area are currently considered equal. Therefore:

$L={\sqrt {A_{s}}}$

$W_{s}=L$

The width of the bottom area is then calculated as:

$W_{b}={\frac {A_{b}}{L}}$

Next, we will calculate the new water level based on $\Delta V$:

$h_{0}=w_{e,0}-z_{b}$

$a={\frac {W_{s}-W_{b}}{4\cdot h_{0}}}$

$A_{t,t}=2\cdot (a\cdot h_{t}^{2}+0.5\cdot W_{b}\cdot h_{t})$

$\Delta A={\frac {\Delta V}{L}}$

$c={\frac {-\Delta A*A_{t,t}}{2}}$

$D=({\frac {W_{b}}{2}})^{2}-4\cdot a\cdot c$

if $D<0$, determine the new water level based on a box shape:

$w_{e,t+1}=w_{e,t}+{\frac {\Delta V}{A_{b}}}$

otherwise, determine the new water level based on the determinant and flow direction:

$w_{-}={\frac {-{\frac {W_{b}}{2}}-{\sqrt {D}}}{2\cdot a}}+z_{b}$

$w_{+}={\frac {-{\frac {W_{b}}{2}}+{\sqrt {D}}}{2\cdot a}}+z_{b}$

$w_{e,t+1}={\begin{cases}w_{+}&{\mbox{ if }}\Delta V<0&{\mbox{ and }}w_{+}\leq w_{e,t}\\w_{-}&{\mbox{ if }}\Delta V<0&{\mbox{ and }}w_{+}>w_{e,t}\\w_{+}&{\mbox{ if }}\Delta V\geq 0&{\mbox{ and }}w_{-}<w_{e,t}\\w_{-}&{\mbox{ if }}\Delta V\geq 0&{\mbox{ and }}w_{-}\geq w_{e,t}\\\end{cases}}$

Where:

$z_{b}$ = Datum height of the bottom of the area water flows from.

$w_{e,t}$ = The water level of the entry area at time $t$; In case of an external area: EXTERNAL_WATER_LEVEL of the breach.

$Q_{w,t}$ = The calculated water flow which takes place at time $t$, based on the weir formula.

$\Delta t$ = The timestep in seconds.

$\Delta V$ = The limited volume of water transported, based on the remaining water.

$A_{s}$ = The size of in the entry area; In case of an external area: EXTERNAL_AREA of the breach.

$A_{b}$ = The (optional) size of the bottom of the entry area; In case of an external area: the second value provided for the EXTERNAL_AREA = of the breach. In case only one value is provided, $A_{b}=A_{s}$

$L$ = Length of the trapezoidal prism

$W_{s}$ = Width of the original top (surface at $t=0$) of the trapezoid

$W_{b}$ = Width of the bottom of the trapezoid

$h_{0}$ = The initial relative water height in the external area.

$a$ = Slope of the trapezoid

$A_{t,t}$ = Area of the trapezoidal slide of the trapezoid prism at time $t$

$c$ = Constant based on new area of the trapezoid $A_{t,t+1}$

$D$ = Polynomial discriminant

Notes

• The water level w_e cannot become lower than the surface height s_e, defined by EXTERNAL_SURFACE_LEVEL in case of an external entry area. In case of an input area within the project, the surface height is determined per grid cell.

Related

The following topics are related to this formula.

Features

Breach

Formulas

Weir formula

Breach growth formula

Models

Surface model

---

Water Module

Water Module

Tutorials

Overview

Theory

Benchmarks

Models

Surface • Ground • Rain • Evaporation • Infiltration • Elevation • Breach • Sewer • Storage • Tracer flow • Border • Order

Model formulas

Timestep • Season • Surface water level • Groundwater level

Flow formulas

Surface flow • Surface infiltration • Ground flow • Ground infiltration • Surface evaporation • Ground evaporation • Ground bottom flow

Hydraulic structure formulas

Culvert • Weir • Breach growth • Breach flow • Pump • Sewer Overflow • Inlet • Drainage

Related

Building attributes • Hydrological features • Hydraulic structures •  Terrain attributes • Model attributes • Results • Simulation data