3Di can simulate the consequences of a breach in a levee. In such case, the dimensions of the breach determine flow through the breach and therefore the flood. The growth of a breach is a very complex process. It is a phenomenon where hydrodynamics, morphology, groundwater and soil mechanics interact. Many aspects of this interaction is still unknown, and therefore it is difficult to model. Moreover, detailed modeling would require very detailed information about the levee. Generally, there is only a limited amount of information about the soil, structure, moisture content etc of the levee. Nevertheless, rules of thumb have been developed that describe the breach growth, which allow us to make an estimate of the breach growth.
13.1. Breaches from 1D channels to 2D grid¶
In a 3Di model, flow may occur between 1D and 2D elements. In certain cases this exchange is limited by a levee. The exchange height is determined by the height of the levee. When a few extra properties of the levee are specified for these connections, a breach can be modelled that can grow over time. More information on levees, can be found in Obstacles and levees
13.1.1. Breach growth formulation¶
In 3Di, the formula of Verheij and van de Knaap (2003) is used to describe the growth of a breach. For this formulation it is expected that the material of the levee (sand or clay), an initial breach width, the maximum breach depth and the period that is needed to reach this depth are known.
This formulation splits the growth of the breach into two phases; in the first phase the breach deepens and in the second phase the breach widens. Mathematically, the first phase can be described as:
Once the minimum breach level is reached, the width of the breach increases according to:
Because the water level gradient in front and behind the breach is included in the formulation, a natural balance may arise. This means that the breach stops growing, when the water levels are equal.
The above formulation also corrects for the presence of different types of materials by using a critical velocity and the parameters and . The growth rate of the breach increases for materials that erode more easily. In principle, this formulation distinguishes only sand and clay.
13.1.2. Exchange formulation¶
The flow between 1D and 2D at the breach is computed based on a simplified momentum balance. A balance is made between the friction and the forcing. Note, that the volume in the breach is neglected.
More details on how to use obstacles, levees and breaches can be found in Flood model.