Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 138430, 8 pages

http://dx.doi.org/10.1155/2013/138430

## Bifurcations in Van der Pol-Like Systems

^{1}Department of Mathematics, Gebze Institute of Technology, Kocaeli, 41400 Gebze, Turkey^{2}Department of Information Systems and Technologies, Yeditepe University, Atasehir, 34755 Istanbul, Turkey^{3}Department of Computer Education and Instructional Technology, Yeditepe University, Atasehir, 34755 Istanbul, Turkey^{4}Department of Physics, Yeditepe University, Atasehir, 34755 Istanbul, Turkey^{5}Department of Physics, Bogazici University, Bebek, 34342 Istanbul, Turkey

Received 2 October 2013; Accepted 25 November 2013

Academic Editor: Panos Liatsis

Copyright © 2013 Orhan Ozgur Aybar et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Generalizations of the Van der Pol system with polynomial interactions involving additional parameters are studied in order to understand qualitative properties such as stability and additional bifurcations in these generalized Van der Pol systems. The generalizations include the one known as the Duffing-Van der Pol system which has properties similar to those of the Maxwell-Bloch system.

#### 1. Introduction

Dynamical systems underlying many oscillatory real life problems that arise in applied sciences can be considered as generalizations of the Van der Pol (henceforth referred to as VP) system [1]. Due to its wide area of applications in a number of scientific areas such as electronics and geology, the VP system stands as an important example of dynamical systems. It was first developed by the Dutch physicist Balthasar Van der Pol (1889–1959) to model power amplifiers in vacuum tube radios exhibiting negative resistance [1]. Methods referred to by names such as averaging by Guckenheimer and Holmes [2] and two-variable expansion by Rand and Armbruster [3] have been applied to this system. Related types of averaging methods have also been applied by Strogatz [4], Verhulst [5, 6], Sanders et al. [7], Kuznetsov [8], Jordan and Smith [9], and Hale and Koçak [10] to similar oscillatory problems. We use standard averaging in the rest of this paper [2, 5, 7]. Oscillatory solutions to dynamical systems using the center manifold and normal form methods, and bifurcation and stability analysis have been extensively studied by Cao and Xiao [11–13]. This paper aims to look at possible generalizations of VP and Duffing-like systems and to study their bifurcation schemes in order to clarify our understanding. The methodology, some of the observed bifurcation scenarios, and attractor structures resemble those reported by Cao and Xiao in other systems [11–13].

Section 2 discusses quadratic generalizations of the VP system including the one known as the Duffing-Van der Pol (henceforth referred to as DVP) system and compares the averaging results with the averaged versions of some commonly known similar dynamical systems. Section 3 introduces a new generalization for the DVP system and investigates results obtained by methods discussed in previous sections. In Section 4, generalizations that introduce higher nonlinearities are studied in the same manner. Overall results are compared and discussed in Section 5.

The VP equation is the second-order ordinary differential equation modeling the VP oscillator [1]. The part of the VP equation including together with its coefficient implies that the friction or the resistance represented by the nonlinear-damping term changes sign and becomes negative as the amplitude increases. This model can be derived from the Rayleigh equation [14] and can be reduced to simple harmonic motion for with the solution .

The generalized models under study can be stated as vector equations given in system (2) by the transformation in order to use dynamical systems approaches for their study [5, 12].

The generalized VP equations [1] under study are special cases of the following system: where and are polynomials in the variable and we can always set by and . System (2) is also known as the Liénard system [15].

Additional parameters in the nonlinearities can be introduced to generate and analyze further bifurcation varieties as follows: so that the generalized VP equation is now

We carry out the stability analysis by looking into the eigenvalues of the linearized matrix of coefficients. In this manner, the system has a trivial equilibrium point at the origin with the eigenvalues . Scaling the time by the undamped angular frequency reveals that the only nontrivial parameter is , so that itself affects directly the stability of the system.

Bifurcation analysis is the key to understanding the qualitative behavior of a dynamical system. Qualitative changes, if any, in the system when any of the parameters are varied, represent bifurcations and help us understand the behavior of the system better. Using the software package MATCONT, bifurcations of the systems of concern are studied numerically [8]. Choosing as the bifurcation parameter, the obvious scenario for the VP system is obtained such as a family of limit cycles growing endlessly from the Hopf bifurcation point at the trivial equilibrium point indicated by a nonnegative first Lyapunov coefficient [2]. Considering that in a physical system all parameters are nonnegative, it can be said that the equilibrium point of system (3) is always unstable. Details of other bifurcation properties for the simple VP system as given in system (3) including transcritical bifurcation are omitted here for brevity [9, 10, 16].

