Verified eigenvalue evaluation for the laplacian over polygonal domains of arbitrary shape

Xuefeng Liu, Shin'ichi Oishi

Research output: Contribution to journalArticle

30 Citations (Scopus)

Abstract

The finite element method (FEM) is applied to bound leading eigenvalues of the Laplace operator over polygonal domains. Compared with classical numerical methods, most of which can only give concrete eigenvalue bounds over special domains of symmetry, our proposed algorithm can provide concrete eigenvalue bounds for domains of arbitrary shape, even when the eigenfunction has a singularity. The problem of eigenvalue estimation is solved in two steps. First, we construct a computable a priori error estimation for the FEM solution of Poisson's problem, which holds even for nonconvex domains with reentrant corners. Second, new computable lower bounds are developed for the eigenvalues. Because the interval arithmetic is implemented throughout the computation, the desired eigenvalue bounds are expected to be mathematically correct. We illustrate several computation examples, such as the cases of an L-shaped domain and a crack domain, to demonstrate the efficiency and flexibility of the proposed method.

Original languageEnglish
Pages (from-to)1634-1654
Number of pages21
JournalSIAM Journal on Numerical Analysis
Volume51
Issue number3
DOIs
Publication statusPublished - 2013 Oct 9

Keywords

  • Eigenvalue problem
  • Elliptic operator
  • Finite element method
  • Min-max principle
  • Prager-Synge's theorem
  • Verified computation

ASJC Scopus subject areas

  • Numerical Analysis
  • Computational Mathematics
  • Applied Mathematics

Fingerprint Dive into the research topics of 'Verified eigenvalue evaluation for the laplacian over polygonal domains of arbitrary shape'. Together they form a unique fingerprint.

  • Cite this