Geometric integrators for piecewise smooth Hamiltonian systems

摘要

In this paper, we consider C-1,C-1 Hamiltonian systems. We prove the existence of a first derivative of the flow with respect to initial values and show that it satisfies the symplecticity condition almost everywhere in the phase-space. In a second step, we present a geometric integrator for such systems (called the SDH method) based on B-splines interpolation and a splitting method introduced by McLachlan and Quispel [Appl. Numer. Math. 45 (2003) 411-418], and we prove it is convergent, and that it preserves the energy and the volume.