Acceleration of nonlinear solvers for natural convection problems

by   Sara Pollock, et al.

This paper develops an efficient and robust solution technique for the steady Boussinesq model of non-isothermal flow using Anderson acceleration applied to a Picard iteration. After analyzing the fixed point operator associated with the nonlinear iteration to prove that certain stability and regularity properties hold, we apply the authors' recently constructed theory for Anderson acceleration, which yields a convergence result for the Anderson accelerated Picard iteration for the Boussinesq system. The result shows that the leading term in the residual is improved by the gain in the optimization problem, but at the cost of additional higher order terms that can be significant when the residual is large. We perform numerical tests that illustrate the theory, and show that a 2-stage choice of Anderson depth can be advantageous. We also consider Anderson acceleration applied to the Newton iteration for the Boussinesq equations, and observe that the acceleration allows the Newton iteration to converge for significantly higher Rayleigh numbers that it could without acceleration, even with a standard line search.



There are no comments yet.


page 15


On the Asymptotic Linear Convergence Speed of Anderson Acceleration, Nesterov Acceleration, and Nonlinear GMRES

We consider nonlinear convergence acceleration methods for fixed-point i...

An efficient nonlinear solver and convergence analysis for a viscoplastic flow model

This paper studies a finite element discretization of the regularized Bi...

Anderson acceleration for contractive and noncontractive operators

A one-step analysis of Anderson acceleration with general algorithmic de...

Anderson Acceleration for Seismic Inversion

The state-of-art seismic imaging techniques treat inversion tasks such a...

Anderson Acceleration as a Krylov Method with Application to Asymptotic Convergence Analysis

Anderson acceleration is widely used for accelerating the convergence of...

Newton-Krylov-BDDC deluxe solvers for non-symmetric fully implicit time discretizations of the Bidomain model

A novel theoretical convergence rate estimate for a Balancing Domain Dec...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.