Numerical methods for solving linearly damped Hamiltonian ordinary differential equations are analyzed and compared. The methods are constructed from the well-known Störmer-Verlet and implicit midpoint methods. The structure preservation properties of each method are shown analytically and numerically. Each method is shown to preserve a symplectic form up to a constant and are therefore conformal symplectic integrators, with each method shown to accurately preserve the rate of momentum dissipation. An analytical linear stability analysis is completed for each method, establishing thresholds between the value of the damping coefficient and the step-size that ensure stability. The methods are all second order and the preservation of the rate of energy dissipation is compared to that of a third order Runge-Kutta method that does not preserve conformal properties. Numerical experiments will include the damped harmonic oscillator and the damped nonlinear pendulum.
| Identifer | oai:union.ndltd.org:ucf.edu/oai:stars.library.ucf.edu:etd-1667 |
| Date | 01 January 2014 |
| Creators | Floyd, Dwayne |
| Publisher | STARS |
| Source Sets | University of Central Florida |
| Language | English |
| Detected Language | English |
| Type | text |
| Format | application/pdf |
| Source | Electronic Theses and Dissertations |
Page generated in 0.0016 seconds