A Dissipation-Preserving Integrator for Damped Oscillatory Hamiltonian Systems

Abstract

In this paper, based on discrete gradient, a dissipation-preserving integrator for weakly dissipative perturbations of oscillatory Hamiltonian system is established. The solution of this system is a damped nonlinear oscillator. Basically, lots of nonlinear oscillatory mechanical systems including frictional forces lend themselves to this approach. The new integrator gives a discrete analogue of the dissipation property of the original system. Meanwhile, since the integrator is based on the variation-of-constants formula for oscillatory systems, it preserves the oscillatory structure of the system. Some properties of the new integrator are derived. The convergence is analyzed for the implicit iterations based on the discrete gradient integrator, and it turns out that the convergence of the implicit iterations based on the new integrator is independent of $\|M\|$, where $M$ governs the main oscillation of the system and usually $\|M\|\gg1$. This significant property shows that a larger stepsize can be chosen for the new schemes than that for the traditional discrete gradient integrators when applied to the oscillatory Hamiltonian system. Numerical experiments are carried out to show the effectiveness and efficiency of the new integrator in comparison with the traditional discrete gradient methods in the scientific literature.