Uniform point-wise error estimates of semi-implicit compact finite difference methods for the nonlinear Schrödinger equation perturbed by wave operator

Abstract

This study is devoted to analysis of semi-implicit compact finite difference (SICFD) methods for the nonlinear Schrödinger equation (NLS) perturbed by the wave operator (NLSW) with a perturbation strength described by a dimensionless parameter ε∈(0,1]. Uniform l∞-norm error bounds of the proposed SICFD schemes are built to give immediate insight on point-wise error occurring as time increases, and the explicit dependence of the mesh size and time step on the parameter ε is also figured out. In the small ε regime, highly oscillations arise in time with O(ε2)-wavelength. This highly oscillatory nature in time as well as the difficulty raised by the compact FD discretization make establishing the l∞-norm error bounds uniformly in ε of the SICFD methods for NLSW to be a very interesting and challenging issue. The uniform l∞-norm error bounds in ε are proved to be of O(h4+τ) and O(h4+τ2/3) with time step τ and mesh size h for well-prepared and ill-prepared initial data. Finally, numerical results are reported to verify the error estimates and show the sharpness of the convergence rates in the respectively parameter regimes.