Composite symmetric general linear methods (COSY-GLMs) for the long-time integration of reversible Hamiltonian systems

Terence J. T. Norton

Research output: Contribution to journalArticle

28 Downloads (Pure)

Abstract

In choosing a numerical method for the long-time integration of reversible Hamiltonian systems one must take into consideration several key factors: order of the method, ability to preserve invariants of the system, and efficiency of the computation. In this paper, 6th-order composite symmetric general linear methods (COSY-GLMs) are constructed using a generalisation of the composition theory associated with Runge–Kutta methods (RKMs). A novel aspect of this approach involves a nonlinear transformation which is used to convert the GLM to a canonical form in which its starting and finishing methods are trivial. Numerical experiments include efficiency comparisons to symmetric diagonally-implicit RKMs, where it is shown that COSY-GLMs of the same order typically require half the number of function evaluations, as well as long-time computations of both separable and non-separable Hamiltonian systems which demonstrate the preservation properties of the new methods.
Original languageEnglish
Pages (from-to)397-421
Number of pages25
JournalBIT Numerical Mathematics
Volume58
Issue number2
Early online date10 Jan 2018
DOIs
Publication statusPublished - 1 Jun 2018

Fingerprint

General Linear Methods
Reversible Systems
Hamiltonians
Runge Kutta methods
Time Integration
Hamiltonian Systems
Composite
Function evaluation
Composite materials
Implicit Runge-Kutta Methods
Numerical methods
Nonlinear Transformation
Nonseparable
Canonical form
Evaluation Function
Runge-Kutta Methods
Preservation
Convert
Trivial
Chemical analysis

Keywords

  • General linear methods
  • Reversible Hamiltonian systems
  • Symmetric composition

ASJC Scopus subject areas

  • Software
  • Computer Networks and Communications
  • Computational Mathematics
  • Applied Mathematics

Cite this

Composite symmetric general linear methods (COSY-GLMs) for the long-time integration of reversible Hamiltonian systems. / Norton, Terence J. T.

In: BIT Numerical Mathematics, Vol. 58, No. 2, 01.06.2018, p. 397-421.

Research output: Contribution to journalArticle

@article{5c564fa043af4923b1fd6e7955e98998,
title = "Composite symmetric general linear methods (COSY-GLMs) for the long-time integration of reversible Hamiltonian systems",
abstract = "In choosing a numerical method for the long-time integration of reversible Hamiltonian systems one must take into consideration several key factors: order of the method, ability to preserve invariants of the system, and efficiency of the computation. In this paper, 6th-order composite symmetric general linear methods (COSY-GLMs) are constructed using a generalisation of the composition theory associated with Runge–Kutta methods (RKMs). A novel aspect of this approach involves a nonlinear transformation which is used to convert the GLM to a canonical form in which its starting and finishing methods are trivial. Numerical experiments include efficiency comparisons to symmetric diagonally-implicit RKMs, where it is shown that COSY-GLMs of the same order typically require half the number of function evaluations, as well as long-time computations of both separable and non-separable Hamiltonian systems which demonstrate the preservation properties of the new methods.",
keywords = "General linear methods, Reversible Hamiltonian systems, Symmetric composition",
author = "Norton, {Terence J. T.}",
year = "2018",
month = "6",
day = "1",
doi = "10.1007/s10543-017-0692-7",
language = "English",
volume = "58",
pages = "397--421",
journal = "BIT Numerical Mathematics",
issn = "0006-3835",
publisher = "Springer Netherlands",
number = "2",

}

TY - JOUR

T1 - Composite symmetric general linear methods (COSY-GLMs) for the long-time integration of reversible Hamiltonian systems

AU - Norton, Terence J. T.

PY - 2018/6/1

Y1 - 2018/6/1

N2 - In choosing a numerical method for the long-time integration of reversible Hamiltonian systems one must take into consideration several key factors: order of the method, ability to preserve invariants of the system, and efficiency of the computation. In this paper, 6th-order composite symmetric general linear methods (COSY-GLMs) are constructed using a generalisation of the composition theory associated with Runge–Kutta methods (RKMs). A novel aspect of this approach involves a nonlinear transformation which is used to convert the GLM to a canonical form in which its starting and finishing methods are trivial. Numerical experiments include efficiency comparisons to symmetric diagonally-implicit RKMs, where it is shown that COSY-GLMs of the same order typically require half the number of function evaluations, as well as long-time computations of both separable and non-separable Hamiltonian systems which demonstrate the preservation properties of the new methods.

AB - In choosing a numerical method for the long-time integration of reversible Hamiltonian systems one must take into consideration several key factors: order of the method, ability to preserve invariants of the system, and efficiency of the computation. In this paper, 6th-order composite symmetric general linear methods (COSY-GLMs) are constructed using a generalisation of the composition theory associated with Runge–Kutta methods (RKMs). A novel aspect of this approach involves a nonlinear transformation which is used to convert the GLM to a canonical form in which its starting and finishing methods are trivial. Numerical experiments include efficiency comparisons to symmetric diagonally-implicit RKMs, where it is shown that COSY-GLMs of the same order typically require half the number of function evaluations, as well as long-time computations of both separable and non-separable Hamiltonian systems which demonstrate the preservation properties of the new methods.

KW - General linear methods

KW - Reversible Hamiltonian systems

KW - Symmetric composition

UR - http://www.scopus.com/inward/record.url?scp=85040322923&partnerID=8YFLogxK

U2 - 10.1007/s10543-017-0692-7

DO - 10.1007/s10543-017-0692-7

M3 - Article

VL - 58

SP - 397

EP - 421

JO - BIT Numerical Mathematics

JF - BIT Numerical Mathematics

SN - 0006-3835

IS - 2

ER -