Energy Transfer and Dissipation in Nonlinear Oscillators

Total Page:16

File Type:pdf, Size:1020Kb

Energy Transfer and Dissipation in Nonlinear Oscillators Energy transfer and dissipation in nonlinear oscillators Maaita Tzamal- Odysseas Physicist, MSc. Computational Physics A thesis submitted for the degree of Doctor of Philosophy in the Department of Physics School of Sciences Aristotle University of Thessaloniki Supervisor: Assistant Professor Efthymia Meletlidou March, 2014 To my grandmothers Mdallal Enayat and Polidou Marika 3 Abstract (In Greek) Sthn παρούσα didaktορική διατριβή μελετάμε èna σύσthma tri¸n suzeugmènwn talantwt¸n me τριβή, dύο grammik¸n me ènan mh γραμμικό. Tètoia susτήμαta talantwt¸n èqoun μεγάλο endiafèron idiaÐtera όtan h μάζα tou mh γραμμικού talantωτή eÐnai poλύ μικρόterh από touc γραμμικούc me sunèpeia o mh γραμμικός talantωτής na leitourgeÐ wc katαβόθρα enèrgeiac. Autού tou eÐdouc ta susτήματa, sta opoÐa συνυπάρχουν ènac αργός kai ènac γρήγορος χρόnoc mpoρούν na μελετηθούν me th βοήθεια thc singularity analysis, twn analloÐwtwn pollaploτήτwn, en¸ σημαντική plhroforÐa gia th δυναμική tou susτήματoc dÐnetai kai από thn δυναμική thc αργής ροής (Slow Flow) tou susτήματoc. Sthn παρούσα διατριβή melετάμε to σύσthma mèsw thc melèthc thc αργής analloÐwthc poλλαπλόthtac (Slow Invariant Manifold- SIM-). Me th βοήθεια tou θεωρήματoc tou Tikhonov kathgoriopoiούμε tic διάφορες peript¸seic thc αργής analloÐwthc poλλαπλόthtac kai orÐzoume αναλυτικά tic συνθήκες me tic opoÐec μπορούμε na οδηγηθούμε sthn κάθε perÐptwsh. Se επόμενο βήμα μελετάμε thn δυναμική thc argήc rοής kai παρατηρούμε όti h δυναμική thc eÐnai πλούσια αφού oi troqièc thc μπορούν na eÐnai kanonikèc, na κάνουν talant¸seic hremÐac (relaxation oscillations), ή na eÐnai qaotikèc. Από thn melèth thc enèrgeiac pou αποθηκεύετai ston mh γραμμικό talantωτή kai από ton ρυθμό απόσβεσης thc sunoλικής enèrgeiac tou susτήματoc παρατηρούμε όti tόσο h ύπαρξη diaklad¸sewn thc argής analloÐwthc poλλαπλόthtac, όσο kai h dunaμική thc αργής ροής paÐzoun kajorisτικό ρόlo sthn metαφορά enèrgeiac apό ton γραμμικό ston mh γραμμικό talantωτή. EpÐshc, stic peript¸seic pou blèpoume metαφορά enèrgeiac παρατηρούμε όti o ρυθμός apόσβεσης thc sunoλικής enèrgeiac tou susτήματoc eÐnai megaλύτερος από ton ρυθμό απόσβεσης όtan den metafèretai enèrgeia ston mh γραμμικό talantωτή. H melèth tou susτήματoc twn tri¸n suzeugmènwn talantwt¸n kleÐnei me thn πρόtash ενός mh γραμμικού hlektrikou kukl¸matoc to opoÐo ulopoieÐ thn mh γραμμική διαφορική exÐswsh δεύτερης τάξης me thn opoÐa proseggÐsame to αρχικό sύsthma. To sugkekrimèno κύκλωμα èqei endiafèron giatÐ mac dÐnei th δυνατόthta na μελετήσουμε kai πειραματικά διάφορα apό ta faiνόμενα pou eÐdame sthn θεωρητική mac ανάλυση. 5 Prologue I was introduced to the subject of this thesis by the late Associate Proffessor of the department of Physics, Simos Ichtiaroglou. Simos Ichtiaroglou participated in the second International Conference on Nonlinear Normal Modes and Localization in Vibrating Systems that took place in Samos at 2006, where he attended a number of paper presentations related to targeted energy transfer. So, when we were discussing the topic of my thesis he proposed to me and to Assistant Professor Efi Meletlidou to work on energy transfer and dissipation because he considered it interesting. The topic was indeed interesting. It included rich theoretical study and practical applications, so we decided to start work on it. The specific system composed of two coupled linear oscillators and a nonlinear oscillator that interacts through an essential nonlinearity with one of the linear oscillators was proposed by proffessor Alex Vakakis in the Euromech Colloquium 503 - Nonlinear Normal Modes, Dimension Reduction and Localization in Vibrating Systems- that took place in Frascati(Rome) at 2009. The work on this thesis was both interesting and challenging. It was interesting because I came into contact with new theories and tools, such as the singular perturbation theory, the multiple scale analysis, the complexification- averaging technique etc... that I was not taught neither in my undergraduate nor graduate studies. It was also interesting because I came into contact with other scientific fields, such as engineering, and I learned a lot of new things, as for example, how we can implement a nonlinear-cubic- term in a mechanical structure. It was challenging, mainly, because there was no financial support and I was obliged to work all the years of my doctoral dissertation outside the university. Our work resulted in four publications in peer reviewed international journals. The first paper is titled “The effect of slow flow dynamics on the oscillations of a singular damped system with an essentially nonlinear attachment” and it was published in the Journal of Applied Nonlinear dynamics in 2013. The second paper is titled “Analytical Homoclinic Solution of a Two-Dimensional Nonlinear System of Differential Equations” and it was published in the 7 Journal of Nonlinear Dynamics in 2013. The third paper is titled “The Study of a Nonlinear Duffing - Type Oscillator Driven by Two Voltage Sources” and it was published in Journal of Engineering Science and Technology Review in 2013. Finaly the fourth paper is titled “The dynamics of the slow flow of a singular damped nonlinear system and its Parametric Study” and it is accepted for publication in the Journal of Applied Nonlinear Dynamics. This doctorate thesis marks the end of a stage and the begining of a new one in my life. It is therefore time for reasoning. I want to thank my supervisor, Assistant Professor, Efi Meletlidou, for helping me to become a researcher, for teaching me many things, especially on mathematics. She was always understanding and patient and despite the difficulties, we managed to cooperate well. I want to thank my advisor, Associate Professor Iannis Kyprianidis, for his support all these years. His help was important for the implementation of the electric circuit, for the numerical work on computing the maximal Lyapunov Characteristic Exponent. I want to thank my advisor, Professor Alex Vakakis, for the proposal of the system and his guidance on its study. His encouraging words always gave me strength and courage to continue. I want to thank, Dr. Christos Volos, for his help on the inplementation and study of the electric circuit and teaching me the Multisim program. I want to thank, Associate Professor Vassilis Rothos, for his advice on books and papers. Professor George Tsaklidis was a mentor on how to write, check and correct a scientific manuscript. He is a good advisor and a friend. I want to thank, Assistant Professor Menios Tsiganis, for his help and advice, especially, for the Lyapunov Characteristic Exponents and the numerical work. We always have good debates and discussions. I want to thank, Professor Harry Varvoglis, for being there from the beginning of my PhD. Although we disagreed on different issues, he was always there, honest and taking an interest. Assistan Professors George Voyatzis and Harry Skokos, and Special Teaching Fellow Fotini Zervaki, thank you for the support on difficult times. I want to thank, Professor Antonis Anagnostopoulos, for his advice and remarks. Special thanks to my “officemates”, Dr. Stella Tzirti and PhD Candidate Kallinikos Nikos, for the good collaboration and the constructive discussion all these years. I want to thank my friend, PhD Candidate Nikos Bastas, for his help on Grid computing. I want to thank, from the depths of my heart, my parents Ifigeneia and Sami Maaita for all the things that they have done to let me become what i am today. My mother for being by my side, working with me on school, supporting me on difficult times. My father for his guidance, rational thinking and encouragment. Also, my brother Tarek Maaita for always being there when I need him. My wife Kiki Moysiadou for her constant support, for her critical remarks and the continuous 8 interaction. For giving birth to my two beloved little girls Tereza- Leila, Ifigeneia- Yasmin and the one that we are waiting to come next June. Finaly, I want to thank all those who contributed in one way or another in the formation of my being. Maaita Jamal- Odysseas, Thessaloniki, Sunday 30th March, 2014. 9 Contents 1 Introduction 15 2 Basic concepts, theorems and tools 21 2.1 Dynamical Systems................................ 21 2.1.1 General definitions ............................ 21 2.1.2 Equilibrium points (Periodic orbits) and stability: definitions and theorems 22 2.1.3 Bifurcation theory............................. 25 2.1.4 Chaos................................... 27 2.1.5 The Averaging method.......................... 29 2.2 Singular Perturbation Theory........................... 30 2.2.1 What is a singular perturbation problem? ................ 30 2.2.2 Multiple Scale Analysis.......................... 30 2.2.3 Geometric Singular Perturbation Theory................. 33 2.3 Electric circuits.................................. 36 2.3.1 Linear and nonlinear electric elements.................. 36 2.3.2 Resistor connections and basic laws ................... 39 2.3.3 Operational Amplifier........................... 40 2.3.4 Examples of nonlinear circuits...................... 41 3 Classification of the SIM 47 3.1 Reduction of the system.............................. 48 3.2 Singular Perturbation Analysis.......................... 51 3.3 Solutions and numerical results.......................... 54 3.4 Conclusions.................................... 57 4 The dynamics of the slow flow 61 4.1 The dynamics of the slow flow.......................... 61 4.2 Conclusions.................................... 71 10 5 Energy transfer and dissipation
Recommended publications
  • Using Standard Syste
    PHYSICAL REVIEW E VOLUME 55, NUMBER 6 JUNE 1997 Cumulative Author Index All authors of papers published in this volume are listed alphabetically. Full titles are included in each ®rst author's entry. For Rapid Communi- cations an R precedes the page number. The letters (C) and (BR) following the page number indicate that a paper is a Comment or a Brief Report, respectively. References with (E) are to Errata. Abarbanel, Henry D. I. — ͑see Huerta, Ramon͒ E 55, R2108 and Jin-Qing Fang — Synchronization of chaos and hyperchaos ͑see Liu, Clif͒ E 55, 6483 using linear and nonlinear feedback functions. E 55, 5285 Abdullaev, F. Kh. and J. G. Caputo — Propagation of an envelope Aliaga, J. — ͑see Gruver, J. L.͒ E 55, 6370 soliton in a medium with spatially varying dispersion. E 55, 6061 Alig, I. — ͑see Mayer, W.͒ E 55, 3102 Abel, T., E. Brener, and H. Mu¨ller-Krumbhaar — Three-dimensional Allahyarov, E. A., L. I. Podloubny, P. P. J. M. Schram, and S. A. growth morphologies in diffusion-controlled channel growth. Trigger — Damping of longitudinal waves in colloidal crystals E 55, 7789͑BR͒ of finite size. E 55, 592 Abramson, Guillermo — Ecological model of extinctions. E 55, 785 Allain, C. — ͑see Senis, D.͒ E 55, 7797͑BR͒ Acharyya, M., J. K. Bhattacharjee, and B. K. Chakrabarti Allegrini, M. — ͑see Bezuglov, N. N.͒ E 55, 3333 — Dynamic response of an Ising system to a pulsed field. E 55, 2392 Allen, Christopher K. and Martin Reiser — Bunched beam envelope Ackerson, Bruce J. — ͑see Paulin, S. E.͒ E 55, 5812 equations including image effects from a cylindrical pipe.
    [Show full text]
  • Fourth SIAM Conference on Applications of Dynamical Systems
    Utah State University DigitalCommons@USU All U.S. Government Documents (Utah Regional U.S. Government Documents (Utah Regional Depository) Depository) 5-1997 Fourth SIAM Conference on Applications of Dynamical Systems SIAM Activity Group on Dynamical Systems Follow this and additional works at: https://digitalcommons.usu.edu/govdocs Part of the Physical Sciences and Mathematics Commons Recommended Citation Final program and abstracts, May 18-22, 1997 This Other is brought to you for free and open access by the U.S. Government Documents (Utah Regional Depository) at DigitalCommons@USU. It has been accepted for inclusion in All U.S. Government Documents (Utah Regional Depository) by an authorized administrator of DigitalCommons@USU. For more information, please contact [email protected]. tI...~ Confers ~'t' '"' \ 1I~c9 ~ 1'-" ~ J' .. c "'. to APPLICAliONS cJ May 18-22, 1997 Snowbird Ski and Summer Resort • Snowbird, Utah Sponsored by SIAM Activity Group on Dynamical Systems Conference Themes The themes of the 1997 conference will include the following topics. Principal Themes: • Dynamics in undergraduate education • Experimental studies of nonlinear phenomena • Hamiltonian systems and transport • Mathematical biology • Noise in dynamical systems • Patterns and spatio-temporal chaos Applications in • Synchronization • Aerospace engineering • Biology • Condensed matter physics • Control • Fluids • Manufacturing • Me;h~~~~nograPhY 19970915 120 • Lasers and o~ • Quantum UldU) • 51a m.@ Society for Industrial and Applied Mathematics http://www.siam.org/meetingslds97/ds97home.htm 2 " DYNAMICAL SYSTEMS Conference Prl Contents A Message from the Conference Chairs ... Get-Togethers 2 Dear Colleagues: Welcoming Message 2 Welcome to Snowbird for the Fourth SIAM Conference on Applications of Dynamica Systems. Organizing Committee 2 This highly interdisciplinary meeting brings together a diverse group of mathematicians Audiovisual Notice 2 scientists, and engineers, all working on dynamical systems and their applications.
    [Show full text]
  • Hopf-Zero Bifurcation of the Ring Unidirectionally Coupled Toda Oscillators with Delay∗
    Nonlinear Analysis: Modelling and Control, Vol.26, No. 3, 375–395 https://doi.org/10.15388/namc.2021.26.23050 Press Hopf-zero bifurcation of the ring unidirectionally coupled Toda oscillators with delay∗ Rina Sua;b, Chunrui Zhangc;1 aCollege of Mechanical and Electrical Engineering, Northeast Forestry University, Harbin, 150040, China [email protected] bCollege of Mathematics and Physics, Inner Mongolia University for Nationalities, Tongliao, 028000, China cDepartment of Mathematics, Northeast Forestry University, Harbin 150040, China [email protected] Received: November 18, 2019 / Revised: November 11, 2020 / Published online: May 1, 2021 Abstract. In this paper, the Hopf-zero bifurcation of the ring unidirectionally coupled Toda oscillators with delay was explored. First, the conditions of the occurrence of Hopf-zero bifurcation were obtained by analyzing the distribution of eigenvalues in correspondence to linearization. Second, the stability of Hopf-zero bifurcation periodic solutions was determined based on the discussion of the normal form of the system, and some numerical simulations were employed to illustrate the results of this study. Lastly, the normal form of the system on the center manifold was derived by using the center manifold theorem and normal form method. Keywords: Toda oscillators, normal form, Hopf-zero. 1 Introduction In more recent decades, the bifurcation theory in dynamical system has become a research hotspot, which is widely used in the fields of physics, chemistry, medicine, finance, biol- ogy, engineering and so on [3,12,13,15,19,20,22,24,30–35]. In particular, the bifurcation phenomenon of the system is analyzed, and the dynamic characteristics of the nonlinear ∗This research was supported by the Science Research Funding of the Inner Mongolia University for the Nationalities under grant No.
    [Show full text]
  • Nonlinear Stochastic Dynamics and Chaos by Numerical Path Integration
    Nonlinear stochastic dynamics and chaos by numerical path integration Eirik Mo January 5, 2008 Acknowledgements The main gratitude goes to Professor Arvid Naess, who has been the advisor for this work, and given many valuable comments. The problem description was written by Naess, the development of path integration for stochastic differential equations is largely devoted to him, and he also applied for and received funding for the project. For four years, I have worked together with Dr. Hans Christian Karlsen, who also was a PhD student advised by Naess. Karlsen's work on running path integration for some specific problems was of great importance for the development of more stable and efficient code. His validation of accuracy and general comments have been very helpful for this work. From August 2004 until February 2005, I visited Professor Mircea Grigoriu at the Department of Structural Engineering at Cornell University, beautifully situated in Ithaca by the Finger Lakes, upstate New York. This was a wonderful stay, where I were warmly welcomed by PhD students at the department, and a lot of people in the area. I would especially like to mention Ms. (Kathleen) Misti Wilcox, who rented out rooms in her house, cooked wonderful meals, proudly presented Ithaca and the surrounding area, and motivated me for working. This visit was a big inspiration, and a turning point in the research. The Research Council of Norway funded the project and the visit to Cornell University, of which I'm very grateful. Motivating words from people around me have been important from time to time.
    [Show full text]
  • Appendix a One-Dimensional Maps
    Appendix A One-Dimensional Maps A.I Introduction The objective of this Appendix is to provide an introduction to the theory of chaos. One-dimensional maps are used as a vehicle to convey the concepts required for study­ ing chaos, firstly because fewer concepts are involved compared to higher dimensional maps and differential equations and secondly because of the fact that one-dimensional maps and their evolutions may easily be visualized, which helps to grasp ideas in­ volved more rapidly. For generalization of these result s to, and additional result s on higher dimensional mappings and differential equations, the reader is referred to the bibliography. Some important theorems are stated in Appendix B without proof. Consider an arbitrary one dimensional, N-th order discrete dynamical system whose time evolution is described by the homogeneous recurrence relation X n+l = f( X n, Xn- l, .•. , X n - N + l ), where f is some function N f: R -t R, and X n is the dynamical system's evolution on the n-th discrete time inst ant , i.e. n = 0,1,2, .... In general f may be nonlinear and is usually at least piecewise differentiable. The function f which describ es some dynamical system is termed a mapping or simply a map. Example A.I Consider the logistic map described by Xn+l = f(x n ), (A.l) where f : I -t R, f(x) := Tx(l - x ), (A.2) 343 344 APPENDIX A. ONE-DIMENSIONAL MAPS 0.9 0.8 0.7 0.6 .-... ~ '-" 0.5 '+-.., 0.4 0.3 0.2 0.1 o x Figure A.I: The logistic map for r = 4 .
    [Show full text]
  • From Hilbert to Dilbert
    From Hilbert to Dilbert A non-orthodox approach to gravitation, psychosynthesis, economics, cosmology, and other issues Victor Christianto & Florentin Smarandache (editors) From Hilbert to Dilbert 1 From Hilbert to Dilbert A non-orthodox approach to gravitation, psychosynthesis, economics, cosmology, and other issues ISBN : Editors: • Victor Christianto, DDiv., Malang Institute of Agriculture (IPM), Indonesia. http://www.sci- 4God.com • Florentin Smarandache, PhD., Dept. Mathematics and Sciences, University of New Mexico, New Mexico, USA. Gallup, New Mexico 87301, USA; e-mail: [email protected] Publisher : Divine Publishing Jl. Sidomakmur No. 75 A Sengkaling Malang, INDONESIA 65151 Email: [email protected] ISBN : 978 - 1 - 59973 - 585 - 6 Cover design : d.v.n. project Layout : d.v.n. project No part of this book is allowed to be copied and distributed by any means, be it in print or electronically, without written permission from the Editors and/or the Publisher, except for citation and review purposes. 2 From Hilbert to Dilbert Dedication ~ A tribute to Mr. Scott Adams the creator of Dilbert’s cartoon series. To many people like us, Dilbert is the real anti-hero of this corporatocracym age… (our wish is that eminent directors such as Oliver Stone will someday create something like “Dilbert the Movie.” — don’t laugh, we’re serious) From Hilbert to Dilbert 3 4 From Hilbert to Dilbert Preface This book took an unconventional theme because we submit an unorthodox theme too. Karl Popper’s epistemology suggests that when the theory is refuted by observation, then it is time to look for a set of new approaches.
    [Show full text]
  • Solving Large Number of Non-Stiff, Low-Dimensional Ordinary
    Solving large number of non-stiff, low-dimensional ordinary differential equation systems on GPUs and CPUs: performance comparisons of MPGOS, ODEINT and DifferentialEquations.jl D´aniel Nagya, Lambert Plavecza, Ferenc Heged˝usa,∗ aDepartment of Hydrodynamic Systems, Faculty of Mechanical Engineering, Budapest University of Technology and Economics, Budapest, Hungary Abstract In this paper, the performance characteristics of different solution techniques and pro- gram packages to solve a large number of independent ordinary differential equation systems is examined. The employed hardware are an Intel Core i7-4820K CPU with 30.4 GFLOPS peak double-precision performance per cores and an Nvidia GeForce Ti- tan Black GPU that has a total of 1707 GFLOPS peak double-precision performance. The tested systems (Lorenz equation, Keller–Miksis equation and a pressure relief valve model) are non-stiff and have low dimension. Thus, the performance of the codes are not limited by memory bandwidth, and Runge–Kutta type solvers are efficient and suitable choices. The tested program packages are MPGOS written in C++ and specialised only for GPUs; ODEINT implemented in C++, which supports execution on both CPUs and GPUs; finally, DifferentialEquations.jl written in Julia that also supports execution on both CPUs and GPUs. Using GPUs, the program package MPGOS is superior. For CPU computations, the ODEINT program package has the best performance. Keywords: ordinary differential equations, non-stiff problems, GPU programming, CPU programming 1. Introduction Low dimensional ordinary differential equation (ODE) systems are still widely used in arXiv:2011.01740v1 [cs.DC] 2 Nov 2020 physics and non-linear dynamical analysis. Among the simplest models (second or third- order systems), one can still find many publications employing the classic Duffng [1–6], Morse [7–10], Toda [11–14], Brusselator [15, 16], van der Pol [17–19] and Lorenz [20–23] equations.
    [Show full text]
  • Appendix Fm =F~M =~(Am +Ibm)
    Appendix A.l Fourier Series [Benedetto 1996 B, Bracewell1986 B, Carslaw 1930 B] A .1.1 General Rules period: x(t) = x(t + T) basic circular frequency: OJ OJ 2n / T T = 2n / OJ = 2n / OJ1 = 1 = complex Fourier series: complex coefficients: m=oo 1 T x(t) = x(t+ T) = IJmexp(-imOJt) Fm = - f x(t) exp(+imOJt)dt m=-oo T o real periodic functions: coefficients: x(t) = x(t+ T) = x*(t) Fo = Fo* = llo = Ao Fm =F~m =~(am +ibm) =.!. ~exp(iam} 2 real Fourier series I: real coefficients: x(t) = x(t+ T) = x*(t) 1 T 00 00 ao = - f x(t)dt T o =ao+ Lamcos(mOJt) + Lbmsin(mOJt) m=1 m=1 2 T am = - f x(t)cos(mOJt)dt T o 2 T bm =-f x(t)sin(mOJt)dt T o real Fourier series TI: real parameters: x(t) = x{t+ T) = x*{t) Ao =ao 00 [2 2]112 = Ao + L~cos(mOJt-am} Am = am +bm ~O m=1 am = arc tan(bm / am} 480 real even periodic functions: coefficients: x{t) = x{t + T) = x * (t) = x{ -t) Fo = Fo* = ao = Ao 1 1 Fm =F-m =Fm *=-a 2 m =-A2 m um =O;bm =0 real odd periodie functions: coefficients: x{t) = x{t+ T) = x*{t) = -x{-t) Fo =ao =Ao = 0 Fm =-F-m =-Fm *=..!..ib 2 m um = ±tr 12;am = 0 bm = Am sign( Um) A .1. 2 Real Periodie Functions Norrnalized parameters: Cü =Cü, = 1 and T =2tr a) Square wave: x{t) = x{t + 2tr) = sign(sin t) = i[sint + ..!..sin 3t + ..!..sin 5t + ..
    [Show full text]
  • Lyapunov Potential Description for Laser Dynamics
    Lyapunov Potential Description for Laser Dynamics Catalina Mayol†,‡, Ra´ul Toral†,‡ and Claudio R. Mirasso† (†) Departament de F´ısica, Universitat de les Illes Balears (‡) Instituto Mediterr´aneo de Estudios Avanzados (IMEDEA, UIB-CSIC) E-07071 Palma de Mallorca, Spain. Web Page: http://www.imedea.uib.es We describe the dynamical behavior of both class A and class B lasers in terms of a Lyapunov potential. For class A lasers we use the potential to analyze both deterministic and stochastic dynamics. In the stochastic case it is found that the phase of the electric field drifts with time in the steady state. For class B lasers, the potential obtained is valid in the absence of noise. In this case, a general expression relating the period of the relaxation oscillations to the potential is found. We have included in this expression the terms corresponding to the gain saturation and the mean value of the spontaneously emitted power, which were not considered previously. The validity of this expression is also discussed and a semi-empirical relation giving the period of the relaxation oscillations far from the stationary state is proposed and checked against numerical simulations. PACs: 42.65.Sf, 42.55.Ah, 42.60.Mi, 42.55.Px I. INTRODUCTION This model for laser dynamics was constructed from the Maxwell-Bloch equations for a single-mode field interact- Even for non-mechanical systems, it is occasionally ing with a two-level medium. The semiclassical laser the- possible to construct a function (called Lyapunov func- ory ignores the quantum-mechanical nature of the elec- tion or Lyapunov potential) that decreases along trajec- tromagnetic field, and the amplifying medium is modeled tories [1].
    [Show full text]