


Chaotic Transitions in Deterministic and Stochastic Dynamical Systems: 
This file is also available in Adobe Acrobat PDF format Chapter One Introduction This work has two main objectives: (1) to present the Melnikov method as a unified theoretical framework for the study of transitions and chaos in a wide class of deterministic and stochastic nonlinear planar systems, and (2) to demonstrate the method's usefulness in applications, particularly for stochastic systems. Our interest in the Melnikov method is motivated by its capability to provide criteria and information on the occurrence of transitions and chaotic behavior in a wide variety of systems in engineering, physics, and the life sciences. To illustrate the type of problem to which the Melnikov method is applicable we consider a celebrated experiment on a system known as the magnetoelastic beam. The experiment demonstrates the remarkable type of dynamic behavior called deterministic chaos (Moon and Holmes, 1979). The system consists of (a) a rigid frame fixed onto a shaking table that may undergo periodic horizontal motions, (b) a beam with a vertical undeformed axis, an upper end fixed onto the frame, and a free lower end, and (c) two identical magnets equidistant from the undeformed position of the beam (Fig. 1.1). The beam experiences nonlinear displacementdependent forces induced by the magnets, linear restoring forces due to its elasticity, dissipative forces due to its internal friction, the viscosity of the surrounding air, and magnetic damping, and periodic excitation forces due to the horizontal motion of the shaking table. Neither the system properties nor the forces acting on the beam vary randomly with time: the system is fully deterministic. In the absence of excitation, and depending upon the initial conditions, the beam settles on one of two possible stable equilibria, that is, with the beam's tip closer to the right magnet or closer to the left magnet. The beam also has an unstable equilibrium positionits vertical undeformed axis. If the excitation is periodic, three distinct types of steadystate dynamic behavior can occur:
The system just described may be modeled as a dynamical systema system that evolves in time in accordance with a specified mathematical expression. In this book we are primarily concerned with dynamical systems capable of exhibiting all three types of behavior illustrated in Fig. 1.2. One basic feature of such systems is that they are multistable, meaning that their unforced counterparts have at least two stable equilibria (the term applied to the case of two stable equilibria is bistable). In the particular case of mechanical systems, the dynamic behavior is modeled by nonlinear differential equations expressing a relationship among terms that represent
Similar terms occur in equations modeling other types of dynamical system, for example, electrical, thermal, or chemical systems. For a large number of systems arising in engineering or physics safe operation requires that steadystate motions occur within a restricted region, called a safe region (in Fig. 1.2(a), the displacement coordinates in the restricted regions are bounded by the vertical line that coincides with the axis of the undeformed beam); transitions to motions visiting another region are undesirable. However, for some systems (e.g., systems that enhance heat transfer, and neurological systems whose activity entails escapes or, in neurological terminology, firings) the occurrence of such transitions is a functional requirement. Although we will also examine systems with a slowly varying third variable, our main focus will be on planar systems, that is, continuous systems with two timedependent variables, for example, displacement and velocity. For planar systems subjected to periodic excitation an analytical condition that guarantees the nonoccurrence of transitions was derived in a seminal paper by Melnikov (1963). That condition involves a functionthe Melnikov functionconsisting of a sum of terms related, through the system's potential, to the system dissipation and excitation. The Melnikov condition for nonoccurrence of transitions states that if the Melnikov function has no zeros or at most a double zero, then transitions cannot occur. The counterpart of that condition, where the Melnikov function has simple zeros, is referred to as the Melnikov necessary condition for the occurrence of transitions. Melnikov theory requires that the perturbationthat is, the excitation and dissipationbe sufficiently small. Numerical calculations show, however, that the theory is useful for perturbations large enough to be of interest in a wide variety of applications. Melnikov was motivated by problems arising in nuclear physics and engineering (e.g., the theory of particle accelerators, controlled thermonuclear reactions). However, his condition has subsequently found application in a wide variety of disciplines. Examples of applications, some of which are reviewed in subsequent chapters, include
The Melnikov method consists of the application, in various forms, of the necessary Melnikov condition for the occurrence of transitions. So far we have mentioned Melnikov's original method which, as indicated earlier, was restricted to periodically excited planar systems. The extension of the Melnikov method to planar systems with nonperiodic excitation is of interest. For deterministic systems with quasiperiodic excitation (i.e., excitation consisting of a sum of harmonics with, in general, incommensurate frequencies) the extension was developed by Wiggins (1988). The extension to stochastic systems was introduced by Frey and Simiu (1993). It is referred to as the stochastic Melnikov methodMelnikov's method applied to stochastic dynamical systems. The stochastic counterpart of the Melnikov function is referred to as a Melnikov process. For systems whose stochastic excitations are bounded the stochastic Melnikov method is a straightforward extension of the original Melnikov method. Systems excited by dichotomous noise, which are of interest in electrical engineering and chemical physics, are one example. However, for systems whose stochastic excitations are modeled by unbounded processes the Melnikov method cannot be used to guarantee the nonoccurrence of escapes. Rather, it can be used to calculate lower bounds for probabilities that transitions will not occur during specified time intervals, or upper bounds for probabilities that motion with transitions can occur. Melnikov (1963) had noted that, depending upon whether it has or does not have simple zeros, the function he developed can be associated with distinct types of dynamic behavior. However, he obtained his original results too early in the development of the theory of chaotic dynamics to describe one of those types of behavior as deterministic chaos. It has been known for several decades that, for deterministic planar multistable systems, the condition that the Melnikov function have simple zeros is necessary for the occurrence not only of transitions but of chaotic behavior as well (see Chapter 3). In fact, for such systems, motions with transitions occur if and only if the motion is chaotic. However, that motions with transitions occurring in stochastic planar multistable systems are also chaotic has only been established relatively recently (see Chapter 5). The following feature of the Melnikov method is especially useful in applications. Through the Melnikov scale factora function of frequency that depends upon the system's potentialthe Melnikov method provides useful information on the relative degree to which the various frequency components of the excitation are effective in promoting transitions. We discuss in this book several applications of this feature, including
Following this Introduction, the book is divided into two parts. Part 1, containing Chapters 2 through 5 is devoted to fundamentals. Part 2, containing Chapters 6 through 12, includes applications. In Chapter 2 we introduce basic elements of the theory of nonlinear dynamical systems needed for the development of the Melnikov method. We then present the derivation of the Melnikov function and the Melnikov necessary condition for the occurrence of transitions. For a particular type of system we interpret this condition from an energy viewpoint. We discuss two useful tools that simplify the study and facilitate the diagnosis of motions in continuous dynamical systems: the Poincaré section, which is applicable to periodically excited systems, and the phase space slice, which is the counterpart of the Poincaré section for quasiperiodically excited systems. We introduce the notion of transport associated with transitions, and a useful measure of the transport, the Melnikov functional known as the phase space flux factor. Finally, for the case of slowly varying planar systems, we present the derivation of the Melnikov function and the associated necessary condition for the occurrence of transitions, and note a shortcoming of the Melnikov approach with respect to the case of planar systems with no slowly varying parameter. In Chapter 3 we establish that motions with transitionsmotions for whose occurrence it is necessary that the Melnikov function have simple zerosare chaotic. As a first step we provide introductory definitions and illustrations of Cantor sets and fractal dimensions. We then study Smale's famous horseshoe map and the shift map, and show the equivalence of the dynamics induced by these two maps. We introduce simple symbolic dynamics techniques and apply them to the shift map to prove that its behavior is chaotic. This establishes that the equivalent dynamics induced by the Smale horseshoe map is chaotic as well. We define chaos mathematically and show that its properties include sensitivity to initial condition and fractal dimensions of geometrical structures associated with the motion, a property that is readily identifiable in certain types of numerical simulations and can therefore serve as a helpful indicator of deterministic chaos. We then state the SmaleBirkhoff theorem, which establishes that the dynamics of motions with transitions in planar multistable systems is chaotic by showing that it is equivalent to the chaotic dynamics induced by the Smale horseshoe map. We provide illustrative numerical examples and examples of chaotic motions of an experimental system, the buckled Stoker column. In Chapter 4 we present simple elements of the theory of stochastic processes needed for the subsequent development of the stochastic Melnikov method. In particular, we define stationary stochastic processes, autocovariance and crosscovariance functions, spectral densities, linear filters, convolutions, and transfer functions. We also derive the relation between the spectral density of the output and the spectral density of the input of a linear filter with known transfer function. (A brief review of basic elements of probability theory is presented in Appendix A5.) In Chapter 5 we develop the stochastic Melnikov method. We show that, like deterministic systems, multistable planar systems with stochastic excitation can exhibit motions with transitions that are in practice indistinguishable from deterministic chaotic motions and, like the latter, are sensitive to initial conditions. We refer to such motions as chaotic stochastic motions, as distinct from stochastic motions that do not exhibit properties associated with chaos. The applicability of the Melnikov method to stochastic systems is demonstrated by approximating the stochastic excitation processes by ensembles of sums of periodic terms, thus reducing the core of the stochastic Melnikov method to the problem, studied in Chapters 3 and 4, of the Melnikov method as applied to systems with periodic or quasiperiodic excitation. By using the stochastic Melnikov method we develop criteria that guarantee the nonoccurrence of transitions in systems excited by dichotomous noise, and simple assessments of the influence of the excitation's spectral shape on the system's escape rate. We now list the applications of the Melnikov method contained in Part 2. Chapter 6 describes an application to a topic in naval architecture: the capsizing problem for a rolling vessel under the action of wave forces. Chapter 7 describes the application of the Melnikov method to openloop control. Chapter 8 uses the Melnikov method to interpret the stochastic resonance phenomenon in chaotic dynamics terms and presents useful results that follow from this interpretation. Chapter 9 describes the use of the Melnikov method for determining the cutoff point for excitation by experimental noise in a onedimensional nonlinear system. Chapter 10 presents a study of a mechanical engineering problem: the snapthrough of a buckled column excited by a transverse stochastic loading. Chapter 11 describes an application to a topic in oceanography: the problem of transitions for windinduced alongshore currents over corrugated bottom topography. Chapter 12 describes an application of the Melnikov method in neurophysiology: the modeling of the behavior of the auditory nerve fiber as a chaotic dynamical system.
File created: 8/7/2007 Questions and comments to: webmaster@pupress.princeton.edu 