Enhancing multigrid solvers for isogeometric analysis of nonlinear problems using polynomial extrapolation techniques

Abstract: When used to accelerate the convergence of fixed-point iterative methods, such as the Picard method, which is a kind of nonlinear fixed-point iteration, polynomial extrapolation techniques can be very effective. The numerical solution of nonlinear problems is further investigated in this study. Particularly, using multigrid with isogeometric analysis as a linear solver of the Picard iterative method, which is accelerated by applying vector extrapolation techniques, is how we address the nonlinear eigenvalue Bratu problem and the Monge-Amp`ere equation. This paper provides quadratic convergence results for polynomial extrapolation methods. Specifically, a new theoretical result on the correlation between the residual norm and the error norm, as well as a new estimation for the generalized residual norm of some extrapolation methods, are given. We perform an investigation between the Picard method, the Picard method accelerated by polynomial extrapolation techniques, and the Anderson accelerated Picard method. Several numerical experiments show that the Picard method accelerated by polynomial extrapolation techniques can solve these nonlinear problems efficiently.