Symmetry preserving discretization of ordinary differential equations. Large symmetry groups and higher order equations

TL;DR

This paper develops symmetry-preserving discretizations of high-order ordinary differential equations invariant under specific Lie groups, demonstrating improved numerical solutions near singularities compared to standard methods.

Contribution

It introduces invariant discretizations of ODEs under large symmetry groups, extending the applicability of symmetry-preserving numerical schemes to higher-order equations.

Findings

Invariant schemes perform as well as Runge-Kutta methods

Invariant schemes excel near solution singularities

Constructs invariant O$ riangle$Ss for complex symmetry groups

Abstract

Ordinary differential equations (ODEs) and ordinary difference systems (O$\Delta$Ss) invariant under the actions of the Lie groups $\mathrm{SL}_x(2)$, $\mathrm{SL}_y(2)$ and $\mathrm{SL}_x(2)\times\mathrm{SL}_y(2)$ of projective transformations of the independent variables $x$ and dependent variables $y$ are constructed. The ODEs are continuous limits of the O$\Delta$Ss, or conversely, the O$\Delta$Ss are invariant discretizations of the ODEs. The invariant O$\Delta$Ss are used to calculate numerical solutions of the invariant ODEs of order up to five. The solutions of the invariant numerical schemes are compared to numerical solutions obtained by standard Runge-Kutta methods and to exact solutions, when available. The invariant method performs at least as well as standard ones and much better in the vicinity of singularities of solutions.