A posteriori error control and adaptivity for Crank–Nicolson finite element approximations for the linear Schrödinger equation

Theodoros Katsaounis, Irene Kyza*

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

7 Scopus citations

Abstract

We derive optimal order a posteriori error estimates for fully discrete approximations of linear Schrödinger-type equations, in the L (L2)-norm. For the discretization in time we use the Crank–Nicolson method, while for the space discretization we use finite element spaces that are allowed to change in time. The derivation of the estimators is based on a novel elliptic reconstruction that leads to estimates which reflect the physical properties of Schrödinger equations. The final estimates are obtained using energy techniques and residual-type estimators. Various numerical experiments for the one-dimensional linear Schrödinger equation in the semiclassical regime, verify and complement our theoretical results. The numerical implementations are performed with both uniform partitions and adaptivity in time and space. For adaptivity, we further develop and analyze an existing time-space adaptive algorithm to the cases of Schrödinger equations. The adaptive algorithm reduces the computational cost substantially and provides efficient error control for the solution and the observables of the problem, especially for small values of the Planck constant.

Original languageEnglish (US)
Pages (from-to)55-90
Number of pages36
JournalNumerische Mathematik
Volume129
Issue number1
DOIs
StatePublished - Jan 1 2014

Keywords

  • 35Q41
  • 65M15
  • 65M60

ASJC Scopus subject areas

  • Computational Mathematics
  • Applied Mathematics

Fingerprint Dive into the research topics of 'A posteriori error control and adaptivity for Crank–Nicolson finite element approximations for the linear Schrödinger equation'. Together they form a unique fingerprint.

Cite this