Composite methods for hyperbolic equations

Abstract

A composite approximation procedure combining the properties of the Lax-Wendroff and leapfrog algorithms is proposed for solving hyperbolic equations. For a one-dimensional equation, a three-step approximation consisting of a two-step Richtmeyer method followed by a leapfrog step is considered. This is a two-level scheme, so all difficulties, including storage requirements, associated with the three-level leapfrog are eliminated. For two-dimensional problems a generalization of the preceding method is used consisting of a rotated Richtmeyer method followed by a modified leapfrog step. It is found that the composite schemes are effective in reducing oscillations and nonlinear instabilities that affect the leapfrog method. The dissipation in the composite schemes is much less than in the Richtmeyer algorithm, and hence can be used for long term integrations.