Stowa 2017 benchmark (Water Module): Difference between revisions
No edit summary |
|||
| (43 intermediate revisions by 4 users not shown) | |||
| Line 1: | Line 1: | ||
The Stowa 2017 benchmark was set up in order to find out which tool currently available to the Dutch Water management industry performs reliable in which circumstances. | The Stowa 2017 [[Water Module benchmarks|benchmark]] was set up in order to find out which tool currently available to the Dutch Water management industry performs reliable in which circumstances. | ||
The goal of the benchmark was to provide insight to waterboards and municipalities which water simulation tool would suit their needs best. | The goal of the benchmark was to provide insight to waterboards and municipalities which water simulation tool would suit their needs best. | ||
The benchmark consist of several topics: | The benchmark consist of several topics: | ||
| Line 7: | Line 7: | ||
* 2D-1D interaction between surface water and water objects such as weirs and culverts. | * 2D-1D interaction between surface water and water objects such as weirs and culverts. | ||
For the complete benchmark, see [https://www.stowa.nl/sites/default/files/assets/PUBLICATIES/Publicaties%202017/STOWA%202017-34.pdf | For the complete benchmark, see [https://www.stowa.nl/sites/default/files/assets/PUBLICATIES/Publicaties%202017/STOWA%202017-34.pdf Stowa benchmark 2017]. | ||
Please note that the referenced document contains results made in 2017, with a model that was not yet using the 2D Saint Venant equations. | Please note that the referenced document contains results made in 2017, with a model that was not yet using the 2D Saint Venant equations. | ||
The Water Module has been much improved since 2017 and is now based on the finite volume method described by Kurganov et al. | The Water Module has been much improved since 2017 and is now based on the finite volume method described by Kurganov et al. | ||
The results presented below have been generated on 04-09-2019. | The results presented below have been generated on '''04-09-2019'''. | ||
==Test case 1== | ==Test case 1== | ||
This test case will validate whether a shallow water calculation tool can calculate sheet flow accurately. An analytic solution of this case is provided by (Di Giammarco, 1996). | |||
The reference grid resolution is 20 by 20 meters. The Manning roughness is 0.015. The surface is expected to be impenetrable, meaning that there is no rain water infiltrating into the underground. Water that reaches the lowest border of this test case is expected to flow out of the system. | |||
[[File:stowa2017-testcase1.png|left|700px|Fig 5.1 Layout of testcase 1 with a sloped ditch and sloped terrain in both directions]] | |||
[[File:stowa2017-testcase1-doel.png|left|700px|Purpose of testcase 1]] | |||
<br style="clear:left"> | |||
===Setup=== | |||
A sloped terrain of 800 by 1000 meters, on which a rain fall is applied of 10.8 mm per hour for a duration of 90 minutes. The total simulation time is 3 hours. | |||
A graphical impression of the height map in this testcase is given below. Again first the original height map and next the height map adjusted for ensuring a constant slope per grid cell x- and y-direction. | |||
<ul> | |||
<li style="display: inline-block;">[[File:stowa2017-testcase1-heightmap.png|thumb|left|200px|Initial]]</li> | |||
<li style="display: inline-block;">[[File:stowa2017-testcase1-heightmap-modified.png|thumb|left|200px|After slope adjustments]]</li> | |||
</ul> | |||
===Results=== | |||
The hydrograph and volumes. | |||
<gallery> | |||
File:stowa2017-testcase1_hydrograph.png|hydrograph | |||
File:stowa2017-testcase1_volumes.png|volumes | |||
</gallery> | |||
The amount of surface water (in meters), color coded per grid cell, presented for several timeframes. | |||
<gallery> | |||
File:stowa2017-testcase1_levels_0h30m.png|0h30m | |||
File:stowa2017-testcase1_levels_1h0m.png|1h00m | |||
File:stowa2017-testcase1_levels_1h30m.png|1h30m | |||
File:stowa2017-testcase1_levels_1h45m.png|1h45m | |||
File:stowa2017-testcase1_levels_2h0m.png|2h00m | |||
File:stowa2017-testcase1_levels_2h30m.png|2h30m | |||
File:stowa2017-testcase1_levels_3h0m.png|3h00m | |||
</gallery> | |||
==Test case 2== | ==Test case 2== | ||
Test case 2 adds the following aspects with respect to test case 1: | |||
* Flow coming from two directions that meet in the middle | |||
* Deeper surface water flow | |||
An analytic solution is available for this case. However, it is more complicated that test case 1. The influence of model boundaries should be less in this test case compared to test case 1. | |||
[[File:stowa2017-testcase2.png|left|700px|Fig 5.9 Layout of testcase 2 with a sloped ditch and sloped terrain in both directions]] | |||
[[File:stowa2017-testcase2-doel.png|left|700px|Purpose of testcase 2]] | |||
<br style="clear:left"> | |||
===Setup=== | |||
The lay-out of the terrain has a slanted V-shape. In the middle is a ditch, which is horizontal in the x-direction but slanted downwards in the y-direction with a value of 0.02. The terrain around the ditch is slanted downwards in the x-direction with a value of 0.05. Furthermore, the ditch is +-0.05m deep at the top. In the y-direction is remains horizontal. For the exact measurements see figure 5.9 above. | |||
A graphical impression of the height map in this testcase is given below. Again first the original height map and next the height map adjusted for ensuring a constant slope per grid cell x- and y-direction. | |||
<ul> | |||
<li style="display: inline-block;">[[File:stowa2017-testcase2-heightmap.png|thumb|left|200px|Initial]]</li> | |||
<li style="display: inline-block;">[[File:stowa2017-testcase2-heightmap-modified.png|thumb|left|200px|After slope adjustments]]</li> | |||
</ul> | |||
===Results=== | |||
The hydrograph and volumes. | |||
<gallery> | |||
File:stowa2017-testcase2_hydrograph.png|hydrograph | |||
File:stowa2017-testcase2_volumes.png|volumes | |||
</gallery> | |||
The amount of surface water (in meters), color coded per grid cell, presented for several timeframes. | |||
<gallery> | |||
File:stowa2017-testcase2_levels_0h30m.png|0h30m | |||
File:stowa2017-testcase2_levels_1h0m.png|1h00m | |||
File:stowa2017-testcase2_levels_1h30m.png|1h30m | |||
File:stowa2017-testcase2_levels_1h45m.png|1h45m | |||
File:stowa2017-testcase2_levels_2h0m.png|2h00m | |||
File:stowa2017-testcase2_levels_2h30m.png|2h30m | |||
File:stowa2017-testcase2_levels_3h0m.png|3h00m | |||
</gallery> | |||
==Test case 3== | ==Test case 3== | ||
This test case adds a slope in the y-direction as well. Due to this extra slope, the flow directions will not be perpendicular anymore to the sides of the grid cells. Instead, the flow direction will vary per grid cell and this is what this test case will validate. | |||
[[File:stowa2017-testcase3.png|left|700px|Fig 5.9 Layout of testcase 3 with a sloped ditch and sloped terrain in both directions]] | |||
[[File:stowa2017-testcase3-doel.png|left|700px|Purpose of testcase 3]] | |||
<br style="clear:left"> | |||
===Setup=== | |||
The lay-out of the terrain has a slanted V-shape. In the middle is a ditch, which is horizontal in the x-direction but slanted downwards in the y-direction. The terrain around the ditch is slanted downwards in the x-direction with a value of 0.05. In the y-direction is is slanted similar to the ditch, with a value of 0.02. For the exact measurements see figure 5.9 above. | |||
A graphical impression of the height map in this testcase is given below. Again first the original height map and next the height map adjusted for ensuring a constant slope per grid cell x- and y-direction. | |||
<ul> | |||
<li style="display: inline-block;">[[File:stowa2017-testcase3-heightmap.png|thumb|left|200px|Initial]]</li> | |||
<li style="display: inline-block;">[[File:stowa2017-testcase3-heightmap-modified.png|thumb|left|200px|After slope adjustments]]</li> | |||
</ul> | |||
===Results=== | |||
The hydrograph and volumes. | |||
<gallery> | |||
File:stowa2017-testcase3_hydrograph.png|hydrograph | |||
File:stowa2017-testcase3_volumes.png|volumes | |||
</gallery> | |||
The amount of surface water (in meters), color coded per grid cell, presented for several timeframes. | |||
<gallery> | |||
File:stowa2017-testcase3_levels_0h30m.png|0h30m | |||
File:stowa2017-testcase3_levels_1h0m.png|1h00m | |||
File:stowa2017-testcase3_levels_1h30m.png|1h30m | |||
File:stowa2017-testcase3_levels_1h45m.png|1h45m | |||
File:stowa2017-testcase3_levels_2h0m.png|2h00m | |||
File:stowa2017-testcase3_levels_2h30m.png|2h30m | |||
File:stowa2017-testcase3_levels_3h0m.png|3h00m | |||
</gallery> | |||
==Test case 4== | ==Test case 4== | ||
This testcase contains a ponding problem. This testcase is subdivided into 3 sub-testcases in order to investigate several aspects of this ponder problem. Each is described below. See image 5.11 for the layout of the problem. Each | This testcase contains a ponding problem. This testcase is subdivided into 3 sub-testcases in order to investigate several aspects of this ponder problem. Each is described below. See image 5.11 for the layout of the problem. Each of the surfaces in this test case have a Manning roughness value of 0.05. A weir is situated at the bottom of the sloped terrain, causing the creation of a pond during the simulation. | ||
[[File:stowa2017-testcase4.png|left|700px|Fig 5.11 Layout of the ponding problem for testcase 4]] | [[File:stowa2017-testcase4.png|left|700px|Fig 5.11 Layout of the ponding problem for testcase 4]] | ||
<br style="clear:left"> | <br style="clear:left"> | ||
==Setup== | ===Setup=== | ||
For each test case the base setup is the same. | For each test case the base setup is the same. A graphical impression of the height map in this testcase is given below. Again first the original height map and next the height map adjusted for ensuring a constant slope per grid cell x- and y-direction. | ||
<ul> | <ul> | ||
<li style="display: inline-block;">[[File:stowa2017-testcase4-heightmap.png|thumb|left|200px|Initial]]</li> | <li style="display: inline-block;">[[File:stowa2017-testcase4-heightmap.png|thumb|left|200px|Initial]]</li> | ||
<li style="display: inline-block;">[[File:stowa2017-testcase4-heightmap-modified.png|thumb|left|200px|After slope adjustments]]</li> | <li style="display: inline-block;">[[File:stowa2017-testcase4-heightmap-modified.png|thumb|left|200px|After slope adjustments]]</li> | ||
</ul> | </ul> | ||
===Test case 4A=== | ===Test case 4A=== | ||
In this testcase, the amount of rainfall is similar to testcase 3. However, it continues to rain until the end of the simulation. The surface and underground are impenetrable. And the weir is considered to be 'perfect'. | In this testcase, the amount of rainfall is similar to testcase 3. However, it continues to rain until the end of the simulation. The surface and underground are impenetrable. And the weir is considered to be 'perfect'. | ||
[[File:stowa2017-testcase4a-doel.png | [[File:stowa2017-testcase4a-doel.png|left|700px]] | ||
<br style="clear:left"> | <br style="clear:left"> | ||
====Results==== | ====Results==== | ||
The hydrograph and | The hydrograph and volumes. | ||
<gallery> | <gallery> | ||
File:stowa2017-testcase4a_hydrograph.png|hydrograph | File:stowa2017-testcase4a_hydrograph.png|hydrograph | ||
| Line 60: | Line 157: | ||
===Test case 4B=== | ===Test case 4B=== | ||
In this testcase the amount of rainfall is 21.6 mm per hour and lasts for 90 minutes. The underground has an infiltration capacity of 10.8 mm per hour. The simulation runs until the water level has lowered to about +/- 2cm below the crest of the weir. The effects of the groundwater table can be ignored. | In this testcase the amount of rainfall is 21.6 mm per hour and lasts for 90 minutes. The underground has an infiltration capacity of 10.8 mm per hour. The simulation runs until the water level has lowered to about +/- 2cm below the crest of the weir. The effects of the groundwater table can be ignored. | ||
[[File:stowa2017-testcase4b-doel.png|left|700px]] | |||
<br style="clear:left"> | |||
====Results==== | ====Results==== | ||
The hydrograph and | The hydrograph and volumes. | ||
<gallery> | <gallery> | ||
File:stowa2017-testcase4b_hydrograph.png|hydrograph | File:stowa2017-testcase4b_hydrograph.png|hydrograph | ||
| Line 81: | Line 181: | ||
===Test case 4C=== | ===Test case 4C=== | ||
In this testcase, the amount of rainfall and the infiltration capacity are equal to [[#testcase 4B|testcase 4B]]. The porosity of the soil is now set to 0.25. The impenetrable layer is strictly horizontal and 10 meters below the lowest surface height. The groundwatertable at t=0 is horizontal, 50 cm below the height of the weir. The initial surface water level is equal to the initial groundwatertable. The hydraulic condictivity is 5x10<sup>-5</sup> for water infiltration from the surface to the underground. The simulation time is equal to testcase 4B. | In this testcase, the amount of rainfall and the infiltration capacity are equal to [[#testcase 4B|testcase 4B]]. The porosity of the soil is now set to 0.25. The impenetrable layer is strictly horizontal and 10 meters below the lowest surface height. The groundwatertable at t=0 is horizontal, 50 cm below the height of the weir. The initial surface water level is equal to the initial groundwatertable. The hydraulic condictivity is 5x10<sup>-5</sup> for water infiltration from the surface to the underground. The simulation time is equal to testcase 4B. | ||
[[File:stowa2017-testcase4c-doel.png|left|700px]] | |||
<br style="clear:left"> | |||
====Results==== | ====Results==== | ||
The hydrograph and | The hydrograph and volumes. | ||
<gallery> | <gallery> | ||
File:stowa2017-testcase4c_hydrograph.png|hydrograph | File:stowa2017-testcase4c_hydrograph.png|hydrograph | ||
| Line 100: | Line 204: | ||
</gallery> | </gallery> | ||
===Test case 5 to 7=== | |||
Currently not implemented. | |||
<!-- | |||
==Complete set of results== | ==Complete set of results== | ||
The complete set of produced results for the Stowa 2017 benchmark can be downloaded as a .zip-file [[File:stowa2017-results2019.zip|here]]. | The complete set of produced results for the Stowa 2017 benchmark can be downloaded as a .zip-file [[File:stowa2017-results2019.zip|here]]. | ||
--> | |||
[[Category:Benchmarks]] | |||
{{Water Module buttons}} | |||
Latest revision as of 15:34, 5 March 2024
The Stowa 2017 benchmark was set up in order to find out which tool currently available to the Dutch Water management industry performs reliable in which circumstances. The goal of the benchmark was to provide insight to waterboards and municipalities which water simulation tool would suit their needs best. The benchmark consist of several topics:
- Rain water discharge for several types of terrain
- Underground water infiltration
- Flow direction and speed
- 2D-1D interaction between surface water and water objects such as weirs and culverts.
For the complete benchmark, see Stowa benchmark 2017. Please note that the referenced document contains results made in 2017, with a model that was not yet using the 2D Saint Venant equations.
The Water Module has been much improved since 2017 and is now based on the finite volume method described by Kurganov et al.
The results presented below have been generated on 04-09-2019.
Test case 1
This test case will validate whether a shallow water calculation tool can calculate sheet flow accurately. An analytic solution of this case is provided by (Di Giammarco, 1996).
The reference grid resolution is 20 by 20 meters. The Manning roughness is 0.015. The surface is expected to be impenetrable, meaning that there is no rain water infiltrating into the underground. Water that reaches the lowest border of this test case is expected to flow out of the system.
Setup
A sloped terrain of 800 by 1000 meters, on which a rain fall is applied of 10.8 mm per hour for a duration of 90 minutes. The total simulation time is 3 hours.
A graphical impression of the height map in this testcase is given below. Again first the original height map and next the height map adjusted for ensuring a constant slope per grid cell x- and y-direction.
Results
The hydrograph and volumes.
hydrograph
volumes
The amount of surface water (in meters), color coded per grid cell, presented for several timeframes.
0h30m
1h00m
1h30m
1h45m
2h00m
2h30m
3h00m
Test case 2
Test case 2 adds the following aspects with respect to test case 1:
- Flow coming from two directions that meet in the middle
- Deeper surface water flow
An analytic solution is available for this case. However, it is more complicated that test case 1. The influence of model boundaries should be less in this test case compared to test case 1.
Setup
The lay-out of the terrain has a slanted V-shape. In the middle is a ditch, which is horizontal in the x-direction but slanted downwards in the y-direction with a value of 0.02. The terrain around the ditch is slanted downwards in the x-direction with a value of 0.05. Furthermore, the ditch is +-0.05m deep at the top. In the y-direction is remains horizontal. For the exact measurements see figure 5.9 above. A graphical impression of the height map in this testcase is given below. Again first the original height map and next the height map adjusted for ensuring a constant slope per grid cell x- and y-direction.
Results
The hydrograph and volumes.
hydrograph
volumes
The amount of surface water (in meters), color coded per grid cell, presented for several timeframes.
0h30m
1h00m
1h30m
1h45m
2h00m
2h30m
3h00m
Test case 3
This test case adds a slope in the y-direction as well. Due to this extra slope, the flow directions will not be perpendicular anymore to the sides of the grid cells. Instead, the flow direction will vary per grid cell and this is what this test case will validate.
Setup
The lay-out of the terrain has a slanted V-shape. In the middle is a ditch, which is horizontal in the x-direction but slanted downwards in the y-direction. The terrain around the ditch is slanted downwards in the x-direction with a value of 0.05. In the y-direction is is slanted similar to the ditch, with a value of 0.02. For the exact measurements see figure 5.9 above. A graphical impression of the height map in this testcase is given below. Again first the original height map and next the height map adjusted for ensuring a constant slope per grid cell x- and y-direction.
Results
The hydrograph and volumes.
hydrograph
volumes
The amount of surface water (in meters), color coded per grid cell, presented for several timeframes.
0h30m
1h00m
1h30m
1h45m
2h00m
2h30m
3h00m
Test case 4
This testcase contains a ponding problem. This testcase is subdivided into 3 sub-testcases in order to investigate several aspects of this ponder problem. Each is described below. See image 5.11 for the layout of the problem. Each of the surfaces in this test case have a Manning roughness value of 0.05. A weir is situated at the bottom of the sloped terrain, causing the creation of a pond during the simulation.
Setup
For each test case the base setup is the same. A graphical impression of the height map in this testcase is given below. Again first the original height map and next the height map adjusted for ensuring a constant slope per grid cell x- and y-direction.
Test case 4A
In this testcase, the amount of rainfall is similar to testcase 3. However, it continues to rain until the end of the simulation. The surface and underground are impenetrable. And the weir is considered to be 'perfect'.
Results
The hydrograph and volumes.
hydrograph
volumes
The amount of surface water (in meters), color coded per grid cell, presented for several timeframes.
0h30m
1h00m
1h30m
1h45m
2h00m
2h30m
3h00m
Test case 4B
In this testcase the amount of rainfall is 21.6 mm per hour and lasts for 90 minutes. The underground has an infiltration capacity of 10.8 mm per hour. The simulation runs until the water level has lowered to about +/- 2cm below the crest of the weir. The effects of the groundwater table can be ignored.
Results
The hydrograph and volumes.
hydrograph
volumes
The amount of surface water (in meters), color coded per grid cell, presented for several timeframes.
0h30m
1h00m
1h30m
1h45m
2h00m
2h30m
3h00m
Test case 4C
In this testcase, the amount of rainfall and the infiltration capacity are equal to testcase 4B. The porosity of the soil is now set to 0.25. The impenetrable layer is strictly horizontal and 10 meters below the lowest surface height. The groundwatertable at t=0 is horizontal, 50 cm below the height of the weir. The initial surface water level is equal to the initial groundwatertable. The hydraulic condictivity is 5x10-5 for water infiltration from the surface to the underground. The simulation time is equal to testcase 4B.
Results
The hydrograph and volumes.
hydrograph
volumes
The amount of surface water (in meters), color coded per grid cell, presented for several timeframes.
0h30m
1h00m
1h30m
1h45m
2h00m
2h30m
3h00m
Test case 5 to 7
Currently not implemented.
