Comparison of Second Order Conformal Symplectic Schemes with Linear Stability Analysis

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.

Author dc:creator

• Floyd, Dwayne

Contributors dc:contributor

• Moore, Brian

• Conformal symplectic; structure preserving; linear stability
• Mathematics
• Dissertations, Academic -- Sciences; Sciences -- Dissertations, Academic

Language dc:language

English

Identifier

CFE0005793

OAI identifier oai:identifier

oai:stars.library.ucf.edu:etd-1667