An iterative adaptive finite element method for elliptic eigenvalue problems

We consider the task of resolving accurately the nth eigenpair of a generalized eigenproblem rooted in some elliptic partial differential equation (PDE), using an adaptive finite element method (FEM). Conventional adaptive FEM algorithms call a generalized eigensolver after each mesh refinement step. This is not practical in our situation since the generalized eigensolver needs to calculate n eigenpairs after each mesh refinement step, it can switch the order of eigenpairs, and for repeated eigenvalues it can return an arbitrary linear combination of eigenfunctions from the corresponding eigenspace. In order to circumvent these problems, we propose a novel adaptive algorithm that only calls a generalized eigensolver once at the beginning of the computation, and then employs an iterative method to pursue a selected eigenvalue-eigenfunction pair on a sequence of locally refined meshes. Both Picard's and Newton's variants of the iterative method are presented. The underlying partial differential equation (PDE) is discretized with higher-order finite elements (hp-FEM) but the algorithm also works for standard low-order FEM. The method is described and accompanied with theoretical analysis and numerical examples. Instructions on how to reproduce the results are provided.

[1]  N. Lalor,et al.  The prediction of low- and mid-frequency internal road vehicle noise: A literature survey , 2007 .

[2]  Chao Yang,et al.  ARPACK users' guide - solution of large-scale eigenvalue problems with implicitly restarted Arnoldi methods , 1998, Software, environments, tools.

[3]  G. Strang,et al.  An Analysis of the Finite Element Method , 1974 .

[4]  David Andrs,et al.  Monolithic discretization of linear thermoelasticity problems via adaptive multimesh hp-FEM , 2010, J. Comput. Appl. Math..

[5]  W. Hackbusch.  Elliptic Differential Equations , 1992 .

[6]  Stefano Giani,et al.  Adaptive finite element methods for computing band gaps in photonic crystals , 2012, Numerische Mathematik.

[7]  Paul Houston,et al.  Adaptive Discontinuous Galerkin Methods for Eigenvalue Problems Arising in Incompressible Fluid Flows , 2009, SIAM J. Sci. Comput..

[8]  I. Babuska,et al.  Finite element-galerkin approximation of the eigenvalues and Eigenvectors of selfadjoint problems , 1989 .

[9]  P. Berini.  Plasmon-polariton waves guided by thin lossy metal films of finite width: Bound modes of asymmetric structures , 2000 .

[10]  Luka Grubisic,et al.  On estimators for eigenvalue/eigenvector approximations , 2009, Math. Comput..

[11]  Donald G. M. Anderson.  Iterative Procedures for Nonlinear Integral Equations , 1965, JACM.

[12]  I. Doležel,et al.  Higher-Order Finite Element Methods , 2003 .

[13]  Stefano Giani,et al.  A Convergent Adaptive Method for Elliptic Eigenvalue Problems , 2009, SIAM J. Numer. Anal..

[14]  Steven G. Johnson,et al.  Photonic Crystals: Molding the Flow of Light , 1995 .