Application of Stability Theory in Study of Local Dynamics of Nonlinear Systems

Investigating local dynamics of equilibrium points of nonlinear systems plays an important role in studying the behavior of dynamical systems. There are many different definitions for stable and unstable solutions in the literature. The main goal to develop stability definitions is exploring the responses or output of a system to perturbation as time approaches infinity. Due to the wide range of application of local dynamical system theory in physics, biology, economics and social science, it still attracts many researchers to play with its definitions to find out the answers for their questions. In this paper, we start with a brief review over continuous time dynamical systems modeling and then we bring useful examples to the playground. We study the local dynamics of some interesting systems and we show the local stable behavior of the system around its critical points. Moreover, we look at local dynamical behavior of famous dynamical systems, Hénon-Heiles system, Duffing oscillator and Van der Pol equation and analyze them. Finally, we discuss about the chaotic behavior of Hamiltonian systems using two different and new examples.


Introduction
A dynamical system describes the evolution of a system over time using a set of mathematical laws. Also, it can be used to predict the interactions between different components of a system [1] [2]. There are two main methods to model the dynamical behaviors of a system, continuous time modeling, discrete-time mod-In this paper we will be concerned with continuous dynamical systems which are defined by differential equations. Indeed, some famous examples of dynamical systems can be written in terms of differential equations: the harmonic oscillator, the pendulum and double pendulum, or the N-body problem [4]- [9]. A dynamical system is a triple ( ) where M is called the phase space and is usually a smooth manifold or a subset of n  , : Stability of a system is one of the most important parts of the studying the dynamical behavior of a system. Generally speaking, an unstable and also a chaotic system are not useful and we like to work with a system with stable and or periodic behavior (although chaos is a known behavior for many systems and sometimes people look for different strategies for chaotification of a system for different purposes [10] [11] [12]).
There are different definitions for stability, however, all they have this common fact that a system is stable if perturbation, external input and or intentionally applied signals cannot make the system get away from the equilibrium point [13] [14]. There are three possibilities for dynamical behavior of a system after applying a perturbation [5]: 1) The system state would return to the equilibrium state.
2) The system state would not return to the equilibrium state but stays near to that state.
3) The system state diverges from the equilibrium state.
Mathematically speaking, the equilibrium state * x is stable if for each initial x , the corresponding trajectory ( ) In this paper, we present some results regarding the study of local dynamics of

Dynamical Systems Playground
When we start to analyze the local dynamics of non-linear systems, the first step is finding the critical points and then exploring how the trajectories of the system evolving in the neighborhood of critical point. This analysis helps us to find out how other solutions or trajectories of the system behave when they get close to the critical points. Another step to analyze a dynamical system is studying the trajectories which trace out a limit cycle or a closed curve. In this case, the solution ( ) x t of the system will go around and create a closed curve C with a certain period T. Therefore, the solutions ( ) x t y t = of the system when it becomes periodic change to be ( ) ( ) ( ) ( ) x t y t T y t + = + = for all t. Any trajectories which are close to the limit cycle C, follow the same behavior as the limit cycle C. For instance, they can get spiral in toward C, or they can spiral away from C, which demonstrates if the closed curve C is stable or unstable. See has the following form: As it is clear from Figure 1, the point ( ) 0, 0 is a stable fixed point and also the maximum of ( ) and it is: In and Figure 2, we see has the following form: The minimum of ( ) Figure 3. The maximum of ( ) and it equals to: Also, the Taylor expansion for ( ) 2 , F X Y at 0 X = has the following form: As we can see in Figure 3, ( ) ( ) are stable. In Figure 4, we can see for are stable. Journal of Applied Mathematics and Physics ( ) ( ) ( ) ( ) ( ) The maximum of ( ) 4 , Also, the Taylor expansion for is a root and Taylor expansion for give the maximum of ( ) Finally, in Figure 6, and Taylor expansion for give the minimum of ( ) ( ) ( ) ( ) T. Azizi, G. Kerr

Hénon-Heiles System
The Hénon-Heiles potential is one of the simplest examples of classical mechanics and Hamiltonian systems [15] [16] [17] [18]. The Hénon-Heiles Hamiltonian demonstrates the emotion of stars around a galactic center. In 1964, Michael Hénon and Carl Heiles simplified the problem of the emotion of stars around a galactic center by using a Hamiltonian to describe the motion of stars near the equilibrium [16]. The Hénon-Heiles system has a vide application is studying chaotic dynamics in a system. If the energy of the motion becomes close to the bounding energy of the potential sink which is surrounding the center of the potential, this system displays chaotic dynamics feit 1984 wave. Consider the following nonlinear system of ordinary differential equations The Hamiltonian function for this system has the form ( ) For any , x y satisfying (13), we have  , H x y demonstrates another kind of dynamics and it is called the Bogdanov-Takens equilibrium point and the Jaobian matrix in this case has two zero eigenvalues. As it can be seen in Figure 7, the critical point ( ) 0, 0 is unstable which is the property of Bogdanov-Takens equilibrium point.

