On the odd-even hopscotch scheme for the numerical integration of time-dependent partial differential equations

This paper is devoted to the odd-even hopscotch scheme for the numerical integration of time-dependent partial differential equations. Attention is focussed on two aspects. Firstly, via the equivalence to the combined leapfrog-Du Fort–Frankel method we derive the explicit expression of the critical time step for von Neumann stability for a class of multi-dimensional convection-diffusion equations. This expression can be derived directly by applying a useful stability theorem due to Hindmarsh, Gresho and Griffiths [9]. The interesting thing on the critical time step is that it is independent of the diffusion parameter and yet smaller than the critical time step for zero diffusion, but only in the multi-dimensional case. This curious phenomenon does not occur for the one-dimensional problem. Secondly, we consider the drawback of the Du Fort–Frankel accuracy deficiency of the hopscotch scheme. To overcome this deficiency we discuss global Richardson extrapolation in time. This simple device can always be used without reducing feasibility. Numerical examples are given to illustrate the outcome of the extrapolation.

[1]  Linear stability of finite difference equations for three-dimensional flow problems , 1975 .

[2]  Ernst Hairer,et al.  Asymptotic expansions of the global error of fixed-stepsize methods , 1984 .

[3]  A. R. Gourlay,et al.  Hopscotch: a Fast Second-order Partial Differential Equation Solver , 1970 .

[4]  A. R. Gourlay.  Some recent methods for the numerical solution of time-dependent partial differential equations , 1971, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[5]  P. Henrici.  Discrete Variable Methods in Ordinary Differential Equations , 1962 .

[6]  H. Stetter.  Analysis of Discretization Methods for Ordinary Differential Equations , 1973 .

[7]  R. D. Richtmyer,et al.  Difference methods for initial-value problems , 1959 .

[8]  John J. H. Miller.  On the Location of Zeros of Certain Classes of Polynomials with Applications to Numerical Analysis , 1971 .

[9]  Paul Gordon,et al.  Nonsymmetric Difference Equations , 1965 .

[10]  David F. Griffiths,et al.  The stability of explicit Euler time‐integration for certain finite difference approximations of the multi‐dimensional advection–diffusion equation , 1984 .

[11]  A. R. Gourlay,et al.  Hopscotch difference methods for nonlinear hyperbolic systems , 1972 .

[12]  J. G. Verwer,et al.  Global extrapolation of a first order splitting method , 1983 .

[13]  L. Chambers.  Linear and Nonlinear Waves , 2000, The Mathematical Gazette.

[14]  Clive A. J. Fletcher,et al.  A comparison of finite element and finite difference solutions of the one- and two-dimensional Burgers' equations , 1983 .

[15]  Lawrence F. Shampine,et al.  Global Error Estimates for Ordinary Differential Equations , 1976, TOMS.