Abstract
New SDIRKN methods specially adapted to the numerical integration of second-order stiff ODE systems with periodic solutions are obtained. Our interest is focused on the dispersion (phase errors) of the dominant components in the numerical oscillations when these methods are applied to the homogeneous linear test model. Based on this homogeneous test model we derive the dispersion and P-stability conditions for SDIRKN methods which are assumed to be zero dissipative. Two four-stage symplectic and P-stable methods with algebraic order 4 and high order of dispersion are obtained. One of the methods is symmetric and sixth-order dispersive whereas the other method is nonsymmetric and eighth-order dispersive. These methods have been applied to a number of test problems (linear as well as nonlinear) and some numerical results are presented to show their efficiency when they are compared with other methods derived by Sharp et al. [IMA J. Numer. Anal. 10 (1990) 489–504].
Similar content being viewed by others
References
L. Brusa and L. Nigro, A one-step method for direct integration of structural dynamic equations, Internat. J. Numer. Methods Engrg. 15 (1980) 685–699.
J.R. Cash, High order, P-stable formulae for the numerical integration of periodic initial value problems, Numer. Math. 37 (1981) 355–370.
M.M. Chawla, P.S. Rao and B. Neta, Two-step fourth order P-stable methods with phase-lag of order six for y″ = f (t,y), J. Comput. Appl. Math. 16 (1986) 233–236.
J.M. Franco, I. Gómez and L. Rández, SDIRK methods for stiff ODEs with oscillating solutions, J. Comput. Appl. Math. 81 (1997) 197–209.
J.D. Lambert and I.A. Watson, Symmetric multistep methods for periodic initial value problems, J. Inst. Math. Appl. 18 (1976) 189–202.
J.M. Sanz-Serna and M.P. Calvo, Numerical Hamiltonian Problems (Chapman and Hall, London, 1994).
P.W. Sharp, J.M. Fine and K. Burrage, Two-stage and three-stage diagonally implicit Runge–Kutta Nyström methods of orders three and four, IMA J. Numer. Anal. 10 (1990) 489–504.
R.M. Thomas, Phase properties of high-order, almost P-stable formulae, BIT 24 (1984) 225–238.
P.J. Van der Houwen and B.P. Sommeijer, Explicit Runge–Kutta (–Nyström) methods with reduced phase errors for computing oscillating solutions, SIAM J. Numer. Anal. 24 (1987) 595–617.
P.J. Van der Houwen and B.P. Sommeijer, Phase-lag analysis of implicit Runge–Kutta methods, SIAM J. Numer. Anal. 26 (1989) 214–229.
P.J. Van der Houwen and B.P. Sommeijer, Diagonally implicit Runge–Kutta–Nyström methods for oscillatory problems, SIAM J. Numer. Anal. 26 (1989) 414–429.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Franco, J., Gómez, I. & Rández, L. Four-stage symplectic and P-stable SDIRKN methods with dispersion of high order. Numerical Algorithms 26, 347–363 (2001). https://doi.org/10.1023/A:1016629706668
Issue Date:
DOI: https://doi.org/10.1023/A:1016629706668