Application of Stability Theory in Study of Local Dynamics of Nonlinear Systems ()
1. 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 modeling [1] [2] [3]. When the time between two measurements is negligible, the continuous time modeling governs the evolution of the system, however, when there is a gap between two measurements, discrete-time system modeling comes to play. Ordinary differential equations are the tool to model a continuous system and iterated maps represent the discrete generations [4] [5].
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
,
, called the evolution, is a smooth action of K in M and K is either a subset of
in the case of a continuous time dynamical system or a subset of
in the case of a discrete time dynamical system. The smooth action
describes the evolution with time
of a point x in the phase space M [4] [5] [6] [7].
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
is stable if for each initial conditions
close enough to
, the corresponding trajectory
remains near
for all
.
In this paper, we present some results regarding the study of local dynamics of non-linear continuous time dynamical systems. We provide different examples to display stable and unstable limit cycles and we demonstrate the numerical results for each case. Also, we study the local dynamics of three well known physical systems, Henon-Heiles system, Duffing oscillator and Van der pol equations and we display the evolution of solutions of the system in time. Finally, we discuss about the chaos in Hamiltonian systems and we provide two examples to show chaos in Hamiltonian systems.
2. 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
of the system will go around and create a closed curve C with a certain period T. Therefore, the solutions
of the system when it becomes periodic change to be
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 Figure 1.
The root point for
is
. The Taylor expansion for
at
has the following form:
(1)
As it is clear from Figure 1, the point
is a stable fixed point and also the maximum of
occurs at
and it is:
(2)
In and Figure 2, we see
. The root for
is
. The Taylor expansion for
at
has the following form:
(3)
The minimum of
happens for
which is equal:
(4)
For
the point
is unstable.
Another example,
which has been displayed in Figure 3. The maximum of
occurs at
and
and it equals to:
(5)
Also, the Taylor expansion for
at
has the following form:
(6)
As we can see in Figure 3,
and
are stable.
In Figure 4, we can see for
,
and
are stable.
Figure 1.
, local dynamics.
Figure 2.
, local dynamics.
Figure 3.
, local dynamics.
Figure 4.
, local dynamics.
The maximum of
occurs at
and
and equals to
(7)
Also, the Taylor expansion for
at
has the form
(8)
For
(Figure 5),
is a root and Taylor expansion for
at
has the form
(9)
here,
and
give the maximum of
which is
(10)
Finally, in Figure 6,
has a root at
and Taylor expansion for
at
has the form
(11)
here,
and
give the minimum of
which is
Figure 5.
, local dynamics.
Figure 6.
, local dynamics.
(12)
3. Application of Continuous Dynamical Systems Modeling
3.1. 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
(13)
The Hamiltonian function for this system has the form
(14)
For any
satisfying (13), we have
. For any solution
of system (13), the Hamiltonian
is constant, it means
. This is a very nice property of Hamiltonian function which is a conserved quantity for a system of ordinary differential equations and it is constant along all solution curves of the system.
The solution curves are given by
. Here, there are two non degenerate critical points
and
. The critical point
is a saddle point and the eigenvectors corresponding to this critical points are
and
. The critical point
is a center. Figure 7 displays the level curves or contours of four different Hamiltonian functions.
Consider the following Hamiltonian functions
(15)
(16)
(17)
(18)
Figure 7. Hénon-Heiles system phase portraits.
These Hamiltonian functions (15)-(18), are corresponding to different system of ordinary differential equations. The Hamiltonian function
has a critical point at
and the Hamiltonian function
has a critical point at
. As we can see in Figure 7, the stable and unstable manifolds from the origin for
and
form a homoclinic orbit which we can not see this property in Hamiltonian functions
and
. This homoclinic loop connects the critical point
to itself and it takes infinite amount of time to make connection. For Hamiltonian functions
and
, the critical point
is called a saddle-node equilibrium and the Jacbian matrix of the system has a zero eigenvalues at this equilibrium point. However, the critical point
for the Hamiltonian functions
and
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
is unstable which is the property of Bogdanov-Takens equilibrium point.
3.2. 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
(19)
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
(20)
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].
Figure 8. Chaotic solutions of Duffing oscillator (20) for
,
,
,
,
, Periodic solutions of Duffing oscillator (20) for
,
,
,
,
.
Figure 9. Periodic solutions of Duffing oscillator (20) for
,
,
,
,
, Chaotic solutions of Duffing oscillator (20) for
,
,
,
,
.
Figure 10. Periodic solutions of Duffing oscillator (20) for
,
,
,
,
, Chaotic solutions of Duffing oscillator (20) for
,
,
,
,
.
3.3. 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
(21)
where constant
is a positive parameter depending on the tube constants. This equation represents describes the current
in a certain type of vacuum tube. We can write (21) as a first order system of differential equations:
(22)
The numerical integration of Equation (22) has been represented in Figures 11-14.
As we can see, 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.
4. 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
. Here, we consider
to be a very small parameter,
an integrable Hamiltonian system and
makes H to be non-integrable. For
and also for
, 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.
Figure 11. Solutions of Van der Pol Equation (22) with
and
.
Figure 12. Solutions of Van der Pol Equation (22) with
and
.
Figure 13. Solutions of Van der Pol Equation (22) with
and
.
Figure 14. Solutions of Van der Pol Equation (22) with
and
.
Figure 15. Dissolving the KAM tori caused by perturbation,
.
Figure 16. Dissolving the KAM tori caused by perturbation,
.
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].
5. 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
is a saddle point and critical point
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
and also for
, 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
.