Lyapunov and reversibility errors for Hamiltonian flows.
Saved in:
| Title: | Lyapunov and reversibility errors for Hamiltonian flows. |
|---|---|
| Authors: | Panichi, F.1 federico.panichi@stud.usz.edu.pl, Turchetti, G.2 turchetti@bo.infn.it |
| Source: | Chaos, Solitons & Fractals. Jul2018, Vol. 112, p83-91. 9p. |
| Subjects: | Lyapunov exponents, Differential equations, Numerical analysis, Hamiltonian systems, Chaos theory |
| 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 e L ( t ) be the Lyapunov error and e R ( t ) be the reversibility error a very simple relation holds e R 2 ( t ) = ∫ 0 t e L 2 ( s ) d s . 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 = H 0 ( 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. [ABSTRACT FROM AUTHOR] |
| Copyright of Chaos, Solitons & Fractals is the property of Pergamon Press - An Imprint of Elsevier Science and its content may not be copied or emailed to multiple sites without the copyright holder's express written permission. Additionally, content may not be used with any artificial intelligence tools or machine learning technologies. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.) | |
| Database: | Engineering Source |
| FullText | Text: Availability: 0 |
|---|---|
| Header | DbId: egs DbLabel: Engineering Source An: 129923470 AccessLevel: 6 PubType: Periodical PubTypeId: serialPeriodical PreciseRelevancyScore: 0 |
| IllustrationInfo | |
| Items | – Name: Title Label: Title Group: Ti Data: Lyapunov and reversibility errors for Hamiltonian flows. – Name: Author Label: Authors Group: Au Data: <searchLink fieldCode="AR" term="%22Panichi%2C+F%2E%22">Panichi, F.</searchLink><relatesTo>1</relatesTo><i> federico.panichi@stud.usz.edu.pl</i><br /><searchLink fieldCode="AR" term="%22Turchetti%2C+G%2E%22">Turchetti, G.</searchLink><relatesTo>2</relatesTo><i> turchetti@bo.infn.it</i> – Name: TitleSource Label: Source Group: Src Data: <searchLink fieldCode="JN" term="%22Chaos%2C+Solitons+%26+Fractals%22">Chaos, Solitons & Fractals</searchLink>. Jul2018, Vol. 112, p83-91. 9p. – Name: Subject Label: Subjects Group: Su Data: <searchLink fieldCode="DE" term="%22Lyapunov+exponents%22">Lyapunov exponents</searchLink><br /><searchLink fieldCode="DE" term="%22Differential+equations%22">Differential equations</searchLink><br /><searchLink fieldCode="DE" term="%22Numerical+analysis%22">Numerical analysis</searchLink><br /><searchLink fieldCode="DE" term="%22Hamiltonian+systems%22">Hamiltonian systems</searchLink><br /><searchLink fieldCode="DE" term="%22Chaos+theory%22">Chaos theory</searchLink> – Name: Abstract Label: Abstract Group: Ab Data: 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 e L ( t ) be the Lyapunov error and e R ( t ) be the reversibility error a very simple relation holds e R 2 ( t ) = ∫ 0 t e L 2 ( s ) d s . 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 = H 0 ( 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. [ABSTRACT FROM AUTHOR] – Name: AbstractSuppliedCopyright Label: Group: Ab Data: <i>Copyright of Chaos, Solitons & Fractals is the property of Pergamon Press - An Imprint of Elsevier Science and its content may not be copied or emailed to multiple sites without the copyright holder's express written permission. Additionally, content may not be used with any artificial intelligence tools or machine learning technologies. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract.</i> (Copyright applies to all Abstracts.) |
| PLink | https://search.ebscohost.com/login.aspx?direct=true&site=eds-live&db=egs&AN=129923470 |
| RecordInfo | BibRecord: BibEntity: Identifiers: – Type: doi Value: 10.1016/j.chaos.2018.03.019 Languages: – Code: eng Text: English PhysicalDescription: Pagination: PageCount: 9 StartPage: 83 Subjects: – SubjectFull: Lyapunov exponents Type: general – SubjectFull: Differential equations Type: general – SubjectFull: Numerical analysis Type: general – SubjectFull: Hamiltonian systems Type: general – SubjectFull: Chaos theory Type: general Titles: – TitleFull: Lyapunov and reversibility errors for Hamiltonian flows. Type: main BibRelationships: HasContributorRelationships: – PersonEntity: Name: NameFull: Panichi, F. – PersonEntity: Name: NameFull: Turchetti, G. IsPartOfRelationships: – BibEntity: Dates: – D: 01 M: 07 Text: Jul2018 Type: published Y: 2018 Identifiers: – Type: issn-print Value: 09600779 Numbering: – Type: volume Value: 112 Titles: – TitleFull: Chaos, Solitons & Fractals Type: main |
| ResultId | 1 |