THE CONTROL OF PARASITISM IN G-SYMPLECTIC METHODS

G-symplectic general linear methods are designed to approximately preserve symplectic invariants for Hamiltonian systems. In this paper, the properties of G-symplectic methods are explored computationally and theoretically. Good preservation properties are observed over long times for many parameter...

Full description

Saved in:
Bibliographic Details
Published inSIAM journal on numerical analysis Vol. 52; no. 5; pp. 2440 - 2465
Main Authors BUTCHER, JOHN C., HABIB, YOUSAF, HILL, ADRIAN T., NORTON, TERENCE J. T.
Format Journal Article
LanguageEnglish
Published Society for Industrial and Applied Mathematics 01.01.2014
Subjects
Online AccessGet full text

Cover

Loading…
More Information
Summary:G-symplectic general linear methods are designed to approximately preserve symplectic invariants for Hamiltonian systems. In this paper, the properties of G-symplectic methods are explored computationally and theoretically. Good preservation properties are observed over long times for many parameter ranges, but, for other parameter values, the parasitic behavior, to which multivalue methods are prone, corrupts the numerical solution by the growth of small perturbations. Two approaches for alleviating this effect are considered. First, compositions of methods with growth parameters of opposite signs can be used to cancel the long-term effect of parasitism. Second, methods can be constructed for which the growth parameters are zero by design. Each of these remedies is found to be successful in eliminating parasitic behavior in long-term simulations using a variety of test problems.
Bibliography:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 23
ISSN:0036-1429
1095-7170
DOI:10.1137/140953277