Numerical Simulations and Parallel Implementation of Some Nonlinear Schrödinger Systems

As it is well known, the nonlinear Schrodinger equation (NLS) $$ \left\{ i \right.{W_t} + {W_{xx}} + a{\left| W \right|^2}W = 0,\quad W\left( {x,0} \right) = {W_0}\left( x \right),$$ (1) where a =const, and W(x, t) being a complex function, can be exactly solved by using the Inverse Scattering Theory. At the same time several numerical schemes have been proposed to solve this equation1–7. In general, the perturbations of the NLS equation as driving forces, dissipations, stochastic potentials, coupled NLS systems are not tractable analytically. The same happens for the NLS systems in two and three space dimensions. As a consequence, large and massive computations are needed and a possible way to make them is by the parallel computing. In this context, we started a project to implement suitable algorithms for the complex nonlinear Schrodinger systems on a distributed array of transputers. Our first result was to propose an alternative finite difference scheme for the unperturbed NLS equation in one space dimension8. This algorithm has been tested in several relevant physical simulations, that is described in Sections 2 and 3. In Sections 4 to 7 we describe some features of the parallelization of the scheme as well as its parallel implementation on an array of transputers.

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

[2]  G. Stegeman,et al.  Soliton switching in fiber nonlinear directional couplers. , 1988, Optics letters.

[3]  Lionel M. Ni,et al.  Scalable Problems and Memory-Bounded Speedup , 1993, J. Parallel Distributed Comput..

[4]  Luigi Brugnano,et al.  A parallel solver for tridiagonal linear systems for distributed memory parallel computers , 1991, Parallel Comput..

[5]  B. Suydam Self-focusing of very powerful laser beams II , 1974 .

[6]  J. M. Sanz-Serna,et al.  Methods for the numerical solution of the nonlinear Schroedinger equation , 1984 .

[7]  Ben M. Herbst,et al.  Numerical Experience with the Nonlinear Schrödinger Equation , 1985 .

[8]  G. C. Fox,et al.  Solving Problems on Concurrent Processors , 1988 .

[9]  J. G. Verwer,et al.  Conerservative and Nonconservative Schemes for the Solution of the Nonlinear Schrödinger Equation , 1986 .

[10]  Yuri S. Kivshar,et al.  Dynamics of Solitons in Nearly Integrable Systems , 1989 .

[11]  W. Strauss,et al.  Numerical solution of a nonlinear Klein-Gordon equation , 1978 .

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

[13]  John Ll Morris,et al.  An investigation into the effect of product approximation in the numerical solution of the cubic nonlinear Schro¨dinger equation , 1988 .

[14]  Govind P. Agrawal,et al.  Nonlinear Fiber Optics , 1989 .

[15]  Curtis R. Menyuk,et al.  Pulse propagation in an elliptically birefringent Kerr medium , 1989 .

[16]  L. Vázquez,et al.  Numerical solution of the sine-Gordon equation , 1986 .

[17]  Luigi A. Lugiato,et al.  Cooperative frequency locking and stationary spatial structures in lasers , 1988 .

[18]  David Potter Computational physics , 1973 .

[19]  J. Ortega Introduction to Parallel and Vector Solution of Linear Systems , 1988, Frontiers of Computer Science.

[20]  John W. Miles,et al.  An Envelope Soliton Problem , 1981 .

[21]  Michel C. Delfour,et al.  Finite-difference solutions of a non-linear Schrödinger equation , 1981 .