DOCSLIB.ORG

Localized Lyapunov Exponents and the Prediction of Predictability

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Localized Lyapunov Exponents and the Prediction of Predictability ¡ ¢ Physics Letters A 9909 2000 xxx ¥ ¤ www.elsevier.nl£ locate pla Localized Lyapunov exponents and the prediction of predictability ¦ Christine Ziehmann a, § , Leonard A. Smith a,b,c, JurgenÈ Kurths a ¨ a Uni© ersitat Potsdam, Institut fur Physik, Nichtlineare Dynamik, Postfach 60 15 53, D-14469, Potsdam, Germany ÈÈ b Mathematical Institute,Uni ersity of Oxford, Oxford, OX1 3LB, UK c London School of Economics, London WC2A 2AE, UK Received 9 December 1999; received in revised form 19 April 2000; accepted 3 May 2000 Communicated by C.R. Doering Abstract Every forecast should include an estimate of its likely accuracy, a current measure of predictability. Two distinct types of localized Lyapunov exponents based on infinitesimal uncertainty dynamics are investigated to reflect this predictability. Regions of high predictability within which any initial uncertainty will decrease are proven to exist in two common chaotic systems; potential implications of these regions are considered. The relevance of these results for finite size uncertainties is discussed and illustrated numerically. 2000 Elsevier Science B.V. All rights reserved. PACS: 05.45.-a; 05.10.-a; 05.45.Tp Keywords: Lyapunov exponents; Predictability; Chaos; Weather forecasting; Ensemble prediction; Nonlinear dynamics 1. Introduction ! difficulty in quantifying predictability: 1 the depen- ! dence of measures of predictability upon the particu- " # $ lar metric adopted 2 the dependence of uncertainty Prediction of predictability refers to the quantita- ! dynamics upon the magnitude of the uncertainty in % & tive attempt to assess the likely error in a particular the initial condition, and 3 the fact that errors in the forecast a priori. There are at least three sources of model are often unknown until after predictions are ' ( ) observed to fail. Lyapunov exponents 1±3 quantify * predictability through globally averaged effective + growth rates of uncertainty in the limits of large time , and small uncertainty; thus by construction they are ) of limited use. To obtain a quantitative estimate of Corresponding author. Tel.: 49-331-977-1302, fax: 49- the accuracy of a particular forecast, the local dy- 331-977-1142. E-mail address: [email protected] C. Ziehmann . - namics of uncertainties about that initial condition / 0 0375-9601. 00 $ - see front matter 2000 Elsevier Science B.V. All rights reserved. 2 PII: S0375-9601 1 00 00336-4 2 C. Ziehmann et al. 3 Physics Letters A 00 ()2000 000±000 4 5 , are more relevant 4±9 . By allowing better risk not considered in this paper as there is no systematic , D D assessment, a prediction of predictability is of value manner to include system H model mismatches, thus it ? in 6 any field from physics to economics; weather is assumed throughout the paper that the perfect forecasting provides a particular example relevant to model is known. 7 both fields. + Effecti 8 e growth rates defined over a fixed dura- tion are also used to quantify predictability; they are I 9 employed daily in the operational weather forecast 2. Localized Lyapunov exponents : ; < centers of Europe and North America 10,11 . In Section 2, the distinction between what will be called The dynamics of infinitesimal uncertainties about , J finite time exponents and finite sample exponents is K L M N = a point x0 in an m-dimensional state space are > J shown to lie in the particular initial orientation of + L O P Q > 9 governed by the linear propagator, x0 , t , which U T the perturbation each considers for a given initial 9 m L R S < evolves any infinitesimal initial uncertainty 0 R > J , condition; this can result in dramatically different , L V + 8 8 about x0 forward for a time t along the system's effecti e growth rates. Both types of exponent are J W < trajectory to x : called `Local Lyapunov exponents', and recognizing t ^ _ ` a > ] J X Y Z [ b the distinction between them resolves some confu- L t x, \ t 0 .1 sion in the literature. In terms of predicting the d forecast accuracy, the finite time exponents are shown c In discrete time maps, the linear propagator over k to be the more relevant quantities in Section 3, where ? iterations is simply the product of Jacobians along e f g h d k l m J J n o L j it is also proven that the mean of the largest finite i the trajectory, that is x0 ,k xk 1 ... p q r s t u v w J time exponent does not provide an unbiased estimate J x10x . x ) of the corresponding global exponent, and similarly For high dimensional systems, interest tends to be for the mean of the smallest finite time exponent. In focused on subspaces which are likely to contain the y z , addition, in regions of state space where the largest fastest growing perturbations 4,10,11 . Two orienta- { | * 6 finite time exponent is less than zero all perturba- tions of particular interest are i that which will tions will shrink independent of their orientation; this } have grown the most under the linearized dynamics ? ~ n d J is investigated in Section 4 where such regions are , after k steps, k x , and ii the local orientation of 1 * proven to exist in two common chaotic maps. the globally fastest growing direction, lx, which 1 Each class of Lyapunov exponent discussed in is sometimes called the Lyapunov vector 13 . The this paper assumes the observational uncertainty is first of these orientations is defined by the singular ? infinitesimal; of course as long as it remains in- B value decomposition 14 of the propagator: the n , finitesimal it cannot limit predictability, and once it k J x are simply the right singular vectors of i d is finite its growth is no longer quantified by Lya- J x,k . Each is associated with a singular value, n * J k punov exponents. Therefore, the rigorous results re- x ; by convention, . The finite time iii 1 stricted to infinitesimal uncertainties are contrasted Lyapuno 8 exponents 15 are @ with numerical demonstrations for finite uncertain- ° ± n ¡ ¢ n £ ¤ ¤ ¥ ¥ ª « ¬ ­ ® ¯ $ ¨ © d ties in Section 4. Part of the popularity of global ¦ 1 § J J k J k Lyapunov exponents stems from the fact that their i x d log 2 x,k i x A k B value does not depend upon the metric or coordinate C ² º » n ³ · ¸ ¹ system used; this is not the case for the exponents 1 ¶ J ´ 7 k ¼ µ d log x .2 based upon a finite length of trajectory, yet in prac- k 2 i tice only the latter are available. We return to this issue in Section 5. Finally, there is the question of ½ Properties of these exponents have been noted by ¾ ¿ À Á Â Ã Ä D model error in nonlinear forecasting, either paramet- Lorenz 4 , Grassberger et al. 16 , Abarbanel 5 and E F Å Æ G ric or structural 12 . Arguably, model error may be references thereof. By Oseledec's Theorem 1 , in Ç n È É Ê d Ë Ì Í < k J more responsible for poor predictions of real nonlin- the limit k the i x converge to a unique set 9 ) ear systems than `chaos.' Model imperfections are of values, the Lyapunov exponents Î i, which are the C. Ziehmann et al. Ï Physics Letters A 00 ()2000 000±000 3 n c d J @ a b , d , k same for almost all x with respect to an ergodic and 11approach , yet for the relatively small k measure. ) over which forecasts are typically made their proper- ties are quite different and neither is constrained by Ú Û Ð Ñ Ò Ó n Ô Õ Ö × Ø $ 1 k Ù e lim d log , the value of . Ý Þ ii Ü 2 i 1 k k à á K i ß 1,2,...,m 3 f 3. Properties of localized Lyapunov exponents If the sum of the â i is negative, a volume 9 element in state space will shrink, on average, as it g 9 evolves along a trajectory, and most initial condi- General constraints on the relative magnitudes of tions will evolve towards an attractor of dimension the largest and the smallest finite time and finite ã ä J ) less than K m 17 . At each point x on such an sample exponents are now derived from the defini- å æ ! , attractor, the orientation lxdenotes the orienta- tions above, and then illustrated below. By construc- 1 tion, maximum growth corresponds to , thus ç n i j k l n m n o tion corresponding to 1, that is the orientation h 1 J q r J J p ? k x k x for each x, and therefore this in- s t è u v towards which almost every uncertainty in the 11 n 9 k x equality also holds for the mean values w { | } ~ n 1 z sufficiently distant past would have evolved, when y ! é ê k J , , where denotes an arithmetic average U the trajectory reaches x. Similarly, define lxas 1 ¼ m ë 2 taken with respect to the natural measure , and U the orientation corresponding to ì for details, see N m í î ï ð ñ ò ó , , < L a numerical approximation with sample size N.Ex- 17 . Numerically, lx10and lxU m 0 can be approx- , amples from two chaotic systems are given in Fig. 1. imated by evolving the singular vectors of ô õ ö ÷ ø ù ù J ú û ý ü The Henon map 24 is x j ,2 j , that is, the 2 j step propagator about ù ù ¼ * J th þ L 6 the j pre-image of x , forward j steps until the 2 0 x 1 ax y , y bx ,6 i 1 iii1 i J L trajectory reaches x .
Load more

