
Browse by Category
 All Solutions
 ALE / Parameterized Geometry (1)
 Error Messages (59)
 Export (1)
 Fluid Dynamics (1)
 General (11)
 Geometry (11)
 Import (9)
 Installation (33)
 Mesh (12)
 Multiphysics (1)
 Physics (10)
 Postprocessing (4)
 Product Information (1)
 Solver (30)
 Structural Mechanics (2)
Solution Number:  103 
Title:  Improving Convergence of Nonlinear Stationary Models 
Platform:  All Platforms 
Applies to:  All Products 
Versions:  All versions 
Categories:  Error Messages, Error Messages 
Keywords: 
Problem Description
Stationary (timeinvariant) models with nonlinearities may converge very slowly. A nonlinearity can be introduced into the model either in the governing equation, or by making any of the material properties, loads, or boundary conditions dependent upon the solution. Multiphysics problems are often nonlinear. If instead the model is linear, see: Knowledgebase 1260: What to do when a linear stationary model is not solving. This guide applies solely to nonlinear stationary models.
Solution
The issue here has do with the iterative algorithm used to solve nonlinear stationary models. The algorithm is, generally speaking, a Newton's method approach. That is, when solving, the software start with the userspecified initial values to evaluate all solutiondependent terms. The software then computes an initial solution and from there it iteratively recomputes the solution, taking into account how the these intermediate solutions affects the nonlinearities. When the difference in the computed solutions between successive iterations is sufficiently small, or when the residual is sufficiently small, the problem is considered converged to within the specified tolerance. Assuming a wellposed problem, the solver may converge slowly (or not at all) if the initial values are poor, if the nonlinear solver is not able to approach the solution via repeated iterations, or if the mesh is not fine enough to resolve the spatial variations in the solution.
Initial Values
The default Initial Values for the unknowns in most physics interfaces are zero. The exception are the Heat Transfer interfaces, which have a default Initial Value of 293.15K, or 20°C, for the temperature fields. Convergence can be poor when the initial values do not provide a good starting point for this iterative approach. If a good estimate to the solution field is known, this can be entered as an an expression in the Initial Value field. With the exception of some thermal problems however, it is often difficult to estimate the solution, so alternative approaches are needed.
Load Ramping
One can say that , in general, if the loads on a nonlinear system are zero the system will be at rest, that is, the solution will be zero. Therefore, an initial value of zero is almost always reasonable if a very small load is applied. Starting from zero initial conditions, the nonlinear solver will converge if a sufficiently small load is applied. That is, start by forst solving a model with a small, but nonzero, load. From there, if an additional small load increment is applied, the previously computed solution is a reasonable initial condition. Extending this logic, if one wants to solve for any arbitrary load on a nonlinear system, it makes sense to solve a sequence of intermediate problems with gradually increasing load values and using the solutions from each previous step as the initial condition for the next step. This is done within the software by using the Continuation method, which is enabled by default when using the Auxiliary sweep study extension, as shown below.
The Auxiliary Sweep can be used to implement load and nonlinearity ramping.
A Global Parameter has to be introduced (in the above screenshot, P) and is ramped from a value near zero up to one. This parameter is used within the physics interfaces to multiply one, some, or all of the applied loads.
The advantages of this technique are twofold. First, it is physically intuitive, often matching how one would perform an experiment. Second, the continuation method will automatically take smaller load increments if a solution cannot be found. For example, if ramping P
over values of: 0.2,0.4,0.6,0.8,1.0
the nonlinear solver may fail to converge for a value of 0.8
. In that case the continuation method will automatically backtrack and try to solve for intermediate values in the range of 0.6
through 0.8
. This is useful since the software will then return an estimation of the maximum possible loadcase for which the solver can converge.
Nonlinearity Ramping
The technique of load ramping is not always reasonable for all problems. In such cases, use the same continuation method, but instead ramp the nonlinearities in the model. Nonlinearities arise as a consequence of the governing equation, as a material nonlinear expression, or as a coupling term between physics. Again, introduce a Global Parameter that gets ramped from zero to one. Use this parameter to modify the nonlinearity expressions in the model. For example, if there is a temperaturedependent material property such as:
k(T) = 10[W/m/K]*exp((T293[K])/100[K])
replace it with the expression:
k(T,P) = 10[W/m/K]*((1P)+P*exp((T293[K])/100[K]))
At a value of P=0
the above expression is linear, and at a value of P=1
the expression is equal to the original nonlinear expression. As P is ramped up, the continuation method uses the previous solutions as initial conditions for the more nonlinear cases. If you define this nonlinearity ramping such that the first case (P=0)
is a purely linear problem, then you are guaranteed to get a solution for this first step in the ramping.
Nonlinearity ramping is an especially useful technique if any of the nonlinear terms in the model are very abrupt. In the extreme case, suppose one wants to model an instantaneous change in properties, such as:
k(T) = 10[W/m/K]+10[W/m/K]*(T>400[K])
That is, the material property changes instantaneously from 10W/m/K
to 20W/m/K
at 400K
. This case is generally difficult, or impossible, to solve since this material property is nonsmooth. Instead, use a nonlinear material property expression that ramps from a very smooth function to a very nearly discontinuous one. The continuation method will again backtrack and try intermediate values of the ramping parameter, thus giving you the nearest approximation to the abrupt transition that is solvable.
Load ramping and nonlinearity ramping can be used in combination, but start with only one or a few of the loads or nonlinearities being ramped. Check the solver log to see if the continuation method is backtracking. If it does so, use a finer increment in that range.
Mesh Refinement
If both load ramping and nonlinearity ramping are still leading to slow convergence, refine the mesh. The finite element mesh must be fine enough to resolve the spatial variations in the solution fields. Ideally, one would use small elements in regions where the solution varies quickly in space, and larger elements elsewhere. However, it is usually not possible to know this ahead of time. Therefore, it is recommended to use Adaptive Mesh Refinement which will uatomatically refine the mesh only in regions where it is needed, and coarsen the mesh elsewhere. It is also possible to manually refine the mesh. For more details, see: Performing a Mesh Refinement Study
Mesh refinement may often needs to be combined with load or nonlinearity ramping and may require a set of studies, first starting with a relatively coarse mesh for nonlinearity ramping, refining the mesh, and the ramping further on the refined mesh. An example model that combines the techniques of nonlinearity ramping and adaptive mesh refinement with multiple study steps is: Cooling and Solidification of Metal
Solution Approach
There are two approaches that can be used when iteratively solving the nonlinear system of equations: a Fully Coupled or a Segregated approach. The former approach solves for all unknowns in the problem at once, and considers all couplings terms between all unknowns within a single iteration. This is relatively expensive to do, but will lead to the most robust convergence. This approach is used by default for most 1D, 2D, and 2Daxisymmetric models.
The segregated approach, on the other hand, solves sets of unknowns separately. The unknowns are segregated into groups, usually according the physics that they represent, and these groups are solved one after another. That is, within each outer Newtontype iteration, the segregated approach solves for each segregated group sequentially. Each physics is thus solved as a standalone problem, using the solution from any previously computed steps as initial values and linearization points. The couplings terms between the different groups are thus neglected. Despite this, the segregated approach can often converge very robustly, unless there are very strong couplings between the physics in the model. The memory requirements will always be lower than with the fully coupled approach, and the overall solution time can often be lower as well. This segregated approach is used by default for most 3D multiphysics models, and the software will automatically segregate the problem into appropriate groups.
To switch between these solver types, go to the Stationary Solver node within the Study sequence. There will always already be either a Segregated or Fully Coupled feature beneath this. Rightclick on the Stationary Solver node and add either the Segregated or Fully Coupled feature. Within either of these features, it can also be helpful to enable the Results While Solving option, as shown in the screenshot below, to visualize the iterations being taken during the solution.
The Fully Coupled solution approach, with the Plot While Solving enabled.
It is sometimes necessary to manually scale the dependent variables. See Knowledge Base 1240: Manually Setting the Scaling of Variables
The other lowlevel default settings within the Stationary Solver are chosen for robustness. That is, they are tuned to achieve convergence in as many cases as possible. It is quite rare that changing these settings is superior to using a combination of the other techniques in this Knowledgebase, although it is possible to tune these settings to reduce solution time and memory requirements, once a model is already converging. Changes to these lowlevel settings from the defaults will usually be quite model and casespecific.
General Approach to Model Checking
If it is not clear that any of the above strategies are working, it is useful to take a more general approach to verifying the general validity of the model. This involves a systematic reduction in the model complexity. Examine the model and identify all terms that introduce nonlinearities, such as multiphysics couplings, nonlinear materials relationships, and nonlinear boundary conditions. With respect to any nonlinearities, replace them by a reasonable linearized term. Repeat this for every nonlinearity of the model. With respect to multiphysics couplings, rather than solving the problem using a fully coupled approach (the default) solve the problem sequentially, with one physics being solved after another. The objective here is to simplify the model to a state where the model will solve, with linear approximations. With sufficient simplification, a model can be reduced to a linear problem, and if this simplified model does not converge, see: What to do when a linear stationary model is not solving.
Once a simplified solvable version of the model has been found, gradually increase the model complexity again, reintroducing nonlinearities and multiphysics couplings. Using this technique systematically, along with the techniques described previously, will usually identify the nonlinearities in the model that are leading to issues. It may also reveal that the model itself is illposed in some way. Sometimes, reducing the model complexity can be quite challenging and it can be better to start from as simple a case as possible and gradually increase the complexity. It is thus always advised to start this procedure with a simplified 2D, or 2Daxisymmetric model.
Converting to a TimeDependent Formulation
If all of the above approaches have been tried and you are certain that the problem itself is wellposed, consider that the nonlinear problem may not, in fact, have a stationary (timeinvariant) solution. An classic example of this is fluid flow around a cylinder with high, but constant, flow rates. At low flow speeds the flow solution will be time invariant, but at higher flow rates there will be vortex shedding, a timevarying change in the flow field behind the cylinder. Such problems must solved in the time domain.
In such cases it will be particularly helpful to ramp the load gradually in time, from consistent initial values. Ramping the nonlinearities over time is not as strongly motivated, but step changes in nonlinearities should be smoothed out throughout the simulation. Use either a very fine mesh throughout the simulation domain or use adaptive mesh refinement.
See also: Knowledge Base 1254: Controlling the Time Dependent solver timesteps
Disclaimer
COMSOL makes every reasonable effort to verify the information you view on this page. Resources and documents are provided for your information only, and COMSOL makes no explicit or implied claims to their validity. COMSOL does not assume any legal liability for the accuracy of the data disclosed. Any trademarks referenced in this document are the property of their respective owners. Consult your product manuals for complete trademark details.