A Posteriori analysis of a multirate numerical method for ordinary differential equations

详细信息	查看全文 | 推荐本文 |

• 作者：D. Estep^a ; ¹ ; ^{estep@stat.colostate.edu} ; V. Ginting^b ; ² ; S. Tavener^c
• 关键词：Adjoint operator ; A posteriori estimates ; Discontinuous Galerkin method ; Iterative method ; Multirate method ; Multiscale integration ; Operator decomposition
• 刊名：Computer Methods in Applied Mechanics and Engineering
• 出版年：2012
• 期刊代码：27_00457825
• 类别：cp
• 出版时间：1 June, 2012
• 卷：223-224
• 期：Complete
• 页码：10-27
• 文件大小：2607 K

摘要

In this paper, we analyze a multirate time integration method for systems of ordinary differential equations that present significantly different scales within the components of the model. The main purpose of this paper is to present a hybrid a priori - a posteriori error analysis that accounts for the effects of projections between the coarse and fine scale discretizations, the use of only a finite number of iterations in the iterative solution of the discrete equations, the numerical error arising in the solution of each component, and the effects on stability arising from the multirate solution. The hybrid estimate has the form of a computable a posteriori leading order expression and a provably-higher order a priori expression. We support this estimate by an a priori convergence analysis. We present several examples illustrating the accuracy of multirate integration schemes and the accuracy of the a posteriori estimate.