View source for HPC Adaptive Timestepping

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 larger, fixed, timestep. Therefore, HPC typically uses a smaller timestep than Classic. The HPC timestep needs to conform to the limiting timestep as determined by three stability criteria or Control Numbers: <br>

*Courant Number, Nu <br>
*Wave Celerity Number, Nc <br>
*Diffusion Number, Nd <br>

Prior to advancing in time, the limiting timestep for all wet 2D cells is determined, and the minimum of the limiting timesteps for all cells is selected as the next timestep. The three stability criteria (Control Numbers) are summarised in Table 1 below.<br>

{|class="wikitable" style="text-align: center;"
!colspan="4" style="background-color:#005581; height:50px; font-weight:bold; color:white;"|Table 1 - 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 explicit solution scheme adopted for HPC requires all three conditions be satisfied under the limits specified. As seen in the above Table, each condition is relative to the timestep and cell size, as well as some hydraulic behaviour (velocity, depth, turbulence). <br>