Numerical solitary wave interaction: The order of the inelastic effect

Solitary wave interaction is examined using an extended Benjamin-Bona-Mahony (eBBM) equation. This equation includes higher-order nonlinear and dispersive effects and is is asymptotically equivalent to the extended Korteweg-de Vries (eKdV) equation. The eBBM formulation is preferable to the eKdV equation for the numerical modelling of solitary wave collisions, due to the stability of its finite-difference scheme. In particular, it allows the interaction of steep waves to be modelled, which due to numerical instability, is not possible using the eKdV equation. Numerical simulations of a number of solitary wave collisions of varying nonlinearity are performed for two special cases corresponding to surface water waves. The mass and energy of the dispersive wavetrain generated by the inelastic collision is tabulated and used to show that the change in solitary wave amplitude after interaction is of O(alpha(4)), verifying previously obtained theoretical predictions.