Homoclinic Orbit and Hidden Attractor in the Lorenz-Like System Describing

Total Page:16

File Type:pdf, Size:1020Kb

Homoclinic Orbit and Hidden Attractor in the Lorenz-Like System Describing Homoclinic orbit and hidden attractor in the Lorenz-like system describing the fluid convection motion in the rotating cavity G. A. Leonov,1 N. V. Kuznetsov,1,2, ∗ and T. N. Mokaev1, 2 1Faculty of Mathematics and Mechanics, St. Petersburg State University, Peterhof, St. Petersburg, Russia 2Department of Mathematical Information Technology, University of Jyv¨askyl¨a, Jyv¨askyl¨a, Finland (Dated: March 20, 2018) In this paper a Lorenz-like system, describing the process of rotating fluid convection, is consid- ered. The present work demonstrates numerically that this system, also like the classical Lorenz system, possesses a homoclinic trajectory and a chaotic self-excited attractor. However, for consid- ered system, unlike the classical Lorenz one, along with self-excited attractor a hidden attractor can be localized. Analytical-numerical localization of hidden attractor is presented. I. INTRODUCTION and λ,µ,β are the coefficients of viscosity, heat conduc- tion, and volume expansion, respectively; q1(t), q2(t), Consider the following physical problem: the convec- and q3(t) are the projections of temperature gradients tion of viscous incompressible fluid flow inside the ellip- on the axes x1, x2 and x3, respectively, in which case q3(t) 0; ω1(t), ω2(t), and ω3(t) are the projections of soid ≡ the vectors of fluid angular velocities on the axis x1, x2, 2 2 2 x1 x2 x3 and x3, in which case + + =1, a1 >a2 >a3 > 0, a1 a2 a3 gβa3 gβa3 ω1 = cos αq1, ω2 = cos αq2. under the condition of stationary inhomogeneous exter- −2a1a2Ω0 −2a1a2Ω0 nal heating. It is assumed that the ellipsoid together The parameters σ, Ta, and Ra are the Prandtl, Taylor, with heat sources rotates with the constant velocity Ω0 and Rayleigh numbers, respectively. around its axis. The axis has a constant angle α with the gravity vector g. This vector is stationary with re- x3 spect to the ellipsoid motion. The value Ω0 is assumed to be such that the centrifugal forces can be neglected Ω0 in comparison with the influence of gravitational field. Consider the case when the ellipsoid rotates around the x2 g axis x3 and the vector is placed in the plane x1x3. The α temperature difference is generated along the axis a1 and ∇Tˆ a constant q0 is a gradient of this temperature. (Fig. 1). Corresponding mathematical model (three-mode g model of convection) was obtained by Glukhovsky and Dolzhansky [1] in the following form x1 x˙ = Ayz + Cz σx, Figure 1: Illustration of the problem setting − y˙ = xz + Ra y, (1) − − After two sequential transformations z˙ = xy z. arXiv:1412.7667v1 [nlin.CD] 24 Dec 2014 − −1 −1 x x, y C y, z C z, Here → → σ → σ x x, y R z, z y, 2 a0R +1 a0R +1 λ Ω0 gβa3qˆ1 → → − → σ = , Ta = 2 , Ra = , µ λ 2a1a2λµ one obtains the following system 2 2 2 a1 a2 2 −1 2a1a2 A = 2 − 2 cos αTa , C = 2 2 σ sin α, x˙ = σ(x y) ayz a + a a3(a + a ) − − − 1 2 1 2 y˙ = rx y xz (2) −1 gβa3 cos α gβa3 − − x = µ ω3 + q3 , y = q1, z˙ = z + xy, 2a1a2Ω0 2a1a2λµ − 2 gβa3 where a0 = A/C , R = Ra C, z = q2, 2a1a2λµ 2 a0σ R a = 2 , r = (a0R + 1). (3) (a0R + 1) σ In the case a =0 system (2) coincides with the classi- ∗ Corresponding author: [email protected] cal Lorenz system [2]. For the first time, system (2) with 2 parameters r,σ > 0 was considered in [3, 1978]. After II. HOMOCLINIC ORBIT IN THE the linear change of variables [4] this system can be re- LORENZ-LIKE SYSTEM duced to the Rabinovich system, describing the waves interaction in plasma [5–8]. As is shown in [4] system Denote by x(t) a trajectory of system (2) starting at a (2) describes the following physical processes: the flow certain initial point. In order to compute numerically a of fluid convection inside the rotating ellipsoid [1], the homoclinic trajectory (limt→+∞ x(t) = limt→−∞ x(t) = rotation of rigid body in viscous fluid [9], the gyrostat S0), one integrates system (2) with the initial data x0 dynamics [10, 11], the convection of a horizontal layer from a δ-vicinity of the saddle point S0 and its one- u of fluid making the harmonic oscillations [12], and the dimensional unstable manifold W (S0) that corresponds model of Kolmogorov’s flow [13]. In [14] for system (2) to a positive eigenvalue of the Jacobian matrix J at in the case σ = ar a detailed analysis of the equilibria ± the saddle point S0. For some parameters of system stability and asymptotic behavior of trajectories is given (2) this trajectory after a certain time intersects a two- and the values of parameters are obtained for which sys- dimensional plane M spanned on the eigenvectors that tem (2) is integrable. Remark also the works [15, 16], in correspond to negative eigenvalues of J. The parameters which the analytical and numerical study of some gener- are chosen in such a way that the point of intersection alizations of system (2) and similar systems is presented. belongs to δ-vicinity of S0. In addition, in [17] system (2) was used to describe a specific scenario of transition to chaos in low-dimensional dynamical systems gluing bifurcations. Note that the Glukhovsky-Dolzhansky system is suffi- ciently different from the classical Lorenz system. In the 12 Lorenz system, the flow of the two-dimensional convec- 10 tion is considered only. In the Glukhovsky-Dolzhansky S S 1 system, the flow of the three-dimensional convection is 8 2 considered which can be interpreted as one of the models Z 6 of ocean flows [1]. 4 10 In what follows system (2) will be considered under the S 5 2 0 condition that the parameter a is positive. In this case 0 0 -5 if r< 1, then (2) has a unique equilibrium S0 = (0, 0, 0), -6 -4 -2 0 2 -10 which is globally asymptotically Lyapunov stable [4, 18]. 4 6 Y If r> 1, then (2) posesses three equilibria: S0 = (0, 0, 0) X and Figure 2: Approximation of the homoclinic butterfly for system (2). S1,2 = ( x1, y1, z1). (4) ± ± Let us fix the parameters: σ =4 and a =0.0052 (such Here values were considered in [1]). For r = 7.44 there is no intersection of the trajectory x(t) with the plane M σ√ξ σξ (see Fig. 3a) and for r = 7.445 the intersection occurs x1 = , y1 = ξ, z1 = , σ + aξ σ + aξ (see Fig. 3c). So, there exists an intermediate value p r∗ [7.44, 7.445] for which one can get the approxima- ∈ ∗ and the number ξ is defined as tion of homoclinic orbit r =7.4430045820796753... (see Fig. 3b). Note that the approximation for the symmet- ric homoclinic orbit can be obtained by the choose in the σ ξ = a(r 2) σ + (σ ar)2 +4aσ . computational procedure the symmetric (with respect to 2a2 − − − h p i S0) initial data (see Fig. 2). From an analytical point of view, the existence of homoclinic trajectory can be justi- The stability of equilibria S1,2 depends on the parameters fied by Fishing principle [20–22]. The Fishing principle σ, r, a (see. Sec. IIIA). is based on the construction of a special two-dimensional For system (2) with the fixed σ, a (or with the fixed manifold such that the separatrix of the saddle point in- σ only) it is possible to observe a classical scenario of tersects or does not intersect the manifold for two differ- transition to the chaos similar to scenario in the Lorenz ent values of a system parameter. The continuity implies system [19]. To demonstrate this, for system (2) with the existence of some intermediate value of the parameter the fixed parameters σ and a and increasing parameter for which the separatrix touches the manifold. Accord- r > 1, a homoclinic trajectory and a self-excited chaotic ing to construction the only possibility for separatrix is to attractor are obtained numerically. Unlike the Lorenz touch the saddle and thus, one can numerically localize system, for system (2) it is also possible to localize a the birth of the homolcinic orbit by changing the variable hidden chaotic attractor. parameter. 3 12 M 12 M 12 M S S 10 1 10 1 10 S S S 2 2 S 1 8 8 8 2 6 6 6 Z Z Z 4 4 S S 4 0 0 S 2 5 2 5 2 0 10 5 0 0 0 0 0 0 -2 -5 -2 -5 -2 -5 -6 -5 -6 -5 -6 -4 -3 -4 -3 -4 -2 -1 -2 -1 -2 0 0 1 -10 0 1 -10 2 -10 2 3 Y 2 3 Y 4 6 Y X X X (a) r = 7.44 (b) r = 7.4430045820796753... (c) r = 7.445 Figure 3: The birth of a homoclinic orbit in system (2) with σ =4, a =0.0052, r [7.44, 7.445], δ =0.01. ∈ xxxxxxxxxxxxxxxx r xxxxxxxxxxxxxxxx III. CHAOTIC ATTRACTOR IN THE xxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxx 1000 xxxxxxxxxxxxxxxx LORENZ-LIKE SYSTEM xxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxx r xxxxxxxxxxxxxxxx 2 xxxxxxxxxxxxxxxx A. Local stability analysis and computation of xxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxr xxxxxxxxxxxxxxxx1 attractors xxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxx 500 xxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxx Let us study the stability of equilibria S1, S2 of system xxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxx rxxxxxxxxxxxxxxxx (2). By the Routh-Hurwitz criterion, one can get the xxxxxxxxxxxxxxxx3 xxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx following xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx Proposition 1 .
Recommended publications
  • Infinitely Many Coexisting Attractors in No-Equilibrium Chaotic System
    Hindawi Complexity Volume 2020, Article ID 8175639, 17 pages https://doi.org/10.1155/2020/8175639 Research Article Infinitely Many Coexisting Attractors in No-Equilibrium Chaotic System Qiang Lai ,1 Paul Didier Kamdem Kuate ,2 Huiqin Pei ,1 and Hilaire Fotsin 2 1School of Electrical and Automation Engineering, East China Jiaotong University, Nanchang 330013, China 2Laboratory of Condensed Matter, Electronics and Signal Processing Department of Physics, University of Dschang, P.O. Box 067, Dschang, Cameroon Correspondence should be addressed to Qiang Lai; [email protected] Received 17 January 2020; Revised 18 February 2020; Accepted 26 February 2020; Published 28 March 2020 Guest Editor: Chun-Lai Li Copyright © 2020 Qiang Lai et al. *is is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. *is paper proposes a new no-equilibrium chaotic system that has the ability to yield infinitely many coexisting hidden attractors. Dynamic behaviors of the system with respect to the parameters and initial conditions are numerically studied. It shows that the system has chaotic, quasiperiodic, and periodic motions for different parameters and coexists with a large number of hidden attractors for different initial conditions. *e circuit and microcontroller implementations of the system are given for illustrating its physical meaning. Also, the synchronization conditions of the system are established based on the adaptive control method. 1. Introduction system may generate hidden attractor. Also, a lot of pre- vious studies have shown that chaotic systems with mul- Encouraging progress has been made on chaos in the past tiple unstable equilibria usually have richer dynamic few decades.
    [Show full text]
  • The Lorenz System : Hidden Boundary of Practical Stability and the Lyapunov Dimension
    This is a self-archived version of an original article. This version may differ from the original in pagination and typographic details. Author(s): Kuznetsov, N. V.; Mokaev, T. N.; Kuznetsova, O. A.; Kudryashova, E. V. Title: The Lorenz system : hidden boundary of practical stability and the Lyapunov dimension Year: 2020 Version: Published version Copyright: © 2020 the Authors Rights: CC BY 4.0 Rights url: https://creativecommons.org/licenses/by/4.0/ Please cite the original version: Kuznetsov, N. V., Mokaev, T. N., Kuznetsova, O. A., & Kudryashova, E. V. (2020). The Lorenz system : hidden boundary of practical stability and the Lyapunov dimension. Nonlinear Dynamics, 102(2), 713-732. https://doi.org/10.1007/s11071-020-05856-4 Nonlinear Dyn https://doi.org/10.1007/s11071-020-05856-4 ORIGINAL PAPER The Lorenz system: hidden boundary of practical stability and the Lyapunov dimension N. V. Kuznetsov · T. N. Mokaev · O. A. Kuznetsova · E. V. Kudryashova Received: 16 March 2020 / Accepted: 29 July 2020 © The Author(s) 2020 Abstract On the example of the famous Lorenz Keywords Global stability · Chaos · Hidden attractor · system, the difficulties and opportunities of reliable Transient set · Lyapunov exponents · Lyapunov numerical analysis of chaotic dynamical systems are dimension · Unstable periodic orbit · Time-delayed discussed in this article. For the Lorenz system, the feedback control boundaries of global stability are estimated and the difficulties of numerically studying the birth of self- excited and hidden attractors, caused by the loss of 1 Introduction global stability, are discussed. The problem of reliable numerical computation of the finite-time Lyapunov In 1963, meteorologist Edward Lorenz suggested an dimension along the trajectories over large time inter- approximate mathematical model (the Lorenz system) vals is discussed.
    [Show full text]
  • Phase-Locked Loops with Active PI Filter: the Lock-In Range Computation JYVÄSKYLÄ STUDIES in COMPUTING 239
    JYVÄSKYLÄ STUDIES IN COMPUTING 239 Konstantin Aleksandrov Phase-Locked Loops with Active PI Filter: the Lock-In Range Computation JYVÄSKYLÄ STUDIES IN COMPUTING 239 Konstantin Aleksandrov Phase-Locked Loops with Active PI Filter: the Lock-In Range Computation Esitetään Jyväskylän yliopiston informaatioteknologian tiedekunnan suostumuksella julkisesti tarkastettavaksi yliopiston Agora-rakennuksen Alfa-salissa kesäkuun 17. päivänä 2016 kello 10. Academic dissertation to be publicly discussed, by permission of the Faculty of Information Technology of the University of Jyväskylä, in building Agora, Alfa hall, on June 17, 2016 at 10 o’clock. UNIVERSITY OF JYVÄSKYLÄ JYVÄSKYLÄ 2016 Phase-Locked Loops with Active PI Filter: the Lock-In Range Computation JYVÄSKYLÄ STUDIES IN COMPUTING 239 Konstantin Aleksandrov Phase-Locked Loops with Active PI Filter: the Lock-In Range Computation UNIVERSITY OF JYVÄSKYLÄ JYVÄSKYLÄ 2016 Editors Timo Männikkö Department of Mathematical Information Technology, University of Jyväskylä Pekka Olsbo, Ville Korkiakangas Publishing Unit, University Library of Jyväskylä URN:ISBN:978-951-39-6688-1 ISBN 978-951-39-6688-1 (PDF) ISBN 978-951-39-6687-4 (nid.) ISSN 1456-5390 Copyright © 2016, by University of Jyväskylä Jyväskylä University Printing House, Jyväskylä 2016 ABSTRACT Aleksandrov, Konstantin Phase-locked loops with active PI filter: the lock-in range computation Jyväskylä: University of Jyväskylä, 2016, 38 p.(+included articles) (Jyväskylä Studies in Computing ISSN 1456-5390; 239) ISBN 978-951-39-6687-4 (nid.) ISBN 978-951-39-6688-1 (PDF) Finnish summary Diss. The present work is devoted to the study of the lock-in range of phase-locked loop (PLL). The PLL concept was originally described by H.
    [Show full text]
  • Hidden Bifurcation in the Multispiral Chua Attractor Menacer Tidjani, René Lozi, Leon Chua
    Hidden bifurcation in the Multispiral Chua Attractor Menacer Tidjani, René Lozi, Leon Chua To cite this version: Menacer Tidjani, René Lozi, Leon Chua. Hidden bifurcation in the Multispiral Chua Attractor. International journal of bifurcation and chaos in applied sciences and engineering , World Scientific Publishing, 2016, 26 (14), pp.1630039-1 1630023-26. 10.1142/S0218127416300391. hal-01488550 HAL Id: hal-01488550 https://hal.archives-ouvertes.fr/hal-01488550 Submitted on 13 Mar 2017 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. January 6, 2017 11:8 WSPC/S0218-1274 1630039 Published in International Journal of Bifurcation and Chaos, Personal file Vol. 26, No. 14, (2016), 1630039-1 to 1630039-26 DOI 10.1142/S0218127416300391 Hidden Bifurcations in the Multispiral Chua Attractor Tidjani Menacer Department of Mathematics, University Mohamed Khider, Biskra, Algeria [email protected] Ren´eLozi Universit´eCˆote d’Azur, CNRS, LJAD, Parc Valrose 06108, Nice C´edex 02, France [email protected] Leon O. Chua Department of Electrical Engineering and Computer Sciences, University of California, Berkeley, CA 94720, USA [email protected] Received August 20, 2016 We introduce a novel method revealing hidden bifurcations in the multispiral Chua attractor in the case where the parameter of bifurcation c which determines the number of spiral is discrete.
    [Show full text]
  • This Is an Electronic Reprint of the Original Article. This Reprint May Differ from the Original in Pagination and Typographic Detail
    This is an electronic reprint of the original article. This reprint may differ from the original in pagination and typographic detail. Author(s): Kuznetsov, Nikolay; Leonov, G. A.; Yuldashev, M. V.; Yuldashev, R. V. Title: Hidden attractors in dynamical models of phase-locked loop circuits : limitations of simulation in MATLAB and SPICE Year: 2017 Version: Please cite the original version: Kuznetsov, N., Leonov, G. A., Yuldashev, M. V., & Yuldashev, R. V. (2017). Hidden attractors in dynamical models of phase-locked loop circuits : limitations of simulation in MATLAB and SPICE. Communications in Nonlinear Science and Numerical Simulation, 51, 39-49. doi:10.1016/j.cnsns.2017.03.010 All material supplied via JYX is protected by copyright and other intellectual property rights, and duplication or sale of all or part of any of the repository collections is not permitted, except that material may be duplicated by you for your research use or educational purposes in electronic or print form. You must obtain permission for any other use. Electronic or print copies may not be offered, whether for sale or otherwise to anyone who is not an authorised user. Accepted Manuscript Hidden attractors in dynamical models of phase-locked loop circuits: limitations of simulation in MATLAB and SPICE Kuznetsov N. V., Leonov G. A., Yuldashev M. V., Yuldashev R. V. PII: S1007-5704(17)30090-4 DOI: 10.1016/j.cnsns.2017.03.010 Reference: CNSNS 4137 To appear in: Communications in Nonlinear Science and Numerical Simulation Received date: 30 August 2016 Revised date: 4 November 2016 Accepted date: 13 March 2017 Please cite this article as: Kuznetsov N.
    [Show full text]
  • A New 4-D Chaotic System with Hidden Attractor and Its Circuit Implementation
    International Journal of Engineering & Technology, 7 (3) (2018) 1245-1250 International Journal of Engineering & Technology Website: www.sciencepubco.com/index.php/IJET doi: 10.14419/ijet.v7i3.9846 Research paper A New 4-D Chaotic System with Hidden Attractor and its Circuit Implementation Aceng Sambas1*, Mustafa Mamat2, Sundarapandian Vaidyanathan3, Mohamad Afendee Mohamed2 and W. S. Mada Sanjaya4 1Department of Mechanical Engineering, Universitas Muhammadiyah Tasikmalaya, Indonesia 2Faculty of Informatics and Computing, Universiti Sultan Zainal Abidin, Malaysia 3Research and Development Centre, Vel Tech University, Avadi, Chennai, India 4Department of Physics, Universitas Islam Negeri Sunan Gunung Djati Bandung, Indonesia *[email protected] Abstract In the chaos literature, there is currently significant interest in the discovery of new chaotic systems with hidden chaotic attractors. A new 4-D chaotic system with only two quadratic nonlinearities is investigated in this work. First, we derive a no-equilibrium chaotic system and show that the new chaotic system exhibits hidden attractor. Properties of the new chaotic system are analyzed by means of phase portraits, Lyapunov chaos exponents, and Kaplan-Yorke dimension. Then an electronic circuit realization is shown to validate the chaotic behavior of the new 4-D chaotic system. Finally, the physical circuit experimental results of the 4-D chaotic system show agreement with numerical simulations. Keywords: Chaos, chaotic systems,circuit simulation, hidden attractors, Lyapunov exponents 1. Introduction chaos exponents of a chaotic system are determined using Wolf’s algorithm [53]. In the literature, there is also good interest shown Chaos theory deals with nonlinear dynamical systems that are highly in building electronic circuit designs of chaotic systems and imple- sensitive to initial conditions.
    [Show full text]
  • Looking More Closely to the Rabinovich-Fabrikant System
    Looking more closely to the Rabinovich-Fabrikant system MARIUS-F. DANCA Dept. of Mathematics and Computer Science, Avram Iancu University of Cluj-Napoca, Romania and Romanian Institute of Science and Technology, Cluj-Napoca, Romania [email protected] MICHAL FECKAN˘ Department of Mathematical Analysis and Numerical Mathematics Comenius University in Bratislava, Slovakia and Mathematical Institute Slovak Academy of Sciences Bratislava, Slovakia [email protected] NIKOLAY KUZNETSOV Department of Applied Cybernetics Saint-Petersburg State University, Russia and University of Jyv¨askyl¨a,Finland [email protected] GUANRONG CHEN Department of Electronic Engineering, City University of Hong Kong, Hong Kong SAR, China [email protected] October 2, 2015 Abstract Recently, we look more closely into the Rabinovich-Fabrikant system, after a decade of the study in [Danca & Chen(2004)], discovering some new characteristics such as cycling chaos, transient chaos, chaotic hidden attractors and a new kind of saddles-like attractor. In addition to extensive and accurate numerical analysis, on the assumptive existence of heteroclinic orbits, we provide a few of their approximations. arXiv:1509.09206v2 [nlin.CD] 1 Oct 2015 Rabinovich-Fabrikant system; cycling chaos; transient chaos; heteroclinic orbit; LIL numerical method 1 Introduction Rabinovich & Fabrikant [1979] introduced and analyzed from physical point of view a model describing the stochasticity arising from the modulation instability in a non-equilibrium dissipative medium. This is a simplification of a complex nonlinear parabolic equation modelling different physical systems, such as the Tollmien-Schlichting waves in hydrodynamic flows, wind waves on water, concentration waves during chemical reactions in a medium where diffusion occur, Langmuir waves in plasma, etc.
    [Show full text]
  • Hidden Attractors on One Path: Glukhovsky-Dolzhansky, Lorenz, and Rabinovich Systems
    Hidden attractors on one path: Glukhovsky-Dolzhansky, Lorenz, and Rabinovich systems G. Chen,1 N.V. Kuznetsov,2, 3 G.A. Leonov,2, 4 and T.N. Mokaev2 1City University of Hong Kong, Hong Kong SAR, China 2Faculty of Mathematics and Mechanics, St. Petersburg State University, Peterhof, St. Petersburg, Russia 3Department of Mathematical Information Technology, University of Jyv¨askyl¨a,Jyv¨askyl¨a,Finland 4Institute of Problems of Mechanical Engineering RAS, Russia (Dated: October 2, 2018) In this report, by the numerical continuation method we visualize and connect hidden chaotic sets in the Glukhovsky-Dolzhansky, Lorenz and Rabinovich systems using a certain path in the parameter space of a Lorenz-like system. I. INTRODUCTION The Glukhovsky-Dolzhansky system describes the con- vective fluid motion inside a rotating ellipsoidal cavity. In 1963, meteorologist Edward Lorenz suggested an ap- If we set proximate mathematical model (the Lorenz system) for a < 0; σ = ar; (7) the Rayleigh-B´enardconvection and discovered numeri- − cally a chaotic attractor in this model [1]. This discovery then after the linear transformation (see, e.g., [8]): stimulated rapid development of the chaos theory, nu- 1 1 1 1 1 merical methods for attractor investigation, and till now (x; y; z) ν1− y; ν1− ν2− h x; ν1− ν2− h z) ; t ν1 t has received a great deal of attention from different fields ! ! [2{7]. The Lorenz system gave rise to various generaliza- with positive ν1; ν2; h, we obtain the Rabinovich system tions, e.g. Lorenz-like systems, some of which are also [10, 11], describing the interaction of three resonantly simplified mathematical models of physical phenomena.
    [Show full text]
  • Unusual Dynamics and Hidden Attractors of the Rabinovich-Fabrikant System
    Nonlinear Dynamics manuscript No. (will be inserted by the editor) Unusual dynamics and hidden attractors of the Rabinovich-Fabrikant system Marius-F. Danca · Nikolay Kuznetsov · Guanrong Chen the date of receipt and acceptance should be inserted later Abstract This paper presents some numerical evidence of unusual dynamics and hidden chaotic attractors of the Rabinovich-Fabrikant system, with some insightful descriptions and discussions. From a generalized Hamiltonian energy perspective, the attractors could be analyzed in more details. Keywords Rabinovich-Fabrikant system · Hidden chaotic attractor · Self-excited attractor · Hamiltonian energy 1 Introduction The Rabinovich-Fabrikant (RF) system, introduced by Rabinovich and Fabrikant [1] and numerically investigated in [2], was initially established as a physical model describing the stochasticity arising from the modulation instability in a non- equilibrium dissipative medium. However, recently in [2] we revealed that besides the physical properties described in [1], the system presents extremely rich com- plex dynamics. More significantly, unlike other familiar chaotic systems containing only second-order nonlinearities (such as the Lorenz system), the RF system with third-order nonlinearities presents some unusual dynamics, namely, some complex dynamics that were rarely or never seen from other systems in the past, such as \virtual" saddles in addition to several chaotic attractors with different shapes as well as hidden chaotic attractors and hidden transient chaotic attractors [3]. In fact, common interest in this system is continuously increasing (see e.g. [4{11]). Marius-F. Danca Dept. of Mathematics and Computer Science, Avram Iancu University of Cluj-Napoca, Ro- mania and Romanian Institute of Science and Technology, Cluj-Napoca, Romania Nikolay Kuznetsov Dept.
    [Show full text]
  • Hidden Chaotic Attractors in a Class of Two-Dimensional Maps
    ORE Open Research Exeter TITLE Hidden chaotic attractors in a class of two-dimensional maps AUTHORS Jiang, H; Liu, Y; Wei, Z; et al. JOURNAL Nonlinear Dynamics DEPOSITED IN ORE 11 October 2016 This version available at http://hdl.handle.net/10871/23849 COPYRIGHT AND REUSE Open Research Exeter makes this work available in accordance with publisher policies. A NOTE ON VERSIONS The version presented here may differ from the published version. If citing, you are advised to consult the published version for pagination, volume/issue and date of publication Hidden chaotic attractors in a class of two-dimensional maps∗ Haibo Jiangay, Yang Liub, Zhouchao Weic, Liping Zhanga a School of Mathematics and Statistics, Yancheng Teachers University, Yancheng 224002, China bSchool of Engineering, Robert Gordon University, Garthdee Road, Aberdeen AB10 7GJ, UK cSchool of Mathematics and Physics, China University of Geosciences, Wuhan 430074, China Abstract This paper studies the hidden dynamics of a class of two-dimensional maps inspired by the H´enonmap. A special consideration is made to the existence of fixed points and their stabilities in these maps. Our concern focuses on three typical scenarios which may generate hidden dynamics, i.e. no fixed point, single fixed point, and two fixed points. A computer search program is employed to explore the strange hidden attractors in the map. Our findings show that the basins of some hidden attractors are tiny, so the standard computational procedure for localization is unavailable. The schematic exploring method proposed in this paper could be generalized for investigating hidden dynamics of high-dimensional maps. Keywords: hidden dynamics; two-dimensional maps; fixed point; stability; coexistence.
    [Show full text]
  • Hidden Periodic and Chaotic Oscillations in Nonlinear Dynamical
    Preprints of the 19th World Congress The International Federation of Automatic Control Cape Town, South Africa. August 24-29, 2014 Hidden attractors in dynamical systems: systems with no equilibria, multistability and coexisting attractors N.V. Kuznetsov, G.A. Leonov Department of Applied Cybernetics, Saint-Petersburg State University, Russia Department of Mathematical Information Technology, University of Jyv¨askyl¨a,Finland (e-mail: [email protected]) Abstract: From a computational point of view it is natural to suggest the classification of attractors, based on the simplicity of finding basin of attraction in the phase space: an attractor is called a hidden attractor if its basin of attraction does not intersect with small neighborhoods of equilibria, otherwise it is called a self-excited attractor. Self-excited attractors can be localized numerically by the standard computational procedure, in which after a transient process a trajectory, started from a point of unstable manifold in a neighborhood of unstable equilibrium, is attracted to the state of oscillation and traces it. Thus one can determine unstable equilibria and check the existence of self-excited attractors. In contrast, for the numerical localization of hidden attractors it is necessary to develop special analytical-numerical procedures in which an initial point can be chosen from the basin of attraction analytically. For example, hidden attractors are attractors in the systems with no-equilibria or with the only stable equilibrium (a special case of multistability and coexistence of attractors); hidden attractors arise in the study of well-known fundamental problems such as 16th Hilbert problem, Aizerman and Kalman conjectures, and in applied research of Chua circuits, phase-locked loop based circuits, aircraft control systems, and others.
    [Show full text]
  • Hidden Oscillations in Dynamical Systems
    Recent Researches in System Science Hidden oscillations in dynamical systems G.A. LEONOV a, N.V. KUZNETSOV b,a, S.M. SELEDZHI a aSt.Petersburg State University, Universitetsky pr. 28, St.Petersburg, 198504, RUSSIA bUniversity of Jyv¨askyl¨a, P.O. Box 35 (Agora), FIN-40014, FINLAND [email protected], [email protected], [email protected] Abstract :- The classical attractors of Lorenz, Rossler, Chua, Chen, and other widely-known attractors are those excited from unstable equilibria. From computational point of view this allows one to use numerical method, in which after transient process a trajectory, started from a point of unstable manifold in the neighborhood of equilibrium, reaches an attractor and identifies it. However there are attractors of another type: hidden attractors, a basin of attraction of which does not contain neighborhoods of equilibria. Study of hidden oscillations and attractors requires the development of new analytical and numerical methods which will be considered in this invited lecture. Key-Words:- Hidden oscillation, attractor localization, hidden attractor, harmonic balance, describing func- tion method, Aizerman conjecture, Kalaman conjecture, Hilbert 16th problem 1 Introduction 5 In establishing and developing the theory of non- 4 linear oscillations in the first half of the twentieth cen- tury [1, 2, 3, 4] most attention was given to analyzing 3 and synthesizing oscillatory systems for which solving 2 the problem of the existence of the oscillation modes did not present any great difficulties. The structure 1 of many mechanical, electromechanical and electronic y 0 systems was such that there were oscillation modes in them, the existence of which was almost obvious −1 — oscillations are excited from unstable equilibria.
    [Show full text]