Cost-optimal adaptive FEM with linearization and algebraic solver for semilinear elliptic PDEs

Abstract

We consider scalar semilinear elliptic PDEs, where the nonlinearity is strongly monotone, but only locally Lipschitz continuous. To linearize the arising discrete nonlinear problem, we employ a damped Zarantonello iteration, which leads to a linear Poisson-type equation that is symmetric and positive definite. The resulting system is solved by a contractive algebraic solver such as a multigrid method with local smoothing. We formulate a fully adaptive algorithm that equibalances the various error components coming from mesh refinement, iterative linearization, and algebraic solver. We prove that the proposed adaptive iteratively linearized finite element method (AILFEM) guarantees convergence with optimal complexity, where the rates are understood with respect to the overall computational cost (i.e., the computational time). Numerical experiments investigate the involved adaptivity parameters.

References in corpus (13)

• Axioms of Adaptivity
• Adaptive FEM with optimal convergence rates for a certain class of non-symmetric and possibly non-linear problems
• On 2D Newest Vertex Bisection: Optimality of Mesh-Closure and H1-Stability of L2-Projection
• Adaptive FEM with coarse initial mesh guarantees optimal convergence rates for compactly perturbed elliptic problems
• Dörfler marking with minimal cardinality is a linear complexity problem
• Rate optimal adaptive FEM with inexact solver for nonlinear operators
• Rate optimality of adaptive finite element methods with respect to the overall computational costs
• Convergence and quasi-optimal cost of adaptive algorithms for nonlinear operators including iterative linearization and algebraic solver
• MooAFEM: An object oriented Matlab code for higher-order adaptive FEM for (nonlinear) elliptic PDEs
• Goal-oriented adaptive finite element method for semilinear elliptic PDEs
• Cost-optimal adaptive iterative linearized FEM for semilinear elliptic PDEs
• Adaptive FEM with quasi-optimal overall cost for nonsymmetric linear elliptic PDEs
• $hp$-robust multigrid solver on locally refined meshes for FEM discretizations of symmetric elliptic PDEs