Before rescaling, the equation has two nontrivial coefficients, one determining the undamped frequency and the other determining the damping term. The latter changes sign as the amplitude is increased. A generalization of the VP equation should involve either or both of these terms. Our object of study will be generalizations where nonlinear functions on will be used instead of the term in equation. Nonlinear functions depending on can be used as generalizations [17].

#### 2. The Generalized Van der Pol System including Quadratic Self-Nonlinearity

Replacing the linear term in the equation of system (3) by the quadratic combination modifies the system to

In addition to the trivial equilibrium point at the origin that always occurs, a second equilibrium point at is observed in this case. The eigenvalues (i.e., the eigenvalues of the coefficients matrix of the linearized system) for the trivial equilibrium point are the same as those of the original one and the eigenvalues for the second equilibrium point are . For , the origin has pure imaginary eigenvalues that indicate Hopf bifurcation. We want to analyze Hopf bifurcation in this system by using the following theorem.

Theorem 1 (Hopf bifurcation, [8]). *A Hopf bifurcation occurs if a planar system , , , has the following eigenvalues at origin for sufficiently small taken as bifurcation parameter:
**
where real part of is zero for and imaginary part of is .**The first Lyapunov coefficient defined as* *is not equal to zero.**Consider (transversality condition).*

By using Theorem 1, we require that the determinant of the Jacobian matrix be greater than zero and trace of the Jacobian matrix be zero at the related equilibrium point. According to criterion in Theorem 1, Hopf point also obeys the transversality condition if the rate of change of real part of eigenvalues at the related equilibrium (codimension-one Hopf) point is greater than zero which means that eigenvalues of the linearized flow cross the imaginary axis with nonzero derivative when bifurcation parameter is zero [8, 9]. We use the first Lyapunov coefficient to express Hopf bifurcation criteria for bifurcation parameter . When the system has pure imaginary eigenvalues at the related equilibrium (Hopf) point, . In the neighborhood of the Hopf point, if the first Lyapunov coefficient is less than zero, supercritical Hopf bifurcation can be observed; if is greater than zero, subcritical Hopf bifurcation can be observed [8, 9]. When these conditions are satisfied, the system has codimension one Hopf bifurcation [8]. If vanishes and becomes zero, this type is called degenerate (generalized) Hopf Bifurcation [8].

We see that trace and determinant of the Jacobian matrix at origin are and , respectively. By applying the conditions, we obtain pure imaginary roots depending on , ; the transversality condition gives at the origin. Under the assumption and sufficiently small bifurcation parameter [18], we have the following conditions at the origin: if , , subcritical Hopf bifurcation is observed with negative first Lyapunov coefficient and if , , supercritical Hopf bifurcation is observed with positive first Lyapunov coefficient [8]. The famous limit cycles of the VP equation are observed to bifurcate from the Hopf bifurcation point (indicated by H) of the system (5) as shown in Figure 1. When , the first Lyapunov coefficient is zero.

The determinant and trace at the other equilibrium point are and zero, respectively. By using the above conditions, is less than zero and rate of change in real part of eigenvalues which is expected to be nonzero vanishes. At this equilibrium point we thus do not observe Hopf bifurcation as shown in Figure 1 according to Theorem 1.

*Definition 2 (normal form expansion). *Let the system of differential equation and be given where and . To consider the system near its equilibrium point and such that and , this point is moved to the origin by and . Taylor expansion near and gives
where , , , are coefficients and and contain terms of degree two and higher. Diagonalizing or bringing the linear part to the Jordan canonical form by a linear transformation produces
where , denote Jordan canonical form terms of the linear part in system (8). Then a sequence of near identity transformations of the form
is applied to simplify the system (8), where and contain terms of degree two and higher. The normal form expansion is given as follows:
where and contain degree two and higher order terms.

