Efficient numerical algorithms for the solution of “good” Boussinesq equation in water wave propagation

This paper proposes three fast and high accuracy numerical methods for solving a nonlinear partial differential equation (PDE) describing water waves and called the Boussinesq (Bq) equation. We numerically solve the Bq equation with fourth-order time-stepping schemes in combination with discrete Fourier transform. We discretize the original PDE with discrete Fourier transform in space and obtain a system of ordinary differential equations (ODEs) which will be solved with fourth-order time-stepping methods. After transforming the equation to a system of ODEs, the linear operator is not diagonal, but we can implement the methods such as diagonal case which reduces the CPU time. Comparing numerical solutions with analytical solutions demonstrates that those methods are accurate and readily implemented. Also we investigate the conservation of mass for Bq equation.

► In this paper the numerical solution of Boussinesq equation is investigated.
► These methods are based on DFT and fourth-order time-stepping schemes.
► Linear operator in the resulting system of ODEs is not diagonal.
► But we can implement the methods such as diagonal case.