Boussinesq-type model with boundary-fitted coordinate system

An efficient numerical method for improved Boussinesq equations is proposed. Because of the use of the boundary-fitted coordinate system, the method facilitates the solution of wave problems with complicated boundaries and topography. The iterative method, combined with an efficient predictor-corrector scheme, is adopted for the numerical solution of the governing differential equations. The proposed numerical scheme is verified by three test cases where laboratory data were available for comparison. The successful simulation of wave runup around a circular cylinder and twin-tandem cylinders and wave propagation over an elliptical shoal shows that the proposed numerical method is both stable and accurate.