Breach growth formula (Water Overlay)

Water can flow through breaches into levee protected areas. These breaches often start small and grow over time.

The water flowing through breaches can originate from an external area outside the project area or an input area within the project area.

This algorithm is based on the incremental timestep formula of Verheij-van der Knaap, as described in section 3.4.4 Implementatie in SOBEK and the conclusion of that paper.

First, the difference in height of the water on either side of the breach is calculated.

$\Delta h_{t}=\left|w_{o,t}-max(w_{i,t},H_{b,t})\right|$

Using the height difference, the breach width increase (m/s) is calculated per computational timestep.

$\Delta W_{b,t}={\frac {f_{1}\cdot f_{2}}{ln(10)}}\cdot {\frac {(g\cdot \Delta h_{t})^{1.5}}{{cs_{b}}^{2}}}\cdot {\frac {1}{1+{\frac {f_{2}\cdot g\cdot \Delta t}{cs_{b}\cdot 3600}}}}$

The current breach width is then equal to the last calculated breach width, plus the calculated breach width increment.

$W_{b,t}=W_{b,t-1}+\Delta W_{b,t}\cdot {\frac {\Delta t}{3600}}$

Where:

$W_{b}$ = The BREACH_WIDTH of the breach.

$H_{b,t}$ = The BREACH_HEIGHT of the breach at time t.

$W_{b,t}$ = The calculated breach width, initially equal to $W_{b}$.

$w_{i,t}$ = Inner water level at breach area at time t.

$w_{o,t}$ = Outer water level at Input Area) (or external area) at time t.

$\Delta h_{t}$ = The difference between the height of the water columns on either side of the breach at time t.

$f_{1}$ = Material factor, set to 1.3 (average for sand and clay levees).

$f_{2}$ = Constant, set to 0.04.

$g$ = the acceleration due to gravity.

$cs_{b}$ = The critical BREACH_SPEED of the breach (e.g. 0.2 for sand and 0.5 for clay).

$\Delta W_{b,t}$ = The calculated width increase of the breach at time t.

$\Delta t$ = Computational timestep.

Notes

• The breach height H_b,t can be defined over time, by configuring the BREACH_HEIGHT as an attribute array.
• The inner water level is to be measured not directly inside the breach area but a few meters away. By default 100m, this can be adjusted via the BREACH_MEASUREMENT_DISTANCE_M attribute. It can also be measured using a manual defined level area.
• For more details on breach growth, we also recommend reading this "Quickscan"-report for Waterschap Rijn en IJssel.

How-to's

• How to add a breach level area to a breach

Related

The following topics are related to this formula.

Features

Breach

Formulas

Breach flow formula

Models

Surface model

See also

For an example of the breach growth, take a look at the Demo Breach Project available in all domains.

References

---

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