Stabilized integration of Hamiltonian systems with hard-sphere inequality constraints

Stephen D. Bond, Benedict J. Leimkuhler

Research output: Contribution to journalArticlepeer-review

Abstract

We consider numerical methods for resolving the dynamics of a Hamiltonian N-body problem subject to hard-sphere inequality constraints. The dynamics of these mixed systems consists of smooth flow of a Hamiltonian system between collisions with an impulsive momentum exchange at the points of collision. The inclusion of these impulses makes traditional backward error analysis inappropriate since the flow is discontinuous and cannot be interpreted using a single modified smooth Hamiltonian system. We introduce two methods which respect the underlying modified smooth Hamiltonian system through the use of a modified map and collision operator at points of collision. In numerical experiments, these new methods show dramatically improved energy conservation over long time intervals.

Original languageEnglish
Pages (from-to)134-147
Number of pages14
JournalSIAM Journal on Scientific Computing
Volume30
Issue number1
DOIs
Publication statusPublished - 28 Nov 2007

Keywords

  • Backward error analysis
  • Collision Verlet
  • Geometric integration
  • Hamiltonian systems
  • Hard spheres
  • Impact dynamics
  • Inequality constraints

Fingerprint

Dive into the research topics of 'Stabilized integration of Hamiltonian systems with hard-sphere inequality constraints'. Together they form a unique fingerprint.

Cite this