Recommended publications
  • What Are Lyapunov Exponents, and Why Are They Interesting?
    BULLETIN (New Series) OF THE AMERICAN MATHEMATICAL SOCIETY Volume 54, Number 1, January 2017, Pages 79–105 http://dx.doi.org/10.1090/bull/1552 Article electronically published on September 6, 2016 WHAT ARE LYAPUNOV EXPONENTS, AND WHY ARE THEY INTERESTING? AMIE WILKINSON Introduction At the 2014 International Congress of Mathematicians in Seoul, South Korea, Franco-Brazilian mathematician Artur Avila was awarded the Fields Medal for “his profound contributions to dynamical systems theory, which have changed the face of the field, using the powerful idea of renormalization as a unifying principle.”1 Although it is not explicitly mentioned in this citation, there is a second unify- ing concept in Avila’s work that is closely tied with renormalization: Lyapunov (or characteristic) exponents. Lyapunov exponents play a key role in three areas of Avila’s research: smooth ergodic theory, billiards and translation surfaces, and the spectral theory of 1-dimensional Schr¨odinger operators. Here we take the op- portunity to explore these areas and reveal some underlying themes connecting exponents, chaotic dynamics and renormalization. But first, what are Lyapunov exponents? Let’s begin by viewing them in one of their natural habitats: the iterated barycentric subdivision of a triangle. When the midpoint of each side of a triangle is connected to its opposite vertex by a line segment, the three resulting segments meet in a point in the interior of the triangle. The barycentric subdivision of a triangle is the collection of 6 smaller triangles determined by these segments and the edges of the original triangle: Figure 1. Barycentric subdivision. Received by the editors August 2, 2016.
    [Show full text]
  • Do Chaos-Based Communication Systems Really Transmit Chaotic Signals?
    Signal Processing 108 (2015) 412–420 Contents lists available at ScienceDirect Signal Processing journal homepage: www.elsevier.com/locate/sigpro Do chaos-based communication systems really transmit chaotic signals? Renato Candido a, Diogo C. Soriano b, Magno T.M. Silva a,n, Marcio Eisencraft a a University of São Paulo, Brazil b Universidade Federal do ABC, Brazil article info abstract Article history: Many communication systems based on the synchronism of chaotic systems have been Received 2 May 2014 proposed as an alternative spread spectrum modulation that improves the level of privacy Received in revised form in data transmission. However, depending on the map and on the encoding function, the 24 September 2014 transmitted signal may cease to be chaotic. Therefore, the sensitive dependence on initial Accepted 2 October 2014 conditions, which is one of the most interesting properties for employing chaos in Available online 14 October 2014 telecommunications, may disappear. In this paper, we numerically analyze the chaotic Keywords: nature of signals modulated using a system that employs the Ikeda map. Additionally, we Chaos propose changes in the communication system in order to guarantee that the modulated Nonlinear systems signals are in fact chaotic. Chaos-based communication systems & 2014 Elsevier B.V. All rights reserved. Ikeda map Attractors 1. Introduction diverges from the original one. These three properties all together are necessary for a signal to be called chaotic Non-linear systems and chaos have been applied in all and are the basis for the alleged advantages of using chaos areas of engineering [1]. This fact is particularly true when in communications, as an improvement in security [15].
    [Show full text]
  • Synthesis of the Advance in and Application of Fractal Characteristics of Traffic Flow
    Synthesis of the Advance in and Application of Fractal Characteristics of Traffic Flow Final Report Contract No. BDK80 977‐25 July 2013 Prepared by: Lehman Center for Transportation Research Florida International University Prepared for: Research Center Florida Department of Transportation Final Report Contract No. BDK80 977-25 Synthesis of the Advance in and Application of Fractal Characteristics of Traffic Flow Prepared by: Kirolos Haleem, Ph.D., P.E., Research Associate Priyanka Alluri, Ph.D., Research Associate Albert Gan, Ph.D., Professor Lehman Center for Transportation Research Department of Civil and Environmental Engineering Florida International University 10555 West Flagler Street, EC 3680 Miami, FL 33174 Phone: (305) 348-3116 Fax: (305) 348-2802 and Hongtai Li, Graduate Research Assistant Tao Li, Ph.D., Associate Professor School of Computer Science Florida International University 11200 SW 8th Street Miami, FL 33199 Phone: (305) 348-6036 Fax: (305) 348-3549 Prepared for: Research Center State of Florida Department of Transportation 605 Suwannee Street, M.S. 30 Tallahassee, FL 32399-0450 July 2013 DISCLAIMER The opinions, findings, and conclusions expressed in this publication are those of the authors and not necessarily those of the State of Florida Department of Transportation. iii METRIC CONVERSION CHART SYMBOL WHEN YOU KNOW MULTIPLY BY TO FIND SYMBOL LENGTH in inches 25.4 millimeters mm ft feet 0.305 meters m yd yards 0.914 meters m mi miles 1.61 kilometers km mm millimeters 0.039 inches in m meters 3.28 feet ft m meters
    [Show full text]
  • Efficient Method for Detection of Periodic Orbits in Chaotic Maps And
    Efficient Method for Detection of Periodic Orbits in Chaotic Maps and Flows A Thesis by Jonathan Crofts in partial fulfillment of the requirements for the degree of Doctor of Philosophy. Department of Mathematics University of Leicester March 2007. c Jonathan J Crofts, 2007. arXiv:0706.1940v2 [nlin.CD] 13 Jun 2007 Acknowledgements I would like to thank Ruslan Davidchack, my supervisor, for his many suggestions and constant support and understanding during this research. I am also thankful to Michael Tretyakov for his support and advice. Further, I would like to acknowledge my gratitude to the people from the Department of Mathematics at the University of Leicester for their help and support. Finally, I would like to thank my family and friends for their patience and sup- port throughout the past few years. In particular, I thank my wife Lisa and my daughter Ellen, without whom I would have completed this research far quicker, but somehow, it just would not have been the same. At this point I would also like to reassure Lisa that I will get a real job soon. Leicester, Leicestershire, UK Jonathan J. Crofts 31 March 2007 i Contents Abstract iv List of figures v List of tables vi 1 Introduction 1 1.1 History,theoryandapplications . 1 1.2 Periodicorbits............................... 3 1.2.1 Periodicorbittheory . 4 1.2.2 Efficient detection of UPOs . 6 1.3 Extendedsystems............................. 9 1.4 Anoteonnumerics............................ 10 1.4.1 Intervalarithmetic . 12 1.5 Overview.................................. 14 1.6 Thesisresults ............................... 16 2 Conventional techniques for detecting periodic orbits 18 2.1 Specialcases................................ 18 2.1.1 One-dimensionalmaps .
    [Show full text]
  • Transformations)
    TRANSFORMACJE (TRANSFORMATIONS) Transformacje (Transformations) is an interdisciplinary refereed, reviewed journal, published since 1992. The journal is devoted to i.a.: civilizational and cultural transformations, information (knowledge) societies, global problematique, sustainable development, political philosophy and values, future studies. The journal's quasi-paradigm is TRANSFORMATION - as a present stage and form of development of technology, society, culture, civilization, values, mindsets etc. Impacts and potentialities of change and transition need new methodological tools, new visions and innovation for theoretical and practical capacity-building. The journal aims to promote inter-, multi- and transdisci- plinary approach, future orientation and strategic and global thinking. Transformacje (Transformations) are internationally available – since 2012 we have a licence agrement with the global database: EBSCO Publishing (Ipswich, MA, USA) We are listed by INDEX COPERNICUS since 2013 I TRANSFORMACJE(TRANSFORMATIONS) 3-4 (78-79) 2013 ISSN 1230-0292 Reviewed journal Published twice a year (double issues) in Polish and English (separate papers) Editorial Staff: Prof. Lech W. ZACHER, Center of Impact Assessment Studies and Forecasting, Kozminski University, Warsaw, Poland ([email protected]) – Editor-in-Chief Prof. Dora MARINOVA, Sustainability Policy Institute, Curtin University, Perth, Australia ([email protected]) – Deputy Editor-in-Chief Prof. Tadeusz MICZKA, Institute of Cultural and Interdisciplinary Studies, University of Silesia, Katowice, Poland ([email protected]) – Deputy Editor-in-Chief Dr Małgorzata SKÓRZEWSKA-AMBERG, School of Law, Kozminski University, Warsaw, Poland ([email protected]) – Coordinator Dr Alina BETLEJ, Institute of Sociology, John Paul II Catholic University of Lublin, Poland Dr Mirosław GEISE, Institute of Political Sciences, Kazimierz Wielki University, Bydgoszcz, Poland (also statistical editor) Prof.
    [Show full text]
  • Chaotic Scattering in Solitary Wave Interactions: a Singular Iterated-Map Description Roy H
    CHAOS 18, 023113 ͑2008͒ Chaotic scattering in solitary wave interactions: A singular iterated-map description Roy H. Goodman Department of Mathematical Sciences, New Jersey Institute of Technology, Newark, New Jersey 07102, USA ͑Received 16 October 2007; accepted 7 March 2008; published online 7 May 2008͒ We derive a family of singular iterated maps—closely related to Poincaré maps—that describe chaotic interactions between colliding solitary waves. The chaotic behavior of such solitary-wave collisions depends on the transfer of energy to a secondary mode of oscillation, often an internal mode of the pulse. This map allows us to go beyond previous analyses and to understand the interactions in the case when this mode is excited prior to the first collision. The map is derived using Melnikov integrals and matched asymptotic expansions and generalizes a “multipulse” Melnikov integral. It allows one to find not only multipulse heteroclinic orbits, but exotic periodic orbits. The maps exhibit singular behavior, including regions of infinite winding. These maps are shown to be singular versions of the conservative Ikeda map from laser physics and connections are made with problems from celestial mechanics and fluid mechanics. © 2008 American Institute of Physics. ͓DOI: 10.1063/1.2904823͔ Solitary waves are solutions to time-dependent partial generally lose their distinct identities and merge into a single differential equations (PDEs) which can be described as larger bump. At the opposite extreme are completely inte- profiles of unchanging spatial shape that move at con- grable systems whose solitary waves ͑solitons͒ survive col- stant velocity. A fundamental question is what happens lisions with their identities intact, due to a linear structure when two of them collide.
    [Show full text]
  • Is Type 1 Diabetes a Chaotic Phenomenon?
    Is type 1 diabetes a chaotic phenomenon? Jean-Marc Ginoux, Heikki Ruskeepää, Matjaž Perc, Roomila Naeck, Véronique Di Costanzo, Moez Bouchouicha, Farhat Fnaiech, Mounir Sayadi, Takoua Hamdi To cite this version: Jean-Marc Ginoux, Heikki Ruskeepää, Matjaž Perc, Roomila Naeck, Véronique Di Costanzo, et al.. Is type 1 diabetes a chaotic phenomenon?. Chaos, Solitons and Fractals, Elsevier, 2018, 111, pp.198-205. 10.1016/j.chaos.2018.03.033. hal-02194779 HAL Id: hal-02194779 https://hal.archives-ouvertes.fr/hal-02194779 Submitted on 29 Jul 2019 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. Is type 1 diabetes a chaotic phenomenon? Jean-Marc Ginoux,1, † Heikki Ruskeep¨a¨a,2, ‡ MatjaˇzPerc,3, 4, 5, § Roomila Naeck,6 V´eronique Di Costanzo,7 Moez Bouchouicha,1 Farhat Fnaiech,8 Mounir Sayadi,8 and Takoua Hamdi8 1Laboratoire d’Informatique et des Syst`emes, UMR CNRS 7020, CS 60584, 83041 Toulon Cedex 9, France 2University of Turku, Department of Mathematics and Statistics, FIN-20014 Turku, Finland 3Faculty of Natural Sciences and Mathematics, University of
    [Show full text]
  • A Novel Measure Inspired by Lyapunov Exponents for the Characterization of Dynamics in State-Transition Networks
    entropy Article A Novel Measure Inspired by Lyapunov Exponents for the Characterization of Dynamics in State-Transition Networks Bulcsú Sándor 1,* , Bence Schneider 1, Zsolt I. Lázár 1 and Mária Ercsey-Ravasz 1,2,* 1 Department of Physics, Babes-Bolyai University, 400084 Cluj-Napoca, Romania; [email protected] (B.S.); [email protected] (Z.I.L.) 2 Network Science Lab, Transylvanian Institute of Neuroscience, 400157 Cluj-Napoca, Romania * Correspondence: [email protected] (B.S.); [email protected] (M.E.-R.) Abstract: The combination of network sciences, nonlinear dynamics and time series analysis provides novel insights and analogies between the different approaches to complex systems. By combining the considerations behind the Lyapunov exponent of dynamical systems and the average entropy of transition probabilities for Markov chains, we introduce a network measure for character- izing the dynamics on state-transition networks with special focus on differentiating between chaotic and cyclic modes. One important property of this Lyapunov measure consists of its non-monotonous dependence on the cylicity of the dynamics. Motivated by providing proper use cases for studying the new measure, we also lay out a method for mapping time series to state transition networks by phase space coarse graining. Using both discrete time and continuous time dynamical systems the Lyapunov measure extracted from the corresponding state-transition networks exhibits similar behavior to that of the Lyapunov exponent. In addition, it demonstrates a strong sensitivity to boundary crisis suggesting applicability in predicting the collapse of chaos. Keywords: Lyapunov exponents; state-transition networks; time series analysis; dynamical systems Citation: Sándor, B.; Schneider, B.; Lázár, Z.I.; Ercsey-Ravasz, M.
    [Show full text]
  • Emergent Chaos in the Verge of Phase Transitions
    UNIVERSIDAD DE LOS ANDES BACHELOR THESIS Emergent chaos in the verge of phase transitions Author: Supervisor: Jerónimo VALENCIA Dr. Gabriel TÉLLEZ A thesis submitted in fulfillment of the requirements for the degree of Bachelor in Physics in the Department of Physics January 14, 2019 iii Declaration of Authorship I, Jerónimo VALENCIA, declare that this thesis titled, “Emergent chaos in the verge of phase transitions” and the work presented in it are my own. I confirm that: • This work was done wholly or mainly while in candidature for a research degree at this University. • Where any part of this thesis has previously been submitted for a degree or any other qualification at this University or any other institution, this has been clearly stated. • Where I have consulted the published work of others, this is always clearly attributed. • Where I have quoted from the work of others, the source is always given. With the exception of such quotations, this thesis is entirely my own work. • I have acknowledged all main sources of help. • Where the thesis is based on work done by myself jointly with others, I have made clear exactly what was done by others and what I have con- tributed myself. Signed: Date: v “The fluttering of a butterfly’s wing in Rio de Janeiro, amplified by atmospheric currents, could cause a tornado in Texas two weeks later.” Edward Norton Lorenz “We avoid the gravest difficulties when, giving up the attempt to frame hypotheses con- cerning the constitution of matter, we pursue statistical inquiries as a branch of rational mechanics” Josiah Willard Gibbs vii UNIVERSIDAD DE LOS ANDES Abstract Faculty of Sciences Department of Physics Bachelor in Physics Emergent chaos in the verge of phase transitions by Jerónimo VALENCIA The description of phase transitions on different physical systems is usually done using Yang and Lee, 1952, theory.
    [Show full text]
  • Application of Lyapunov Exponents to Strange Attractors and Intact & Damaged Ship Stability
    Application of Lyapunov Exponents to Strange Attractors and Intact & Damaged Ship Stability William R. Story Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University in partial fulfillment of the requirements for the degree of Master of Science In Ocean Engineering Leigh McCue, Chair Alan Brown Wayne Neu April 29, 2009 Blacksburg, Virginia Tech Keywords: Stability, Capsize, Lyapunov, Attractor, Lorenz Application of Lyapunov Exponents to Strange Attractors and Intact & Damaged Ship Stability William R. Story (ABSTRACT) The threat of capsize in unpredictable seas has been a risk to vessels, sailors, and cargo since the beginning of a seafaring culture. The event is a nonlinear, chaotic phenomenon that is highly sensitive to initial conditions and difficult to repeatedly predict. In extreme sea states most ships depend on an operating envelope, relying on the operator’s detailed knowledge of headings and maneuvers to reduce the risk of capsize. While in some cases this mitigates this risk, the nonlinear nature of the event precludes any certainty of dynamic vessel stability. This research presents the use of Lyapunov exponents, a quantity that measures the rate of trajectory separation in phase space, to predict capsize events for both intact and damaged stability cases. The algorithm searches backwards in ship motion time histories to gather neighboring points for each instant in time, and then calculates the exponent to measure the stretching of nearby orbits. By measuring the periods between exponent maxima, the lead‐ time between period spike and extreme motion event can be calculated. The neighbor‐ searching algorithm is also used to predict these events, and in many cases proves to be the superior method for prediction.
    [Show full text]
  • Fdlyapu: Estimating the Lyapunov Exponents Spectrum Sylvain Mangiarotti & Mireille Huc 2019-08-29
    FDLyapu: Estimating the Lyapunov exponents spectrum Sylvain Mangiarotti & Mireille Huc 2019-08-29 The FDLyapu package is an extention of the GPoM package1. Its aim is to assess the spectrum of the Lyapunov exponents for dynamical systems of polynomial form. The fractal dimension can also be deduced from this spectrum. The algorithms available in the package are based on the algebraic formulation of the equations. The GPoM package, which enables the numerical formulation of algebraic equations in a polynomial form, is thus directly required for this purpose. A user-friendly interface is also provided with FDLyapu, although the codes can be used in a blind mode. The connexion with GPoM is very natural here since this package aims to obtain Ordinary Differential Equations (ODE) from observational time series. The present package can thus be applied afterwards to characterize the global models obtained with GPoM but also to any dynamical system of ODEs expressed in polynomial form. Algorithms The Lyapunov exponents quantify the rate of separation of infinitesimally close trajectories. Depending on the initial conditions, this rate can be positive (divergence) or negative (convergence). A n-dimensionnal system is characterized by n characteristic Lyapunov exponents. For continuous dynamical systems, one Lyapunov exponent must correspond to the direction of the flow and should thus equal zero in average. To estimate the Lyapunov exponents spectrum, two algorithms are made available to the user, both based on the formal derivation of the Jacobian matrix (derivation is actually semi-formal only since the numeric coefficients are used in the derivation process). ¤ Method 1 (Wolf) The first algorithm made available in the package was introduced by Wolf et al.
    [Show full text]
  • Intermittent Large Deviation of Chaotic Trajectory in Ikeda Map: Signature of Extreme Events Arnob Ray,1 Sarbendu Rakshit,1 Dibakar Ghosh,1, A) and Syamal K
    Intermittent large deviation of chaotic trajectory in Ikeda map: Signature of extreme events Arnob Ray,1 Sarbendu Rakshit,1 Dibakar Ghosh,1, a) and Syamal K. Dana2, 3 1)Physics and Applied Mathematics Unit, Indian Statistical Institute, 203 B. T. Road, Kolkata 700108, India 2)Department of Mathematics, Jadavpur University, Kolkata 700032, India 3)Division of Dynamics, Technical University of Lodz, Stefanowskiego 1/15, 90-924 Lodz, Poland We notice signatures of extreme events-like behavior in a laser based Ikeda map. The trajectory of the system occasionally travels a large distance away from the bounded chaotic region, which appears as intermittent spiking events in the temporal dynamics. The large spiking events satisfy the conditions of extreme events as usually observed in dynamical systems. The probability density function of the large spiking events shows a long-tail distribution con- sistent with the characteristics of rare events. The inter-event intervals obey a Poisson-like distribution. We locate the parameter regions of extreme events in phase diagrams. Furthermore, we study two Ikeda maps to explore how and when extreme events terminate via mutual interaction. A pure diffusion of in- formation exchange is unable to terminate extreme events where synchronous occurrence of extreme events is only possible even for large interaction. On the other hand, a threshold-activated coupling can terminate extreme events above a critical value of mutual interaction. PACS numbers: 89.75.Fb, 05.45.a Rare and recurrent large amplitude deviations of normally bounded dynamics are seen in many systems. Such occasional large amplitude spiking events are larger than a nominal value and their statistical distribution of occurrence shows qualitative similarities, in dynamical sense, with data records of natural disas- ters, rogue waves, tsunami, flood and share market crashes.
    [Show full text]