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...
Saved in:
Published in | SIAM journal on numerical analysis Vol. 52; no. 5; pp. 2440 - 2465 |
---|---|
Main Authors | , , , |
Format | Journal Article |
Language | English |
Published |
Society for Industrial and Applied Mathematics
01.01.2014
|
Subjects | |
Online Access | Get full text |
Cover
Loading…
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 |