Nonlinear Solver
Default Nonlinear Settings for Heat Transfer Interfaces
Nonlinear solver settings depend on the heat transfer model and on the study type.
Fully Coupled Solver Attribute
Heat transfer models use a fully coupled nonlinear solver attribute by default. The Jacobian update is set to On every iteration for Stationary studies, and to On first iteration for Time-Dependent studies. A Newton nonlinear method is set by default with
•
Automatic damping factor computation for stationary studies
•
Constant damping factor for time-dependent studies, with a default damping factor set to 0.9.
Segregated Solver Attribute
The segregated solver attribute is set by default in the following cases:
•
Another physics interface is solved together with heat transfer. The dependent variables of the Heat Transfer interface are placed in a separate segregated group.
•
Radiation in participating media using the Discrete ordinates method defines a large number of dependent variables (up to 168 for a single wavelength), which are placed in segregated groups. The number of dependent variables per segregated group and the nonlinear method settings depend on the Performance index parameter available in the heat transfer interface settings in the Participating Media Settings section. When multiple wavelength are considered, the variables relative to distinct wavelengths are not mixed together in the segregated groups.
•
The Thermal Damage subfeature (added under Biological Tissue feature) defines an additional variable alpha that is placed in a dedicated segregated group.
Default Nonlinear Settings for Moisture Transport Interfaces
Nonlinear solver settings depend on the study type.
Fully Coupled Solver Attribute
A fully coupled solver attribute is set by default. A Newton nonlinear method with a constant damping factor is set by default, with the following settings:
•
The Damping factor is set to 0.7.
•
The Jacobian update is set to On first iteration for Time-Dependent studies, and to On every iteration for Stationary studies or when the liquid concentration on surfaces, cl_evap, is solved.
•
The Maximum number of iterations is set to 5 for Time-Dependent studies, and to 50 for Stationary studies.
•
The Tolerance factor is set to 0.01 for Time-Dependent studies, and to 0.1 for Stationary studies.
Segregated Solver Attribute
The segregated solver attribute is set by default when another physics interface is solved together with moisture transport. The dependent variables of the Moisture Transport interface are placed in a separate segregated group.
A Newton nonlinear method with a constant damping factor is set by default, with the following settings:
•
The Damping factor is set to 0.7.
•
The Jacobian update is set to On first iteration for Time-Dependent studies, and to On every iteration for Stationary studies or when the liquid concentration on surfaces, cl_evap, is solved.
•
The termination technique is Iterations, with Number of iterations set to 2.
Default Nonlinear Settings for The Lumped Thermal System Interface
A fully coupled solver attribute is set by default. A Newton nonlinear method with a constant damping factor is set by default, with the following settings:
•
The Damping factor is set to 1.
•
The Jacobian update is set to On every iteration.
•
The Maximum number of iterations is set to 5 for Time-Dependent studies, and to 50 for Stationary studies.
Tuning the nonlinear solver
Default solver settings are defined to handle efficiently classical configurations. For particular applications, the default settings may need modifications to improve the robustness and performance of the solver.
Optimize Nonlinear Solver for Robustness
When the nonlinear solver fails or converges erratically, different options can be considered:
•
Using the Automatic highly nonlinear (Newton) option forces to start the computation with a very low damping factor and increases it carefully. Alternatively a low constant damping factor can be used. The damping factor ranges between 0 and 1. A constant damping factor equal to 0.1 is a very low value and should be robust but slow to converge. For low values of the damping factor, it is thus usually needed to increase the number of nonlinear iterations. If the nonlinear solver is unstable with such a damping factor then the automatic option should be used because it makes it possible to start with a lower damping factor and gradually increases it.
•
A good initial value, as close as possible from the expected solution and consistent with the boundary conditions, helps to guide the nonlinear solver to a stable physical solution. To do that:
-
Try to ramp the temperature on the boundary from the initial to the desired value by using a auxiliary sweep — for stationary problems — or a time-dependent step function — for time-dependent problems.
-
Use results from a simplified problem, for instance with no temperature dependency, or using a one-way multiphysics coupling, as initial value.
Note that it is sometimes easier to update the boundary conditions than the initial condition to get consistent initial settings (see the Heat Conduction in a Finite Slab model).
•
When it is not possible to provide a good initial value, the segregated solver associated with low damping factors in each segregated step helps to achieve convergence.
•
Forcing the Jacobian update at every iteration ensures that the nonlinear solver iterates using optimal information from the equation system. This is needed when nonlinearities are due to the temperature itself — for example, in case of strong temperature dependency of material properties — or to another variable solved in the same segregated group as the temperature — for example, in natural convection models.
Optimize Convergence Speed
Low convergence can be improved by following ways:
•
Using a constant damping factor equal to 1 for linear problems. The linearity is determined at the beginning of the resolution and indicated in the Log section of the solver window.
•
Providing a good initial value is an asset for computational speed.
•
In the convergence area, the fully coupled solver has a better convergence rate than the segregated solver.
•
Using minimal Jacobian update option avoid to spend time in Jacobian computation. This is suited for linear models and models with mild nonlinearities.