TUFLOW 2D Cell Size Selection

Introduction

This page of the TUFLOW Wiki discusses 2D cell size convergence. In well design modelling software such as TUFLOW, cell size convergence refers to the tendency for model results to trend towards a common answer as cell size decreases. This behaviour occurs due to topographic features that influence the hydraulic flow behaviour better approximating reality as resolution increases. The series of creek cross-section images below demonstrate this. As model resolution increases from 20m to 5m the modelled topography progressively matches the real-world geometry more closely.

Whilst high resolution modelling is always desirable, unfortunately it isn't practical in all situations due to the impact increasing resolution has on simulation speed. Increasing resolution, increases the computation load and therfore makes a model run slower. Even though TUFLOW with GPU hardware technology is arguably the fastest simulation software available, the time it takes for a model to run is still an important consideration. The challenge for modellers is knowing what resolution is necessary to achieve results that are fit for purpose with sufficient accuracy whilst also having a run time that's practical (i.e. hours, not days).

Australian Rainfall and Runoff - Two Dimensional Modelling in Urban and Rural Floodplain Guideline provides some recommendations on this topic. It states:

The resolution of a 2D model grid/mesh determines the scale of physical features and flow behaviour that can be modelled for a given study area. Selection of an appropriate resolution is generally driven by a combination of the following factors:

* The scale of topographic and/or flow phenomena to be modelled

* The desired level of detail to be achieved in the model outputs

* The length of event time and consequent run time

* The size of the area of interest

Details of the model schematisation process including resolution aspects are described in Chapter 6. Chapter 7 also highlights the importance of grid/mesh resolutions in achieving manageable run times to maximise calibration outcomes. Table 10-2 (below) provides guidance on levels of model resolution that may be appropriate in certain typical situations.

Modelling Case	Typical 2D Cell Resolution
Flow in Channel	In order to adequately resolve flow in a channel it is desirable to provide at least 5 grid/mesh elements laterally across the channel
Urban Overland	Most urban flood models employ grid/mesh resolutions of 2m to 5m.
Flow in Floodplain	Rural floodplain models typically employ grid/mesh resolutions of 10m and 50m (although resolutions up to 200m have been used) depending on the size of the area to be analysed, the characteristics/dimensions of the floodplain and the desired level of output detail.
Lakes and Estuaries	These situations often include areas of open water where less detail is required than along the water body boundary. Such situations are well suited to a flexible mesh rather than fixed grid based model as the mesh is able to incorporate a change of resolution across the model domain. Element resolutions for these models can span the full range as described above depending on project requirements.
Flow Over and Embankment	Embankments effectively act as weirs in the floodplain context. Many 2D modelling packages have automatic or manually activated corrections that compensate for the error in head loss typically associated with modelling broad-crested weir flow with a 2D scheme. For practical purposes, a single 2D element is generally adequate to represent the impact of a levee, road or railway embankment. The resolution of these elements is generally not a significant limitation on the schematisation of most domains.

This Wiki page uses two test cases to discuss this topic. Quantitative results are presented demonstrating how and when resolution assumptions have a tangible impact on model results. TUFLOW HPC and it's GPU Module have been used for all simulations documented in the following sections. The computer used for the modelling has a NVIDIA GeForce GTX 1080 Ti GPU card.

Test Case 1 - Rural Dam Break

This test case has been sourced from the UK Environment Agency 2D Hydraulic Model Benchmark Test dataset. It is referred to as Test 5 in the Environment Agency (EA) document. The EA designed this test to simulate flood wave propagation down a river valley following the failure of a dam. The valley DEM is ~0.8 km by ~17 km and the valley slopes downstream on a slope of ~0.01 in its upper region, easing to ~0.001 at lower elevations. The model uses a single manning's n value of 0.04 across the entire domain. The model topography and EA reporting points are shown below.

The model has a single inflow at the top of the catchment. The inflow hydrograph that has been used is shown below. Although the inflow only introduces water into the model for 100 minutes the model has been run for a simulation period of 30 hours. This was a requirement in the original EA model benchmark documentation.

The EA benchmark testing originally assumed a 50m cell resolution. For the purpose of this assessment a range of cell sizes have been used to determine the impact changing grid resolution has on the model results. The following grid resolutions were used:

• 10m (33ft) resolution - 188,240 cell count
• 20m (66ft) resolution - 47,080 cell count
• 50m (164ft) resolution - 7,540 cell count
• 100m (328ft) resolution - 1,880 cell count
• 150m (492ft) resolution - 840 cell count
• 200m (656ft) resolution - 480 cell count
• 250m (820ft) resolution - 300 cell count

Test Case 1 - Results

In the interest or keeping this page brief we have focused our result presentation on EA reporting points 4 and 5. Their location furthest downstream in the catchment makes them the most sensitive of all EA 7 reporting points. This is due to any divergence in result caused by poor representation of the upstream topography accumulating, amplifying the result difference associated with a change in cell resolution at these selected locations.

Results are presented below. Although the figures are shown in sequence (from fine to coarse resolution), the results are overlayed on one another moving through the dataset so it is obvious if poor convergence occurs.

Cell Size	Location 4 Result	Location 5 Result	Mesh Resolution Figure
10m
20m
50m
100m
150m
200m
250m

Test Case 1 - Discussion

The graphed results indicate convergence is observed for all cases with a cell resolution equal to or less than 100m. The grid figures presented to the right of the graphs suggest the 100m cell resolution case has approximately (4) four cells laterally across the valley that the dam break flow is contained within. Although the 400-600m wide valley isn't a traditional creek or river channel, the scale of flow associated with the dam break means the valley is behaving like one. Minor topographic features within the valley are too small to have a significant impact on the flood behaviour for flow conditions of this large magnitude. This result trend is consistent with the Australian Rainfall and Runoff (ARR) cell size recommendation for the flow in a channel, "In order to adequately resolve flow in a channel it is desirable to provide at least 5 grid/mesh elements laterally across the channel".

Simulation time can influence model resolution selection. This is not the case in this situation. All models have run in under 5 minutes using TUFLOW HPC's GPU hardware module.

Cell Size	Model Size (Cell Count)	Simulation Run Time (seconds)	Judgement of Convergence (Yes / No)
10m (33ft)	188,240	284	Yes
20m (66ft)	47,080	98	Yes
50m (164ft)	7,540	32	Yes
100m (328ft)	1,880	15	Yes
150m (492ft)	840	10	No
200m (656ft)	480	9	No
250m (820ft)	300	7	No

Since simulation time is not a limiting factor in this case selection of cell size only needs to consider the result accuracy. What is appropriate in terms of accuracy is influenced by the intended use of the model. For example, if lot scale assessment of inundation and structural damage risk is required, selection of the 10m cell size may be appropriate. Alternatively, if the model is only intended to inform broad-scale risk, the 50m or 100m cell size may be sufficient.

Test Case 2 - Urbanised Catchment

This test case uses a hypothetical urban model. The model encompasses the entire catchment and covers an area of approximately 21miles2 (54km2). Development within the catchment ranges from rural undeveloped in the mountainous upper catchment to dense urban development in the lower catchment.

A direct rainfall approach has been used, applying rainfall directly to every cell within the model. The event hyetograph is a hypothetical 24 hour extreme event. The simulation duration has also been set to 24 hours.
Note, if you are unfamiliar with the direct rainfall modelling approach, this Hydrology and Water Symposium paper introduces the concept.

The following grid resolutions were used to test the impact of cell size on the assessment results and simulation time:

• 10ft (3.0m) cell resolution – 5,821,533 total cell count
• 12ft (3.7m) cell resolution – 4,042,960 total cell count
• 15ft (4.6m) cell resolution – 2,587,629 total cell count
• 20ft (6.1m) cell resolution – 1,455,869 total cell count
• 30ft (9.1m) cell resolution – 647,003 total cell count
• 50ft (15.2m) cell resolution – 232,957 total cell count
• 75ft (22.9m) cell resolution – 103,525 total cell count

