Numerical simulation of solitary waves of Rosenau–KdV equation by Crank–Nicolson meshless spectral interpolation method

Abstract

An efficient and accurate Crank–Nicolson meshless spectral radial point interpolation (CN-MSRPI) method is proposed for the numerical solution of nonlinear Rosenau–KdV equation. The proposed method uses meshless shape functions, owing Kronecker delta property, for approximation of spatial operator. Crank–Nicolson difference scheme is used for temporal operator approximation. Single solitary wave motion, interaction of double and triple solitary waves as well as generation of train of solitary waves from initial data are numerically simulated. Error analysis is made via computation of discrete $$L_{\infty }$$, $$L_{2}$$ and $$L_{\mathrm{rms}}$$ error norms. Efficiency of the proposed numerical scheme is assessed via variation of number of nodes N and time step-size $$\tau $$. Two invariant quantities correspond to mass and energy are computed using the proposed method for further validation. Stability of the proposed method is discussed and verified computationally. Comparison of obtained results made with exact and existing results in the literature revealed the proposed CN-MSRPI method superiority.