Conservative finite difference schemes for dynamical systems

Abstract

The article presents the implementation of one of the approaches to the integration of dynamical systems, which preserves algebraic integrals in the original fdm for Sage system. This approach, which goes back to the paper by del Buono and Mastroserio, makes it possible, based on any two explicit difference schemes, including any two explicit Runge-Kutta schemes, to construct a new numerical algorithm for integrating a dynamical system that preserves the given integral. This approach has been implemented and tested in the original fdm for Sage system. Details and implementation difficulties are discussed. For testing, two Runge-Kutta schemes were taken having the same order, but different Butcher tables, which does not complicate the method due to paralleling. Two examples are considered - a linear oscillator and a Jacobi oscillator with two quadratic integrals. The second example shows that the preservation of one integral of motion does not lead to the conservation of the other. Moreover, this method allows us to propose a practical application of the well-known ambiguity in the definition of Butcher tables.