Abstract
In this work, we consider the numerical integration of the nonlinear Dirac equation and the Dirac-Poisson system (NDEs) under rough initial data. We propose an ultra low-regularity integrator (ULI) for solving the NDEs which enables optimal first-order time convergence in Hr for solutions in Hr, i.e., without requiring any additional regularity on the solution. In contrast to classical methods, a ULI overcomes the numerical loss of derivatives and is therefore more efficient and accurate for approximating low regular solutions. Convergence theorems and the extension of a ULI to second order are established. Numerical experiments confirm the theoretical results and underline the favourable error behaviour of the new method at low regularity compared to classical integration schemes.
| Original language | English |
|---|---|
| Pages (from-to) | 189-214 |
| Number of pages | 26 |
| Journal | Mathematics of Computation |
| Volume | 90 |
| Issue number | 327 |
| Early online date | 7 Aug 2020 |
| DOIs | |
| Publication status | Published - Jan 2021 |
Keywords
- Dirac-Poisson system
- exponential-type integrator
- low regularity
- Nonlinear Dirac equation
- optimal convergence
- splitting schemes.
ASJC Scopus subject areas
- Algebra and Number Theory
- Computational Mathematics
- Applied Mathematics