Jacobi-Gauss-Lobatto collocation method for the numerical solution of 1+1 nonlinear Schrödinger equations

Authors:
E. H. Doha;A. H. Bhrawy;M. A. Abdelkawy;Robert A. Van Gorder
Affiliations:
Department of Mathematics, Faculty of Science, Cairo University, Giza, Egypt;Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah, Saudi Arabia and Department of Mathematics, Faculty of Science, Beni-Suef University, Beni-Suef, Egypt;Department of Mathematics, Faculty of Science, Beni-Suef University, Beni-Suef, Egypt;Department of Mathematics, University of Central Florida, Orlando, FL 32816, USA
Venue:
Journal of Computational Physics
Year:
2014

Citing 6
Cited 0

Efficient spectral-Galerkin algorithms for direct solution of fourth-order differential equations using Jacobi polynomials

Applied Numerical Mathematics
Numerical solution of coupled nonlinear Schrödinger equation by Galerkin method

Mathematics and Computers in Simulation
The use of compact boundary value method for the solution of two-dimensional Schrödinger equation

Journal of Computational and Applied Mathematics
The spectral collocation method with three different bases for solving a nonlinear partial differential equation arising in modeling of nonlinear waves

Mathematical and Computer Modelling: An International Journal
An efficient direct solver for multidimensional elliptic Robin boundary value problems using a Legendre spectral-Galerkin method

Computers & Mathematics with Applications
New algorithms for solving high even-order differential equations using third and fourth Chebyshev-Galerkin methods

Journal of Computational Physics

Quantified Score

Hi-index	31.45

Visualization

Abstract

A Jacobi-Gauss-Lobatto collocation (J-GL-C) method, used in combination with the implicit Runge-Kutta method of fourth order, is proposed as a numerical algorithm for the approximation of solutions to nonlinear Schrodinger equations (NLSE) with initial-boundary data in 1+1 dimensions. Our procedure is implemented in two successive steps. In the first one, the J-GL-C is employed for approximating the functional dependence on the spatial variable, using (N-1) nodes of the Jacobi-Gauss-Lobatto interpolation which depends upon two general Jacobi parameters. The resulting equations together with the two-point boundary conditions induce a system of 2(N-1) first-order ordinary differential equations (ODEs) in time. In the second step, the implicit Runge-Kutta method of fourth order is applied to solve this temporal system. The proposed J-GL-C method, used in combination with the implicit Runge-Kutta method of fourth order, is employed to obtain highly accurate numerical approximations to four types of NLSE, including the attractive and repulsive NLSE and a Gross-Pitaevskii equation with space-periodic potential. The numerical results obtained by this algorithm have been compared with various exact solutions in order to demonstrate the accuracy and efficiency of the proposed method. Indeed, for relatively few nodes used, the absolute error in our numerical solutions is sufficiently small.