•
|
Automatic (Newton) to let the solver automatically determine a damping factor in each iteration of Newton’s method. Go to Automatic (Newton) and Automatic Highly Nonlinear (Newton) for settings.
|
•
|
Constant (Newton) to manually specify a constant damping factor that is used in all iterations of Newton’s method. Go to Constant (Newton) for settings.
|
•
|
Automatic highly nonlinear (Newton) if the solver does not converge with Automatic (Newton) first. It is similar to Automatic (Newton) but this method can make the solver more careful when solving highly nonlinear problems. Go to Automatic (Newton) and Automatic Highly Nonlinear (Newton) for settings.
|
•
|
For stationary problems, choose Double dogleg to use the double dogleg nonlinear solver. Go to Double Dogleg for settings.
|
•
|
Initial damping factor, to specify a damping factor for the first Newton iteration. The default value is 1 for Automatic (Newton) and 1.0·10−4 for Automatic highly nonlinear (Newton).
|
•
|
Minimum damping factor, to specify the smallest allowed damping factor. The default value is 1.0·10−4 for Automatic (Newton) and 1.0·10−8 for Automatic highly nonlinear (Newton).
|
•
|
Restriction for step-size update, to specify a factor that limits how much the damping factor is allowed to change in a Newton iteration. The damping factor can change up or down by at most this factor. The default is 10.
|
•
|
Restriction for step-size increase, to specify the maximum for the allowed absolute increase in the damping factor for a Newton iteration, as opposed to the Restriction for step-size update value, which limits how much the damping factor can change relative to the current damping factor. The default value of 1 means that this setting is not active. A smaller value imposes a restriction (0–1 is the valid range for this value).
|
The default, Automatic, is equivalent to On for stationary problems and Off for time-dependent problems. For stationary parametric continuation problems, Automatic corresponds to On when solving for the first parameter value and Off when solving for subsequent parameter values.
|
•
|
Enter a value for the Damping factor to specify a constant damping factor for Newton’s method. The default is 1.
|
•
|
With a Time-Dependent Solver, select the Limit on nonlinear convergence rate check box to force the nonlinear solver to terminate as soon as the convergence is estimated to be too slow. The default is 0.9. Enter a limit on the convergence rate in the field as required.
|
•
|
With a Time-Dependent Solver, choose a Jacobian update: Minimal (the default), On every iteration, or Once per time step:
|
-
|
Minimal reuses the Jacobian for several nonlinear systems whenever deemed possible.
|
-
|
On every iteration computes a new Jacobian for all iterations of Newton’s method.
|
-
|
Once per time step computes a new Jacobian on the first iteration of each time step.
|
•
|
With a Stationary Solver or a parametric solver, choose a Jacobian update: Minimal, On every iteration (the default), or On first iteration:
|
-
|
Minimal reuses the Jacobian for several nonlinear systems whenever deemed possible.
|
-
|
On every iteration computes a new Jacobian for all iterations of Newton’s method.
|
-
|
On first iteration computes a new Jacobian on the first iteration of each parameter step.
|
•
|
The Update weights for automatic scales check box, available for stationary and parametric studies, is selected by default. When selected, the solver updates the weights if they change by two orders of magnitude, and then the Newton solver is restarted from the current solution. Clear this check box to turn off the update of weights for automatic scales.
|
•
|
Enter a value for the Initial damping factor, to specify a damping factor for the first Newton iteration. The default value is 1.0·10−4.
|
•
|
Choose a Residual scaling: Field-wise or Uniform. Field-wise scales the equations based on the field-wise sizes of the initial residual. If Uniform is selected, the algorithm terminates on the relative residual based on the initial residual.
|
•
|
The double dogleg solver can restart, if the solver iterations are considered as stagnated. The stagnation is defined by a given number of consecutive iterations that all have a Newton damping factor less than 0.1 and a step size smaller than 0.1 times the tolerance. By default, the restart is active with 7 iterations before the restart. To turn off the restart, clear the Number of iterations before restart check box. The number of iterations before restart must be a positive integer.
|
•
|
Tolerance to terminate the Newton iterations when the estimated relative error is smaller than a specified tolerance. Then enter the Maximum number of iterations to limit the number of Newton iterations. When the maximum number of iterations have been performed, Newton’s method is terminated even if the tolerance is not fulfilled.
|
•
|
Iterations or tolerance to terminate the Newton iterations when the estimated tolerance is smaller than a specified tolerance or after a specified number of iterations, whichever comes first. Then enter the Number of iterations to specify a fixed number of iterations to perform.
|
•
|
If Tolerance or Iterations or tolerance are set as the Termination technique, then enter a Tolerance factor to modify the tolerance used for termination of the Newton iterations. The actual tolerance used is this factor times the value specified in the Relative tolerance field in the General sections of the Stationary Solver and Time-Dependent Solver.
|
•
|
Solution to terminate the Newton iterations on a solution-based estimated relative error.
|
•
|
Residual to terminate the Newton iterations on a residual-based estimated relative error.
|
•
|
Solution or residual to terminate the Newton iterations on the minimum of the solution-based and residual-based estimated relative errors. Enter a scalar Residual factor multiplying the residual error estimate. The default is 1000.
|
•
|
Solution and residual to terminate the Newton iterations on a solution-based estimated relative error and a residual-based estimated relative error (that is, a more stringent termination criterion). Enter a scalar Residual factor multiplying the residual error estimate. The default is 1000.
|
•
|
None (the default) to not use any stabilization or acceleration method.
|
•
|
Pseudo time-stepping to use a pseudo time-stepping method to stabilize convergence toward steady state for a stationary solver. Pseudo time stepping is not available for time-dependent solvers. See Pseudo Time Stepping for more information. For the pseudo time-stepping method, specify the following parameters:
|
-
|
Initial CFL number. The default is 5.
|
-
|
PID controller - proportional. The default is 0.65.
|
-
|
PID controller - integral. The default is 0.05.
|
-
|
PID controller - derivative. The default is 0.05.
|
-
|
Target error estimate. The default is 0.1.
|
-
|
Select the Anderson acceleration check box to activate Anderson acceleration for the pseudo time-stepping method. See below.
|
-
|
Select the Override Jacobian update check box to override updates of the Jacobian. The CFL threshold value (default: 100), which is the value of the CFL number where overriding of the Jacobian update becomes active. From the Jacobian update list, choose On first iteration (the default) or Minimal, which updates the Jacobian at least once and then only when the nonlinear solver fails during parameter stepping. It reuses the Jacobian for several nonlinear systems whenever deemed possible.
|
•
|
Anderson acceleration, which is a nonlinear convergence acceleration method that uses information from previous Newton iterations to accelerate convergence. The Anderson acceleration method is primarily intended for acceleration of nonlinear iterations in transport problems involving, for example, crosswind diffusion stabilization. It is useful is for solving linear or almost linear problems using the segregated solver, where convergence can be improved and the performance increased. You can control the number of iteration increments to store using the Dimension of iteration space field (default: 10) and the mixing factor as a value between 0 and 1 using the Mixing factor field (default: 1.0). The Iteration delay field (default 0) contains the number of iterations between pseudo time stepping becomes inactive and Anderson acceleration becomes active. When used for pseudo time stepping, you can also enter a CFL threshold value (default: 100), which is the value of the CFL number where Anderson acceleration becomes active and pseudo time stepping becomes inactive.
|
For an example using a Stationary Solver, see Tubular Reactor with Nonisothermal Cooling Jacket: Application Library path COMSOL_Multiphysics/Chemical_Engineering/tubular_reactor.
For an example using a Time-Dependent Solver, see Effective Diffusivity in Porous Materials: Application Library path COMSOL_Multiphysics/Diffusion/effective_diffusivity.
|