The performance of phase-lag enhanced explicit Runge-Kutta Nyström pairs on \(N\)-body problems
type: Item
instance of: scholarly article
title: The performance of phase-lag enhanced explicit Runge-Kutta Nyström pairs on \(N\)-body problems
zbMATH Open document ID: 1251.70001
author: J. C. Castillo-Rogez
author: P. W. Sharp
author: Kevin R. Grazier
published in: Journal of Computational and Applied Mathematics
publication date: 2012-03-19T00:00:00Z
review text: Explicit Runge-Kutta Nyström methods (RKN\(p\)(\(q\)), where \(p\) is the order of the formulae to propagate the solution, and \(q\) is the order of the formulae to estimate the local error of integration) with enhanced phase-lag order (i.e., with free parameters selected to ensure the highest possible order for the integration of \(\ddot{y}=-\omega^2 y\) with \(\omega\) constant) are compared to general purpose RNK\(p\)(\(q\)) methods for long-time integrations of \(N\)-body problems. The aim of the authors is to see gains in the efficiency for long-time integrations, attending to the fact that the solutions of systems describing the motion of planets or satellites are someway close to oscillatory solutions of an harmonic oscillator. In the paper, the emphasis is made on the RKN8(6) pairs suitable for long-time integrations, but RKN4(3) and RKN6(4) have also been included just to gain insight in what way the efficiency depends on the order. Also integrations have been performed and compared using double and quadruple precision. The \(N\)-body problems are considered, and some of the main results achieved are the following: \begin{itemize}\item[1)] For Kepler's two-body problem used as a basic reference problem, it is detected that the pairs with enhanced phase-lag order are usually more efficient than the general purpose pairs for low eccentricities, and less efficient for high eccentricities. \item[2)] For the motion of the Sun and the planets Jupiter, Saturn, Uranus and Neptune the authors show that the enhanced RKN8(6) pair is more efficient than the general RKN8(6) pair for local error tolerances greater than 10\(^{-11}\), and less efficient for more severe ones. \item[3)] The same question as before, but including the motion of the Hellin-Roman-Crocket comet which has close encounters with Jupiter. Here the enhanced RKN8(6) appears to be slightly more efficient for lax accuracy requirements, and of similar efficiency otherwise. \item[4)] Finally, for the orbit…
Mathematics Subject Classification ID: 70-08
Mathematics Subject Classification ID: 70F10
Mathematics Subject Classification ID: 70F15
zbMATH DE Number: 6015751
zbMATH Keywords: long-time integration
zbMATH Keywords: local error
zbMATH Keywords: oscillatory solutions
reviewed by: Josep J. Masdemont
describes a project that uses: mctoolbox
MaRDI profile type: Publication
full work available at URL: j.cam.2011.11.024
OpenAlex ID: W2083707802
cites work: Explicit Runge–Kutta (–Nyström) Methods with Reduced Phase Errors for Computing Oscillating Solutions
cites work: Families of Runge-Kutta-Nystrom Formulae
cites work: High-Order Embedded Runge-Kutta-Nystrom Formulae
cites work: High Phase-Lag-Order Runge--Kutta and Nyström Pairs
cites work: High order Runge--Kutta--Nyström codes for the integration of oscillatory problems.
cites work: Algorithm 670: a Runge-Kutta-Nyström code
cites work: Accuracy and Stability of Numerical Algorithms
cites work: N-body simulations
DOI: 10.1016/J.CAM.2011.11.024
description: scientific article; zbMATH DE number 6015751