Second Order Conformal Symplectic Schemes For Damped Hamiltonian Systems

Keywords

Birkhoffian method; Conformal symplectic; Dissipation preservation; Linear damping; Structure-preserving algorithm

Abstract

Numerical methods for solving linearly damped Hamiltonian systems are constructed using the popular Störmer–Verlet and implicit midpoint methods. Each method is shown to preserve dissipation of symplecticity and dissipation of angular momentum of an N-body system with pairwise distance dependent interactions. Necessary and sufficient conditions for second order accuracy are derived. Analysis for linear equations gives explicit relationships between the damping parameter and the step size to reveal when the methods are most advantageous; essentially, the damping rate of the numerical solution is exactly preserved under these conditions. The methods are applied to several model problems, both ODEs and PDEs. Additional structure preservation is discovered for the discretized PDEs, in one case dissipation in total linear momentum and in another dissipation in mass are preserved by the methods. The numerical results, along with comparisons to standard Runge–Kutta methods and another structure-preserving method, demonstrate the usefulness and strengths of the methods.

Publication Date

3-1-2016

Publication Title

Journal of Scientific Computing

Volume

66

Issue

3

Number of Pages

1234-1259

Document Type

Article

Personal Identifier

scopus

DOI Link

https://doi.org/10.1007/s10915-015-0062-z

Socpus ID

84957428627 (Scopus)

Source API URL

https://api.elsevier.com/content/abstract/scopus_id/84957428627

This document is currently not available here.

Share

COinS