Duffing Oscillator
The Duffing oscillator is a single ordinary differential equation which represents a nonlinear damped driven oscillator. This simple nonlinear system displays different kinds of dynamical behaviors from periodic and regular behaviors to chaos. When we add a driving force and friction, we can see this simple equation exhibit chaotic behavior [19] [20] [21]. The Duffing oscillator equation has the following form The law of energy conservation mentions that this is impossible to see chaotic motion in a single degree of freedom. Therefore, with adding a driving force and damping, the energy conservation would be eliminated. Then, the equations of motion has the form We have demonstrated different dynamical behaviors of (20) in Figures 8-10. The Duffing oscillator can be used to model different physical phenomenon such as stiffening springs, beam buckling, nonlinear electronic circuits, superconducting Josephson parametric amplifiers, and ionization waves in plasmas [22].   T. Azizi, G. Kerr

The Van Der Pol Equation
Van der Pol equation which is a well known second order ordinary differential equation with cubic nonlinearity has attracted many researchers in different field of science. This self oscillatory system, Van der Pol oscillator, has been considered as very useful mathematical model for many complicated systems [23] [24] [25]. Mathematical representation of the Van der Pol system has the form ( ) where constant µ is a positive parameter depending on the tube constants.
This equation represents describes the current ( ) x t in a certain type of vacuum tube. We can write (21) as a first order system of differential equations: The numerical integration of Equation (22)

Chaos in Continuous Dynamical Systems
In this section, we assume that there is a Hamiltonian function with two degrees of freedom and it is given by H makes H to be non-integrable. For 0 ε = and also for 0 1 ε <  , there exist quasi periodic cycles which are known as KAM tori. However, under perturbation, these quasi periodic cycles will be deformed and KAM tori will be dissolved gradually as we increase ε . This phenomenon can be observed in Figure 15 and Figure 16.        According to, KAM theory when x is irrational, then the torus is preserved for small perturbation ε . But, proportional tori and adjacent irrational tori would be destroyed. Also, the stable manifold and unstable manifold of the saddle point which are intersecting transversely, appear to be Smale horseshoe and chaotic motion. As ε increases gradually, these chaotic layers grow and they envelope larger area in phase space [26] [27].

Conclusion
Dynamic systems modeling have been frequently used to describe different physical systems and have a very important role in predicting the interactions between multiple components of a system over time. In the present study, we explored different dynamical behaviors of some continuous dynamical systems, from stable and regular motions to periodic and limit cycles, and then chaotic and irregular oscillations. We started with studying the local dynamics of some vector fields and we demonstrated the local stable behavior of the system around its critical points. We continued this paper with studying the well known problems which have been used a lot for different physical purposes. Hénon-Heiles system, Duffing oscillator and Van der Pol equation are three important dynamical examples which have been widely studied numerically. We demonstrated the stable and unstable manifolds from the origin form a homoclinic orbit in Hénon-Heiles system and we discussed about the local dynamical behaviors of its critical points. We showed that the critical point ( ) 0, 0 is a saddle point and critical point ( ) 1, 0 − is a center. For Duffing oscillator, which can be used to model different physical phenomenon, we showed the periodic and chaotic motions of the system using time series. Also, for Van der Pol equation, we presented the limit cycle solutions and periodic behavior of the system. We concluded that depending on different values for µ , solutions look like periodic motion. When µ small values, this motion is nearly sinusoidal, however for larger values of µ , the solutions seem to be relaxation oscillations, which means solutions are similar to a series of step functions and jump twice per cycle between the positive and the negative values. Finally, we discussed about the chaos in Hamiltonian systems and we provided two interesting and different examples which exhibit chaotic behaviors. We assume a Hamiltonian function with two degrees of freedom and it can be obtained by adding an integrable Hamiltonian system and a non-integrable Hamiltonian system. We showed that for 0 ε = and also for 0 1 ε <  , there exist quasi periodic cycles which are known as KAM tori. However, under perturbation, these quasi periodic cycles will be deformed and KAM tori will be dissolved gradually as we increase ε .