For advanced models, the automatically damped Newton method might not be robust enough. A pseudo time-stepping algorithm can then be invoked. See Pseudo Time Stepping for Laminar Flow Models.

For small 2D and 3D models, the default solver suggestion is a direct solver. Direct solvers can handle most nonsingular systems and are very robust and also very fast for small models. Unfortunately, they become slow for large models and their memory requirement scales as somewhere between N1.5 and N2, where N is the number of degrees of freedom in the model. The default suggestion for large 2D and 3D models is therefore the iterative GMRES solver. The memory requirement for an iterative solver optimally scales as N.

Multigrid methods need smoothers, but the saddle point character of the linear system restricts the number of applicable smoothers. The choices are further restricted by the anisotropic meshes frequently encountered in fluid-flow problems. The efficiency of the smoothers is highly dependent on the numerical stabilization. Iterative solvers perform at their best when both Streamline Diffusion and Crosswind Diffusion are active.

The default smoother for P1+P1 elements is SCGS. This is an efficient and robust smoother specially designed to solve saddle point systems on meshes that contain anisotropic elements. The SCGS smoother works well even without crosswind diffusion. SCGS can sometimes work for higher-order elements, especially if Method in the SCGS settings is set to Mesh element lines. But there is no guarantee for this, so the default smoother for higher order elements is an SOR Line smoother. SOR Line handles mesh anisotropy but does not formally address the saddle point character. It does, however, function in practice provided that streamline diffusion and crosswind diffusion are both active.

A different kind of saddle point character can arise if the equation system contains ODE variables. Some advanced boundary conditions, for example, Fully Developed Flow (Inlet), can add equations with such variables. These variables must be treated with the Vanka algorithm. SCGS includes an option to invoke Vanka. Models with higher-order elements must apply SCGS or use the Vanka smoother. The latter is the default suggestion for higher-order elements, but it does not work optimally for anisotropic meshes.

BDF methods have been used for a long time and are known for their stability. However, they can have severe damping effects, especially the lower-order methods. Hence, if robustness is not an issue, a model can benefit from using the generalized-α method instead. Generalized-α is a solver that has properties similar to those of the second-order BDF solver but it is much less diffusive.

Both BDF and generalized-α are per default set to automatically adjust the time step. While this works well for many models, extra efficiency and accuracy can often be gained by specifying a maximum time step. It is also often beneficial to specify an initial time step to make the solver progress smoothly in the beginning of the time series.