Do you want to publish a course? Click here

Comparative study of variational chaos indicators and ODEs numerical integrators

196   0   0.0 ( 0 )
 Added by Luciano Darriba
 Publication date 2012
  fields Physics
and research's language is English




Ask ChatGPT about the research

The reader can find in the literature a lot of different techniques to study the dynamics of a given system and also, many suitable numerical integrators to compute them. Notwithstanding the recent work of Maffione et al. (2011a) for mappings, a detailed comparison among the widespread indicators of chaos in a general system is still lacking. Such a comparison could lead to select the most efficient algorithms given a certain dynamical problem. Furthermore, in order to choose the appropriate numerical integrators to compute them, more comparative studies among numerical integrators are also needed. This work deals with both problems. We first extend the work of Maffione et al. (2011) for mappings to the 2D Henon & Heiles (1964) potential, and compare several variational indicators of chaos: the Lyapunov Indicator (LI); the Mean Exponential Growth Factor of Nearby Orbits (MEGNO); the Smaller Alignment Index (SALI) and its generalized version, the Generalized Alignment Index (GALI); the Fast Lyapunov Indicator (FLI) and its variant, the Orthogonal Fast Lyapunov Indicator (OFLI); the Spectral Distance (D) and the Dynamical Spectras of Stretching Numbers (SSNs). We also include in the record the Relative Lyapunov Indicator (RLI), which is not a variational indicator as the others. Then, we test a numerical technique to integrate Ordinary Differential Equations (ODEs) based on the Taylor method implemented by Jorba & Zou (2005) (called taylor), and we compare its performance with other two well-known efficient integrators: the Prince & Dormand (1981) implementation of a Runge-Kutta of order 7-8 (DOPRI8) and a Bulirsch-Stoer implementation. These tests are run under two very different systems from the complexity of their equations point of view: a triaxial galactic potential model and a perturbed 3D quartic oscillator.



rate research

Read More

161 - Ch. Skokos , E. Gerlach 2010
We present and compare different numerical schemes for the integration of the variational equations of autonomous Hamiltonian systems whose kinetic energy is quadratic in the generalized momenta and whose potential is a function of the generalized positions. We apply these techniques to Hamiltonian systems of various degrees of freedom, and investigate their efficiency in accurately reproducing well-known properties of chaos indicators like the Lyapunov Characteristic Exponents (LCEs) and the Generalized Alignment Indices (GALIs). We find that the best numerical performance is exhibited by the textit{`tangent map (TM) method}, a scheme based on symplectic integration techniques which proves to be optimal in speed and accuracy. According to this method, a symplectic integrator is used to approximate the solution of the Hamiltons equations of motion by the repeated action of a symplectic map $S$, while the corresponding tangent map $TS$, is used for the integration of the variational equations. A simple and systematic technique to construct $TS$ is also presented.
It is known that the asymptotic invariant manifolds around an unstable periodic orbit in conservative systems can be represented by convergent series (Cherry 1926, Moser 1956, 1958, Giorgilli 2001). The unstable and stable manifolds intersect at an infinity of homoclinic points, generating a complicated homoclinic tangle. In the case of simple mappings it was found (Da Silva Ritter et al. 1987) that the domain of convergence of the formal series extends to infinity along the invariant manifolds. This allows in practice to study the homoclinic tangle using only series. However in the case of Hamiltonian systems, or mappings with a finite analyticity domain,the convergence of the series along the asymptotic manifolds is also finite. Here, we provide numerical indications that the convergence does not reach any homoclinic points. We discuss in detail the convergence problem in various cases and we find the degree of approximation of the analytical invariant manifolds to the real (numerical) manifolds as i) the order of truncation of the series increases, and ii) we use higher numerical precision in computing the coefficients of the series. Then we introduce a new method of series composition, by using action-angle variables, that allows the calculation of the asymptotic manifolds up to an a arbitrarily large extent. This is the first case of an analytic development that allows the computation of the invariant manifolds and their intersections in a Hamiltonian system for an extent long enough to allow the study of homoclinic chaos by analytical means.
120 - D. J. Albers , J. C. Sprott 2004
This paper examines the most probable route to chaos in high-dimensional dynamical systems in a very general computational setting. The most probable route to chaos in high-dimensional, discrete-time maps is observed to be a sequence of Neimark-Sacker bifurcations into chaos. A means for determining and understanding the degree to which the Landau-Hopf route to turbulence is non-generic in the space of $C^r$ mappings is outlined. The results comment on previous results of Newhouse, Ruelle, Takens, Broer, Chenciner, and Iooss.
We summarize various cases where chaotic orbits can be described analytically. First we consider the case of a magnetic bottle where we have non-resonant and resonant ordered and chaotic orbits. In the sequence we consider the hyperbolic Henon map, where chaos appears mainly around the origin, which is an unstable periodic orbit. In this case the chaotic orbits around the origin are represented by analytic series (Moser series). We find the domain of convergence of these Moser series and of similar series around other unstable periodic orbits. The asymptotic manifolds from the various unstable periodic orbits intersect at homoclinic and heteroclinic orbits that are given analytically. Then we consider some Hamiltonian systems and we find their homoclinic orbits by using a new method of analytic prolongation. An application of astronomical interest is the domain of convergence of the analytical series that determine the spiral structure of barred-spiral galaxies.
The mechanism for transitions from phase to defect chaos in the one-dimensional complex Ginzburg-Landau equation (CGLE) is presented. We introduce and describe periodic coherent structures of the CGLE, called Modulated Amplitude Waves (MAWs). MAWs of various period P occur naturally in phase chaotic states. A bifurcation study of the MAWs reveals that for sufficiently large period P, pairs of MAWs cease to exist via a saddle-node bifurcation. For periods beyond this bifurcation, incoherent near-MAW structures occur which evolve toward defects. This leads to our main result: the transition from phase to defect chaos takes place when the periods of MAWs in phase chaos are driven beyond their saddle-node bifurcation.
comments
Fetching comments Fetching comments
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا