A First-Order Fourier Integrator for the Nonlinear Schrödinger Equation on T Without Loss of Regularity.

Yifei Wu,Fangyan Yao
DOI: https://doi.org/10.1090/mcom/3705
IF: 2.118
2021-01-01
Mathematics of Computation
Abstract:In this paper, we propose a first-order Fourier integrator for solving the cubic nonlinear Schrödinger equation in one dimension. The scheme is explicit and can be implemented using the fast Fourier transform. By a rigorous analysis, we prove that the new scheme provides the first-order accuracy in H γ H^\gamma for any initial data belonging to H γ H^\gamma , for any γ > 3 2 \gamma >\frac 32 . That is, up to some fixed time T T , there exists some constant C = C ( ‖ u ‖ L ∞ ( [ 0 , T ] ; H γ ) ) > 0 C=C(\|u\|_{L^\infty ([0,T]; H^{\gamma })})>0 , such that ‖ u n − u ( t n ) ‖ H γ ( T ) ≤ C τ , \begin{equation*} \|u^n-u(t_n)\|_{H^\gamma (\mathbb {T})}\le C \tau , \end{equation*} where u n u^n denotes the numerical solution at t n = n τ t_n=n\tau . Moreover, the mass of the numerical solution M ( u n ) M(u^n) verifies | M ( u n ) − M ( u 0 ) | ≤ C τ 5 . \begin{equation*} \left |M(u^n)-M(u_0)\right |\le C\tau ^5. \end{equation*} In particular, our scheme does not cost any additional derivative for the first-order convergence and the numerical solution obeys the almost mass conservation law. Furthermore, if u 0 ∈ H 1 ( T ) u_0\in H^1(\mathbb {T}) , we rigorously prove that ‖ u n − u ( t n ) ‖ H 1 ( T ) ≤ C τ 1 2 − , \begin{equation*} \|u^n-u(t_n)\|_{H^1(\mathbb {T})}\le C\tau ^{\frac 12-}, \end{equation*} where C = C ( ‖ u 0 ‖ H 1 ( T ) ) > 0 C= C(\|u_0\|_{H^1(\mathbb {T})})>0 .
What problem does this paper attempt to address?