The reduction to invariant cones for nonsmooth systems. Tassilo Küpper Mathematical Institute University of Cologne
|
|
- Kristopher Horn
- 5 years ago
- Views:
Transcription
1 The reduction to invariant cones for nonsmooth systems Tassilo Küpper Mathematical Institute University of Cologne International Workshop on Resonance Oscillations and Stability of Nonsmooth Systems 16 June to 25 June, 2009, Imperial College London, United Kingdom.
2 Aims Determine invariant manifold for non-smooth system. Extend concept of center manifold approach to non-smooth system. Bifurcation analysis of reduced system. Ex. Alteration of stable/unstable systems(damping/ excitation) Ex. Mechanics - Cell interaction: Neural dynamics - Electrical circuits. Results Stable+Stable Unstable [Carmona,et al (2005)] Generalized Hopf Bifurcation: Generation of Periodic orbits n=2 : Küpper and Moritz [2002], Beyn, Küpper and Zou (2006) n > 2 : Küpper (2008), Küpper and Hosham (2009). Reduction techniques (center manifold approach).
3 1. Smooth systems: Center manifold approach. Review ξ = f(ξ, λ), (ξ R n, λ R p ), (1) ξ = 0, stationary solution linearization ξ = A(λ)ξ + g(ξ, λ) A = f (0, λ) = ξ ( ) A1 0 ξ = 0 A 2 λ Ev ofa 1 : Reλ < 0. λ Ev ofa 2 : Reλ = 0. ( ) x. y
4 ẋ = A 1 x + g 1 (x, y, λ) (2) ẏ = A 2 y + g 2 (x, y, λ). Center manifold reduction h = h(y), h(0) = 0, h (0) = 0, x = h(y). (3) y Reduced system ( locally equivalent to (2)) ẏ = A 2 (λ)y + g 2 (h(y), y, λ), (4) Remark: dim y n, once (4) has been established the dynamics, stability and bifurcation behavior of (1) can be obtained by studying (4). First Prev Next Last Go Back Full Screen Close Quit
5 2. Non-smooth system: general analysis R n = N M i, disjoint (5) i=1 f i = R n R p R n. ξ = f i (ξ, λ), (ξ M i ). (6) ξ = 0, (7) stationary solution in some M i, but no interior point. Local reduction to equivalent system?
6 3. Simplified setting ξ = f + (ξ, λ), e T 1 ξ > 0, (R n +) f 0 (ξ, λ), e T 1 ξ = 0, (R n 0) f (ξ, λ), e T 1 ξ < 0, (R n ) (8)
7
8 Special cases on M 0 Transversal crossing [e T 1 f + (ξ, λ)][e T 1 f (ξ, λ)] > 0. Sliding [e T 1 f + (ξ, λ)][e T 1 f (ξ, λ)] < 0. Higher order conditions [Dieci and Lopez (2008,2009)]. Jumps in phase space ξ(t + ) = Rξ(t ). Ex. Impact Pendulum Ex. Bells S. Köker: Zur Dynamik des Glockenläutens, Diploma Thesis (2009).
9 Approximation by PWLS, (0 = f + (0, λ) = f (0, λ)) f ± = A ± (λ)ξ + g ± (ξ, λ), g ± smooth, h.o.t. (9) 4. PWLS (transversal) a. Def. of PWLS For fixed λ set, ξ = { A + ξ, e T 1 ξ > 0, A ξ, e T 1 ξ < 0, W ± > = {ξ M e T 1 f ± ξ > 0}, W ± < = {ξ M e T 1 f ± ξ < 0}, W ± 0 = {ξ M e T 1 f ± ξ = 0}. transverse crossing if ξ W < W + < or ξ W > W + >, (10)
10 b. Set up of Poincarè map Flow: e t (ξ)a (λ) ξ e t+ A + (λ) η For ξ W < W + < : t (ξ) = inf{t > 0 e T 1 ϕ (t, ξ) = 0}, P (ξ) := e t (ξ)a ξ For η W > W + >, t + (η) = inf{t > 0 e T 1 ϕ + (t, η) = 0}, P + (η) := e t+ (η)a + η. Return map P (ξ) := P + (P (ξ)) = e t+ (η)a + e t (ξ)a ξ. P (ξ) nonlinear t + and t nonlinear constant on half-rays t (rξ) = t (ξ)(0 < r < ) and t + (rξ) = t + (ξ)(0 < r < ).
11
12 Theorem 1. (T. Küpper (2008)) Let ξ M be an eigenvector of the nonlinear eigenvalue problem, P ( ξ) = µ ξ with some real and positive eigenvalue µ. Then there is an invariant cone for (10) If µ > 1, then the stationary solution 0 is unstable. If µ = 1, then the cone consists of periodic orbits. If µ < 1, then the stability of 0 depends on the stability of P with respect to the complementary directions. Remark (i) ξ is an eigenvector of P ξ ( ξ) with eigenvalue µ. If we assume that the remaining (n 2) eigenvalues λ 1,...,λ n 2 of P ξ ( ξ) satisfy λ j < min{1, µ}, (j = 1,..., n 2) then the invariant cone is attractive under the flow of (10).
13 Remark (ii); Interpretation Cone: 2-dim invariant surface. Smooth system: flat cone 5. Numerical Illustration 3D PWLS [+] system [ - ]system ẋ = λ + x w + y, : ẋ = (λ w )x 2w y + α(λ w mu )z ẏ = w + x + λ + y, ẏ = w x + (λ + w )y + αw z, ż = µ + z. ż = µ z
14 Figure 1: Attractive/ periodic / repulsive invariant cone for case ii (1), µ = 1.0, w + = w = 1, λ = , t + = π, t = 3π/4, δ = 0.0, α = 0.5. (a)λ + = , µ + = (b) λ + = , µ + = (c) λ + = , µ + = First Prev Next Last Go Back Full Screen Close Quit
15 6. Nonlinear perturbation of PWLS f ± = A ± ξ }{{} Basic P W LS + g ± (ξ) }{{} h.o.t. consider truncated version of (11); i.e for same r > 0 (11) g ± (ξ) 0, ( ξ r). Theorem 2 (Küpper,Hosham,Weiß (2009)). Reduction to invariant manifold Assume Ex. invariant cone for (PWLS), i.e. µ, ξ P ( ξ) = µ ξ. Ev of P ξ, µ, λ 1, λ 2,...λ n 2 λ j < min{1, µ}, (j = 1,..., n 2) Then Poincarè map P for (11) h : [0, ) R n 1, L := (y, h(y)) y 0, invariant cone under P h(y) = ξy + O(y 2 ) tangent to cone. First Prev Next Last Go Back Full Screen Close Quit
16 Remark Flow of (11) applied to L gives invariant cone-like manifold C. Attractive if λ j < min{1, µ}, j = 1,..., n 1. Bifurcation analysis reduced to C res. L Numerical example (Nonlinear system) [+] system [ - ]system ẋ = λ + x w + y, ẋ = (λ w )x 2w y ẏ = w + x + λ + y + υy 3, ẏ = w x + (λ + w )y + υy 3 ż = (µ + υ)y. ż = δw x δ(λ +w µ )y+(µ w δ)z.
17 Figure 2: Branch of instable periodic orbits for nonlinear example, µ = 1.0, w + = w = 1, λ = 1.0, δ = 0.5, α = 0.0. λ + = 1.0, µ + = 1.5. First Prev Next Last Go Back Full Screen Close Quit
18 Remark Existence of multiple cones Example: 3D PWLS [+] system [ - ]system ẋ = λ + x w + y, : ẋ = (λ w )x 2w y + α(λ w mu )z ẏ = w + x + λ + y, ẏ = w x + (λ + w )y + αw z, ż = µ + z. ż = µ z
19 Figure 3: 2 attractive invariant cones,λ + = 0.5, µ + = w + = w = 1, λ = 0.5, µ = 1.0, t + = π, δ = 0.0, α = 0.5. where t = π for a 1 and t = for a 2. First Prev Next Last Go Back Full Screen Close Quit
20 7. Six-dimensional non-smooth Brake-system A brake pad 1 on a rigid frame acts on a brake disc 2. Between brake pad and brake disc there is a relative displacement with constant velocity v, thus the frictional forces depend only on the normal force and the constant friction µ 1. The brake pad is equipped with three mechanical degrees of freedom modelled. Vertical movement x 1. Horizontal movement x 2. Rotation φ.
21 The Mathematical Model mẍ 1 + (d 1 + d 2 ) x 1 + b 2 (d1 d2) φ + (c 1 + c 2 )x 1 + b 2 (c 1 c 2 )φ = µ 2 sgn( x 1 a φ)c 3 x 2, (12) mẍ 2 (d 1 + d 2 )µ 1 x 1 µ 1b 2 (d1 d2) φ (c 1 + c 2 )µ 1 x 1 +c 3 x 2 + µ 1b 2 (c 1 c 2 )φ = 0, (13) j φ + ( b 2 (d 1 d 2 ) (d 1 + d 2 )hµ 1 ) x 1 + ( b2 4 (d1 + d2) + bhµ 1 2 (d 2 d 1 )) φ +( b 2 (c 1 c 2 ) (c 1 + c 2 )hµ 1 )x 1 c 3 sx 2 + ( b2 4 (c 1 + c 2 ) + bhµ 1 2 (c 2 c 1 ))φ = µ 2 sgn( x 1 a φ)c 3 ax 2, (14) First Prev Next Last Go Back Full Screen Close Quit
22 The Covering is connected with a brake holder through a velocity depending friction force R(v rel ) (v rel is the relative velocity between the contacting surfaces) acting on the contact surfaces. Initially this dependance is modelled by simple Coulomb friction characteristic. This means R(v rel ) = Nµ 2 sgn(v rel ), with v rel = x 1 a φ. If a more realistic form of the friction force with a typically nonlinear characteristic will be taken, nonlinear terms have to be added.
23 8. Simplification and reduction of the problem by setting some parameters and transformation The equations of motion (1-3), contain six unknown variables (x 1, x 1, x 2, x 2, φ, φ) and 13 parameters. It is clear, the exact analytic solutions are unavailable. That s why we fix some parameters to simplify the problem for a start. We take the following: (i)- c := c 1 = c 2, (ii)- d := d 1 = d 2. The Dynamical system (1-3) reduces to the following mẍ 1 + 2dx 1 + 2cx 1 = µ 2 sgn( x 1 a φ)c 3 x 2, (15) mẍ 2 2dµ 1 x 1 2cµ 1 x 1 + c 3 x 2 = 0, (16) j φ 2dhµ 1 x 1 + db2 2 φ 2chµ 1 x 1 c 3 sx 2 + cb2 2 φ = µ 2sgn( x 1 a φ)c 3 ax 2, (17) First Prev Next Last Go Back Full Screen Close Quit
24 After the transformation z 1 := x 1 z 2 := x 2 z 3 := x 1 aφ and the scaling z 4 := µ 1 x 1 z 5 := x 2 z 6 := x 2 a φ t maµ 1 t.
25 We Obtain ż = { A + z, z 6 > 0, A z, z 6 < 0, with the simple form of the matrices ã b 0 A ± = b c α 0 d 0 0 c γ 0 d 0 0 ẽ β g h 0 f (18) (19) where A ± are constant matrices containing various parameters, (Note that the elements of matrices are constant functions of eleven parameters i.e. ã = ma,...). The general structure is given by the observation that A + and A only differ in two entrances due to the simple (piecewise constant) form of the friction force.
26 Summary of analysis (α = 0) We can determine the eigenvalues and eigenvectors in simple form as: Eigenvalues Eigenvectors, λ 1,2 = l ± im, ν 12 = u 1 iw 1, λ 3,4 = ±in, ν 34 = u 2 iw 2, λ 5,6 = k ± ip, ν 56 = u 3 iw 3, where the λ 1,...,6 and ν 1,...,6 are to be represented in terms of 11 parameters in the Brake-system. So, we can write the general solution as: z = e lt (c 1 (cos(mt)u 1 + sin(mt)w 1 ) + c 2 ( sin(mt)u 1 + cos(mt)w 1 )) +c 3 (cos(nt)u 2 + sin(nt)w 2 ) + c 4 ( sin(nt)u 2 + cos(nt)w 2 ) + e kt (c 5 (cos(pt)u 3 + sin(pt)w 3 ) + c 6 ( sin(pt)u 3 + cos(pt)w 3 )), (20) where this solution corresponds to the system ż = A + z, z > 0 β = β + and ż = A z, z < 0 t = t with β = β. t = t + with
27 The constants c 1,...,6 are determined form the initial data z(0) = c 1 u 1 + c 2 w 1 + c 3 u 2 + c 4 w 2 + c 5 u 3 + c 6 w 3. The Poincarè map is of the form. c c 14 0 c P = P + P 21 c 22 0 c 24 c 25 = c 31 c 32 c 33 a 34 c 35 c c (21) c 51 c 52 0 c 54 c 55 Note that, c 11,..., c 55 are functions of the parameters of t + and t, so the Poincare map depends on the quantities t + and t in a non-linear way. Both times are determined by the following equations: H 1 (t + ) = {e T 6 z = 0, z > 0} := a 16 z 1 + a 26 z 2 + a 36 z 3 + a 46 z 4 + a 56 z 5 = 0,(22) H 2 (t ) = {e T 6 z = 0, z < 0} := b 16 z 1 + b 26 z 2 + b 36 z 3 + b 46 z 4 + b 56 z 5 = 0,(23) The equations (22) and (23) can only be solved numerically.
28 The eigenvalues of P as: µ 1 = c 33, µ 2,3 = 1 2 (c 11 + c 44 ± c c 14 c 41 2c 11 c 44 + c 2 44) µ 4,5 = 1 2 (c 22 + c 55 ± c c 25 c 52 2c 22 c 55 + c 2 55) Some Results In order to obtain an invariant cone consisting of periodic orbits, we solve equations (22) and (23) for µ = 1, we determine the eigenvector (corresponding to µ = 1) of the Poincare map (21), which defines the cone, Fig. (5); The periodic orbit is stable when all eigenvalues µ 1,...,5 1 In Fig. (6) and (7) we see the influence of the friction parameter for various values of µ 1 (µ 1 = 0.01 and 0.9).
29 Figure 4: sets of stable periodic orbits where all eigenvalues of P µ 1,...,5 1 First Prev Next Last Go Back Full Screen Close Quit
30 Figure 5: sets of stable periodic orbits where all eigenvalues of P µ 1,...,5 1 First Prev Next Last Go Back Full Screen Close Quit
31 Figure 6: sets of stable periodic orbits where all eigenvalues of P µ 1,...,5 1 First Prev Next Last Go Back Full Screen Close Quit
32 References [1] M. di Bernardo, C. Budd, A. R. Champneys, P. Kowalczyk, A. B. Nordmark, G. Olivar and P.T. Piiroinen, Bifurcations in Nonsmooth Dynamical Systems, SIAM Review, 50(4), (2008) [2] M. di Bernardo, C. Budd, A. R. Champneys, P. Kowalczyk, Piecewise-smooth Dynamical Systems: Theory and Applications (2008). [3] M. di Bernardo, C. J. Budd, and A. R. Champneys, Corner collision implies bordercollision bifurcation. Phys. D. 154 (3-4)(2001) [4] V. Carmona, E. Freire, E. Ponce, F. Torres, Bifurcation of invariant cones in piecwise linear homogeneous system. Int. J. Bifur. Chaos 15 (8) (2005) [5] V. Carmona, E. Freire, J. Ros, F. Torres, Limit cycle bifurcation in 3d continuous piecewise linear systems with two zones. Application to Chua s circuit, Int. J. Bifur. Chaos 15(10)(2005) [6] H. A. Hosham, T. Küpper and D. Weiß, Invariant manifold in nonsmooth systems, in preparation. [7] M. Kunze, Non-smooth Dynamical Systems, Springer-Verlag, Berlin, 2000 (Lecture Notes in Mathematics 1744). [8] T. Küpper, Invariant cones for non-smooth systems, Mathematics and Computers in Simulation, 79 (2008) First Prev Next Last Go Back Full Screen Close Quit
33 [9] T. Küpper, Concepts of non-smooth dynamical systems in: Mathematics and 21st Century, World Scientific Publishing Co, 2001, pp [10] T. Küpper and H. A. Hosham, Reduction to Invariant Cones for Non-smooth Systems (2009), It will appear in Special Issue of Mathematics and Computers in Simulation (2009). [11] T. Küpper, S. Moritz, Generalized Hopf bifurcation of nonsmooth planar systems. R. Soc. Lond. Philos, Trans. Ser. A Math. Phys. Eng. Sci. 359 (1789) (2001) [12] R. I. Leine, and H. Nijmeijer, Dynamics and Bifurcations of Nonsmooth Mechanical Systems, vol. 18 of Lecture Notes in Applied and Computational Mechanics. Springer- Verlag, Berlin, [13] H. E. Nusse, E. Ott and J. A. Yorke, Border-collision bifurcations: an explanation for observed difurcation phenomena. Phys. Rev. E(3) 49(2) (1994) [14] Y. Zou. T. Küpper, Generalized Hopf bifurcation for nonsmooth planar systems: the coner case, Notheast. Math. J. 17(4)(2001) [15] Y. Zou, T. Küpper, Hopf bifurcation for nonsmooth planar dynamical systems, Notheast. Math. J. 17(3)(2001) [16] Y. Zou, T. Küpper, Generalized Hopf bifurcation emanated from a corner, Nonlin. Anal. TAM 62 (1)(2005) [17] Y. Zou. T. Küpper, W. J. Beyn, Generalized Hopf bifurcation for planar Filippov systems continuous at the origin. J. Nonlin. Sci. 16(2)(2006) First Prev Next Last Go Back Full Screen Close Quit
Reduction to invariant cones for non-smooth systems
Available online at www.sciencedirect.com Mathematics and Computers in Simulation 8 () 98 995 Reduction to invariant cones for non-smooth systems T. Küpper, H.A. Hosham Mathematisches Institut, Universität
More informationDerivation of border-collision maps from limit cycle bifurcations
Derivation of border-collision maps from limit cycle bifurcations Alan Champneys Department of Engineering Mathematics, University of Bristol Mario di Bernardo, Chris Budd, Piotr Kowalczyk Gabor Licsko,...
More informationPIECEWISE LINEAR DIFFERENTIAL SYSTEMS WITH TWO REAL SADDLES. 1. Introduction
PIECEWISE LINEAR DIFFERENTIAL SYSTEMS WITH TWO REAL SADDLES JOAN C. ARTÉS1, JAUME LLIBRE 1, JOAO C. MEDRADO 2 AND MARCO A. TEIXEIRA 3 Abstract. In this paper we study piecewise linear differential systems
More informationModelling Research Group
Modelling Research Group The Importance of Noise in Dynamical Systems E. Staunton, P.T. Piiroinen November 20, 2015 Eoghan Staunton Modelling Research Group 2015 1 / 12 Introduction Historically mathematicians
More informationSzalai, R., & Osinga, H. M. (2007). Unstable manifolds of a limit cycle near grazing.
Szalai, R., & Osinga, H. M. (2007). Unstable manifolds of a limit cycle near grazing. Early version, also known as pre-print Link to publication record in Explore Bristol Research PDF-document University
More informationarxiv: v1 [nlin.cd] 7 Jun 2010
arxiv:006.236v [nlin.cd] 7 Jun 200 Simple models of bouncing ball dynamics and their comparison Andrzej Okninski, Bogus law Radziszewski 2 Physics Division, Department of Mechatronics and Mechanical Engineering
More informationBifurcations in Nonsmooth Dynamical Systems
SIAM REVIEW Vol. 5, No. 4, pp. 629 7 c 28 Society for Industrial and Applied Mathematics Bifurcations in Nonsmooth Dynamical Systems Mario di Bernardo Chris J. Budd Alan R. Champneys Piotr Kowalczyk Arne
More informationBIFURCATION PHENOMENA Lecture 4: Bifurcations in n-dimensional ODEs
BIFURCATION PHENOMENA Lecture 4: Bifurcations in n-dimensional ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits: equilibria cycles connecting orbits compact invariant manifolds strange
More informationAnalysis of the Takens-Bogdanov bifurcation on m parameterized vector fields
Analysis of the Takens-Bogdanov bifurcation on m parameterized vector fields Francisco Armando Carrillo Navarro, Fernando Verduzco G., Joaquín Delgado F. Programa de Doctorado en Ciencias (Matemáticas),
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 informationThe instantaneous local transition of a stable equilibrium to a chaotic attractor in piecewise-smooth systems of differential equations.
arxiv:1606.00073v1 [math.ds] 31 May 2016 The instantaneous local transition of a stable equilibrium to a chaotic attractor in piecewise-smooth systems of differential equations. D.J.W. Simpson Institute
More informationDynamics of Non-Smooth Systems
Dynamics of Non-Smooth Systems P Chandramouli Indian Institute of Technology Madras October 5, 2013 Mouli, IIT Madras Non-Smooth System Dynamics 1/39 Introduction Filippov Systems Behaviour around hyper-surface
More informationSliding Modes in Filippov Systems.
in Filippov Systems. School of Mathematics, Georgia Tech Joint work with Luciano Lopez (Univ. Bari) Geneva Conference In honor of Ernst Hairer 60th birthday References Filippov [DEs with discontinuous
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 informationME8230 Nonlinear Dynamics
ME8230 Nonlinear Dynamics Lecture 1, part 1 Introduction, some basic math background, and some random examples Prof. Manoj Srinivasan Mechanical and Aerospace Engineering srinivasan.88@osu.edu Spring mass
More informationNonsmooth systems: synchronization, sliding and other open problems
John Hogan Bristol Centre for Applied Nonlinear Mathematics, University of Bristol, England Nonsmooth systems: synchronization, sliding and other open problems 2 Nonsmooth Systems 3 What is a nonsmooth
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre
More informationDynamical behaviour of a controlled vibro-impact system
Vol 17 No 7, July 2008 c 2008 Chin. Phys. Soc. 1674-1056/2008/17(07)/2446-05 Chinese Physics B and IOP Publishing Ltd Dynamical behaviour of a controlled vibro-impact system Wang Liang( ), Xu Wei( ), and
More informationThe Two-Fold Singularity of Discontinuous Vector Fields
SIAM J. APPLIED DYNAMICAL SYSTEMS Vol. 8, No. 2, pp. 624 640 c 2009 Society for Industrial and Applied Mathematics The Two-Fold Singularity of Discontinuous Vector Fields M. R. Jeffrey and A. Colombo Abstract.
More informationControlling Attractivity of Friction-Induced Equilibrium Sets
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 2005 Seville, Spain, December 12-15, 2005 TuA16.5 Controlling Attractivity of Friction-Induced Equilibrium
More information8.385 MIT (Rosales) Hopf Bifurcations. 2 Contents Hopf bifurcation for second order scalar equations. 3. Reduction of general phase plane case to seco
8.385 MIT Hopf Bifurcations. Rodolfo R. Rosales Department of Mathematics Massachusetts Institute of Technology Cambridge, Massachusetts MA 239 September 25, 999 Abstract In two dimensions a Hopf bifurcation
More informationOn rigorous integration of piece-wise linear systems
On rigorous integration of piece-wise linear systems Zbigniew Galias Department of Electrical Engineering AGH University of Science and Technology 06.2011, Trieste, Italy On rigorous integration of PWL
More informationDynamical Analysis of a Hydraulic Pressure Relief Valve
Proceedings of the World Congress on Engineering 9 Vol II WCE 9, July -, 9, London, U.K. Dynamical Analysis of a Hydraulic Pressure Relief Valve Gábor Licskó Alan Champneys Csaba Hős Abstract A mathematical
More informationCommunications in Nonlinear Science and Numerical Simulation
Accepted Manuscript Periodic motions generated from non-autonomous grazing dynamics M.U. Akhmet, A. Kıvılcım PII: S17-574(17)342-4 DOI: 1.116/j.cnsns.217.2.2 Reference: CNSNS 419 To appear in: Communications
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 informationEntrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.
Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the
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 informationA conjecture on sustained oscillations for a closed-loop heat equation
A conjecture on sustained oscillations for a closed-loop heat equation C.I. Byrnes, D.S. Gilliam Abstract The conjecture in this paper represents an initial step aimed toward understanding and shaping
More informationHopf Bifurcation and Limit Cycle Analysis of the Rikitake System
ISSN 749-3889 (print), 749-3897 (online) International Journal of Nonlinear Science Vol.4(0) No.,pp.-5 Hopf Bifurcation and Limit Cycle Analysis of the Rikitake System Xuedi Wang, Tianyu Yang, Wei Xu Nonlinear
More informationLecture 5. Numerical continuation of connecting orbits of iterated maps and ODEs. Yu.A. Kuznetsov (Utrecht University, NL)
Lecture 5 Numerical continuation of connecting orbits of iterated maps and ODEs Yu.A. Kuznetsov (Utrecht University, NL) May 26, 2009 1 Contents 1. Point-to-point connections. 2. Continuation of homoclinic
More informationNBA Lecture 1. Simplest bifurcations in n-dimensional ODEs. Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011
NBA Lecture 1 Simplest bifurcations in n-dimensional ODEs Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011 Contents 1. Solutions and orbits: equilibria cycles connecting orbits other invariant sets
More informationStability analysis of sliding-grazing phenomenon in dry-friction oscillator using TS fuzzy approach
Stability analysis of sliding-grazing phenomenon in dry-friction oscillator using TS fuzzy approach Kamyar Mehran School of Engineering University of Warwick Coventry, CV4 7LW, UK Email: k.mehran@warwick.ac.uk
More informationContinuation of cycle-to-cycle connections in 3D ODEs
HET p. 1/2 Continuation of cycle-to-cycle connections in 3D ODEs Yuri A. Kuznetsov joint work with E.J. Doedel, B.W. Kooi, and G.A.K. van Voorn HET p. 2/2 Contents Previous works Truncated BVP s with projection
More informationCenter Bifurcation of a point on the Poincaré Equator
ITERATION THEORY (ECIT 08) A.N. Sharkovsky, I.M. Sushko (Eds.) Grazer Math. Ber., ISSN 1016 7692 Bericht Nr. 354 (2009), 254-279 Center Bifurcation of a point on the Poincaré Equator Iryna Sushko and Laura
More informationVERSAL DEFORMATIONS OF BILINEAR SYSTEMS UNDER OUTPUT-INJECTION EQUIVALENCE
PHYSCON 2013 San Luis Potosí México August 26 August 29 2013 VERSAL DEFORMATIONS OF BILINEAR SYSTEMS UNDER OUTPUT-INJECTION EQUIVALENCE M Isabel García-Planas Departamento de Matemàtica Aplicada I Universitat
More informationSIMPLE CHAOTIC FLOWS WITH ONE STABLE EQUILIBRIUM
International Journal of Bifurcation and Chaos, Vol. 23, No. 11 (2013) 1350188 (7 pages) c World Scientific Publishing Company DOI: 10.1142/S0218127413501885 SIMPLE CHAOTIC FLOWS WITH ONE STABLE EQUILIBRIUM
More information1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system:
1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system: ẋ = y x 2, ẏ = z + xy, ż = y z + x 2 xy + y 2 + z 2 x 4. (ii) Determine
More informationThe problem of singularity in impacting systems
The problem of singularity in impacting systems Soumitro Banerjee, Department of Physics Indian Institute of Science Education & Research, Kolkata, India The problem of singularity in impacting systems
More informationLectures on Dynamical Systems. Anatoly Neishtadt
Lectures on Dynamical Systems Anatoly Neishtadt Lectures for Mathematics Access Grid Instruction and Collaboration (MAGIC) consortium, Loughborough University, 2007 Part 3 LECTURE 14 NORMAL FORMS Resonances
More informationOn the Takens-Bogdanov Bifurcation in the Chua s Equation
1722 PAPER Special Section on Nonlinear Theory and Its Applications On the Takens-Bogdanov Bifurcation in the Chua s Equation Antonio ALGABA, Emilio FREIRE, Estanislao GAMERO, and Alejandro J. RODRÍGUEZ-LUIS,
More informationComputational Methods in Dynamical Systems and Advanced Examples
and Advanced Examples Obverse and reverse of the same coin [head and tails] Jorge Galán Vioque and Emilio Freire Macías Universidad de Sevilla July 2015 Outline Lecture 1. Simulation vs Continuation. How
More informationChaotifying 2-D piecewise linear maps via a piecewise linear controller function
Chaotifying 2-D piecewise linear maps via a piecewise linear controller function Zeraoulia Elhadj 1,J.C.Sprott 2 1 Department of Mathematics, University of Tébéssa, (12000), Algeria. E-mail: zeraoulia@mail.univ-tebessa.dz
More informationDRIVEN and COUPLED OSCILLATORS. I Parametric forcing The pendulum in 1:2 resonance Santiago de Compostela
DRIVEN and COUPLED OSCILLATORS I Parametric forcing The pendulum in 1:2 resonance Santiago de Compostela II Coupled oscillators Resonance tongues Huygens s synchronisation III Coupled cell system with
More informationDESYNCHRONIZATION TRANSITIONS IN RINGS OF COUPLED CHAOTIC OSCILLATORS
Letters International Journal of Bifurcation and Chaos, Vol. 8, No. 8 (1998) 1733 1738 c World Scientific Publishing Company DESYNCHRONIZATION TRANSITIONS IN RINGS OF COUPLED CHAOTIC OSCILLATORS I. P.
More informationPeriodic solutions of the perturbed symmetric Euler top
Periodic solutions of the perturbed symmetric Euler top Universitatea Babeş-Bolyai (Cluj-Napoca, Romania) abuica@math.ubbcluj.ro and Universitat de Lleida (Lleida, Spain) Plan of the talk 1 Our problem
More informationAn Introduction to Numerical Continuation Methods. with Application to some Problems from Physics. Eusebius Doedel. Cuzco, Peru, May 2013
An Introduction to Numerical Continuation Methods with Application to some Problems from Physics Eusebius Doedel Cuzco, Peru, May 2013 Persistence of Solutions Newton s method for solving a nonlinear equation
More informationHOPF BIFURCATION FOR A 3D FILIPPOV SYSTEM
Dynamics of Continuous, Discrete and Impulsive Systems Series A: Mathematical Analysis 16 (2009) 759-775 Copyright c 2009 Watam Press http://www.watam.org HOPF BIFURCATION FOR A 3D FILIPPOV SYSTEM M. U.
More informationConnecting orbits and invariant manifolds in the spatial three-body problem
C C Dynamical A L T E C S H Connecting orbits and invariant manifolds in the spatial three-body problem Shane D. Ross Control and Dynamical Systems, Caltech Work with G. Gómez, W. Koon, M. Lo, J. Marsden,
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 informationApril 13, We now extend the structure of the horseshoe to more general kinds of invariant. x (v) λ n v.
April 3, 005 - Hyperbolic Sets We now extend the structure of the horseshoe to more general kinds of invariant sets. Let r, and let f D r (M) where M is a Riemannian manifold. A compact f invariant set
More information27. Topological classification of complex linear foliations
27. Topological classification of complex linear foliations 545 H. Find the expression of the corresponding element [Γ ε ] H 1 (L ε, Z) through [Γ 1 ε], [Γ 2 ε], [δ ε ]. Problem 26.24. Prove that for any
More informationCo-dimension-two Grazing Bifurcations in Single-degree-of-freedom Impact Oscillators
Co-dimension-two razing Bifurcations in Single-degree-of-freedom Impact Oscillators Phanikrishna Thota a Xiaopeng Zhao b and Harry Dankowicz c a Department of Engineering Science and Mechanics MC 219 Virginia
More informationHomoclinic bifurcations in Chua s circuit
Physica A 262 (1999) 144 152 Homoclinic bifurcations in Chua s circuit Sandra Kahan, Anibal C. Sicardi-Schino Instituto de Fsica, Universidad de la Republica, C.C. 30, C.P. 11 000, Montevideo, Uruguay
More informationParameter-Dependent Eigencomputations and MEMS Applications
Parameter-Dependent Eigencomputations and MEMS Applications David Bindel UC Berkeley, CS Division Parameter-Dependent Eigencomputations and MEMS Applications p.1/36 Outline Some applications Mathematical
More informationTHREE SCENARIOS FOR CHANGING OF STABILITY IN THE DYNAMIC MODEL OF NERVE CONDUCTION
THREE SCENARIOS FOR CHANGING OF STABILITY IN THE DYNAMIC MODEL OF NERVE CONDUCTION E.A. Shchepakina Samara National Research University, Samara, Russia Abstract. The paper deals with the specific cases
More informationBifurcation of Fixed Points
Bifurcation of Fixed Points CDS140B Lecturer: Wang Sang Koon Winter, 2005 1 Introduction ẏ = g(y, λ). where y R n, λ R p. Suppose it has a fixed point at (y 0, λ 0 ), i.e., g(y 0, λ 0 ) = 0. Two Questions:
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationConnecting Orbits with Bifurcating End Points
Connecting Orbits with Bifurcating End Points Thorsten Hüls Fakultät für Mathematik Universität Bielefeld Postfach 100131, 33501 Bielefeld Germany huels@mathematik.uni-bielefeld.de June 24, 2004 Abstract
More informationThe Stick-Slip Vibration and Bifurcation of a Vibro-Impact System with Dry Friction
Send Orders for Reprints to reprints@benthamscience.ae 308 The Open Mechanical Engineering Journal, 2014, 8, 308-313 Open Access The Stick-Slip Vibration and Bifurcation of a Vibro-Impact System with Dry
More information2D-Volterra-Lotka Modeling For 2 Species
Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose
More informationAN ELECTRIC circuit containing a switch controlled by
878 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING, VOL. 46, NO. 7, JULY 1999 Bifurcation of Switched Nonlinear Dynamical Systems Takuji Kousaka, Member, IEEE, Tetsushi
More information2 Lecture 2: Amplitude equations and Hopf bifurcations
Lecture : Amplitude equations and Hopf bifurcations This lecture completes the brief discussion of steady-state bifurcations by discussing vector fields that describe the dynamics near a bifurcation. From
More informationBoundary value problem with integral condition for a Blasius type equation
ISSN 1392-5113 Nonlinear Analysis: Modelling and Control, 216, Vol. 21, No. 1, 114 12 http://dx.doi.org/1.15388/na.216.1.8 Boundary value problem with integral condition for a Blasius type equation Sergey
More information7 Two-dimensional bifurcations
7 Two-dimensional bifurcations As in one-dimensional systems: fixed points may be created, destroyed, or change stability as parameters are varied (change of topological equivalence ). In addition closed
More informationPersonal notes on renormalization
Personal notes on renormalization Pau Rabassa Sans April 17, 2009 1 Dynamics of the logistic map This will be a fast (and selective) review of the dynamics of the logistic map. Let us consider the logistic
More informationThe Infinitesimal Phase Response Curves of Oscillators in Piecewise Smooth Dynamical Systems
arxiv:1603.03503v1 [math.ds] 11 Mar 2016 The Infinitesimal Phase Response Curves of Oscillators in Piecewise Smooth Dynamical Systems Youngmin Park Department of Mathematics, Applied Mathematics, and Statistics.
More information(8.51) ẋ = A(λ)x + F(x, λ), where λ lr, the matrix A(λ) and function F(x, λ) are C k -functions with k 1,
2.8.7. Poincaré-Andronov-Hopf Bifurcation. In the previous section, we have given a rather detailed method for determining the periodic orbits of a two dimensional system which is the perturbation of a
More informationComplex Dynamic Systems: Qualitative vs Quantitative analysis
Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems Chiara Mocenni Department of Information Engineering and Mathematics University of Siena (mocenni@diism.unisi.it) Dynamic
More informationThe Conley Index and Rigorous Numerics of Attracting Periodic Orbits
The Conley Index and Rigorous Numerics of Attracting Periodic Orbits Marian Mrozek Pawe l Pilarczyk Conference on Variational and Topological Methods in the Study of Nonlinear Phenomena (Pisa, 2000) 1
More informationSWITCHED dynamical systems are useful in a variety of engineering
1184 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS I: REGULAR PAPERS, VOL. 51, NO. 6, JUNE 2004 Bifurcation Analysis of Switched Dynamical Systems With Periodically Moving Borders Yue Ma, Hiroshi Kawakami,
More informationOn the Dynamics of a n-d Piecewise Linear Map
EJTP 4, No. 14 2007 1 8 Electronic Journal of Theoretical Physics On the Dynamics of a n-d Piecewise Linear Map Zeraoulia Elhadj Department of Mathematics, University of Tébéssa, 12000, Algeria. Received
More informationApplication demonstration. BifTools. Maple Package for Bifurcation Analysis in Dynamical Systems
Application demonstration BifTools Maple Package for Bifurcation Analysis in Dynamical Systems Introduction Milen Borisov, Neli Dimitrova Department of Biomathematics Institute of Mathematics and Informatics
More informationHomoclinic trajectories of non-autonomous maps
Homoclinic trajectories of non-autonomous maps Thorsten Hüls Department of Mathematics Bielefeld University Germany huels@math.uni-bielefeld.de www.math.uni-bielefeld.de:~/huels Non-autonomous dynamical
More informationON THE ROBUSTNESS OF PERIODIC SOLUTIONS IN RELAY FEEDBACK SYSTEMS. Mario di Bernardo Λ Karl Henrik Johansson ΛΛ;1 Ulf Jönsson ΛΛΛ Francesco Vasca Λ
Copyright IFAC 1th Triennial World Congress, Barcelona, Spain ON THE ROBUSTNESS OF PERIODIC SOLUTIONS IN RELAY FEEDBACK SYSTEMS Mario di Bernardo Λ Karl Henrik Johansson ΛΛ;1 Ulf Jönsson ΛΛΛ Francesco
More informationarxiv: v2 [math.ds] 8 Feb 2016
Canard explosion and relaxation oscillation in planar, piecewise-smooth, continuous systems Andrew Roberts Department of Mathematics, Cornell University (Dated: February 9, 2016) arxiv:1412.0263v2 [math.ds]
More informationDynamical Systems with Applications
Stephen Lynch Dynamical Systems with Applications using MATLAB Birkhauser Boston Basel Berlin Preface xi 0 A Tutorial Introduction to MATLAB and the Symbolic Math Toolbox 1 0.1 Tutorial One: The Basics
More informationSlow Manifold of a Neuronal Bursting Model
Slow Manifold of a Neuronal Bursting Model Jean-Marc Ginoux 1 and Bruno Rossetto 2 1 PROTEE Laboratory, Université du Sud, B.P. 2132, 83957, La Garde Cedex, France, ginoux@univ-tln.fr 2 PROTEE Laboratory,
More informationTWO DIMENSIONAL FLOWS. Lecture 5: Limit Cycles and Bifurcations
TWO DIMENSIONAL FLOWS Lecture 5: Limit Cycles and Bifurcations 5. Limit cycles A limit cycle is an isolated closed trajectory [ isolated means that neighbouring trajectories are not closed] Fig. 5.1.1
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 informationCONTROLLING CHAOTIC DYNAMICS USING BACKSTEPPING DESIGN WITH APPLICATION TO THE LORENZ SYSTEM AND CHUA S CIRCUIT
Letters International Journal of Bifurcation and Chaos, Vol. 9, No. 7 (1999) 1425 1434 c World Scientific Publishing Company CONTROLLING CHAOTIC DYNAMICS USING BACKSTEPPING DESIGN WITH APPLICATION TO THE
More informationA Model of Evolutionary Dynamics with Quasiperiodic Forcing
paper presented at Society for Experimental Mechanics (SEM) IMAC XXXIII Conference on Structural Dynamics February 2-5, 205, Orlando FL A Model of Evolutionary Dynamics with Quasiperiodic Forcing Elizabeth
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 informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 1. Linear systems Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Linear systems 1.1 Differential Equations 1.2 Linear flows 1.3 Linear maps
More informationSelf-Excited Vibration
Wenjing Ding Self-Excited Vibration Theory, Paradigms, and Research Methods With 228 figures Ö Springer Contents Chapter 1 Introduction 1 1.1 Main Features of Self-Excited Vibration 1 1.1.1 Natural Vibration
More informationChaotic behavior analysis based on corner bifurcations
Author manuscript, published in "Nonlinear Analysis: Hybrid Systems 3 (29) 543-55" Chaotic behavior analysis based on corner bifurcations D. Benmerzouk a,, J.P. Barbot b, a University of Tlemcen,Department
More informationAttractors near grazing-sliding bifurcations. Glendinning, Paul and Kowalczyk, Piotr and Nordmark, Arne. MIMS EPrint:
Attractors near grazing-sliding bifurcations Glendinning, Paul and Kowalczyk, Piotr and Nordmark, Arne 11 MIMS EPrint: 11.63 Manchester Institute for Mathematical Sciences School of Mathematics The University
More informationSt. Petersburg Math. J. 13 (2001),.1 Vol. 13 (2002), No. 4
St Petersburg Math J 13 (21), 1 Vol 13 (22), No 4 BROCKETT S PROBLEM IN THE THEORY OF STABILITY OF LINEAR DIFFERENTIAL EQUATIONS G A Leonov Abstract Algorithms for nonstationary linear stabilization are
More informationAMADEU DELSHAMS AND RAFAEL RAMíREZ-ROS
POINCARÉ-MELNIKOV-ARNOLD METHOD FOR TWIST MAPS AMADEU DELSHAMS AND RAFAEL RAMíREZ-ROS 1. Introduction A general theory for perturbations of an integrable planar map with a separatrix to a hyperbolic fixed
More informationStability of Direct Spring Operated Pressure Relief Valves from CFD to spreadsheet
Stability of Direct Spring Operated Pressure Relief Valves from CFD to spreadsheet 12th September 217 HDR ... with major contributions by Prof. Alan Champneys (University of Bristol, Dept. of Eng. Mathematics)
More informationA BENDIXSON DULAC THEOREM FOR SOME PIECEWISE SYSTEMS
A BENDIXSON DULAC THEOREM FOR SOME PIECEWISE SYSTEMS LEONARDO P. C. DA CRUZ AND JOAN TORREGROSA Abstract. The Bendixson Dulac Theorem provides a criterion to find upper bounds for the number of limit cycles
More informationInvariant manifolds of the Bonhoeffer-van der Pol oscillator
Invariant manifolds of the Bonhoeffer-van 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 informationAnti-synchronization of a new hyperchaotic system via small-gain theorem
Anti-synchronization of a new hyperchaotic system via small-gain theorem Xiao Jian( ) College of Mathematics and Statistics, Chongqing University, Chongqing 400044, China (Received 8 February 2010; revised
More informationBifurcations of phase portraits of pendulum with vibrating suspension point
Bifurcations of phase portraits of pendulum with vibrating suspension point arxiv:1605.09448v [math.ds] 9 Sep 016 A.I. Neishtadt 1,,, K. Sheng 1 1 Loughborough University, Loughborough, LE11 3TU, UK Space
More informationPeriodic Orbits and Transport: Some Interesting Dynamics in the Three-Body Problem
Periodic Orbits and Transport: Some Interesting Dynamics in the Three-Body Problem Shane Ross Martin Lo (JPL), Wang Sang Koon and Jerrold Marsden (Caltech) CDS 280, January 8, 2001 shane@cds.caltech.edu
More informationSTABILITY OF 2D SWITCHED LINEAR SYSTEMS WITH CONIC SWITCHING. 1. Introduction
STABILITY OF D SWITCHED LINEAR SYSTEMS WITH CONIC SWITCHING H LENS, M E BROUCKE, AND A ARAPOSTATHIS 1 Introduction The obective of this paper is to explore necessary and sufficient conditions for stability
More informationDynamical Systems Analysis Using Differential Geometry
Dynamical Systems Analysis Using Differential Geometry Jean-Marc GINOUX and Bruno ROSSETTO Laboratoire P.R.O.T.E.E., équipe EBMA / ISO Université du Sud Toulon-Var, B.P. 3, 83957, La Garde, ginoux@univ-tln.fr,
More informationCHARACTERIZATION OF SADDLE-NODE EQUILIBRIUM POINTS ON THE STABILITY BOUNDARY OF NONLINEAR AUTONOMOUS DYNAMICAL SYSTEM
CHARACTERIZATION OF SADDLE-NODE EQUILIBRIUM POINTS ON THE STABILITY BOUNDARY OF NONLINEAR AUTONOMOUS DYNAMICAL SYSTEM Fabíolo Moraes Amaral, Josaphat Ricardo Ribeiro Gouveia Júnior, Luís Fernando Costa
More informationStudy of Chaos and Dynamics of DC-DC Converters BY SAI RAKSHIT VINNAKOTA ANUROOP KAKKIRALA VIVEK PRAYAKARAO
Study of Chaos and Dynamics of DC-DC Converters BY SAI RAKSHIT VINNAKOTA ANUROOP KAKKIRALA VIVEK PRAYAKARAO What are DC-DC Converters?? A DC-to-DC converter is an electronic circuit which converts a source
More informationDynamical Systems with Applications using Mathematica
Stephen Lynch Dynamical Systems with Applications using Mathematica Birkhäuser Boston Basel Berlin Contents Preface xi 0 A Tutorial Introduction to Mathematica 1 0.1 A Quick Tour of Mathematica 2 0.2 Tutorial
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 information