Surface model (Water Overlay): Difference between revisions

Revision as of 07:19, 17 June 2019

The Water Module's primary function is the simulation of the 2-dimensional flow of water on the surface. In order to simulate flowing water, the project area is discretized into x by y cells, based on the configured grid cell size.

Secondly a model is required which describes the rules that need to be followed. For this we use the two-dimensional Saint-Venant system, which reads:

where:

h_t = height of water column on time t

u_t = flow speed in the x direction on time t

v_t = flow speed in the y direction on time t

B_x = slope in the x direction

B_y = slope in the y direction

In the Tygron Platform this model is implemented using the Well Balanced Positivity Preserving Central-Upwind Scheme described in Kurganov and Petrova (2007)^[1]. For more information about the implementation, see Surface flow formula.

B = surface elevation of grid cell
h = water height of grid cell
w = water level of grid cell.

The implemented method rewrites the h from the original Saint-Venant system to w = h + B. See the image on the right for clarification on the terms water level (w), surface (or bottom) elevation (B) and water height (h).

Flowing water

Based on variations in the surface elevation and water levels, which may cause unbalance, water will start flowing, until it eventually is balanced in terms of water level and fluxes hu and hv.

Water level initialization

In theory, each grid cell can have a unique water level and accompanying water height. In practice though, water levels are often initialized for large groups of cells, since it is assumed that a particular area within the project area has a given water level. Therefore, the surface water level is initialized based on hydrological features. We make the following distinction:

Water levels of water terrains.
Potential water levels for inundated land.

For all water terrains, the amount of water set on the grid cell is such that the resulting water level in that location is equal to the WATER_LEVEL attribute provided by the Water level area. If no water level area is overlapping that grid cell, or the water level is below the surface elevation of the grid cell, the water height is assumed to be 0 and the water level equal to the surface elevation.

For the second case, inundation areas have been added to the model. Water is placed in all grid cells which are covered by an inundated area, (regardless of the terrain type in that location, in contrast to the water level areas), such that the resulting water level is again equal to the inundation area's INUNDATION_LEVEL attribute. If no inundation area is overlapping that grid cell, or the inundation level is below the surface elevation of the grid cell, the water height is assumed to be 0 and the water level equal to the surface elevation.

The terrain height is defined by the terrain height in the project. The surface height is further influenced by the height of constructions present in the project (bounded by DESIGN_FLOOD_ELEVATION_M) and by the BREACH_HEIGHT of breaches.

Notes

Water can be added to and removed from the described surface system by the rain, evaporation and infiltration model, as well as certain hydraulic structures, sewers and breaches.
Underground flow uses a different flow system and formula's, described in underground model.
In addition to water flowing between adjacent grid cells, water can also flow through hydraulic constructions. This is described in Hydraulic structures.

Related formulas

Related models

References

↑ Kurganov A, Petrova G (2007) ∙ A Second-Order Well-Balanced Positivity Preserving Central-Upwind Scheme for the Saint-Venant System ∙ found at: http://www.math.tamu.edu/~gpetrova/KPSV.pdf (last visited 2019-04-11)

Template:WaterOverlay nav

Water Module

Tutorials

Overview

Theory

Benchmarks

[Kurganov2-1] Kurganov A, Petrova G (2007) ∙ A Second-Order Well-Balanced Positivity Preserving Central-Upwind Scheme for the Saint-Venant System ∙ found at: http://www.math.tamu.edu/~gpetrova/KPSV.pdf (last visited 2019-04-11)

[1]

@@ Line 33: / Line 33: @@
 For the second case, [[Inundation area (Water_Overlay)|inundation area]]s have been added to the model. Water is placed in all grid cells which are covered by an inundated area, (regardless of the terrain type in that location, in contrast to the water level areas), such that the resulting water level is again equal to the inundation area's [[Inundation level (Water Overlay)|INUNDATION_LEVEL]] attribute. If no inundation area is overlapping that grid cell, or the inundation level is below the surface elevation of the grid cell, the water height is assumed to be 0 and the water level equal to the surface elevation.
-The terrain height is defined by the [[Terrain height|terrain height]] in the project. The surface height is further influenced by the height of constructions present in the project (bounded by [[Design flood elevation m model attribute (Water Overlay)|DESIGN_FLOOD_ELEVATION_M]]) and by the [[Breach height (Water Overlay)|BREACH_HEIGHT]] of [[Breach (Water Overlay)|breaches]].
+The terrain height is defined by the [[terrain height]] in the project. The surface height is further influenced by the height of constructions present in the project (bounded by [[Design flood elevation m model attribute (Water Overlay)|DESIGN_FLOOD_ELEVATION_M]]) and by the [[Breach height (Water Overlay)|BREACH_HEIGHT]] of [[Breach (Water Overlay)|breaches]].
 ==Notes==