We use normal form expansion around the trivial equilibrium point (origin) to yield Hopf bifurcation as defined in Definition 2. In the rest of this work, and will refer to the variables in the near identity transformation [19]. For example, system (5) has the following normal form expansion up to the fifth order which leads to degenerate Hopf bifurcation at origin for which is also consistent with the Hopf bifurcation observed at origin in Figure 1.

Figure 2 gives the bifurcation curves with lines indicating pitchfork bifurcation points (indicated by BP). The bifurcation analysis while relaxing the bifurcation parameter shows a pitchfork and a degenerate limit point bifurcation (indicated by LP) point with limit point coefficient equal to zero as shown in Figure 2. For , , the axis contains equilibrium points with corresponding eigenvalues that satisfy the limit point (saddle-node) bifurcation criteria [5, 8].

##### 2.1. Averaging of the Generalized Van der Pol System including Quadratic Self-Nonlinearity

The method of averaging is a powerful perturbation method which gives a normal form that can be compared to more commonly known dynamical systems. Numerically obtained results concerning stability and bifurcation can also be analytically confirmed using the method of averaging. In order to understand the behavior of the system, we investigate the averaged version of system (5) and observe its similarity to averaged versions of some better known systems.

The averaged equations of system (3) and (5) are calculated by the usual method discussed in [2, 5, 7]. We introduce a small parameter , and replace and by , , respectively, to analyze this system for small parameter values while averaging is carried out over a period in . We use to obtain According to the theorem given by Verhulst [5, 6], the origin is expected to be the only equilibrium point. We calculate and and introduce and to get in the standard averaging sense [2, 5, 7]. This result shows a typical Hopf bifurcation form [4, 7, 8]. The steady state of leads either to the trivial case or to the limit cycle condition where or implying the existence of a limit cycle of radius . Numerical and averaged solutions are given in Figure 3 that also gives unstable (stable) equilibrium point for . For small , as is increased, the averaged solution does not tend to the numerical solution of the original equation as shown in Figure 4. The determinant of the Jacobian matrix at the related equilibrium point is greater than zero . The trace becomes zero if is zero. With the aid of Theorem 1, subcritical (supercritical) case is observed for sufficiently small , , with the first Lyapunov coefficient being less (greater) than zero.

System (5) is a special case of the following averaged version of the general form of generalized VP systems: where we introduce parameters , and as coefficients of the generalization terms in parentheses. The averaged form of system (16) calculated by the usual method discussed in [5, 7] becomes by introducing , and by replacing , by , in that given order.

A number of other systems with linear periodic parts also give similar results upon averaging. For example, the cubic generalization of the Lotka-Volterra system in [19] upon averaging (replacing , , , by , , , ) by the usual method discussed in [5, 7] gives where for values and for .

#### 3. The Duffing-Van der Pol System

A further generalization of the VP system that introduces the cubic property of the Duffing system in the potential function is the DVP system given by where parameters and are introduced in the potential function [14]. Depending on the sign of the parameter , single or double well potential cases exist. There are three equilibrium points, that is, the trivial equilibrium point at the origin and two equilibrium points at . The eigenvalues at the equilibrium points are respectively. When a sufficiently small is chosen as bifurcation parameter, trace of the Jacobian is zero and determinant of the Jacobian matrix is at the origin. By using Theorem 1, eigenvalues at the origin are that give Hopf bifurcation at the origin as shown in Figure 5. For and , the first Lyapunov coefficient is less than zero that gives supercritical Hopf bifurcation and, for and , the first Lyapunov coefficient is greater than zero that gives subcritical Hopf bifurcation.

The Jacobian matrix of system (20) at is where eigenvalues of the Jacobian matrix are purely imaginary if the determinant of the Jacobian matrix is and the trace of the Jacobian matrix is . These two conditions are used to show Hopf bifurcation. We see that the trace becomes zero when or . By assuming the latter one, the second and third equilibrium points become . Another necessary condition is that the rate of change of real part of eigenvalues with respect to bifurcation parameter should be greater than zero where . Combining these results, we see that Hopf bifurcation occurs if , , , or , , , given in Figures 6 and 7 for the bifurcation parameter chosen as .

