Parallel implementation of the split-step and the pseudospectral methods for solving higher KdV equation

Numerical simulations show that higher order KdV equation under certain conditions has a self-focusing singularity, which means that the solution of the equation blows up in finite time. In this paper, two numerical schemes: the split-step Fourier transform and the pseudospectral methods are used to investigate this self-focusing singularity problem. Parallel algorithms for the proposed schemes are designed and implemented. FFTW-MPI algorithm designed by Matteo Frigo and Steven Johnson is used for parallel implementation of the discrete Fourier transform (DFT). The parallel algorithms are implemented on an SGI Origin 2000 multiprocessor computer and experiments show that a considerable speedup is attained.

[1]  P. Souganidis,et al.  Stability and instability of solitary waves of Korteweg-de Vries type , 1987, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

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

[3]  Bengt Fornberg,et al.  A numerical and theoretical study of certain nonlinear wave phenomena , 1978, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[4]  J. Bona,et al.  Fully discrete galerkin methods for the korteweg-de vries equation☆ , 1986 .

[5]  D. Korteweg,et al.  XLI. On the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary waves , 1895 .

[6]  Steven G. Johnson,et al.  FFTW: an adaptive software architecture for the FFT , 1998, Proceedings of the 1998 IEEE International Conference on Acoustics, Speech and Signal Processing, ICASSP '98 (Cat. No.98CH36181).

[7]  Thiab R. Taha,et al.  Analytical and numerical aspects of certain nonlinear evolution equations. I. Analytical , 1984 .