Abstract

In this paper, we construct an integrator that conserves volume in phase space. We compare the results obtained using this method and a symplectic integrator. The results of our experiments do not reveal any superiority of the symplectic over strictly volume-preserving integrators. We also investigate the effect of numerically conserving energy in a numerical process by rescaling velocities to keep energy constant at every step. Our results for Henon-Heiles problem show that keeping energy constant in this way destroys ergodicity and forces the solution onto a periodic orbit.

Department(s)

Computer Science

Comments

This work was supported in part by NSF Grant DMS 90 15533 while the author was at the University of Illinois at Champaign-Urbana

Keywords and Phrases

Hamiltonian Systems; Energy Conservation; Symplectic Integrators; Louville Integrators

Report Number

CSC-93-26

Document Type

Technical Report

Document Version

Final Version

File Type

text

Language(s)

English

Rights

© 1993 University of Missouri--Rolla, All rights reserved.

Publication Date

29 Sep 1993

Share

 
COinS