Optimal error bounds on an exponential wave integrator Fourier spectral method for the logarithmic Schrödinger equation
Fuente:
arXiv
Guardado en:
| Autores principales: | , , |
|---|---|
| Formato: | Preprint |
| Publicado: |
2024
|
| Materias: | |
| Acceso en línea: | |
| Etiquetas: |
Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
|
| _version_ | 1866911128749080576 |
|---|---|
| author | Bao, Weizhu Ma, Ying Wang, Chushan |
| author_facet | Bao, Weizhu Ma, Ying Wang, Chushan |
| contents | We prove a nearly optimal error bound on the exponential wave integrator Fourier spectral (EWI-FS) method for the logarithmic Schrödinger equation (LogSE) under the assumption of $H^2$-solution, which is theoretically guaranteed. Subject to a CFL-type time step size restriction $τ|\ln τ| \lesssim h^2/|\ln h|$ for obtaining the stability of the numerical scheme affected by the singularity of the logarithmic nonlinearity, an $L^2$-norm error bound of order $O(τ|\ln τ|^2 + h^2 |\ln h|)$ is established, where $τ$ is the time step size and $h$ is the mesh size. Compared to the error estimates of the LogSE in the literature, our error bound either greatly improves the convergence rate under the same regularity assumptions or significantly weakens the regularity requirement to obtain the same convergence rate. Moreover, our result can be directly applied to the LogSE with low regularity $L^\infty$-potential, which is not allowed in the existing error estimates. Two main ingredients are adopted in the proof: (i) an $H^2$-conditional $L^2$-stability estimate, which is established using the energy method to avoid singularity of the logarithmic nonlinearity, and (ii) mathematical induction with inverse inequalities to control the $H^2$-norm of the numerical solution. Numerical results are reported to confirm our error estimates and demonstrate the necessity of the time step size restriction imposed. We also apply the EWI-FS method to investigate soliton collisions in one dimension and vortex dipole dynamics in two dimensions. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2412_16902 |
| institution | arXiv |
| publishDate | 2024 |
| record_format | arxiv |
| spellingShingle | Optimal error bounds on an exponential wave integrator Fourier spectral method for the logarithmic Schrödinger equation Bao, Weizhu Ma, Ying Wang, Chushan Numerical Analysis 35Q55, 65M15, 65M70, 81Q05 We prove a nearly optimal error bound on the exponential wave integrator Fourier spectral (EWI-FS) method for the logarithmic Schrödinger equation (LogSE) under the assumption of $H^2$-solution, which is theoretically guaranteed. Subject to a CFL-type time step size restriction $τ|\ln τ| \lesssim h^2/|\ln h|$ for obtaining the stability of the numerical scheme affected by the singularity of the logarithmic nonlinearity, an $L^2$-norm error bound of order $O(τ|\ln τ|^2 + h^2 |\ln h|)$ is established, where $τ$ is the time step size and $h$ is the mesh size. Compared to the error estimates of the LogSE in the literature, our error bound either greatly improves the convergence rate under the same regularity assumptions or significantly weakens the regularity requirement to obtain the same convergence rate. Moreover, our result can be directly applied to the LogSE with low regularity $L^\infty$-potential, which is not allowed in the existing error estimates. Two main ingredients are adopted in the proof: (i) an $H^2$-conditional $L^2$-stability estimate, which is established using the energy method to avoid singularity of the logarithmic nonlinearity, and (ii) mathematical induction with inverse inequalities to control the $H^2$-norm of the numerical solution. Numerical results are reported to confirm our error estimates and demonstrate the necessity of the time step size restriction imposed. We also apply the EWI-FS method to investigate soliton collisions in one dimension and vortex dipole dynamics in two dimensions. |
| title | Optimal error bounds on an exponential wave integrator Fourier spectral method for the logarithmic Schrödinger equation |
| topic | Numerical Analysis 35Q55, 65M15, 65M70, 81Q05 |
| url | https://arxiv.org/abs/2412.16902 |