As an example, we set the determinant of the Jacobian matrix at the origin to and obtain . Since the pair of eigenvalues cannot have real part according to the Hopf bifurcation condition, we can assume that where gives pure imaginary eigenvalues as shown in Figure 5. The second and third equilibrium points do not give pure imaginary eigenvalues; hence, we do not observe Hopf bifurcation at these equilibrium points. When , the second and third equilibrium points possess positive determinants and by using and , system has pure imaginary eigenvalues as shown in Figures 6 and 7. This behavior observed for the single well potential case of the DVP system with parameter free is similar to that of the Lorenz system in which an unstable equilibrium point is between two stable equilibrium points. Two families of limit cycles bifurcating from the two unstable Hopf points intersect on the bifurcation curve (line) of the system. The pitchfork bifurcation point forks a quadratic curve with a negative slope which passes through the two subcritical Hopf bifurcation points so that the projection of the Lorenz attractor created would seem to intersect at the pitchfork bifurcation point BP in the two-dimensional graph. In this case, the equilibrium points of the system are again , , and and the eigenvalues are , , and , respectively. The limit cycles go around two subcritical Hopf bifurcation points indicated by two positive first Lyapunov coefficients. The bifurcation analysis when is varied shows a pitchfork bifurcation curve with numerous limit point bifurcation points on it. To understand the bifurcation properties in system (20), we calculate the averaged equations by the usual methods [2, 5, 7] by introducing a small parameter , and replacing , by , , respectively.

Numerical and averaged solutions are given in Figure 8 where the origin is a stable (unstable) equilibrium point for , , , .

##### 3.1. The Maxwell-Bloch System

It has been of considerable interest to relate models for laser dynamics to VP and DVP systems. DVP model is known to give oscillatory solutions which includes several Hopf bifurcations [20]. The Maxwell-Bloch equations (to be referred to as MB henceforth) are that contain one unstable equilibrium point at , with one positive and two negative eigenvalues and pair of equilibrium points at either side with eigenvalues . When , , , and , MB system (24) can be transformed into the Lorenz system about the equilibrium point by setting , , [20]. It was reported that operating point shows chaotic behavior that is similar to that of the Lorenz system [20]. The behavior of the MB system around the point is compared to the behavior of the DVP system around each of the equilibrium points , , and in Figure 9. For both systems, Figure 9 shows trajectories involving similar structures.

#### 4. The Generalized Van der Pol System including Quartic Self Nonlinearity

Now we introduce a further generalized system including a quartic term in the parentheses. The generalization with this additional quartic term is The averaged version of system (25) is by introducing a small parameter , and by replacing , by , in that given order.

According to the sign of , the system exhibits different bifurcation scenarios where some of them are of interest and others are trivial.

A bifurcation analysis selecting as the bifurcation parameter gives a pitchfork bifurcation point which forks a curve with numerous limit point bifurcation points as in the case for the DVP system. Selecting as the bifurcation parameter, a pitchfork bifurcation is observed, and selecting as the bifurcation parameter, the usual subcritical Hopf bifurcation point scenario is observed. Other parameters do not give any bifurcation in this generalization.

The averaged version of the following further generalized system: is where , , , and scenarios depend on the characteristic of the numerous roots of .

#### 5. Conclusion

The subcritical Hopf bifurcation point of the VP system at the origin is inherited through all polynomial generalizations in the VP system. Generally a pitchfork bifurcation is observed in all bifurcation graphs which forks into higher bifurcations. Many limit point bifurcations are observed for all generalizations. Pitchfork and limit point bifurcations are observed for generalizations of VP system including quadratic self-coupling. Hopf bifurcation has also been studied using the normal form method as given in Definition 2, for the eigenvalues of the linearized system at equilibrium points. For generalized VP systems considered in this study, averaging over the period of the linear part, a limit cycle property is observed for all possible cases. To this approximation, a change in the nonlinearity of the system does not change the occurrence of limit cycle phenomena. However in polar coordinates does not remain constant. Also the case where the bifurcation curves cycle and oscillate around a pitchfork bifurcation point between two Hopf points only occurs for the original DVP system which is nonintegrable. A relation between the DVP and MB systems where both exhibit Lorenz-like behavior has been shown.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The authors would also like to thank the editor and referees for their valuable and constructive comments. The authors would like to thank Professor Dr. Ferdinand Verhulst for his valuable suggestions and comments that helped to improve the paper throughout the process.

#### References

