EconPapers    
Economics at your fingertips  
 

Lyapunov and reversibility errors for Hamiltonian flows

F. Panichi and G. Turchetti

Chaos, Solitons & Fractals, 2018, vol. 112, issue C, 83-91

Abstract: We discuss the stability of a Hamiltonian system by comparing the standard Lyapunov error (LE) with the forward error (FE) due to a small random perturbation. We introduce also the reversibility error (RE) where the evolution is computed forward up to time t and backwards to t=0 in presence of noise. This procedure has been investigated in the case of symplectic maps, but it turns out that the results are simpler in the case of a noisy flow, in the limit of zero noise amplitude. Indeed the stochastic processes defined by the displacement of the noisy orbit at time t for FE, or at time 0 for RE after the evolution up to time t, satisfy linear Langevin equations, are Gaussian processes, and the errors are just their root mean square deviations. All the errors are expressed in terms of the fundamental matrix L(t) of the tangent flow and can be evaluated numerically using a symplectic integrator. Letting eL(t) be the Lyapunov error and eR(t) be the reversibility error a very simple relation holds eR2(t)=∫0teL2(s)ds. The integral relation is quite natural since the local errors due to a random perturbations accumulate during the evolution whereas for the Lyapunov case the error is introduced only at time zero and propagated. The plot of errors for initial conditions in a Poincaré section reflects the phase portrait, whereas in the action plane it allows to single out the resonance strips. We have applied the method to a 3D Hamiltonian model H=H0(J)+λV(Θ), where analytic estimates can be obtained for the single resonances from perturbation theory. This allows to inspect the double resonance structure where the single resonance strips intersect. We have also considered the Hénon–Heiles Hamiltonian to show numerically the equivalence of the errors apart from a shift of 1/2 in the power law exponent in the case of regular orbits. The reversibility error method (REM), previously introduced as the error due to round off in the symplectic integration, appears to be comparable with RE also for the models considered here.

Keywords: Hamiltonian systems; Stability analysis; Reversibility; Resonances (search for similar items in EconPapers)
Date: 2018
References: View references in EconPapers View complete reference list from CitEc
Citations:

Downloads: (external link)
http://www.sciencedirect.com/science/article/pii/S0960077918301152
Full text for ScienceDirect subscribers only

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:eee:chsofr:v:112:y:2018:i:c:p:83-91

DOI: 10.1016/j.chaos.2018.03.019

Access Statistics for this article

Chaos, Solitons & Fractals is currently edited by Stefano Boccaletti and Stelios Bekiros

More articles in Chaos, Solitons & Fractals from Elsevier
Bibliographic data for series maintained by Thayer, Thomas R. ().

 
Page updated 2025-03-19
Handle: RePEc:eee:chsofr:v:112:y:2018:i:c:p:83-91