Conserving first integrals under discretization with variable step size integration procedures
conservation of first integralsflowJacobi integralstabilizationthree-body problemvariable step size discretization methods
Dynamics induced by flows and semiflows (37C10) Stability theory for smooth dynamical systems (37C75) Computational methods for ergodic theory (approximation of invariant measures, computation of Lyapunov exponents, entropy, etc.) (37M25) Dynamical systems in classical and celestial mechanics (37N05) Numerical nonlinear stabilities in dynamical systems (65P40) Three-body problems (70F07)
This paper is concerned with the conservation of first integrals of systems of ordinary differential equations when they are integrated by means of variable stepsize discretization methods. The technique proposed by the author in order to achieve the desired stabilization consists in modifying the flow of the given differential system so that the first integrals become exponentially attractive invariant sets in the new flow. Then the author proves that under suitable assumptions on the stepsizes \(h_j\) of a discretization method of order \(p\), such invariant sets are also \({\mathcal O}(h_{\max}^{p+1})\)-positive invariant for the method, where \(h_{\max} = \max h_j\). This implies the \({\mathcal O}(h_{\max}^{p+1})\)-conservation of first integrals for arbitrarily large integration intervals in contrast with the behaviour of standard integrators in which the errors grow along the numerical integration. The author illustrates the stabilization technique with the restricted three-body problem by using the well-known Jacobi integral. It is seen that, with a particular variable stepsize NAG routine, the numerical solution computed without stabilization introduces an error in the Jacobi integral that grows linearly with the time whereas with the stabilized technique it remains almost constant. Further the portrait phase of some trajectories are presented to show the long term effect of the proposed stabilization.
- Solving ODEs numerically while preserving a first integral
- Solving ODEs numerically while preserving a first integral
- Variable step implementation of geometric integrators
- Numerical stabilization of all laws of conservation in the many body problem
- Discretization of dynamical systems with first integrals
- A note on multistep methods and attracting sets of dynamical systems
- Differential Topology
- scientific article; zbMATH DE number 88931 (Why is no real title available?)
- scientific article; zbMATH DE number 979834 (Why is no real title available?)
- scientific article; zbMATH DE number 750193 (Why is no real title available?)
- scientific article; zbMATH DE number 1436449 (Why is no real title available?)
- scientific article; zbMATH DE number 3247443 (Why is no real title available?)
- Stabilization of DAEs and invariant manifolds
- Stable Attracting Sets in Dynamical Systems and in Their One-Step Discretizations
- Variable time step integration with symplectic methods
- One- and multistep discretizations of index 2 differential algebraic systems and their use in optimization
- Time-dependent first integrals, nonlinear dynamical systems, and numerical integration
- scientific article; zbMATH DE number 1014969 (Why is no real title available?)
- scientific article; zbMATH DE number 1016808 (Why is no real title available?)
- Numerical stabilization of all laws of conservation in the many body problem
This page was built for publication: Conserving first integrals under discretization with variable step size integration procedures
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1971850)
