Stabilized integration of Hamiltonian systems with hard-sphere inequality constraints

Stephen D. Bond, Benedict J. Leimkuhler

Research output: Contribution to journalArticlepeer-review

Abstract / Description of output

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
Issue number1
Publication statusPublished - 28 Nov 2007

Keywords / Materials (for Non-textual outputs)

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


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

Cite this