UK Benchmark 1 (Water Module): Difference between revisions
Jump to navigation
Jump to search
No edit summary |
|||
| (19 intermediate revisions by 3 users not shown) | |||
| Line 1: | Line 1: | ||
This | This [[Water Module benchmarks|benchmark]] performs test case 1 of the UK benchmark <ref name="Néelz"/> | ||
named ''Test 1 – Uncovering of a beach'' as well as its result generated by the [[Water Module]] in the {{software}}. | |||
''This test consists of a sloping topography with a depression. An inflow boundary condition is applied at the low end, causing the water to rise to the level indicated by a thick blue line. The inflow is then replaced by a sink term until the water level becomes as indicated by the thin blue line. A similar test has been used by EDF (2000) for validation of the TELEMAC package. | ''This test consists of a sloping topography with a depression. An inflow boundary condition is applied at the low end, causing the water to rise to the level indicated by a thick blue line. The inflow is then replaced by a sink term until the water level becomes as indicated by the thin blue line. A similar test has been used by EDF (2000) for validation of the TELEMAC package. | ||
| Line 5: | Line 6: | ||
==Description== | ==Description== | ||
[[File:Sloping_topography_ukbm.png| | [[File:Sloping_topography_ukbm.png|left|frame|Fig. a: Sloping topography with depression]] | ||
''This test consists of a sloping topography with a depression as illustrated in Figure (a). The modeled domain is a perfect 700m x 100m rectangle. A varying water level, see Figure (b), is applied as a boundary condition along the entire length of the left-hand side of the rectangle, causing the water to rise to level 10.35m. This elevation is maintained for long enough for the water to fill the depression and become horizontal over the entire domain. It is then lowered back to its initial state, causing the water level in the pond to become horizontal at the same elevation as the sill, 10.25m.'' | |||
This test consists of a sloping topography with a depression as illustrated in Figure (a). The modeled domain is a perfect 700m x 100m rectangle. A varying water level, see Figure (b), is applied as a boundary condition along the entire length of the left-hand side of the rectangle, causing the water to rise to level 10.35m. This elevation is maintained for long enough for the water to fill the depression and become horizontal over the entire domain. It is then lowered back to its initial state, causing the water level in the pond to become horizontal at the same elevation as the sill, 10.25m. | |||
[[File:Topdown_1_ukbm.png|left|frame|Fig. b: Top down situation]] | [[File:Topdown_1_ukbm.png|left|frame|Fig. b: Top down situation]] | ||
[[File:Waterlevel_1_ukbm.png|left|frame|Fig. c: Water level rise at the left boundary]] | [[File:Waterlevel_1_ukbm.png|left|frame|Fig. c: Water level rise at the left boundary]] | ||
<br style="clear:both"> | <br style="clear:both"> | ||
==Boundary and initial conditions== | ==Boundary and initial conditions== | ||
Varying water level along the dashed red line in Figure (a). Table provided as part of dataset. All other boundaries are closed. | *''Varying water level along the dashed red line in Figure (a). Table provided as part of dataset.'' | ||
Initial condition: Water level elevation = 9.7m. | *''All other boundaries are closed.'' | ||
*''Initial condition: Water level elevation = 9.7m.'' | |||
==Parameter values== | ==Parameter values== | ||
*Manning’s n: 0.03 (uniform) | *Manning’s n: 0.03 (uniform) | ||
*Model grid resolution: 10m (or 700 nodes in the area | *Model grid resolution: 10m (or 700 nodes in the area modeled) | ||
* Time of end: the model is to be run until time t = 20 hours | * Time of end: the model is to be run until time t = 20 hours | ||
==Technical setup== | ==Technical setup== | ||
====Heightmap==== | |||
The provided ascii height file named test1DEM.asc is first imported. It has a cell size of 2m, while the test is expected to run on a 10m grid. Therefore, it will be automatically rescaled by the [[Grid Rasterizer|grid rasterizer]]. The resulting rescaled asc file is packed in the {{anchor|test case zip}} down below. | The provided ascii height file named test1DEM.asc is first imported. It has a cell size of 2m, while the test is expected to run on a 10m grid. Therefore, it will be automatically rescaled by the [[Grid Rasterizer|grid rasterizer]]. The resulting rescaled asc file is packed in the {{anchor|test case zip}} down below. | ||
====Overlay settings==== | |||
General parameters that have been used: | |||
* Type of overlay: [[Flooding_(Overlay)|Flooding Overlay]] | |||
* [[Weather]]: | |||
** Evaporation: 0 mm | |||
** Rain: 0 mm | |||
** Simulate for: 20 hours. | |||
* Groundwater: 0 (off) | |||
* Calculation mode: Accuracy | |||
* Quad cell: 0 | |||
====Inlets==== | |||
In order to regulate the water level according to the water level graph, we used the following setup: | In order to regulate the water level according to the water level graph, we used the following setup: | ||
On gridcells with x = 1 and x = 2 Inlet objects were placed. Each inlet had its own grid cell. The inlets were configured as: | On gridcells with x = 1 and x = 2 Inlet objects were placed. Each inlet had its own grid cell. The inlets were configured as: | ||
[[File:inletpositions_case1_ukbm.png| | [[File:inletpositions_case1_ukbm.png|left]] | ||
<br/><br/><br/><br/><br/><br/><br/><br/> | |||
* External area (m2): 1 000 000 000; | * External area (m2): 1 000 000 000; | ||
* Water level (m): 1; | * Water level (m): 1; | ||
* Inlet Q (m) : | * Inlet Q (m) : | ||
{|class="wikitable" | {|class="wikitable" | ||
! time (hrs) | ! time (hrs) | ||
| Line 47: | Line 64: | ||
* Threshold (m): dynamic; | * Threshold (m): dynamic; | ||
The threshold array was set using a .csv-file that approximated there water level over time. It is approximated since our water levels remain the same for a given period of time (no linear interpolation between time key - values). However, this setup proved accurate enough. The ''inletthresholds.csv''-file can be found in the zip-file below. | The threshold array was set using a .csv-file that approximated there water level over time. It is approximated since our water levels remain the same for a given period of time (no linear interpolation between time key - values). However, this setup proved accurate enough. The ''inletthresholds.csv''-file can be found in the zip-file below. | ||
[[File:thresholds_case1_ukbm.png]] | [[File:thresholds_case1_ukbm.png|left]] | ||
<br style="clear:both"> | <br style="clear:both"> | ||
| Line 61: | Line 78: | ||
*Grid resolution: 10 m. | *Grid resolution: 10 m. | ||
*Simulation time: 28s. | *Simulation time: 28s. | ||
*Remaining volume water: | *Remaining volume water: 10050 m3 | ||
Measured point graphs are displayed below: | Measured point graphs are displayed below: | ||
| Line 77: | Line 94: | ||
==Notes== | ==Notes== | ||
* Tests are run with multi gpu setup. For small cases like this, running it on a single gpu is actually faster: 17 seconds, which is +- 40% compared to 2. | * Tests are run with multi gpu setup. For small cases like this, running it on a single gpu is actually faster: 17 seconds, which is +- 40% compared to 2. | ||
==References== | |||
<references> | |||
<ref name="Néelz">Néelz, S., & Pender, G. (2013). Benchmarking the latest generation of 2D hydraulic modelling packages. Report: SC120002. Environment Agency, Horison House, Deanery Road, Bristol, BS1 9AH. ISBN: 978-1-84911-306-9. Retrieved from: https://www.gov.uk/government/publications/benchmarking-the-latest-generation-of-2d-hydraulicflood-modelling-packages </ref> | |||
</references> | |||
[[Category:Benchmarks]] | |||
{{Water Module buttons}} | |||
Latest revision as of 15:08, 5 March 2024
This benchmark performs test case 1 of the UK benchmark [1] named Test 1 – Uncovering of a beach as well as its result generated by the Water Module in the Tygron Platform.
This test consists of a sloping topography with a depression. An inflow boundary condition is applied at the low end, causing the water to rise to the level indicated by a thick blue line. The inflow is then replaced by a sink term until the water level becomes as indicated by the thin blue line. A similar test has been used by EDF (2000) for validation of the TELEMAC package. The aim of the test is to assess basic package capabilities such as handling disconnected water bodies and wetting and drying of floodplains.
Description
Fig. a: Sloping topography with depression
This test consists of a sloping topography with a depression as illustrated in Figure (a). The modeled domain is a perfect 700m x 100m rectangle. A varying water level, see Figure (b), is applied as a boundary condition along the entire length of the left-hand side of the rectangle, causing the water to rise to level 10.35m. This elevation is maintained for long enough for the water to fill the depression and become horizontal over the entire domain. It is then lowered back to its initial state, causing the water level in the pond to become horizontal at the same elevation as the sill, 10.25m.
Fig. b: Top down situation
Fig. c: Water level rise at the left boundary
Boundary and initial conditions
- Varying water level along the dashed red line in Figure (a). Table provided as part of dataset.
- All other boundaries are closed.
- Initial condition: Water level elevation = 9.7m.
Parameter values
- Manning’s n: 0.03 (uniform)
- Model grid resolution: 10m (or 700 nodes in the area modeled)
- Time of end: the model is to be run until time t = 20 hours
Technical setup
Heightmap
The provided ascii height file named test1DEM.asc is first imported. It has a cell size of 2m, while the test is expected to run on a 10m grid. Therefore, it will be automatically rescaled by the grid rasterizer. The resulting rescaled asc file is packed in the down below.
Overlay settings
General parameters that have been used:
- Type of overlay: Flooding Overlay
- Weather:
- Evaporation: 0 mm
- Rain: 0 mm
- Simulate for: 20 hours.
- Groundwater: 0 (off)
- Calculation mode: Accuracy
- Quad cell: 0
Inlets
In order to regulate the water level according to the water level graph, we used the following setup: On gridcells with x = 1 and x = 2 Inlet objects were placed. Each inlet had its own grid cell. The inlets were configured as:
- External area (m2): 1 000 000 000;
- Water level (m): 1;
- Inlet Q (m) :
| time (hrs) | inlet Q (m3/s) |
|---|---|
| 0, | 1000, |
| 11 | 1000 |
| 11 + 1 sec | -1000 |
| 22 | -1000 |
- Threshold (m): dynamic;
The threshold array was set using a .csv-file that approximated there water level over time. It is approximated since our water levels remain the same for a given period of time (no linear interpolation between time key - values). However, this setup proved accurate enough. The inletthresholds.csv-file can be found in the zip-file below.
Output as required
- Software package used: Tygron Platform
- Numerical scheme: FV (Kurganov, Bollerman, Horvath)*
- Specification of hardware used to undertake the simulation:
- Processor: Intel Xeon @2.10GHz x 8,
- RAM 62.8 GiB,
- GPU: 2x NVidia 1080
- Operating system: Linux 4.13
- Time increment used: adaptive: 0.25 - 1.13 s.
- Grid resolution: 10 m.
- Simulation time: 28s.
- Remaining volume water: 10050 m3
Measured point graphs are displayed below:
- Graph of point 1 generated by the Tygron Platform on the left and generated by others on the right.
- Graph of point 1 generated by the Tygron Platform on the left and generated by others on the right.
- Slice of the situation in the x -direction, generated by the Tygron Platform.
Notes
- Tests are run with multi gpu setup. For small cases like this, running it on a single gpu is actually faster: 17 seconds, which is +- 40% compared to 2.
References
- ↑ Néelz, S., & Pender, G. (2013). Benchmarking the latest generation of 2D hydraulic modelling packages. Report: SC120002. Environment Agency, Horison House, Deanery Road, Bristol, BS1 9AH. ISBN: 978-1-84911-306-9. Retrieved from: https://www.gov.uk/government/publications/benchmarking-the-latest-generation-of-2d-hydraulicflood-modelling-packages