- B. van der Pol, “On relaxation-oscillations,”
*The London and Edinburgh Philosophical Magazine and Journal of Science*, vol. 7, pp. 978–992, 1927. View at Google Scholar - J. Guckenheimer and P. Holmes,
*Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields*, Springer, New York, NY, USA, 2002. View at MathSciNet - R. H. Rand and D. Armbruster,
*Perturbation Methods, Bifurcation Theory and Computer Algebra*, vol. 65 of*Applied Mathematical Sciences*, Springer, New York, NY, USA, 1987. View at Publisher · View at Google Scholar · View at MathSciNet - S. H. Strogatz,
*Nonlinear Dynamics and Chaos: With Applications to Physics, Biology, Chemistry and Engineering*, Westview Press, 2001. - F. Verhulst,
*Nonlinear Differential Equations and Dynamical Systems*, Springer, New York, NY, USA, 2nd edition, 2006. View at Publisher · View at Google Scholar · View at MathSciNet - F. Veerman and F. Verhulst, “Quasiperiodic phenomena in the Van der Pol-Mathieu equation,”
*Journal of Sound and Vibration*, vol. 326, no. 1-2, pp. 314–320, 2009. View at Publisher · View at Google Scholar · View at Scopus - J. A. Sanders, F. Verhulst, and J. Murdock,
*Averaging Methods in Nonlinear Dynamical Systems*, vol. 59 of*Applied Mathematical Sciences*, Springer, New York, NY, USA, 2007. View at MathSciNet - Y. A. Kuznetsov,
*Elements of Applied Bifurcation Theory*, vol. 112 of*Applied Mathematical Sciences*, Springer, New York, NY, USA, 1995. View at MathSciNet - D. W. Jordan and P. Smith,
*Nonlinear Ordinary Differential Equations: An Introduction to Dynamical Systems*, vol. 2 of*Oxford Texts in Applied and Engineering Mathematics*, Oxford University Press, Oxford, UK, 1999. View at MathSciNet - J. K. Hale and H. Koçak,
*Dynamics and Bifurcations*, vol. 3 of*Texts in Applied Mathematics*, Springer, New York, NY, USA, 1991. View at Publisher · View at Google Scholar · View at MathSciNet - M. Xiao and J. Cao, “Genetic oscillation deduced from Hopf bifurcation in a genetic regulatory network with delays,”
*Mathematical Biosciences*, vol. 215, no. 1, pp. 55–63, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Cao and M. Xiao, “Stability and Hopf bifurcation in a simplified BAM neural network with two time delays,”
*IEEE Transactions on Neural Networks*, vol. 18, no. 2, pp. 416–430, 2007. View at Publisher · View at Google Scholar · View at Scopus - M. Xiao and J. Cao, “Delayed feedback-based bifurcation control in an Internet congestion model,”
*Journal of Mathematical Analysis and Applications*, vol. 332, no. 2, pp. 1010–1027, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. J. W. S. Rayleigh,
*The Theory of Sound*, Dover, New York, NY, USA, 1945. View at MathSciNet - A. Liénard, “Etude des oscillations entretenues,”
*Revue Générale de l'Électricité*, vol. 23, pp. 901–912, 946–954, 1928. View at Google Scholar - F. C. Moon and P. J. Holmes, “A magnetoelastic strange attractor,”
*Journal of Sound and Vibration*, vol. 65, no. 2, pp. 275–296, 1979. View at Google Scholar · View at Scopus - W. Yu and J. Cao, “Hopf bifurcation and stability of periodic solutions for van der Pol equation with time delay,”
*Nonlinear Analysis: Theory, Methods & Applications*, vol. 62, no. 1, pp. 141–165, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Doelman and F. Verhulst, “Bifurcations of strongly nonlinear self-excited oscillations,”
*Mathematical Methods in the Applied Sciences*, vol. 17, no. 3, pp. 189–207, 1994. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - I. Kusbeyzi, O. O. Aybar, and A. Hacinliyan, “Stability and bifurcation in two species predator-prey models,”
*Nonlinear Analysis: Real World Applications*, vol. 12, no. 1, pp. 377–387, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - I. Kusbeyzi, O. O. Aybar, and A. S. Hacinliyan, “Approximate solutions of Maxwell Bloch equations and possible Lotka Volterra type behavior,”
*Nonlinear Dynamics*, vol. 62, no. 1-2, pp. 17–26, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet