Stability of the Richardson Extrapolation combined with some implicit Runge-Kutta methods

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

卷:310

Article

Zlatev, Zahari¹  Dimov, Ivan²  Farago, Istvan^3,4  Georgiev, Krassimir²  Havasi, Agnes^3,4 

[1] Aarhus Univ, Dept Environm Sci, Roskilde, Denmark

[2] Bulgarian Acad Sci, Inst Informat & Commun Technol, Sofia, Bulgaria

[3] MTA ELTE Numer Anal & Large Networks Res Grp, Pazmany Ps 1-C, H-1117 Budapest, Hungary

[4] Eotvos Lorand Univ, Dept Appl Anal & Computat Math, Pazmany Ps 1-C, H-1117 Budapest, Hungary

关键词: Systems of ODEs;    Implicit Runge-Kutta methods;    Richardson Extrapolation;    Stability properties;    Atmospheric chemical scheme;   

DOI  :  10.1016/j.cam.2016.03.018

来源: Elsevier

【 摘 要 】

The implementation of the Richardson Extrapolation in combination with different numerical methods for solving systems of ordinary differential equations (ODEs) is relatively simple, but the important requirement for stability of the computational process may cause serious difficulties. For example, the commonly used by scientists and engineers Trapezoidal Rule has good stability properties, but its combination with the Richardson Extrapolation is unstable. Therefore, it is necessary to study in advance and very carefully the stability of the new numerical methods arising when the scientists and the engineers use this computational device in combination with different algorithms for solving systems of ODEs. We are presenting a systematic investigation of the implementation of Richardson Extrapolation for two implicit Runge-Kutta methods. Three numerical examples, including an atmospheric chemical scheme used successfully in several extensive environmental studies and described mathematically by a very stiff and badly scaled nonlinear system of ODEs, are presented to illustrate the advantages of the presented approach. The numerical results show that not only are the computations stable, but also the achieved accuracy is higher when the Richardson Extrapolation is additionally applied. It will be possible to derive similar stability and accuracy results for other implicit Runge-Kutta methods. (C) 2016 Elsevier B.V. All rights reserved.

【 授权许可】

Free