Projection methods and discrete gradient methods for preserving first integrals of ODEs
Richard A. Norton David I. McLaren G. R. W. Quispel Ari Stern Antonella Zanna
Discrete & Continuous Dynamical Systems - A 2015, 35(5): 2079-2098 doi: 10.3934/dcds.2015.35.2079
In this paper we study linear projection methods for approximating the solution and simultaneously preserving first integrals of autonomous ordinary differential equations. We show that each (linear) projection method is equivalent to a class of discrete gradient methods, in both single and multiple first integral cases, and known results for discrete gradient methods also apply to projection methods. Thus we prove that in the single first integral case, under certain mild conditions, the numerical solution for a projection method exists and is locally unique, and preserves the order of accuracy of the underlying method. Our results allow considerable freedom for the choice of projection direction and do not have a time step restriction close to critical points.
keywords: discrete gradients Hamiltonian systems. projection energy preserving integrators Geometric integration

Year of publication

Related Authors

Related Keywords

[Back to Top]