Arxiv:1911.08207V1 [Nlin.CD] 19 Nov 2019
Total Page:16
File Type:pdf, Size:1020Kb
Power-law trapping in the volume-preserving Arnold-Beltrami-Childress map Swetamber Das1 and Arnd B¨acker2, 1 1Max-Planck-Institut f¨urPhysik komplexer Systeme, N¨othnitzerStraße 38, 01187 Dresden, Germany 2Technische Universit¨atDresden, Institut f¨urTheoretische Physik and Center for Dynamics, 01062 Dresden, Germany (Dated: November 20, 2019) Understanding stickiness and power-law behavior of Poincar´erecurrence statistics is an open problem for higher-dimensional systems, in contrast to the well-understood case of systems with two degrees-of-freedom. We study such intermittent behavior of chaotic orbits in three-dimensional volume-preserving systems using the example of the Arnold-Beltrami-Childress map. The map has a mixed phase space with a cylindrical regular region surrounded by a chaotic sea for the considered parameters. We observe a characteristic overall power-law decay of the cumulative Poincar´erecur- rence statistics with significant oscillations superimposed. This slow decay is caused by orbits which spend long times close to the surface of the regular region. Representing such long-trapped orbits in frequency space shows clear signatures of partial barriers and reveals that coupled resonances play an essential role. Using a small number of the most relevant resonances allows for classifying long-trapped orbits. From this the Poincar´erecurrence statistics can be divided into different ex- ponentially decaying contributions which very accurately explains the overall power-law behavior including the oscillations. I. INTRODUCTION Thus already purely for topological reasons additional types of transport are possible, leading for example to Hamiltonian systems generically have a mixed phase the famous Arnold diffusion [47{49]. space in which regular and chaotic motion coexist. A Still, power-law trapping is observed in higher- typical chaotic trajectory is often trapped intermittently dimensional systems [20, 28, 29, 33, 35], but not due to a close to regular regions of phase space for arbitrary long generalized hierarchy [33, 35], so that the mechanism of but finite times. The existence of such sticky behavior is stickiness is different from that in two-dimensional area- known to have important consequences for chaotic trans- preserving maps. While some generalizations have been port, e.g. the existence of anomalous kinetics, L´evypro- proposed [29, 50{52], a profound understanding of power- cesses, and L´evyflights and finds interesting applications law trapping in higher-dimensional systems remains a sig- in a variety of physical contexts such as chaotic advec- nificant open question [46]. tion in fluids [1], in the three-body problem [2], driven This motivates to investigate power-law trapping for coupled Morse oscillators [3], confinement in particle ac- 3d volume-preserving maps, whether features of 2d maps celerators [4], and plasma physics [5], to name a few. The still apply, while those present for 4d maps might already behavior of stickiness in a mixed phase-space can be char- become relevant. For example in a 3d volume-preserving acterized by the cumulative Poincar´erecurrence statistics map there can be 2d regular tori which therefore are char- P (t), which is the probability that a chaotic orbit has not acterized by two frequencies, as in 4d symplectic maps. returned to an initial region until time t [6]. Fully chaotic On the other hand, the 2d tori can provide absolute bar- systems typically show an exponential decay of P (t) [7{ riers to the motion in a 3d volume-preserving map as in 10]. However, systems with a mixed phase space usually 2d area-preserving maps. The dynamical properties of display a much slower decay, usually well-described by a 3d volume-preserving maps have been studied in much power-law [2, 11{35]. This so-called power-law trapping detail, see e.g. [53{69]. However, there are only a few has striking implications for transport in many physical studies of power-law trapping in 3d volume-preserving systems [2, 3, 36{39]. maps, see Refs. [70{72]. In Ref. [70] it is argued that arXiv:1911.08207v1 [nlin.CD] 19 Nov 2019 Power-law trapping is well-understood for Hamilto- stickiness in their model originates around the hyperbolic nian systems and maps with two degrees-of-freedom invariant sets embedded in the sticky domain. In Ref. [71] [15, 17, 32, 40{46]. For example in area-preserving maps, recurrence times statistics of an extended standard map regular tori are one-dimensional and therefore act as com- have been analyzed. Power-law decay in the presence of plete barriers to transport in phase space. If such tori accelerator modes is found in Ref. [72]. break up they may turn into so-called cantori which pro- In this paper, we investigate power-law trapping in vide partial barriers with slow transport across them. three-dimensional volume-preserving maps using the spe- Usually there is a whole hierarchy of such partial bar- cific example of the Arnold-Beltrami-Childress (ABC) riers which is the origin of the power-law trapping. In map. The ABC map is a discretized version of the ABC contrast, in higher-dimensional systems regular tori have flows with time-periodic forcing [73, 74]. The ABC flows an insufficient dimension to provide complete barriers to are a class of non-turbulent flows which appear as a so- transport. For example for a 4d symplectic map, the reg- lution of the three-dimensional Euler equation and is ular tori are two-dimensional and therefore cannot divide used as a simple model for the study of chaos in three- the phase space into dynamically disconnected regions. dimensional steady incompressible flows. The 3d phase 2 space of the flow and the corresponding map is known choice of the real parameters A, B, and C. The dynam- to have invariant surfaces and chaotic streamlines. The ics of the system can be very different, showing a mixed flow was first introduced by Arnold [75] and later Chil- phase space or chaotic motion [55, 74]. For the rest of dress [76] demonstrated that these are important models the paper, we fix the parameters (A; B) = (2:0; 1:5) and for the kinematic dynamo effect in astrophysical plasma, consider different values of C. i.e., a magnetic field can be generated and sustained For C = 0 the dynamics in (x; z) becomes independent by the motion of an electrically conducting fluid. The of the motion in the y-coordinate and therefore reduces flow also serves as a prototypical model of a force-free to a 2d area-preserving subsystem, magnetic field to study spatial diffusion of charged par- ticles [77]. The ABC map exhibits most of the basic features of a typical three-dimensional, time-periodic, xn+1 = xn + A sin(zn) volume-preserving ABC flow. The map has a mixed (2) z = z + B cos(x ): phase space with periodic, quasi-periodic, and chaotic n+1 n n+1 behavior in the 3d phase space. For specific parameters of the map, the phenomenon of resonance-induced diffu- Figure 1 shows the 2d phase space for this map for sion has been observed which leads to global transport in (A; B) = (2:0; 1:5), which has the typical structure of phase space [74]. The map has been employed as a rep- a mixed phase space with regular and chaotic motion. resentative of volume-preserving flows to study chaotic The fixed point at (x; z) = (3π=2; π) is elliptic and there- advection of inertial particles in bailout embedding mod- fore surrounded by regular orbits forming invariant tori els [78{80]. (full curves). Embedded within this regular region one The outline of the paper is as follows. Section II de- has a prominent period-4 regular island. In addition to scribes the ABC map, its basic properties, and illustrates this central regular island one has a symmetry-related is- regular orbits in the 3d phase space in Sec. II A. The land around the elliptic fixed point at (x; z) = (π=2; 0). variation of the regular region is illustrated in Sec. II B Outside of these regular islands one has a large area with using a sequence of 2d plots of the finite-time Lyapunov predominantly chaotic motion, as illustrated by a longer exponent. The frequency space of the system is intro- orbit in the figure. Due to the symmetry, dynamical duced in Sec. II C and a boundary cylinder, separating properties, such as stickiness at the islands are identi- regular and chaotic motion, is identified. In Sec. III, cal, so that we can restrict to the central region around the stickiness near the regular region is studied using the fixed point at (x; z) = (3π=2; π). the Poincar´erecurrence statistics and power-law trapping has been explained in Sec. III A. As illustration, one long- trapped orbit is displayed in the 3d phase space, on a 2d Poincar´esection, and in frequency-time representation in 2π Sec. III B. A quantitative analysis of the approach to the boundary cylinder is done in Sec. III C, demonstrating z the significant role played by coupled resonances. Fur- 3 thermore, in Sec. III D we use the frequency analysis to π identify a small number of the most relevant resonances. 2 This allows for classifying long-trapped orbits and by this to divide the Poincar´erecurrence statistics into different exponentially decaying contributions. This very accu- rately explains the overall power-law behavior including π the oscillations. In Sec. III E possible partial barriers are briefly discussed. Section IV gives a summary and out- look. 1 π 2 II. ARNOLD-BELTRAMI-CHILDRESS MAP The ABC map [74] is the discrete-time dynamical sys- tem (xn; yn; zn) (xn+1; yn+1; zn+1) given by 0 7! 1 3 x = x + A sin(z ) + C cos(y ) 0 π π π 2π n+1 n n n 2 2 x yn+1 = yn + B sin(xn+1) + A cos(zn) (1) zn+1 = zn + C sin(yn+1) + B cos(xn+1); FIG.