Since the direct rainfall modelling approach means the entire model is "wet", output has been filtered to only show results in locations where the flood depth exceeds 0.3ft. This is done in TUFLOW using the TCF command "Map Cutoff Depth == ".

Test Case 2 - Results

The cell convergence test results are summarised in the table below. They have been presented in a histogram form. The histograms are calculated by subtracting the peak (maximum) flood level result raster grid for the larger cell resolution simulation from the 10ft cell resolution result. The result difference has been sampled at every pixel in the raster datasets where flood information is present in both result files (ie. they overlap).

Cell Resolution, Model Size and Simulation Run Time	Difference in Water Level Result (ft) (Relative to 10ft Resolution Model)	Sample Peak (Maximum) Water Level Result
10ft (3.0m) cell size 5,821,533 total cell count 20.3hr simulation run time	N/A The histogram graphs below are calculated by subtracting the peak (maximum) flood level result from the 10ft cell resolution peak (maximum) flood level result (shown to the right)
12ft (3.7m) cell size 4,042,960 total cell count 15.0hr simulation run time
15ft (4.6m) cell size 2,587,629 total cell count 9.2hr simulation run time
20ft (6.1m) cell size 1,455,869 total cell count 3.9hr simulation run time
30ft (9.1m) cell size 647,003 total cell count 1.3hr simulation run time
50ft (15.2m) cell size 232,957 total cell count 0.5hr simulation run time
75ft (22.9m) cell size 103,525 total cell count 0.3hr simulation run time

Test Case 2 - Discussion

The urban model testing indicates result convergence occurs when the cell resolution is less than 20ft. This observation agrees with the ARR guideline recommendation for urban models to use a cell resolution of 2m (6.6ft) to 5m (16.4ft).

Cell Size	Model Size (Cell Count)	Simulation Run Time (hours)	Judgement of Convergence (Yes / No)
10ft (3.0m)	5,821,533	20.3hr	Yes
12ft (3.7m)	4,042,960	15.0hr	Yes
15ft (4.6m)	2,587,629	9.2hr	Yes
20ft (6.1m)	1,455,869	3.9hr	Yes
30ft (9.1m)	647,003	1.3hr	Mostly
50ft (15.2m)	232,957	0.5hr	No
75ft (22.9m)	103,525	0.3hr	No

As per Test Case 1, selection of the appropriate cell resolution during a real-world project requires consideration of the intended use of the results and also how the simulation run time impacts the overall project timeline.

• Detailed integrated 1D/2D urban modelling may adopt the 15ft (4.6m) or slightly finer 12ft (3.7m) cell resolution if the cumulative simulation time for all model events/scenarios does not adversely impact the project timeline. This timeframe consideration is particularly relevant in locations such as Australia where an ensemble methodology is currently replacing single event approaches as a way to address uncertainty associated with hydrologic assumptions.
• Broad-scale quick run simulation for urban flood forecasting purposes may adopt a coarser 20ft cell resolution or investigate the suitability of a 25ft cell resolution to achieve a quicker simulation turnaround time.
  

The run time results from this test case highlight a useful relationship used by modellers who wish to estimate likely simulations times for models that use a higher resolution than that presently simulated. The simulation time trend is as follows:

• A factor 1.5 increase in model resolution (eg. increasing the resolution from 30ft to 20ft) will increase the simulation time by a factor of 4.5.
• A factor 2.0 increase in model resolution (eg. increasing the resolution from 30ft to 15ft) will increase the simulation time by a factor of 8.0.
• A factor 3.0 increase in model resolution (eg. increasing the resolution from 30ft to 10ft) will increase the simulation time by a factor of 16.0.

This relationship is due to doubling the resolution (halving the cell size) increasing the number of cells by a factor (4) four and requiring half the timestep to maintain stability (i.e. 2 two times the number of calculations) (4x2=8). The ratios for this test case don't match the above generalised trend exactly, though are close enough to be useful when optimising the cell size to meet project requirements.
700px