Comparison of leapfrog, Smolarkiewicz, and MacCormack schemes applied to nonlinear equations

Abstract

A detailed comparison of these schemes is presented using the nonlinear advection equation and the Euler equations for a vareity of conditions at different Courant numbers. The unifiltered leapfrog is the least acceptable for the solution of nonlinear equations. Although it is numerically stable for linear problems, when used to solve nonlinear equations (without using any explicit filtering) it becomes numerically unstable or nonlinearly unstable. Furthermore, it introduces large phase errors, and produces better results with small Courant numbers. The MacCormack scheme is nonlinearly stable, produces modest amounts of numerical dispersion and diffusion, has no phase speed error, and works best at Courant numbers close to 1. When applied to nonlinear equations, the Smolarkiewicz scheme exhibits the least amount of numerical diffusion but more numerical diffusion but more numerical dispersion than the MacCormack scheme. For stability it requires Courant numbers equal to or smaller than 0.5. In practical applications, the MacCormack scheme is recommended for the solution of the nonlinear equations, and either the Smolarkiewicz or the MacCormack scheme for equations involving conservation of a passive scalar. (from Authors)