Error estimates of a regularized finite difference method for the Logarithmic Schrödinger equation with Dirac delta potential
Abstract
In this paper, we introduce a conservative Crank-Nicolson-type finite difference schemes for the regularized logarithmic Schrödinger equation (RLSE) with Dirac delta potential in 1D. The regularized logarithmic Schrödinger equation with a small regularized parameter $0<\eps \ll 1$ is adopted to approximate the logarithmic Schrödinger equation (LSE) with linear convergence rate $O(\eps)$. The numerical method can be used to avoid numerical blow-up and/or to suppress round-off error due to the logarithmic nonlinearity in LSE. Then, by using domain-decomposition technique, we can transform the original problem into an interface problem. Different treatments on the interface conditions lead to different discrete schemes and it turns out that a simple discrete approximation of the Dirac potential coincides with one of the conservative finite difference schemes. The optimal $H^1$ error estimates and the conservative properties of the finite difference schemes are investigated. The Crank-Nicolson finite difference methods enjoy the second-order convergence rate in time and space. Numerical examples are provided to support our analysis and show the accuracy and efficiency of the numerical method.
- Publication:
-
arXiv e-prints
- Pub Date:
- April 2024
- DOI:
- 10.48550/arXiv.2404.15791
- arXiv:
- arXiv:2404.15791
- Bibcode:
- 2024arXiv240415791Z
- Keywords:
-
- Mathematics - Numerical Analysis;
- 35Q41;
- 65M70;
- 65N35
- E-Print:
- 20 pages, 3 figures