Efficient implementation of stable Richardson Extrapolation algorithms

Istvn Farag, Gnes Havasi, Zahari Zlatev

Research output: Contribution to journalArticle

18 Citations (Scopus)

Abstract

Richardson Extrapolation is a powerful computational tool which can successfully be used in the efforts to improve the accuracy of the approximate solutions of systems of ordinary differential equations (ODEs) obtained by different numerical methods (including here combined numerical methods consisting of appropriately chosen splitting procedures and classical numerical methods). Some stability results related to two implementations of the Richardson Extrapolation (Active Richardson Extrapolation and Passive Richardson Extrapolation) are formulated and proved in this paper. An advanced atmospheric chemistry scheme, which is commonly used in many well-known operational environmental models, is applied in a long sequence of experiments in order to demonstrate the fact that it is indeed possible to improve the accuracy of the numerical results when the Richardson Extrapolation is used (also when very difficult, badly scaled and stiff non-linear systems of ODEs are to be treated),the computations can become unstable when the combination of the Trapezoidal Rule and the Active Richardson Extrapolation is used,the application of the Active Richardson Extrapolation with the Backward Euler Formula is leading to a stable computational process,experiments with different algorithms for solving linear systems of algebraic equations are very useful in the efforts to select the most suitable approach for the particular problems solved andthe computational cost of the Richardson Extrapolation is much less than that of the underlying numerical method when a prescribed accuracy has to be achieved.

Original languageEnglish
Pages (from-to)2309-2325
Number of pages17
JournalComputers and Mathematics with Applications
Volume60
Issue number8
DOIs
Publication statusPublished - Oct 2010

Keywords

  • Atmospheric chemistry scheme
  • Backward Euler Formula
  • Richardson Extrapolation
  • Sparse matrix technique
  • Stability
  • Trapezoidal rule

ASJC Scopus subject areas

  • Modelling and Simulation
  • Computational Theory and Mathematics
  • Computational Mathematics

Fingerprint Dive into the research topics of 'Efficient implementation of stable Richardson Extrapolation algorithms'. Together they form a unique fingerprint.

  • Cite this