Extended local fourier analysis for multigrid optimal smoothing, coarse grid correction, and preconditioning

Multigrid methods are fast iterative solvers for partial di erential equations. Especially for elliptic equations they have been proven to be highly e cient. For problems with nonelliptic and nonsymmetric features--as they often occur in typical real life applications--a rigorous mathematical theory is generally not available. For such situations Fourier smoothing and two-grid analysis can be considered as the main analysis tools to obtain quantitative convergence estimates and to optimize different multigrid components like smoothers or inter-grid transfer operators. In general, it is difficult to choose the correct multigrid components for large classes of problems. A popular alternative to construct a robust solver is the use of multigrid as a preconditioner for a Krylov subspace acceleration method like GMRES. Our contributions to the Fourier analysis for multigrid are two-fold. Firstly we extend the range of situations for which the Fourier analysis can be applied. More precisely, the Fourier analysis is generalized to k-grid cycles and to multigrid as a preconditioner. With a k-grid analysis it is possible to investigate real multigrid effects which cannot be captured by the classical two-grid analysis. Moreover, the k-grid analysis allows for a more detailed investigation of possible coarse grid correction difficulties. Additional valuable insight is obtained by evaluating multigrid as a preconditioner for GMRES. Secondly we extend the range of discretizations and multigrid components for which detailed Fourier analysis results exist. We consider four well-known singularly perturbed model problems to demonstrate the usefulness of the above generalizations: The anisotropic Poisson equation, the rotated anisotropic diffusion equation, the convection diffusion equation with dominant convection, and the driven cavity problem governed by the incompressible Navier Stokes equations. Each of these equations represents a larger class of problems with similar features and complications which are of practical relevance. With the help of the newly developed Fourier analysis methods, a comprehensive study of characteristic difficulties for singular perturbation problems can be performed. Based on the insights from this analysis it is possible to identify remedies resulting in an improved multigrid convergence. The theoretical considerations are validated by numerical test calculations

Practical Fourier analysis for multigrid methods

Before applying multigrid methods to a project, mathematicians, scientists, and engineers need to answer questions related to the quality of convergence, whether a development will pay out, whether multigrid will work for a particular application, and what the numerical properties are. Practical Fourier Analysis for Multigrid Methods uses a detailed and systematic description of local Fourier k-grid (k=1,2,3) analysis for general systems of partial differential equations to provide a framework that answers these questions. This volume contains software that confirms written statements about convergence and efficiency of algorithms and is easily adapted to new applications. Providing theoretical background and the linkage between theory and practice, the text and software quickly combine learning by reading and learning by doing. The book enables understanding of basic principles of multigrid and local Fourier analysis, and also describes the theory important to those who need to delve deeper into the details of the subject