HydroGeoSphere/Newton Iteration Parameters
The following parameters can be used to control the Newton-Raphson iteration scheme for solution of the variably-saturated flow problem as described in Section 3.13.2.
Newton maximum iterationsEdit
- maxnewt Maximum number of Newton iterations.
Assigns a new value for the maximum number of Newton iterations, which defaults to 15. If this number is exceeded during a time step, the current time step value is reduced by half and a new solution is attempted.
- • • •
- • • •
Jacobian epsilonEdit
- epsilon Jacobian epsilon.
Assigns a new value for the Jacobian epsilon, which defaults to 1 × 10^{−4}. The Jacobian epsilon is the shift in pressure head used to numerically compute the derivatives in the Jacobian matrix. As a rule of thumb, a value equal to 10^{−5} times the average pressure head in the domain is recommended.
- • • •
- • • •
Newton absolute convergence criteriaEdit
- delnewt Newton absolute convergence criteria.
Assigns a new value for the Newton absolute convergence criteria, which defaults to 1 × 10^{−5}. Convergence of the solution occurs when the maximum absolute nodal change in pressure head over the domain for one Newton iteration is less than this value.
- • • •
- • • •
Newton residual convergence criteriaEdit
- resnewt Newton residual convergence criteria.
Assigns a new value for the Newton residual convergence criteria, which defaults to 1 × 10^{−8}. Convergence of the solution occurs when the maximum absolute nodal residual (see Section 5.5.3) in the domain for one Newton iteration exceeds this value.
- • • •
- • • •
Newton maximum update for headEdit
- NR_dhtol Newton maximum update for head.
Assigns a new value for the Newton maximum update for head, which defaults to 1.0. This is used to calculate the underrelaxation factor in such a way that:
- (Equation 5.6)
- (Equation 5.7)
where is the computed maximum update for head in the Newton iteration, and is the head flow solution after iterations. As becomes smaller, the Newton solution becomes more stable but with possibly more iterations. For highly nonlinear problems for which Newton linearization easily fails to converge, it is recommended to set this value smaller.
- • • •
- • • •
Newton maximum update for depthEdit
- NR_ddtol Newton maximum update for depth.
Assigns a new value for the Newton maximum update for water depth, which defaults to 1 × 10^{−2}. The same as above but applies only to water depth.
- • • •
- • • •
Newton maximum residual increaseEdit
- NR_resnorm_fac Newton maximum residual increase.
Assigns a new value for the Newton maximum residual increase, which defaults to 1 × 10^{30}.
If the Newton maximum residual increases by more than , the Newton loop is restarted with a smaller time step.
- • • •
- • • •
Remove negative coefficientsEdit
Forces negative inter-nodal conductances to zero. Negative inter-nodal conductances result in inter-nodal flow from lower to higher heads and can cause oscillatory behavior during Newton iterations [Letniowski and Forsyth, 1991].
- • • •
- • • •
Nodal flow check toleranceEdit
- n_flow_check_tol Assigns a new value for the nodal flow check tolerance, which defaults to 1 × 10^{−2}.
This checks that the relative local (nodal) fluid mass balance is acceptable for all nodes:
Absolute fluid mass balance can always be satisfied against the global convergence criteria, if local inflows and outflows, and mass accumulation are very small.
and where , , and (mass accumulation) are much less than 1.0. This can deteriorate the transport solution, as the concentration is defined as solute mass per unit fluid mass.
- • • •
- • • •
No nodal flow checkEdit
Turns off the nodal flow check feature. In cases where a transport solution is not required, the nodal flow check is not necessary.
- • • •
- • • •
Underrelaxation factorEdit
- under_rel Underrelaxation factor.
Assigns a new value for the underrelaxation factor for the Newton iteration, which defaults to 1. This value can range from 0 (full underrelaxation) to 1 (no underrelaxation).
- • • •
- • • •
Compute underrelaxation factorEdit
Causes the underrelaxation factor ω to be computed according to the following method described by Cooley [1983]:
- (Equation 5.8)
where:
- (Equation 5.9)
In the equations presented above, and represent the previous and current iteration level, and represent the underrelaxation factor for the previous and current iteration levels, and represents the maximum value of the largest difference between head values for 2 successive iterations, .
- • • •
- • • •
Compute underrelaxation factor limitEdit
- dellim Upper limit on the computed underrelaxation factor.
Assigns a new value for the upper limit on the computed underrelaxation factor, which defaults to 1000. A suggested value is 10 times the system domain thickness.
- • • •
- • • •
Minimum relaxation factor allowedEdit
- min_relfac_allowed Minimum value allowed for computed underrelaxation factor.
Assigns a new value for the lower limit on the computed underrelaxation factor, which defaults to 0.001. If not the first timestep, and the computed underrelaxation factor is less than this value, the current timestep is cut in half and the Newton Raphson iteration loop is re-started.
- • • •
- • • •
Newton informationEdit
Causes HydroGeoSphere to write more detailed information to the listing file about the performance of the Newton iteration process.
- • • •
- • • •