Pryce pre-analysis adapted to some DAE solvers |
| |
Authors: | T Nguyen-Ba H Yagoub H Hao R Vaillancourt |
| |
Institution: | Department of Mathematics and Statistics, University of Ottawa, Ottawa, Ontario, Canada K1N 6N5 |
| |
Abstract: | The ODE solver HBT(12)5 of order 12 (T. Nguyen-Ba, H. Hao, H. Yagoub, R. Vaillancourt, One-step 5-stage Hermite-Birkho-Taylor ODE solver of order 12, Appl. Math. Comput. 211 (2009) 313-328. doi:10.1016/j.amc.2009.01.043), which combines a Taylor series method of order 9 with a Runge-Kutta method of order 4, is expanded into the DAE solver HBT(12)5DAE of order 12. Dormand-Prince’s DP(8, 7)13M is also expanded into the DAE solver DP(8, 7)DAE. Pryce structural pre-analysis, extended ODEs and ODE first-order forms are adapted to these DAE solvers with a stepsize control based on local error estimators and a modified Pryce algorithm to advance integration. HBT(12)5DAE uses only the first nine derivatives of the unknown variables as opposed to the first 12 derivatives used by the Taylor series method T12DAE of order 12. Numerical results show the advantage of HBT(12)5DAE over T12DAE, DP(8, 7)DAE and other known DAE solvers. |
| |
Keywords: | General linear method Non-stiff DAEs Hermite-Birkhoff method Taylor method Pryce structural analysis |
本文献已被 ScienceDirect 等数据库收录! |
|