Abstract:In this article, we study numerical approximation of eigenvalue problems of the Schrödinger operator $-\Delta u + \frac{c^2}{|x|^2}u$. There are three stages in our investigation: We start from a ball of any dimension, in which case the exact solution in the radial direction can be expressed by Bessel functions of fractional degrees. This knowledge helps us to design two novel spectral methods by modifying the polynomial basis to fit the singularities of the eigenfunctions. At the second stage, we move to circular sectors in the two-dimensional setting. Again the radial direction can be expressed by Bessel functions of fractional degrees. Only in the tangential direction some modifications are needed from stage one. At the final stage, we extend the idea to arbitrary polygonal domains. We propose a mortar spectral element approach: a polygonal domain is decomposed into several subdomains with each singular corner including the origin covered by a circular sector, in which origin and corner singularities are handled similarly as in the former stages, and the remaining domains are either a standard quadrilateral/triangle or a quadrilateral/triangle with a circular edge, in which the traditional polynomial based spectral method is applied. All subdomains are linked by mortar elements (note that we may have hanging nodes). In all three stages, exponential convergence rates are achieved. Numerical experiments indicate that our new methods are superior to standard polynomial based spectral (or spectral element) methods and $hp$-adaptive methods. Our study offers a new and effective way to handle eigenvalue problems of the Schrödinger operator including the Laplacian operator on polygonal domains with reentrant corners.

Tailored Finite Point Methods for Solving Singularly Perturbed Eigenvalue Problems with Higher Eigenvalues.

Tailored Finite Point Method for Numerical Solutions of Singular Perturbed Eigenvalue Problems

Asymptotic Analysis and Numerical Method for Singularly Perturbed Eigenvalue Problems.

Tailored Finite Point Method for a Singular Perturbation Problem with Variable Coefficients in Two Dimensions

The Tailored Finite Point Method

A Tailored Finite Point Method for a Singular Perturbation Problem on an Unbounded Domain

Multiscale Tailored Finite Point Method for Second Order Elliptic Equations with Rough or Highly Oscillatory Coefficients

A Parameter-Uniform Tailored Finite Point Method for Singularly Perturbed Linear Ode Systems

Tailored Finite Point Method for Parabolic Problems

Tailored Finite Point Method for Steady-State Reaction-Diffusion Equations

Tailored finite point method based on exponential bases for convection-diffusion-reaction equation.

An Equation Decomposition Method for the Numerical Solution of a Fourth‐order Elliptic Singular Perturbation Problem

A Tailored Finite Point Method for the Helmholtz Equation with High Wave Numbers in Heterogeneous Medium

Fitted finite element methods for singularly perturbed elliptic problems of convection-diffusion type

A Semi-Discrete Tailored Finite Point Method for a Class of Anisotropic Diffusion Problems.

Tailored Finite Point Method for Diffusion Equations with Interfaces on Distorted Meshes

Tailored Finite Point Method for First Order Wave Equation

Local and Parallel Finite Element Discretizations for Eigenvalue Problems

Tailored Finite Point Method for the Interface Problem

Finite Element Analysis of the Oseen eigenvalue problem

Efficient Spectral and Spectral Element Methods for Eigenvalue Problems of Schrödinger Equations with an Inverse Square Potential