On localized solutions of chains of oscillators with cubic nonlinearity


 Bernadette Gardner
 11 months ago
 Views:
Transcription
1 On localized solutions of chains of oscillators with cubic nonlinearity Francesco Romeo, Giuseppe Rega Dipartimento di Ingegneria Strutturale e Geotecnica, SAPIENZA Università di Roma, Italia Keywords: localization, wave propagation, nonlinear maps. SUMMARY. Stationary localized solutions arising from the dynamics of chains of linearly coupled mechanical oscillators characterized by on site cubic nonlinearity are addressed aiming to explore their relationship with the underlying nonlinear wave propagation regions. 1 INTRODUCTION In this study the dynamics of chains of linearly coupled mechanical oscillators characterized by on site cubic nonlinearity is addressed. It is well known that such nonlinear periodic systems possess propagation and attenuation regions in the frequency domain. Nonlinear travelling wave solutions exist inside the propagation zones the boundaries of which can be analytically determined through a map approach when dealing with a discrete lattice structure [1]. Given the nature of the considered maps, namely reversible area preserving maps, periodic orbits can be also identified by resorting to their symmetry lines, i.e. fixed points of the involutions with which the map can be factorized. The analytical description of both, propagation regions boundaries and periodic orbits, enables to predict the nature of the periodic orbits and thereby to distinguish between hyperbolic and elliptic fixed points. Such a distinction plays a crucial role for the onset of stationary localized solutions in the nonlinear discrete lattice under investigation. A number of studies have indeed recently confirmed the ability of nonlinear lattices to support timeperiodic and spatially localized solutions (discrete breathers) by linking them to the map homoclinic and heteroclinic orbits [2]. According to this geometric approach discrete breathers are identified with sequences of homoclinic or heteroclinic intersection points of the mapping so that the unstable fixed points govern their existence since stable and unstable manifolds emanate from them. The cases in which the trivial solution is hyperbolic or elliptic, namely outside or inside the linear passing band, respectively, were extensively investigated to show the existence of a variety of breather solutions. Simple, solitonlike solutions were identified with sequences of homoclinic or heteroclinic primary intersection points. Less insight is available for the other unstable fixed points belonging to the symmetry lines; in this respect, information drawn by both nonlinear propagation regions and internal bifurcation thresholds, merged with the higher order symmetry lines, are useful to locate and classify the relevant localized solutions, along with their features. In Section 1 the equations of motion of the mechanical system are introduced and the main aspects of its Hamiltonian dynamics described by a planar map are recalled. In particular, the emphasis is placed on the chain propagation regions showing how they affect the nature of the fixed point of the underlying map. In Section 2 the symmetry lines of the nonlinear map are obtained by means of recurrence formulas; these lines are used to identify the periodic orbits and to give more insight on the map phase portraits. In Section 3, the nonlinear propagation regions are used to locate and classify different types of localized solutions given by sequences of homoclinic or heteroclinic intersection points of the mapping. Some remarks and future developments are discussed in the concluding section. 1
2 2 THE OSCILLATORY CHAIN MODEL A simple mass and spring chain is considered to analyze the regions of existence of nonlinear localization phenomenon and their connection with the underlying nonlinear propagation regions. Discrete models similar to the studied one are often used to describe the nonlinear dynamics of micromechanical arrays [3, 4]. In this context the oscillation of a single cantilever is modeled following the EulerBernoulli beam theory, characterized by the cubic hardening nonlinearity. In essence through a Galerkin approach, in which the beam displacement is expanded in terms of the first linear mode, a sdof Duffing type equation is obtained, where the oscillator displacements describe the tip displacement amplitude of the nth cantilever. Moreover, in the simplest model, the driven cantilever array accounts for nearest neighbor interaction only. In this study we consider a chain of linearly coupled nonlinear oscillators described by the following equations of motion m q n + k 1 q n + k 3 q 3 n + k(2q n q n 1 q n+1 )=0 (1) Time periodic solutions of equation (1) are sought for by assuming the harmonic solution q n = a n cos(ωt). Equating coefficients of cos(ωt) and setting α = mω2 k 1 k 2 and β = 3k3 4k the following secondorder difference equation for the stationary amplitude is obtained (α + βa 2 n )a n + a n+1 + a n 1 =0 (2) In the frequency dependent parameter α in (2), k is the linear coupling stiffness while m and k 1 represent mass and linear stiffness of the single oscillator, respectively. In the nonlinear parameter β, k 3 represents the softening (β >0) or hardening (β <0) oscillator cubic stiffness. Equation (2), relating the amplitudes a in adjacent chain sites n 1, n and n +1, can be rewritten in matrix form as a n+1 = T(a n )a n, where a n+1 =(a n+1,a n ) and a n =(a n,a n 1 ); T(a n ) is the frequency dependent nonlinear transfer matrix belonging to the class of area preserving maps such that det(dt(a n )) = 1, where DT is the Jacobian or tangent map with reciprocal eigenvalues. Therefore, in order to study the stationary wave transmission properties of the onedimensional nonlinear chain (1) of length N, it is convenient to rely on the eigenvalues of the linearized map equations in the neighborhood of an orbit ranging from (a 0,a 1 ) to (a N 1,a N ). Indeed, according to [5], by considering the introduced map as a dynamical system where the chain index n plays the role of discrete time n, the analysis of the transmission properties is equivalent to the stability analysis of the orbits. The interest lies in the linear stability of spatially periodic orbits a n+m = a n with cycle length m. The eigenvalues of DT are determined by its trace. If the eigenvalues lie on the unit circle, then stable elliptic periodic cycles or oscillating solutions occur; if the eigenvalues are real, then unstable hyperbolic periodic cycles or exponentially increasing solutions occur. By setting a n+1 = x n+1 and a n = y n+1, the map T can be written as T : x n+1 = y n x n (βx 2 n + α) y n+1 = x n (3) The stability of periodm orbits is governed by the condition tr(dt m ) =2; so, for m =1, the stability condition gives r := {(x, α, β) x 2 + α +2=0} s := {(x, α, β) x 2 + α 2=0} (4) 2
3 Table 1: Hyperbolic (r) and reflection hyperbolic (s) thresholds obtained from tr(dt m ) =2and mperiodic orbits p. m r s p 1 ± 2 α x = y = ± β 2 2 α ±i α x = y = 2 α β 3 4 ± 2 α 1 α ± 1 α ±i α 2 α ± 2 α 2 α x = y = 1 α β x = ±y = α β and the curves r and s represent hyperbolic (λ 1 = λ 2 =1) and reflection hyperbolic (λ 1 = λ 2 = 1) thresholds, respectively [1]. For periodm orbits, with m =1, 2, 3, 4, the latter curves are reported in Table 1. 3 SYMMETRY LINES AND PERIODIC ORBITS In this section the special properties of the map (3) are exploited to define and compute its symmetry lines with a twofold goal: to identify the periodic orbits and to give more insight on the map phase portraits. We begin by considering that the planar map (3) is reversible since it satisfies the condition R T R = T 1 where denotes the composition of two maps and R is called a reversing operator of T [6]; reversible maps have a number of special properties and occur often in applications. In particular, T is an area preserving map (4) for which the reversor R is an involution I 0, namely I 2 0 = 1. The higher order involutions are defined as I m := T m I 0, m = ±1, ±2,... The fixed points of the involutions with which the map can be factorized provide the symmetry lines, defined as γ m = Fix(I m )={ξ I m (ξ) =(ξ), ξ R} (5) The set γ m of I m is called the mth order symmetry line satisfying the recursion relation (l = 0, 1, 2,...) [6] The main symmetry lines of the map (3) are summarized in Table 2, where the change of notation = n+1, = n is introduced. By defining the spatially mperiodic orbits by p m := {ξ T m (ξ) =(ξ)} the following relation 3
4 T l (γ m )=γ 2l+m m =0 T l (γ 0 )=γ 2l even symmetry line m =1 T l (γ 1 )=γ 2l+1 odd symmetry line Table 2: Main symmetry lines of the map (3) for β>0 and β<0. β>0 I 0 : I 1 : x = y y = x x = x y = y x(α + βx 2 ) γ 0 : y = x γ 1 : y = 1 2 x(βx2 + α) β<0 I 0 I 2 : x = y y = x x = x(α + βx 2 )+y y = x x(α + βx 2 ) 2 y(α + βx 2 ) γ 0 : y = x γ 2 : y = x( 1+α + βx 2 ) between periodic orbits and symmetry lines holds γ m γ l p m l (6) The symmetry lines and the relation (6) can be used to find spatially periodic solution of any order; in the last column of Table 1 they are reported for m =1, 2, 3, 4. In Figure 1 and Figure 2 the first 10 symmetry lines, generated by the involutions reported in Table 2, are superimposed on the phase portraits of the map (3) for two values of α. The periodic solutions in the figures are discussed against the T propagation region in Fig. 3a. In Figure 1, where α = 0.5 and β =1.0, the elliptic fixed point at (0, 0) as well as the period4 orbits, located at x = ±y = ± 0.5 (both inside the pass band in Fig. 3a), are clearly highlighted by the relevant symmetry lines intersections obeying to (6); the quasiperiodic orbits scenario surrounding these points is also clearly related to the symmetry lines (Figure 1a). In Figure 1b, the symmetry lines are shown on a wider x, y range showing a large number of intersections (i.e. periodic solutions) confined in the range 2.5 x 2.5 (light gray); these boundaries correspond to hyperbolic period2 orbits. Figure 2 refers to α = 3.0 and β =1.0; in Figure 2a the phase portrait refers to orbits corresponding to positive initial conditions x 0,y 0. The period1 orbit at x = y =1undergoes a quadrupling bifurcation and the arrangement of the associated quasiperiodic orbits is governed by the main symmetry lines. The global view in Figure 2b shows that the scenario is mirrored with respect to both the main diagonals and that the symmetry lines intersections (periodic solutions) occur only within two separated regions given by 5 x 1 and 1 x 5. These regions are bounded below by the elliptic period1 orbits and above by the hyperbolic period2 ones. The analytical description of both, propagation regions boundaries and periodic orbits, enables to predict the nature of the periodic orbits and thereby to distinguish between hyperbolic and elliptic 4
5 Figure 1: Phase portrait of the map for α = 0.5, β =1.0 and the first 10 symmetry lines generated by the involutions I 0, I 1 in Table 2. solutions. As it will be shown in the next section, such a distinction plays a crucial role for the onset of stationary localized solutions in the nonlinear discrete lattice under investigation. 4 LOCALIZED SOLUTIONS AGAINST PROPAGATION REGIONS According to [7] the discrete breathers, i.e. timeperiodic spatially localized solutions, are solutions of (1) with period T b =2πΩ b that can be expanded in a Fourier series as follows q n (t) = k A kn e ikω bt with A k, n 0 (7) where the assumption of Fourier coefficients localized in space is introduced. Following the approach proposed in [8] and [9] for hardening and softening type of nonlinearity, respectively, we look for localized states by relating them to the intersections of the invariant manifolds emanating from unstable fixed points of the map (3). As known, hyperbolic fixed points P possess stable W s (P ) and unstable W u (P ) manifolds such that a point belonging to W s (P )(W u (P )) approaches P under map iteration T n for n (n ). A homoclinic point Q is defined as Q W s W u ; thus, any Q belongs to both the stable and unstable manifold of P so, when iterated either forward or backward, it will converge to P and it corresponds to a breather solution. Due to the entangled pattern of the manifolds, additional intersection points are generated implying an infinite number of different breathers and multibreathers, where more than one site exhibits large amplitude oscillations [7]. In Figure 3 the regions of existence of localized solutions are shown on the (α, x) plane for β =1.0 (Figure 3a) and β = 1.0 (Figure 3b); these regions are classified according to the nature of the map T fixed points, namely the trivial solution (x = y =0) and the period1 solutions (x = y = β ). For the sake of completeness, the period2 orbits are also reported as they play a twofold role. As shown in Figures 1b and 2b, they represent the upper boundary of all the periodic solutions; moreover, their location on the x = y symmetry line affects the topology of the invariant manifolds. Regions I +,I correspond to discrete breathers associated to the homoclinic connections emanating 5
6 Figure 2: Phase portrait of the map for α = 3.0, β =1.0 and the first 10 symmetry lines generated by the involutions I 0, I 1 in Table 2. from the trivial solutions for β>0 and β<0, respectively. Region II corresponds to nested homoclinic connections arising from the period doubling bifurcations series of period1 orbits beginning at the intersection with the propagation region boundary s; for β =1, the latter bifurcation series occurs in the range α 4.0. As shown in [10], each period doubling bifurcation can be identified by the intersections of the newborn elliptic solutions and the relevant stype higher order internal thresholds. Region III refers to the heteroclinic connections between the hyperbolic period1 orbits. To validate the analytical predictions, numerical simulations pertaining to localized solutions belonging to regions I +, II and III are reported in the following. Figure 3: Projection on the (α, x) plane of the regions of existence of localized solutions against the T propagation regions. The dashed and solid lines correspond to hyperbolic and elliptic points, respectively; a) β = 1.0 (softening nonlinearity); b) β = 1.0 (hardening nonlinearity). In Figure 4 the localization scenario corresponding to the parameters α = 2.8, β = 1 be 6
7 longing to region I + is illustrated. In Figure 4a stable (blue) and unstable (red) manifolds of (0, 0) surrounding the elliptic period1 points (green) located at x = y = ±0.894 are shown; in Figure 4b a closer view of the manifolds is displayed with two families of homoclinic orbits; the ensuing discrete breathers are also shown in Figures 4c,d. Figure 4: Stable (blue) and unstable (red) manifolds of (0, 0) at α = 2.8, β =1(region I + ) and two homoclinic orbits; a) global view; b) closer view with intersections forming two families of homoclinic orbits; c) discrete breather corresponding to black dots in b); discrete breather corresponding to white dots in b). Figure 5 refers to region II localization, characterized by nested homoclinic orbits at α = 4.05). As shown by the global view of the invariant manifolds (Figure 5a), besides the global homoclinic connection emanating from the trivial solution, two more homoclinic tangles arise. They stem from new born period1 hyperbolic branch following the period doubling bifurcations occurring at (± 2.0, ± 2.0) for α = 4.0. The closer view in Figure 5b shows one of the two homoclinic connections in which the newborn elliptic points are located along the γ 1 (see Table 2) symmetry line (dashed line). The associated homoclinic orbit (black dots) gives rise to the multibreather displayed in Figure 5d, whereas Figure 5c refers to a breather corresponding to a global homoclinic orbit. The numerical results presented in Figure 6 are related to the localized solutions associated with the heteroclinic connections occurring in region III. In particular, the invariant manifolds of ( 2, 2) and ( 2, 2) shown in Figure 6a are obtained for the parameters α =0.0, β = 1. The black 7
8 Figure 5: Stable (blue) and unstable (red) manifolds of (0, 0) and (± 2.05, ± 2.05) at α = 4.05, β =1and homoclinic orbits; a) global view; b) closer view of the manifolds of ( 2.05, 2.05); c) discrete breather corresponding to a homoclinic orbit of (0, 0); discrete multibreather corresponding to a homoclinic orbit of ( 2.05, 2.05). dots correspond to the discrete breather depicted in Figure 6b. 5 CONCLUSIONS The nonlinear propagation regions of chains of oscillators with cubic nonlinearity are used to identify the regions of existence of discrete breathers, i.e. time periodic and spatially localized solutions, and to guide their analysis. By describing the map symmetry lines through recurrence formulas, the phase portrait topological features and the periodic orbits arrangement have been discussed. Moreover, by merging the analytical descriptions of the map symmetry lines and propagation bands, the regions where the spatially periodic orbits are hyperbolic (saddle) or elliptic are determined. Since discrete breathers correspond to sequences of homoclinic or heteroclinic intersection points, the unstable fixed points govern their existence. Besides the known regions corresponding to localized solutions stemming from the trivial fixed points, located either outside (homoclinic connections) or inside (heteroclinic connections) the linear passing band, a further region has been identified. It corresponds to nested homoclinic orbits, emanating from period1 unstable solutions, giving rise to multibreathers. In agreement with recent literature findings, the larger region of existence of 8
9 Figure 6: Stable (blue) and unstable (red) manifolds of ( 2, 2) and ( 2, 2) at α =0.0, β = 1; a) heteroclinic orbits; b) discrete breather corresponding to the heteroclinic orbit highlighted in a). breathers for the hardening cubic nonlinearity with respect to the softening one has been confirmed. Further studies will address the stability analysis in time of the different families of localized solutions. References [1] Romeo, F., Rega, G., Wave propagation properties of chains of oscillators with cubic nonlinearities via nonlinear map approach, Chaos, Solitons and Fractals, 27, (2006). [2] Hennig, D., Rasmussen, K.., Gabriel, H. and Blow, A., Solitonlike solutions of the generalized discrete nonlinear Schrödinger equation, Phys. Rev. E, 54, (1996). [3] Sato, M., Hubbard, B.E., Sievers, A.J., Nonlinear energy localization and its manipulation in micromechanical oscillator arrays, Rev. Modern Phys., 78, , (2006). [4] Sato, M., Sievers, A.J., Driven localized excitations in the acoustic spectrum of small nonlinear macroscopic and microscopic lattices, Phys. Rev. Lett., 98, , (2007). [5] Hennig, D., Tsironis, G.P., Wave transmission in nonlinear lattices, Physics Reports, 307, , (1999). [6] Shi, Y.G., Chen, L., Reversible maps and their symmetry lines, Commun Nonlinear Sci Numer Simulat., 16, (2011). [7] Flach, S., Gorbach, A.V., Discrete breathers Advances in theory and applications, Phys. Reports, 467, (2008). [8] Bountis, T., Capel, H.W., Kollmann, M., Ross, J.C., Bergamin, J.M., van der Weele J.P., Multibreathers and homoclinic orbits in 1dimensional nonlinear lattices, Phys. Lett. A, 268, 5060, (2000). 9
10 [9] Panagopoulos, P., Bountis, T., Skokos, C., Existence and stability of localized oscillations in 1dimensional lattices with softspring and hardspring potentials, J. of Vib. and Acoustics, 126, , (2004). [10] Romeo, F., Rega, G., Propagation properties of bicoupled nonlinear oscillatory chains: analytical prediction and numerical validation, Int. J. of Bif. and Chaos, 18, , (2008). 10
arxiv:nlin/ v1 [nlin.cd] 17 Jan 2002
Symbolic Dynamics of Homoclinic Orbits in a Symmetric Map arxiv:nlin/0201031v1 [nlin.cd] 17 Jan 2002 ZaiQiao Bai and WeiMou Zheng Institute of Theoretical Physics, Academia Sinica, Beijing 100080, China
More informationHamiltonian Chaos and the standard map
Hamiltonian Chaos and the standard map Outline: What happens for small perturbation? Questions of long time stability? Poincare section and twist maps. Area preserving mappings. Standard map as time sections
More informationNONSTATIONARY RESONANCE DYNAMICS OF THE HARMONICALLY FORCED PENDULUM
CYBERNETICS AND PHYSICS, VOL. 5, NO. 3, 016, 91 95 NONSTATIONARY RESONANCE DYNAMICS OF THE HARMONICALLY FORCED PENDULUM Leonid I. Manevitch Polymer and Composite Materials Department N. N. Semenov Institute
More informationTitle Coupled Cantilever Arrays( Disserta. https://doi.org/ /doctor.k14
Title Studies on the Manipulation of Intr Coupled Cantilever Arrays( Disserta Author(s) Kimura, Masayuki Citation Kyoto University ( 京都大学 ) Issue Date 29323 URL https://doi.org/1.14989/doctor.k14 Right
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A
More informationChapter 23. Predicting Chaos The Shift Map and Symbolic Dynamics
Chapter 23 Predicting Chaos We have discussed methods for diagnosing chaos, but what about predicting the existence of chaos in a dynamical system. This is a much harder problem, and it seems that the
More informationBifurcations of Traveling Wave Solutions for a Generalized CamassaHolm Equation
Computational and Applied Mathematics Journal 2017; 3(6): 5259 http://www.aascit.org/journal/camj ISSN: 23811218 (Print); ISSN: 23811226 (Online) Bifurcations of Traveling Wave Solutions for a Generalized
More informationInvariant manifolds of the Bonhoeffervan der Pol oscillator
Invariant manifolds of the Bonhoeffervan der Pol oscillator R. Benítez 1, V. J. Bolós 2 1 Dpto. Matemáticas, Centro Universitario de Plasencia, Universidad de Extremadura. Avda. Virgen del Puerto 2. 10600,
More informationSuppression of the primary resonance vibrations of a forced nonlinear system using a dynamic vibration absorber
Suppression of the primary resonance vibrations of a forced nonlinear system using a dynamic vibration absorber J.C. Ji, N. Zhang Faculty of Engineering, University of Technology, Sydney PO Box, Broadway,
More informationSymbolic dynamics and chaos in plane Couette flow
Dynamics of PDE, Vol.14, No.1, 7985, 2017 Symbolic dynamics and chaos in plane Couette flow Y. Charles Li Communicated by Y. Charles Li, received December 25, 2016. Abstract. According to a recent theory
More informationBarcelona, Spain. RTBP, collinear points, periodic orbits, homoclinic orbits. Resumen
XX Congreso de Ecuaciones Diferenciales y Aplicaciones X Congreso de Matemática Aplicada Sevilla, 2428 septiembre 27 (pp. 1 8) The dynamics around the collinear point L 3 of the RTBP E. Barrabés 1, J.M.
More informationAdditive resonances of a controlled van der PolDuffing oscillator
Additive resonances of a controlled van der PolDuffing oscillator This paper has been published in Journal of Sound and Vibration vol. 5 issue  8 pp.. J.C. Ji N. Zhang Faculty of Engineering University
More informationChaotic Synchronization of Symbolic Information. in the Discrete Nonlinear Schrödinger Equation. Abstract
Chaotic Synchronization of Symbolic Information in the Discrete Nonlinear Schrödinger Equation. C. L. Pando L. IFUAP, Universidad Autónoma de Puebla, arxiv:nlin/3511v1 [nlin.cd] 8 May 23 Apdo. Postal J48.
More informationExample of a Blue Sky Catastrophe
PUB:[SXG.TEMP]TRANS2913EL.PS 16OCT2001 11:08:53.21 SXG Page: 99 (1) Amer. Math. Soc. Transl. (2) Vol. 200, 2000 Example of a Blue Sky Catastrophe Nikolaĭ Gavrilov and Andrey Shilnikov To the memory of
More informationApproximation of Solitons in the Discrete NLS Equation
University of Massachusetts Amherst ScholarWorks@UMass Amherst Mathematics and Statistics Department Faculty Publication Series Mathematics and Statistics 8 Approximation of Solitons in the Discrete NLS
More informationSelfReplication, SelfDestruction, and SpatioTemporal Chaos in the GrayScott Model
Letter Forma, 15, 281 289, 2000 SelfReplication, SelfDestruction, and SpatioTemporal Chaos in the GrayScott Model Yasumasa NISHIURA 1 * and Daishin UEYAMA 2 1 Laboratory of Nonlinear Studies and Computations,
More informationarxiv: v1 [nlin.cd] 20 Jul 2010
Invariant manifolds of the Bonhoeffervan der Pol oscillator arxiv:1007.3375v1 [nlin.cd] 20 Jul 2010 R. Benítez 1, V. J. Bolós 2 1 Departamento de Matemáticas, Centro Universitario de Plasencia, Universidad
More informationOne Dimensional Dynamical Systems
16 CHAPTER 2 One Dimensional Dynamical Systems We begin by analyzing some dynamical systems with onedimensional phase spaces, and in particular their bifurcations. All equations in this Chapter are scalar
More informationChapter 3. GumowskiMira Map. 3.1 Introduction
Chapter 3 GumowskiMira Map 3.1 Introduction Non linear recurrence relations model many real world systems and help in analysing their possible asymptotic behaviour as the parameters are varied [17]. Here
More informationSection 9.3 Phase Plane Portraits (for Planar Systems)
Section 9.3 Phase Plane Portraits (for Planar Systems) Key Terms: Equilibrium point of planer system yꞌ = Ay o Equilibrium solution Exponential solutions o Halfline solutions Unstable solution Stable
More informationA Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term
ETASR  Engineering, Technology & Applied Science Research Vol., o.,, 95 9 A Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term Fei Yu College of Information Science
More informationYulin, A. V., & Champneys, A. R. (2009). Discrete snaking: multiple cavity solitons in saturable media.
Yulin, A. V., & Champneys, A. R. (2009). Discrete snaking: multiple cavity solitons in saturable media. Early version, also known as preprint Link to publication record in Explore Bristol Research PDFdocument
More informationEffect of various periodic forces on Duffing oscillator
PRAMANA c Indian Academy of Sciences Vol. 67, No. 2 journal of August 2006 physics pp. 351 356 Effect of various periodic forces on Duffing oscillator V RAVICHANDRAN 1, V CHINNATHAMBI 1, and S RAJASEKAR
More informationTORSION PENDULUM: THE MECHANICAL NONLINEAR OSCILLATOR
TORSION PENDULUM: THE MECHANICAL NONLINEAR OSCILLATOR Samo Lasič, Gorazd Planinšič,, Faculty of Mathematics and Physics University of Ljubljana, Slovenija Giacomo Torzo, Department of Physics, University
More informationHopf Bifurcation Analysis and Approximation of Limit Cycle in Coupled Van Der Pol and Duffing Oscillators
The Open Acoustics Journal 8 93 9 Open Access Hopf ifurcation Analysis and Approximation of Limit Cycle in Coupled Van Der Pol and Duffing Oscillators Jianping Cai *a and Jianhe Shen b a Department of
More information8.1 Bifurcations of Equilibria
1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integrodifferential equations
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More informationPhonons and lattice dynamics
Chapter Phonons and lattice dynamics. Vibration modes of a cluster Consider a cluster or a molecule formed of an assembly of atoms bound due to a specific potential. First, the structure must be relaxed
More informationDifference Resonances in a controlled van der PolDuffing oscillator involving time. delay
Difference Resonances in a controlled van der PolDuffing oscillator involving time delay This paper was published in the journal Chaos, Solitions & Fractals, vol.4, no., pp.97598, Oct 9 J.C. Ji, N. Zhang,
More informationINVARIANT CURVES AND FOCAL POINTS IN A LYNESS ITERATIVE PROCESS
International Journal of Bifurcation and Chaos, Vol. 13, No. 7 (2003) 1841 1852 c World Scientific Publishing Company INVARIANT CURVES AND FOCAL POINTS IN A LYNESS ITERATIVE PROCESS L. GARDINI and G. I.
More informationCOMPLEX DYNAMICS AND CHAOS CONTROL IN DUFFINGVAN DER POL EQUATION WITH TWO EXTERNAL PERIODIC FORCING TERMS
International J. of Math. Sci. & Engg. Appls. (IJMSEA) ISSN 09739424, Vol. 9 No. III (September, 2015), pp. 197210 COMPLEX DYNAMICS AND CHAOS CONTROL IN DUFFINGVAN DER POL EQUATION WITH TWO EXTERNAL
More informationStrange dynamics of bilinear oscillator close to grazing
Strange dynamics of bilinear oscillator close to grazing Ekaterina Pavlovskaia, James Ing, Soumitro Banerjee and Marian Wiercigroch Centre for Applied Dynamics Research, School of Engineering, King s College,
More informationDaba Meshesha Gusu and O.Chandra Sekhara Reddy 1
International Journal of Basic and Applied Sciences Vol. 4. No. 1 2015. Pp.2227 Copyright by CRDEEP. All Rights Reserved. Full Length Research Paper Solutions of Non Linear Ordinary Differential Equations
More informationHomoclinic bifurcations in Chua s circuit
Physica A 262 (1999) 144 152 Homoclinic bifurcations in Chua s circuit Sandra Kahan, Anibal C. SicardiSchino Instituto de Fsica, Universidad de la Republica, C.C. 30, C.P. 11 000, Montevideo, Uruguay
More informationGrowth of Sobolev norms for the cubic defocusing NLS Lecture 1
Growth of Sobolev norms for the cubic defocusing NLS Lecture 1 Marcel Guardia February 8, 2013 M. Guardia (UMD) Growth of Sobolev norms February 8, 2013 1 / 40 The equations We consider two equations:
More informationIntroduction LECTURE 1
LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in
More informationIntroduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams.
Outline of Continuous Systems. Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Vibrations of Flexible Strings. Torsional Vibration of Rods. BernoulliEuler Beams.
More informationClearly the passage of an eigenvalue through to the positive real half plane leads to a qualitative change in the phase portrait, i.e.
Bifurcations We have already seen how the loss of stiffness in a linear oscillator leads to instability. In a practical situation the stiffness may not degrade in a linear fashion, and instability may
More informationKronigPenney model. 2m dx. Solutions can be found in region I and region II Match boundary conditions
KronigPenney model I II d V( x) m dx E Solutions can be found in region I and region II Match boundary conditions Linear differential equations with periodic coefficients Have exponentially decaying solutions,
More informationPseudoChaotic Orbits of Kicked Oscillators
Dynamical Chaos and NonEquilibrium Statistical Mechanics: From Rigorous Results to Applications in NanoSystems August, 006 PseudoChaotic Orbits of Kicked Oscillators J. H. Lowenstein, New York University
More informationSimulation of twolevel quantum system using classical coupled oscillators. Abstract
Simulation of twolevel quantum system using classical coupled oscillators Namit Anand and Ashok Mohapatra School of Physical Sciences, National Institute of Science Education and Research, Bhubaneswar
More informationSymmetry Properties of Confined Convective States
Symmetry Properties of Confined Convective States John Guckenheimer Cornell University 1 Introduction This paper is a commentary on the experimental observation observations of Bensimon et al. [1] of convection
More informationIntroduction to Applied Nonlinear Dynamical Systems and Chaos
Stephen Wiggins Introduction to Applied Nonlinear Dynamical Systems and Chaos Second Edition With 250 Figures 4jj Springer I Series Preface v L I Preface to the Second Edition vii Introduction 1 1 Equilibrium
More informationTorus Maps from Weak Coupling of Strong Resonances
Torus Maps from Weak Coupling of Strong Resonances John Guckenheimer Alexander I. Khibnik October 5, 1999 Abstract This paper investigates a family of diffeomorphisms of the two dimensional torus derived
More informationAMADEU DELSHAMS AND RAFAEL RAMíREZROS
POINCARÉMELNIKOVARNOLD METHOD FOR TWIST MAPS AMADEU DELSHAMS AND RAFAEL RAMíREZROS 1. Introduction A general theory for perturbations of an integrable planar map with a separatrix to a hyperbolic fixed
More informationTWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 SYMMETRY. P. Yu 1,2 and M. Han 1
COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 3, Number 3, September 2004 pp. 515 526 TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 SYMMETRY P. Yu 1,2
More informationMathematical Foundations of Neuroscience  Lecture 7. Bifurcations II.
Mathematical Foundations of Neuroscience  Lecture 7. Bifurcations II. Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland Winter 2009/2010 Filip
More information2.034: Nonlinear Dynamics and Waves. Term Project: Nonlinear dynamics of piecewise linear oscillators Mostafa Momen
2.034: Nonlinear Dynamics and Waves Term Project: Nonlinear dynamics of piecewise linear oscillators Mostafa Momen May 2015 Massachusetts Institute of Technology 1 Nonlinear dynamics of piecewise linear
More informationModels Involving Interactions between Predator and Prey Populations
Models Involving Interactions between Predator and Prey Populations Matthew Mitchell Georgia College and State University December 30, 2015 Abstract Predatorprey models are used to show the intricate
More informationStabilization of Hyperbolic Chaos by the Pyragas Method
Journal of Mathematics and System Science 4 (014) 75576 D DAVID PUBLISHING Stabilization of Hyperbolic Chaos by the Pyragas Method Sergey Belyakin, Arsen Dzanoev, Sergey Kuznetsov Physics Faculty, Moscow
More information= F ( x; µ) (1) where x is a 2dimensional vector, µ is a parameter, and F :
1 Bifurcations Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 A bifurcation is a qualitative change
More informationIgor A. Khovanov,Vadim S. Anishchenko, Astrakhanskaya str. 83, Saratov, Russia
Noise Induced Escape from Dierent Types of Chaotic Attractor Igor A. Khovanov,Vadim S. Anishchenko, Dmitri G. Luchinsky y,andpeter V.E. McClintock y Department of Physics, Saratov State University, Astrakhanskaya
More informationDynamics of four coupled phaseonly oscillators
Communications in Nonlinear Science and Numerical Simulation 13 (2008) 501 507 www.elsevier.com/locate/cnsns Short communication Dynamics of four coupled phaseonly oscillators Richard Rand *, Jeffrey
More informationMAPPING CLASS ACTIONS ON MODULI SPACES. Int. J. Pure Appl. Math 9 (2003),
MAPPING CLASS ACTIONS ON MODULI SPACES RICHARD J. BROWN Abstract. It is known that the mapping class group of a compact surface S, MCG(S), acts ergodically with respect to symplectic measure on each symplectic
More informationHamiltonian Dynamics
Hamiltonian Dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS Feb. 10, 2009 Joris Vankerschaver (CDS) Hamiltonian Dynamics Feb. 10, 2009 1 / 31 Outline 1. Introductory concepts; 2. Poisson brackets;
More informationSubharmonic Oscillations and Chaos in Dynamic Atomic Force Microscopy
Subharmonic Oscillations and Chaos in Dynamic Atomic Force Microscopy John H. CANTRELL 1, Sean A. CANTRELL 2 1 NASA Langley Research Center, Hampton, Virginia 23681, USA 2 NLS Analytics, LLC, Glencoe,
More informationNBA Lecture 1. Simplest bifurcations in ndimensional ODEs. Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011
NBA Lecture 1 Simplest bifurcations in ndimensional ODEs Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011 Contents 1. Solutions and orbits: equilibria cycles connecting orbits other invariant sets
More informationAPPLIED SYMBOLIC DYNAMICS AND CHAOS
DIRECTIONS IN CHAOS VOL. 7 APPLIED SYMBOLIC DYNAMICS AND CHAOS BaiLin Hao WeiMou Zheng The Institute of Theoretical Physics Academia Sinica, China Vfö World Scientific wl Singapore Sinaaoore NewJersev
More information7 Twodimensional bifurcations
7 Twodimensional bifurcations As in onedimensional systems: fixed points may be created, destroyed, or change stability as parameters are varied (change of topological equivalence ). In addition closed
More informationAn Exactly Solvable 3 Body Problem
An Exactly Solvable 3 Body Problem The most famous nbody problem is one where particles interact by an inverse squarelaw force. However, there is a class of exactly solvable nbody problems in which
More informationDynamics of Modified LeslieGower PredatorPrey Model with Predator Harvesting
International Journal of Basic & Applied Sciences IJBASIJENS Vol:13 No:05 55 Dynamics of Modified LeslieGower PredatorPrey Model with Predator Harvesting K. Saleh Department of Mathematics, King Fahd
More informationExploiting pattern transformation to tune phononic band gaps in a twodimensional granular crystal
Exploiting pattern transformation to tune phononic band gaps in a twodimensional granular crystal The Harvard community has made this article openly available. Please share how this access benefits you.
More information16 Period doubling route to chaos
16 Period doubling route to chaos We now study the routes or scenarios towards chaos. We ask: How does the transition from periodic to strange attractor occur? The question is analogous to the study of
More informationSolutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10)
Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10) Mason A. Porter 15/05/2010 1 Question 1 i. (6 points) Define a saddlenode bifurcation and show that the first order system dx dt = r x e x
More informationNonlinear dynamics & chaos BECS
Nonlinear dynamics & chaos BECS114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes
More informationSolitons. Nonlinear pulses and beams
Solitons Nonlinear pulses and beams Nail N. Akhmediev and Adrian Ankiewicz Optical Sciences Centre The Australian National University Canberra Australia m CHAPMAN & HALL London Weinheim New York Tokyo
More informationSufficient conditions for a period incrementing big bang bifurcation in onedimensional maps.
Sufficient conditions for a period incrementing big bang bifurcation in onedimensional maps. V. Avrutin, A. Granados and M. Schanz Abstract Typically, big bang bifurcation occur for one (or higher)dimensional
More informationSTABILITY OF ORBITS NEAR LARGE MASS RATIO BINARY SYSTEMS
IAAAASDyCoSS20508 STABILITY OF ORBITS NEAR LARGE MASS RATIO BINARY SYSTEMS Natasha Bosanac, Kathleen C. Howell and Ephraim Fischbach INTRODUCTION With recent scientific interest into the composition,
More informationA classification of explosions in dimension one
Ergod. Th. & Dynam. Sys. (29), 29, 715 731 doi:1.117/s143385788486 c 28 Cambridge University Press Printed in the United Kingdom A classification of explosions in dimension one E. SANDER and J. A. YORKE
More informationExtremely slow edge waves in mechanical graphene with rotating grains
Phononic Crystals and Acoustic Metamaterials: Paper ICA2016104 Extremely slow edge waves in mechanical graphene with rotating grains LiYang Zheng (a), Vincent Tournat (b), Georgios Theocharis (c),vitalyi
More informationOn Riddled Sets and Bifurcations of Chaotic Attractors
Applied Mathematical Sciences, Vol. 1, 2007, no. 13, 603614 On Riddled Sets and Bifurcations of Chaotic Attractors I. Djellit Department of Mathematics University of Annaba B.P. 12, 23000 Annaba, Algeria
More information2.10 Saddles, Nodes, Foci and Centers
2.10 Saddles, Nodes, Foci and Centers In Section 1.5, a linear system (1 where x R 2 was said to have a saddle, node, focus or center at the origin if its phase portrait was linearly equivalent to one
More informationThe SD oscillator and its attractors
The SD oscillator and its attractors Qingjie Cao, Marian Wiercigroch, Ekaterina Pavlovskaia Celso Grebogi and J Michael T Thompson Centre for Applied Dynamics Research, Department of Engineering, University
More informationCalculus and Differential Equations II
MATH 250 B Second order autonomous linear systems We are mostly interested with 2 2 first order autonomous systems of the form { x = a x + b y y = c x + d y where x and y are functions of t and a, b, c,
More informationGeometry of Crystal Lattice
0 Geometry of Crystal Lattice 0.1 Translational Symmetry The crystalline state of substances is different from other states (gaseous, liquid, amorphous) in that the atoms are in an ordered and symmetrical
More informationME 680 Spring Geometrical Analysis of 1D Dynamical Systems
ME 680 Spring 2014 Geometrical Analysis of 1D Dynamical Systems 1 Geometrical Analysis of 1D Dynamical Systems Logistic equation: n = rn(1 n) velocity function Equilibria or fied points : initial conditions
More informationBIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs
BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.
More informationComplex Behavior in Coupled Nonlinear Waveguides. Roy Goodman, New Jersey Institute of Technology
Complex Behavior in Coupled Nonlinear Waveguides Roy Goodman, New Jersey Institute of Technology Nonlinear Schrödinger/GrossPitaevskii Equation i t = r + V (r) ± Two contexts for today: Propagation of
More informationPhysics 211B : Problem Set #0
Physics 211B : Problem Set #0 These problems provide a cross section of the sort of exercises I would have assigned had I taught 211A. Please take a look at all the problems, and turn in problems 1, 4,
More informationHardening nonlinearity effects on forced vibration of viscoelastic dynamical systems to external step perturbation field
Int. J. Adv. Appl. Math. and Mech. 3(2) (215) 16 32 (ISSN: 23472529) IJAAMM Journal homepage: www.ijaamm.com International Journal of Advances in Applied Mathematics and Mechanics Hardening nonlinearity
More informationInfluence of moving breathers on vacancies migration
Physics Letters A 315 (2003) 364 371 www.elsevier.com/locate/pla Influence of moving breathers on vacancies migration J. Cuevas a,, C. Katerji a, J.F.R. Archilla a,j.c.eilbeck b,f.m.russell b a Grupo de
More informationTWO VAN DER POLDUFFING OSCILLATORS WITH HUYGENS COUPLING
ENOC8, Saint Petersburg, Russia, June, 3 July, 8 TWO VAN DER POLDUFFING OSCILLATORS WITH HUYGENS COUPLING V.N. Belykh, E.V. Pankratova Mathematics Department Volga State Academy Nizhny Novgorod, 3, Russia.
More informationTwo models for the parametric forcing of a nonlinear oscillator
Nonlinear Dyn (007) 50:147 160 DOI 10.1007/s1107100691483 ORIGINAL ARTICLE Two models for the parametric forcing of a nonlinear oscillator Nazha Abouhazim Mohamed Belhaq Richard H. Rand Received: 3
More informationUNIFORM SUBHARMONIC ORBITS FOR SITNIKOV PROBLEM
Manuscript submitted to Website: http://aimsciences.org AIMS Journals Volume 00, Number 0, Xxxx XXXX pp. 000 000 UNIFORM SUBHARMONIC ORBITS FOR SITNIKOV PROBLEM CLARK ROBINSON Abstract. We highlight the
More informationOn low speed travelling waves of the KuramotoSivashinsky equation.
On low speed travelling waves of the KuramotoSivashinsky equation. Jeroen S.W. Lamb Joint with Jürgen Knobloch (Ilmenau, Germany) MarcoAntonio Teixeira (Campinas, Brazil) Kevin Webster (Imperial College
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: August 22, 2018, at 08 30 12 30 Johanneberg Jan Meibohm,
More informationComputational Neuroscience. Session 42
Computational Neuroscience. Session 42 Dr. Marco A Roque Sol 06/21/2018 TwoDimensional TwoDimensional System In this section we will introduce methods of phase plane analysis of twodimensional systems.
More informationLecture 2. Contents. 1 Fermi s Method 2. 2 Lattice Oscillators 3. 3 The SineGordon Equation 8. Wednesday, August 28
Lecture 2 Wednesday, August 28 Contents 1 Fermi s Method 2 2 Lattice Oscillators 3 3 The SineGordon Equation 8 1 1 Fermi s Method Feynman s Quantum Electrodynamics refers on the first page of the first
More informationdynamical zeta functions: what, why and what are the good for?
dynamical zeta functions: what, why and what are the good for? Predrag Cvitanović Georgia Institute of Technology November 2 2011 life is intractable in physics, no problem is tractable I accept chaos
More informationDelay Coordinate Embedding
Chapter 7 Delay Coordinate Embedding Up to this point, we have known our state space explicitly. But what if we do not know it? How can we then study the dynamics is phase space? A typical case is when
More informationBifurcation Trees of Periodic Motions to Chaos in a Parametric, Quadratic Nonlinear Oscillator
International Journal of Bifurcation and Chaos, Vol. 24, No. 5 (2014) 1450075 (28 pages) c World Scientific Publishing Company DOI: 10.1142/S0218127414500758 Bifurcation Trees of Periodic Motions to Chaos
More informationHamiltonian dynamics of breathers with thirdorder dispersion
1150 J. Opt. Soc. Am. B/ Vol. 18, No. 8/ August 2001 S. Mookherjea and A. Yariv Hamiltonian dynamics of breathers with thirdorder dispersion Shayan Mookherjea* and Amnon Yariv Department of Applied Physics
More informationHow Does the Smaller Alignment Index (SALI) Distinguish Order from Chaos?
Progress of Theoretical Physics Supplement No. 150, 2003 439 How Does the Smaller Alignment Index (SALI) Distinguish Order from Chaos? Charalampos Skokos, 1,2, ) Chris Antonopoulos, 1, ) Tassos C. Bountis
More informationBreather Modes Induced by Localized RF Radiation: Analytical and Numerical Approaches
Proceedings of the 5th International Conference on Nonlinear Dynamics NDKhPI2016 September 2730, 2016, Kharkov, Ukraine Breather Modes Induced by Localized RF Radiation: Analytical and Numerical Approaches
More informationOn dynamical properties of multidimensional diffeomorphisms from Newhouse regions: I
IOP PUBLISHING Nonlinearity 2 (28) 923 972 NONLINEARITY doi:.88/95775/2/5/3 On dynamical properties of multidimensional diffeomorphisms from Newhouse regions: I S V Gonchenko, L P Shilnikov and D V Turaev
More informationEXACT DARK SOLITON, PERIODIC SOLUTIONS AND CHAOTIC DYNAMICS IN A PERTURBED GENERALIZED NONLINEAR SCHRODINGER EQUATION
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 17, Number 1, Spring 9 EXACT DARK SOLITON, PERIODIC SOLUTIONS AND CHAOTIC DYNAMICS IN A PERTURBED GENERALIZED NONLINEAR SCHRODINGER EQUATION JIBIN LI ABSTRACT.
More informationDIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS
DIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS Morris W. Hirsch University of California, Berkeley Stephen Smale University of California, Berkeley Robert L. Devaney Boston University
More information1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point
Solving a Linear System τ = trace(a) = a + d = λ 1 + λ 2 λ 1,2 = τ± = det(a) = ad bc = λ 1 λ 2 Classification of Fixed Points τ 2 4 1. < 0: the eigenvalues are real and have opposite signs; the fixed point
More informationPart 1: Fano resonances Part 2: Airy beams Part 3: Paritytime symmetric systems
Lecture 3 Part 1: Fano resonances Part 2: Airy beams Part 3: Paritytime symmetric systems Yuri S. Kivshar Nonlinear Physics Centre, Australian National University, Canberra, Australia http://wwwrsphysse.anu.edu.au/nonlinear/
More informationSimplest Chaotic Flows with Involutional Symmetries
International Journal of Bifurcation and Chaos, Vol. 24, No. 1 (2014) 1450009 (9 pages) c World Scientific Publishing Company DOI: 10.1142/S0218127414500096 Simplest Chaotic Flows with Involutional Symmetries
More information