Exactly Conservative Integrators

DOI10.1137/S0036139995289313MaRDI QIDQ4258445zbMATH OpenOpenAlexFDO

Authors B. A. Shadwick, John C. Bowman, Philip J. Morrison

Publication date 13 September 1999

Published in SIAM Journal on Applied Mathematics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/chao-dyn/9507012, https://digital.library.unt.edu/ark:/67531/metadc691320/

zbMATH Keywords

algorithmsconservative systemsEuler schemeHamiltonianKepler problemLotka-Volterra predator-prey modelnonlinear invariantspredictor-corrector schemesymplecticthree-waves problem

Mathematics Subject Classification ID

Topological and differentiable equivalence, conjugacy, moduli, classification of dynamical systems (37C15) Discretization methods and integrators (symplectic, variational, geometric, etc.) for dynamical systems (37M15) Numerical methods for Hamiltonian systems including symplectic integrators (65P10)

Abstract: Traditional numerical discretizations of conservative systems generically yield an artificial secular drift of any nonlinear invariants. In this work we present an explicit nontraditional algorithm that exactly conserves these invariants. We illustrate the general method by applying it to the three-wave truncation of the Euler equations, the Lotka--Volterra predator--prey model, and the Kepler problem. This method is discussed in the context of symplectic (phase space conserving) integration methods as well as nonsymplectic conservative methods. We comment on the application of our method to general conservative systems.

Review text

The paper deals with numerical discretization methods for conservative systems. Traditional schemes generically predict artificial secular drifts of nonlinear invariants. The authors present, for three examples, explicit algorithms which preserve exactly such invariants. The idea is to find a right perturbation or correction of the former system which both vanishes when the time discretization step tends to 0 and preserves the invariants of the system. The choice of such a perturbation derives from the very scheme of approximation and depends on the time discretization step. They apply this technique to the three-waves problem. They point out that for the Lotka-Volterra predator-prey model one can modify the classical predictor-corrector scheme so as to obtain a new conservative predictor-corrector scheme. Looking at the Kepler problem, they obtain the suitable correction, which keeps the Hamiltonian invariant, by expressing the problem in new variables to which they apply a modified Euler scheme.

Recommendations

• Conservative integrators for many-body problems
• Integral-preserving integrators
• Line integral methods for conservative problems
• Exponential integrators preserving first integrals or Lyapunov functions for conservative or dissipative systems
• scientific article; zbMATH DE number 3972116
• Conservative numerical methods for integrating differential equations in nonlinear mechanics
• scientific article; zbMATH DE number 1302556
• On integrable conservation laws
• Conservative methods for dynamical systems
• Line integral methods which preserve all invariants of conservative problems

Cited in

(18)

• A survey of shadowing methods for numerical solutions of ordinary differential equations
• An exact energy conservation property of the quantum lattice Boltzmann algorithm
• Explicit exactly energy-conserving methods for Hamiltonian systems
• Conservative integrators for many-body problems
• Invariantization of numerical schemes using moving frames
• A new conservative numerical integration algorithm for the three-dimensional Kepler motion based on the Kustaanheimo-Stiefel regularization theory
• Links between dissipation, intermittency, and helicity in the GOY model revisited
• Locally exact modifications of numerical schemes
• Conservative integrators for vortex blob methods on the plane
• Numerical analysis of nonlinear conservative oscillators
• Numerical methods for non conservative perturbations of conservative problems
• Integral-preserving integrators
• Conservative methods for dynamical systems
• On the arbitrarily long-term stability of conservative methods
• Volume-preserving maps with an invariant
• Conservative integrators for piecewise smooth systems with transversal dynamics
• Improving the accuracy of the Newmark method through backward error analysis
• From chaos of lines to Lagrangian structures in flux conservative fields

This page was built for publication: Exactly Conservative Integrators

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4258445)