ABSTRACT
A form of Rosenbrock-type methods optimal in terms of the number of non-zero parameters and computational costs per step is considered. A technique of obtaining \((m, k)\)-methods from some well-known Rosenbrock-type methods is justified. Formulas for transforming the parameters of \((m,k)\)-schemes and for obtaining a stability function are given for two canonical representations of the schemes. An \(L\)-stable \((3, 2)\)-method of order 3 is proposed, which requires two evaluations of the function: one evaluation of the Jacobian matrix and one \(LU\)-decomposition per step. A variable step size integration algorithm based on the \((3,2)\)-method is formulated. It provides a numerical solution for both explicit and implicit systems of ODEs. Numerical results are presented to show the efficiency of the new algorithm.
Similar content being viewed by others
REFERENCES
Hairer, E. and Wanner, G., Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems, Berlin: Springer-Verlag, 1996.
Novikov, E.A. and Shornikov, Yu.V., Komp’yuternoe modelirovanie zhestkikh gibridnykh sistem (Computer Simulation of Stiff Hybrid Systems), Novosibirsk: NSTU, 2012.
Rosenbrock, H.H., Some General Implicit Processes for the Numerical Solution of Differential Equations,Computer, 1963, vol. 5, pp. 329–330.
Artem’ev, S.S. and Demidov, G.V., An Algorithm of Variable Order and Step for the Numerical Solution of Stiff Systems of Ordinary Differential Equations, Dokl. Akad. Nauk SSSR, 1978, vol. 238, no. 3, pp. 517–520.
Kaps, P. and Wanner, G., A Study of Rosenbrock-Type Methods of High Order, Numerische Math., 1981, vol. 38, pp. 279–298.
Steihaug, T. and Wolfbrandt, A., An Attempt to Avoid Exact Jacobian and Nonlinear Equations in the Numerical Solution of Stiff Differential Equations, Math. Computat., 1979, vol. 33, pp. 521–534.
Werwer, J.G., Scholz, S., Blom, J.G., and Louter-Nool, M., A Class of Runge–Kutta–Rosenbrock Methods for Solving Stiff Differential Equations, ZAAM, 1983, vol. 63, pp. 13–20.
Voevodin, V.V. and Kuznetsov, Yu.A., Matritsy i vychisleniya (Matrices and Computations), Novosibirsk: Nauka, 1984.
Norsett, S.P., Restricted Padè Approximations to the Exponential Function, SIAM J. Num. An., 1978, vol. 15, no. 5, pp. 1008–1029.
Levykin, A.I. and Novikov, E.A., A Study of \((m,k)\)-Methods for Solving Differential-Algebraic Systems of Index 1, Comm. Comp. Informat. Sci., 2015, vol. 549, pp. 94–107.
Roche, M., Rosenbrock Methods for Differential Algebraic Equations, Numerische Math., 1988, vol. 52, pp. 45–63.
Levykin, A.I., Novikov, A.E., and Novikov, E.A., Third Order \((m, k)\)-Method for Solving Stiff Systems of ODEs and DAEs, Proc. XIV Int. Scientific-Technical Conf. APEIE-2018, vol. 1, part 4, Novosibirsk: NSTU, 2018, pp. 158–163.
Funding
This work was supported by the Russian Foundation for Basic Research (project no. 17-07-01513 A). The work of the first author was supported by ICM&MG SB RAS (state assignment no. 0315-2019-0002).
Author information
Authors and Affiliations
Corresponding authors
Rights and permissions
About this article
Cite this article
Levykin, A.I., Novikov, A.E. & Novikov, E.A. Schemes of (m, k)-Type for Solving Differential-Algebraic and Stiff Systems. Numer. Analys. Appl. 13, 34–44 (2020). https://doi.org/10.1134/S1995423920010036
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1134/S1995423920010036