The analytical evolution of NLS solitons due to the numerical discretization error

• Published in: Journal of physics. A, Mathematical and theoretical Vol. 44; no. 50; pp. 505205 - 17
• Main Authors: Hoseini, S M, Marchant, T R
• Format: Journal Article
• Language: English
• Published: Bristol IOP Publishing 16.12.2011; IOP
• Subjects: Discretization; Error analysis; Exact sciences and technology; Mathematical analysis; Mathematical models; Nonlinearity; Physics; Schroedinger equation; Shelves; Solitary waves; Solitons; Solitary wave; Perturbation theory; Second order; Singularity; Inverse scattering; Discretization; First order; Digital simulation; Solitons; Schroedinger equation; Analytical solution; Corrections
• ISSN: 1751-8121; 1751-8113; 1751-8121
• DOI: 10.1088/1751-8113/44/50/505205

Soliton perturbation theory is used to obtain analytical solutions describing solitary wave tails or shelves, due to numerical discretization error, for soliton solutions of the nonlinear Schrodinger equation. Two important implicit numerical schemes for the nonlinear Schrodinger equation, with second-order temporal and spatial discretization errors, are considered. These are the Crank-Nicolson scheme and a scheme, due to Taha [1], based on the inverse scattering transform. The first-order correction for the solitary wave tail, or shelf, is in integral form and an explicit expression is found for large time. The shelf decays slowly, at a rate of t super(-1/2), which is characteristic of the nonlinear Schrodinger equation. Singularity theory, usually used for combustion problems, is applied to the explicit large-time expression for the solitary wave tail. Analytical results are then obtained, such as the parameter regions in which qualitatively different types of solitary wave tails occur, the location of zeros and the location and amplitude of peaks. It is found that three different types of tail occur for the Crank-Nicolson and Taha schemes and that the Taha scheme exhibits some unusual symmetry properties, as the tails for left and right moving solitary waves are different. Optimal choices of the discretization parameters for the numerical schemes are also found, which minimize the amplitude of the solitary wave tail. The analytical solutions are compared with numerical simulations, and an excellent comparison is found.