View source for HPC Adaptive Timestepping

=Introduction=

The HPC solver, by default, uses adaptive timestepping to progress through the simulation. The timestep is adjusted so that it complies with the mathematical stability criteria of a 2D SWE explicit solution. In contrast, TUFLOW Classic, being an implicit solution, is not numerically bound by these criteria, and can progress through the simulation using a fixed timestep. Due to the underlying solution scheme, HPC typically uses a smaller timestep than Classic. the following sections outline how TUFLOW HPC automatically selects its adaptive timestep. 

=HPC 2D Timestep=
TUFLOW HPC uses adaptive timestepping to maintain stability. There are three primary processes that determine the maximum timestep that an explicit solution to the Shallow Water Equations can use: 
* Courant Number, Nu <br>
* Wave Celerity Number, Nc <br>
* Diffusion Number, Nd <br>

{|class="wikitable" style="text-align: center;"
!colspan="4" style="background-color:#005581; height:50px; font-weight:bold; color:white;"|HPC 2D Adaptive Timestep Controlling Numbers
|-
|style="width: 11%; height:50px;" | <b>Control Number</b>||style="width: 40%;" | <b>Description</b>||style="width: 18%;" | <b>Expression</b>||style="width: 9%;" | <b>Control Number Limit</b>
|-
|style="height:80px;" | Courant Number (Nu)
|This condition ensures that water entering one side of a 2D cell does not pass through the other side within one timestep. For this to be satisfied, the product of the water velocity (𝑢) and model timestep (∆𝑡) must be less than the cell size (∆𝑥).
|[[FILE: Courant Number Equation.PNG |250px]]
|<1.0
|-
|style="height:110px;" | The Shallow Wave Celerity <br> Number (Nc)
|This numerical condition relates to the shallow water wave celerity (wave speed) and is derived from the fluid flow equations to represent long waves (i.e. wave length is substantially longer than the water depth). The product of the model timestep (∆𝑡) and the long wave speed (square root of the gravity (g) and water depth (h)) must be less than the cell size (∆𝑥), for the condition to be satisfied.
|[[FILE: SWC Number Equation.PNG |250px]]
|<1.0
|-
|style="height:80px;" | Diffusion Number (Nd)
|This numerical condition relates to the sub-grid scale eddy viscosity term which causes diffusion of momentum. To maintain stability the product of the eddy viscosity coefficient (ν_T) and the timestep (∆𝑡) divided by the square of the grid spacing (∆𝑥2) must remain below 0.3.
|[[FILE: Diffusion Number Equation.PNG |250px]]
|<0.3
|}

The method TUFLOW HPC uses to calculate a timestep and achieve unconditional stability is as follows:
* The <font color="blue"><tt>Timestep </tt></font> command included in the TUFLOW Control File (TCF) is only used for the first calculation timestep. The specified value should be consistent with what would be appropriate for a TUFLOW Classic model (i.e. 1/2 to 1/5 the 2D cell size). Internally, TUFLOW HPC divides this value by 10 to apply a value that is suitable for an explicit solution scheme. All subsequent calculations are completed using the adaptive timestep approach outlined in the following bulletpoints.
* The HPC timestep is calculated using the hydraulic conditions from the end of the previous timestep.
* If the hydraulic conditions have changed significantly it is possible for one or more of the Nu, Nc, Nd control number criteria to be violated at one or more locations within the model.  For example, a sudden change in rainfall from one timestep to the next (which occurs with stepped rainfall boundaries) would potentially cause a violation.  The HPC solver, by default, treats a 20% exceedance of a control number as a violation and will implement a repeat timestep feature. HPC uses a repeat timestep feature to maintain unconditional stability. The repeat timestep feature involves retaining the complete hydraulic solution from the previous (good) timestep.  Should a control number anywhere within the model be exceeded by more than 20%, the solution reverts to the retained timestep, the timestep is reduced and then repeated. 
* Each timestep is also tested for the occurrence of NaNs.  A NaN is “Not a Number” and occurs due to undefined mathematical calculations such as a divide by zero or square root of a negative number.  The occurrence of a NaN is also indicative of a sudden instability.  Should a NaN occur, the repeat timestep feature is implemented.<br>
* Should a timestep need to be repeated more than ten times consecutively, the solution stops.  
* The simulation will also stop if the default minimum permissible timestep of 0.1 seconds has been reached. This value can be manually adjusted using <font color="blue"><tt>Timestep Minimum</tt></font>. <br>
Repeated timesteps are displayed to the Console Window and the number of them for a time interval are provided in the nRS_NaNs and nRS_HCNs columns in the _HPC.csv file output in the results folder.  They are also reported to the _messages layer.<br>
Repeated timesteps are an indication the 2D HPC solution is numerically “on-the-edge”.  Models that have a high number of repeated timesteps should be sensitivity tested by reducing the control number limits using <font color="blue"><tt>Control Number Factor</tt></font>.  For example, repeat the simulation using <font color="blue"><tt>Control Number Factor</tt></font><font color="red"><tt>==</tt></font><tt> 0.8 </tt> and compare the results.  If there are acceptably immeasurable changes in the results, then running at the default control number limits can be considered satisfactory.

''NOTE: Models that require an excessively small 2D timestep to remain stable are symptomatic of models with input topography or boundary condition errors. Please refer to [Tutorial_Module03 | Tutorial Module 03] for discussion on this topic.''

=HPC 1D Timesteping=
HPC links with the 1D solver, ESTRY. When run with HPC (instead of Classic) ESTRY has been reconfigured to automatically act as an adaptive/varying timestep solution, and can step at different multiples of steps to the HPC 2D solver.  Both 1D and 2D solutions are always synchronising at the 2D target timestep, or a multiple of the 2D target timestep if the 1D timestep is sufficiently greater for the 2D to perform more than one step.  If the 1D limiting timestep is less than half the 2D target timestep, the 1D proceeds in two or more steps eventually synchronising with the 2D timestep.  Where there is not a one to one synchronisation of the 1D and 2D timesteps, a usually negligible mass error may occur and can be checked by reviewing the CME% values shown on the Console Window, the .tlf file or the _MB.csv file in the same manner as Classic.
The ESTRY 1D <font color="blue"><tt>Timestep </tt></font> for a HPC 1D/2D linked model is the maximum limiting timestep the 1D solver can use.