Perturbation solution for the 2D shallow-water waves

The Boussinesq model of shallow water flow is considered, which contains nonlinearity and fourth-order dispersion. Boussinesq equation has been one of the main soliton models in 1D. To find its 2D solutions, a perturbation series with respect to the small parameter ε:=c² is developed in the present work, where c is the phase speed of the localized wave. Within the order O(ε²)=O(c⁴), a hierarchy is derived consisting of fourth-order ordinary differential equations (ODEs). The Bessel operators involved are reformulated to facilitate the creation of difference schemes for the ODEs from the hierarchy. The numerical scheme uses a special approximation for the behavioral condition in the singularity point (the origin). The results of this work show that at infinity the 2D wave shape decays algebraically, rather than exponentially as in the 1D cases. The new result can be instrumental for understanding the interaction of 2D Boussinesq solitons, and for creating more efficient numerical algorithms explicitly acknowledging the asymptotic behavior of the solution.