skip to: onlinetools | mainnavigation | content | footer

Stephen D. Bond

Citation:

Stephen D. Bond and Benedict J. Leimkuhler, Stabilized Integration of Hamiltonian Systems with Hard-Sphere Inequality Constraints, SIAM Journal on Scientific Computing 30:1 (2007) 134-147.

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.

Full Text:

DOI: 10.1137/06066552X

AIP: SCE/30/134

Preprint:

Stabilized Integration of Hamiltonian Systems with Hard-Sphere Inequality Constraints (248K, PDF)

Bibtex:
@article{BoLe2007,
  author  = {Stephen D. Bond and Benedict J. Leimkuhler},
  title   = {Stabilized Integration of Hamiltonian Systems with Hard-Sphere
             Inequality Constraints},
  journal = {SIAM Journal on Scientific Computing},
  volume  = 30,
  number  = 1,
  year    = 2007,
  pages   = {134--147},
  url     = {http://link.aip.org/link/?SCE/30/134},
  doi     = {10.1137/06066552X}
}

Stephen Bond