FEM for elliptic eigenvalue problems: how coarse can the coarsest mesh be chosen? An experimental study

Abstract

In this paper, we consider the numerical discretization of elliptic eigenvalue problems by Finite Element Methods and its solution by a multigrid method. From the general theory of finite element and multigrid methods, it is well known that the asymptotic convergence rates become visible only if the mesh width h is sufficiently small, h ≤ h 0. We investigate the dependence of the maximal mesh width h 0 on various problem parameters such as the size of the eigenvalue and its isolation distance. In a recent paper (Sauter in Finite elements for elliptic eigenvalue problems in the preasymptotic regime. Technical Report. Math. Inst., Univ. Zürich, 2007), the dependence of h 0 on these and other parameters has been investigated theoretically. The main focus of this paper is to perform systematic experimental studies to validate the sharpness of the theoretical estimates and to get more insights in the convergence of the eigenfunctions and -values in the preasymptotic regime.