Efficient numerical methods for solving nonlinear wave equations and studying the propagation
and stability properties of their solitary waves (solitons) are applied to a Boussinesq type equation in
one space dimension. These methods use a pseudospectral (Fourier transform) treatment of the space
dependence, together with (a) finite differences, or (b) a fourth-order Runge-Kutta scheme (RK4), for
the time evolution. Our schemes follow very accurately single solitons, which are given by simple
closed formulas and are known to be stable for all allowed velocities. However, as a parameter of
the problem tends to the critical value b = 0.5, where the velocity of the exact soliton vanishes, our
solutions destabilize due to numerical errors, producing two small solitons in the place of the exact
one. On the other hand, when we study the interaction of two such solitons, starting far apart from
each other, we find in the b1, b2 parameter plane a curve beyond which the solution becomes unstable
by exponential blow-up of the amplitudes, independently of our space and time discretization. We
claim that this is due to a dynamical resonance rather than the accumulation of numerical errors,
in agreement with theoretical predictions. Our implementation relies on the fast Fourier transform
(FFT) algorithm and no major differences are observed, when either scheme (a) or (b) is used for the
evolution of time. |
