Comparative Review on Computational Performance of Multistep Schemes in Solving One-Dimensional Linear Wave Equation

Among several numerical methods used to solve the hyperbolic model of the linear wave equation, single-step algorithms can be the more popular ones. However, these algorithms are time-consuming while incurring numerical inaccuracy. Thus, multistep methods can be a suitable option as it has a high order of accuracy. This study aims to investigate and compare the computational performance of these multistep schemes in solving hyperbolic model based on one-dimensional linear wave equation. The techniques studied in this paper comprise the two-step Lax-Wendroff method, MacCormack method, second-order upwind method, Rusanov-Burstein-Mirin method, Warming-Kutler-Lomax method, and fourth-order Runge-Kutta method. Finite difference method is applied in discretisation. Our simulation found that although higher-order multistep methods are more stable than single-step algorithm, they suffer numerical diffusion. The two-step Lax-Wendroff method outperforms other schemes, although it is relatively simple compared with the other three and four steps schemes. The second-order upwind method is attractive as well because it is executable even with a high Courant number.