SMS:ADH Model Control Iterations

From XMS Wiki
Revision as of 20:55, 16 November 2012 by Jcreer (talk | contribs)
Jump to navigationJump to search

This window is part of the ADH model control dialog (ADH|Model Control... command). It is accessed from the Model Control window by clicking on the Iterations tab.

Iterations

The iterations tab allows the user to specify the level of precision for the conservation of mass and momentum as ADH performs its calculations.

ADH solves conservation statements concerning water volume, momentum, and constituent mass. ADH is written in conservative form and thus can be regarded as both a finite element method and also a finite volume method. As such ADH can be shown to represent a sum of fluxes around the edge of an element to be in balance with the mass or volume change within the element. In the case of the momentum equation it will be a sum of fluxes and forces balancing the momentum change within an element.

ADH computes the water levels and velocities at each computation point in the domain for each time step in the simulation in an iterative process. At the end of each iteration, calculations are made to determine how much the solution is changing and what the errors in conservation are. If the errors are small, and the solution is not changing, the process is said to be converged and calculations can proceed to the next time step. If the solution is changing and/or the errors in conservation are large, ADH can attempt another iteration.

There are several iteration parameter controls in ADH. These are described in detail in the ADH users manual. The Solver options tab controls the linear iterations. This dialog controls the non-linear iteration parameters.

At the end of each iteration ADH computes the conservation properties for the current approximations of water depth and velocity. If the errors are less than specified tolerance, the solver moves on to the next time step. If the errors are greater than the specified tolerance, but still reasonable, and the maximum number of iterations specified on the IP NIT card has not been reached, ADH will attempt another iteration. If the errors have become unreasonable (divergence), ADH will exit the loop with a failure notification. The controls in this dialog allow the user to define how ADH should determine convergence and how to proceed when convergence is not reached.

Maximum number of iterations per time step (IP NIT and IP FNI cards)

In this edit field the user specifies the maximum number of iterations that ADH will try for a single time step. By examining the output for successive iterations, a user can determine if the process is converging in general (trending towards zero), or oscillating. If the trend is towards convergence, this maximum number of iterations can be increased. As the value increases, ADH can spend more time processing individual time steps which may decrease the overall efficiency of the solver. For this reason, a user should not simply increase this parameter without verifying that it will be beneficial.

SMS writes either the IP NIT or IP FNI card based on the setting of the If tolerance is not satisfied radio group.

  • IP NIT card is written if the Reduce timestep option is selected. (Recommended)
  • IP FNI card is written if the Accept solution option is selected.

The IP FNI card instructs ADH to continue with the solution as if it had converged. This is not recommended.

When the Reduce timestep option is selected and the solution has diverged or the maximum number of iterations have been attempted, ADH will reduce the time step size to ¼ of the previous and continue the calculations.

When the Accept solution option is selected, and the maximum number of iterations have been attempted, ADH proceeds as if convergence has been reached to the next time step. (This option is not recommended.)

NTL/ITL/Both combo box

NTL (IP NTL card)

If the NTL option is selected, ADH will check for non-linear convergence on the "maximum residual norm". An initial estimate for an appropriate tolerance value for a problem can be computed as:

This is the maximum allowance, so you will likely want to set the tolerance higher than this estimate so your simulation will progress.

ITL (IP ITL card)

If the ITL option is selected, ADH will check for non-linear convergence on the "maximum increment norm" or the change in the solution (velocity, depth and concentration).

Both

If the Both option is selected, SMS will write both cards to the BC file and ADH will check for convergence on both terms. Both must be satisfied for the solution to proceed.


Go to SMS:ADH page