Share This Article:

Filtered Leapfrog Time Integration with Enhanced Stability Properties

Full-Text HTML XML Download Download as PDF (Size:1057KB) PP. 1354-1370
DOI: 10.4236/jamp.2016.47145    501 Downloads   672 Views  


The Leapfrog method for the solution of Ordinary Differential Equation initial value problems has been historically popular for several reasons. The method has second order accuracy, requires only one function evaluation per time step, and is non-dissipative. Despite the mentioned attractive properties, the method has some unfavorable stability properties. The absolute stability region of the method is only an interval located on the imaginary axis, rather than a region in the complex plane. The method is only weakly stable and thus exhibits computational instability in long time integrations over intervals of finite length. In this work, the use of filters is examined for the purposes of both controlling the weak instability and also enlarging the size of the absolute stability region of the method.

Cite this paper

Aluthge, A. , A. Sarra, S. and Estep, R. (2016) Filtered Leapfrog Time Integration with Enhanced Stability Properties. Journal of Applied Mathematics and Physics, 4, 1354-1370. doi: 10.4236/jamp.2016.47145.


[1] Richardson, L.F. (1922) Weather Prediction by Numerical Process. Cambridge University Press, Cambridge.
[2] Williams, P. (2009) A Proposed Modification to the Robert-Asselin Time Filter. Monthly Weather Review, 137, 2538-2546.
[3] Coiffier, J. (2011) Fundamentals of Numerical Weather Prediction. Cambridge University Press, Cam-bridge.
[4] Butcher, J. (2003) Numerical Methods for Ordinary Differential Equations. Wiley.
[5] Milne, W. and Reynolds, R. (1959) Stability of a Numerical Solutions of Differential Equations. Journal of the ACM, 6, 196-203.
[6] Asselin, R. (1972) Frequency Filter for Time Integrations. Monthly Weather Review, 100, 487-490.<0487:FFFTI>2.3.CO;2
[7] Gragg, W.B. (1963) Re-peated Extrapolation to the Limit in the Numerical Solution of Ordinary Differential Equations. PhD Thesis, UCLA.
[8] Bulirsch, R. and Stoer, J. (1966) Numerical Treatment of Ordinary Differential Equations by Extrapolation Methods. Numerische Mathematik, 8, 1-13.
[9] Shampine, L. and Baca, L. (1983) Smoothing the Extrapolated Midpoint Rule. Numerische Mathematik, 41, 165-175.
[10] Iri, M. (1963) A Stabilizing Device for Unstable Numerical Solutions of Ordinary Differential Equations—Design Principle and Application of a Filter. Journal of the Information Processing Society of Japan, 4, 249-260.
[11] Li, Y. and Trenchea, C. (2014) A Higher-Order Robert-Asselin Type Time Filter. Journal of Computational Physics, 259, 23-32.
[12] Norton, T. and Hill, A. (2015) An Iterative Starting Method to Control Parasitism for the Leapfrog Method. Applied Numerical Mathematics, 87, 145-156.
[13] Dahlquist, G. (1985) 33 Years of Numerical Instability, Part I. BIT, 188-204.
[14] Todd, J. (1950) Solution of Differential Equations by Recurrence Relations. Mathematical Tables and Other Aids to Computation, 4, 39-44.
[15] Hairer, E., Norsett, S. and Wanner, G. (2000) Solving Ordinary Differential Equations I: Nonstiff Problems. Springer.
[16] Hairer, E. and Wanner, G. (2000) Solving Ordinary Differential Equations II: Stiff and Differential—Algebraic Problems. 2nd Edition, Springer.
[17] Henrici, P. (1962) Discrete Variable Methods in Ordinary Differential Equations. Wiley.
[18] Iserles, A. (1973) A First Course in the Analysis of Differential Equations. Wiley.
[19] Lambert, J. (1973) Computational Methods in Ordinary Differential Equations. Wiley.
[20] Aluthge, A. (1985) Filtering and Extrapolation Techniques in the Numerical Solution of Ordinary Differential Equations. Master’s Thesis, University of Ottawa, Ottawa.
[21] Boyd, J.P. (2000) Chebyshev and Fourier Spectral Methods. 2nd Edition, Dover Publications, Inc., New York.
[22] Canuto, C., Hussaini, M., Quarteroni, A. and Zang, T. (1988) Spectral Methods in Fluid Dynamics. Springer.
[23] Weideman, J.A.C. and Reddy, S. (2000) A MATLAB Differentiation Matrix Suite. ACM Transactions on Mathematical Software, 26, 465-519.

comments powered by Disqus

Copyright © 2017 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.