Higher-order split-step Fourier schemes for the generalized nonlinear Schrödinger equation

The generalized nonlinear Schrodinger (GNLS) equation is solved numerically by a split-step Fourier method. The first, second and fourth-order versions of the method are presented. A classical problem concerning the motion of a single solitary wave is used to compare the first, second and fourth-order schemes in terms of the accuracy and the computational cost. This numerical experiment shows that the split-step Fourier method provides highly accurate solutions for the GNLS equation and that the fourth-order scheme is computationally more efficient than the first-order and second-order schemes. Furthermore, two test problems concerning the interaction of two solitary waves and an exact solution that blows up in finite time, respectively, are investigated by using the fourth-order split-step scheme and particular attention is paid to the conserved quantities as an indicator of the accuracy. The question how the present numerical results are related to those obtained in the literature is discussed.

[1]  J. A. C. Weideman,et al.  A numerical study of the nonlinear Schro¨dinger equation involving quintic terms , 1990 .

[2]  M. Suzuki,et al.  General theory of higher-order decomposition of exponential operators and symplectic integrators , 1992 .

[3]  Gulcin M. Muslu,et al.  A split-step Fourier method for the complex modified Korteweg-de Vries equation☆ , 2003 .

[4]  Spectral methods and mappings for evolution equations on the infinite line , 1990 .

[5]  J. M. Sanz-Serna,et al.  Numerical Hamiltonian Problems , 1994 .

[6]  T. Taha,et al.  Analytical and numerical aspects of certain nonlinear evolution equations. II. Numerical, nonlinear Schrödinger equation , 1984 .

[7]  John P. Boyd,et al.  Spectral Modeling of Nonlinear Dispersive Waves , 1998 .

[8]  Qianshun Chang,et al.  Difference Schemes for Solving the Generalized Nonlinear Schrödinger Equation , 1999 .

[9]  H. Yoshida Construction of higher order symplectic integrators , 1990 .

[10]  T. Driscoll,et al.  Regular Article: A Fast Spectral Algorithm for Nonlinear Wave Equations with Linear Dispersion , 1999 .

[11]  G. Quispel,et al.  Acta Numerica 2002: Splitting methods , 2002 .

[12]  D. Pathria,et al.  Pseudo-spectral solution of nonlinear Schro¨dinger equations , 1990 .

[13]  B. Herbst,et al.  Split-step methods for the solution of the nonlinear Schro¨dinger equation , 1986 .

[14]  D. Pathria,et al.  Exact solutions for a generalized nonlinear Schrödinger equation , 1989 .

[15]  R. McLachlan Symplectic integration of Hamiltonian wave equations , 1993 .

[16]  Mark P. Robinson,et al.  The solution of nonlinear Schrödinger equations using orthogonal spline collocation , 1997 .

[17]  J. Boyd The Rate of Convergence of Fourier Coefficients for Entire Functions of Infinite Order with Application to the Weideman-Cloot Sinh-Mapping for Pseudospectral Computations on an Infinite Interval , 1994 .

[18]  Qin Sheng,et al.  Solving the Generalized Nonlinear Schrödinger Equation via Quartic Spline Approximation , 2001 .