Stability of the Richardson Extrapolation applied together with the θ-method
作者:
Highlights:
•
摘要
Consider a system of ordinary differential equations (ODEs) dy/dt=f(t,y) where (a) t∈[a,b] with b>a, (b) y is a vector containing s components and (c) y(a) is given. The θ-method is applied to solve approximately the system of ODEs on a set of prescribed grid points. If N is the number of time steps that are to be carried out, then this numerical method can be defined using the following set of relationships: yn=yn−1+h(1−θ)f(tn−1,yn−1)+hθf(tn,yn),θ∈[0.5,1.0],n=1,2,…,N,h=(b−a)/N,tn=tn−1+h=t0+nh,t0=a,tN=b. As a rule, the accuracy of approximations {yn∣n=1,2,…,N} can be improved by applying the Richardson Extrapolation under the assumption that the stability of the computational process is preserved. Therefore, it is natural to require that the combined numerical method (Richardson Extrapolation + the θ-method) is in some sense stable. It is proved in this paper that the combined method is strongly A-stable when θ∈[2/3,1.0]. It is furthermore shown that some theorems proved in a previous paper by the same authors, Faragó et al. (2009) [1], are simple corollaries of the main result obtained in the present work.The usefulness of the main result in the solution of many problems arising in different scientific and engineering areas is demonstrated by performing a series of tests with an extremely badly scaled and very stiff atmospheric chemistry scheme which is actually used in several well-known large-scale air pollution models.
论文关键词:Systems of ordinary differential equations,Numerical methods,Richardson Extrapolation,Stability,Atmospheric chemistry scheme
论文评审过程:Received 31 October 2009, Revised 26 May 2010, Available online 8 June 2010.
论文官网地址:https://doi.org/10.1016/j.cam.2010.05.052