Abstract and Applied AnalysisVolume 2012 (2012), Article ID 974293, 10 pagesdoi:10.1155/2012/974293

Research Article

Application of Multistage Homotopy Perturbation Method to the Chaotic Genesio System

M. S. H. Chowdhury,¹ I. Hashim,² S. Momani,³ and M. M. Rahman⁴

¹Deptartment of Science in Engineering, Faculty of Engineering, International Islamic University Malaysia, Jalan Gombak, 53100 Kuala Lumpur, Malaysia²School of Mathematical Sciences, Universiti Kebangsaan Malaysia, 43600 Bangi Selangor, Malaysia³Department of Mathematics, The University of Jordan, Amman 11942, Jordan⁴Department of Physics, Faculty of Science, Universiti Putra Malaysia, 43400 Selangor, Malaysia

Received 31 December 2011; Accepted 9 March 2012

Academic Editor: Muhammad Aslam Noor

Copyright © 2012 M. S. H. Chowdhury 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.

Finding accurate solution of chaotic system by using efficient existing numerical methods is very hard for its complex dynamical behaviors. In this paper, the multistage homotopy-perturbation method (MHPM) is applied to the Chaotic Genesio system. The MHPM is a simple reliable modification based on an adaptation of the standard homotopy-perturbation method (HPM). The HPM is treated as an algorithm in a sequence of intervals for finding accurate approximate solutions to the Chaotic Genesio system. Numerical comparisons between the MHPM and the classical fourth-order Runge-Kutta (RK4) solutions are made. The results reveal that the new technique is a promising tool for the nonlinear chaotic systems of ordinary differential equations.

1. Introduction

Chaos is very interesting nonlinear phenomenon and has been intensively studied in the last three decades. The dynamical systems that exhibit chaotic behavior are sensitive to initial conditions. Chaotic behavior can be found in a variety of systems such as electrical circuits, lasers, fluid dynamics, mechanical devices, time evolution of the magnetic field of celestial bodies, population growth in ecology, the dynamics of molecular vibrations, and not forgetting the weather. The history of chaos theory has come a long way since Jacques Hadamard who in 1898 published a significant study of a free particle gliding frictionlessly on a surface of constant negative curvature which exhibits chaotic motion. The Genesio-Tesi system, proposed by Genesio and Tesi [1], is one of paradigms of chaos since it captures many features of chaotic systems. It includes a simple part and three simple ordinary differential equations (ODEs) that depend on three positive real parameters. The dynamic equation of the system is as follows:

where , , are state variables, and , , are positive real constants satisfying . For instance, the system is chaotic for the parameter , , .

In recent years, much attention has been devoted to the application of the HPM [2], to the solutions of various scientific models [3–6]. HPM yields rapidly convergent series solutions [7, 8]. Now, the application of HPM has been extended to the thermal problems [9]. The homotopy perturbation method admits some unknown parameters in the obtained series solutions. In [10], Ji-Huan suggested the least square method can be used to identify the unknown parameters involved in the series solutions.

Recently, Chowdhury et al. [11] was the first to successfully apply the multistage homotopy-perturbation method (MHPM) to the chaotic and nonchaotic Lorenz system. The MHPM is also applied to solve the Chen system and a class of systems of ODEs [12–14].

In this paper we are again interested in the accuracy of the MHPM for nonlinear systems of ODEs capable of exhibiting chaotic behavior. The system which is of interest to us is the chaotic Genesio system (1.1). Numerical comparisons with the HPM and fourth-order Runge-Kutta (RK4) solutions show that MHPM is accurate and efficient.

2. Solution Approaches

We consider a general system of first-order ODEs: subject to the initial conditions:

First we write system (2.1) in the operator form:subject to the initial conditions (2.2), where is a linear operator and are the nonlinear operators. We will next present the solution approaches for (2.3) based on the standard HPM and MHPM separately.

2.1. Solution by HPM

According to HPM, we construct a homotopy for (2.3) which satisfies the following relations:

where is an embedding parameter and are initial approximations satisfying the given conditions. It is obvious that when the perturbation parameter , (2.4) become a linear system and when we get the original nonlinear system.

Let us take the initial approximations as follows:

where , are functions yet to be determined. Substituting (2.5) into (2.4) and arranging the coefficients of the same powers of , we get

and so forth. We solve the above systems of equations for the unknowns by applying the inverse operator

Therefore, according to HPM the -term approximations for the solutions of (2.3) can be expressed as

2.2. Solution by MHPM

The approximate solutions (2.8) are generally, as will be shown in the numerical experiments of this paper, not valid for large . A simple way of ensuring validity of the approximations for large is to treat (2.6) as an algorithm for approximating the solutions of (2.1) in a sequence of intervals choosing the initial approximations as

where is the left-end point of each subinterval.

Now we solve (2.6) for the unknowns by applying the inverse linear operator

In order to carry out the iterations in every subinterval of equal length , , , , we need to know the values of the following:

But, in general, we do not have these information at our clearance except at the initial point . A simple way for obtaining the necessary values could be by means of the previous -term approximations of the preceding subinterval given by (2.8), that is,

3. Application

In this section,we will study the Genesio system (1.1) subject to the initial conditions:

According to the HPM, we can construct a homotopy which satisfies the following relation:

We take the initial approximations as:

where ,   are functions yet to be determined. Substituting (3.3) into (3.2) and collecting terms the same powers of , we have

In order to obtain the unknowns ,  , we solve the above system taking the initial conditions ,  , we obtain,

Hence, the solution to the Genesio system (1.1) is: To carry out the iterations on every subinterval of equal length , we need to know the values of the following initial conditions:

In general, we do not have these information at our clearance except at the initial point , but we can obtain these values following the MHPM as given in Section 2.2. We note that the 10-term approximations of , , and are denoted as , , and .

4. Results and Discussions

The MHPM algorithm is coded in the computer algebra package Maple and we employ the Maple’s built-in fourth-order Runge-Kutta procedure (RK4). The Maple environment variable digits controlling the number of significant digits are set to 16 in all the calculations done in this paper. For the comparison, we set the parameter , , and where the system exhibits chaotic behavior, alongside with its initial conditions , , and . The time range studied in this work is to 20. In general, there is not a known exact solution for the Genesio system. Thus, the accuracy of the present method will be determined by its comparison to the numerical solution by Runge-kutta method. We determine the accuracy of RK4 for the solution of (1.1) for different time steps. From the results presented in Table 1 we see that the maximum difference between the RK4 solutions on time steps and is of the order of magnitude of . So based on these observations we choose the RK4 solutions on the time step as the benchmark for our comparison purposes. So now we can compare the accuracy of the HPM and MHPM with the RK4 method on the chosen time step . We choose this time step since a smaller one is computationally costly. We note that increasing the number of terms in the series solutions (1.1) improves the accuracy of the MHPM solutions but at the expense of increased computational efforts. In this work we fix the number of terms used to be ten. The details of the differences between the 10-term HPM solutions and 10-term MHPM solutions () and the RK4 solutions on are given in Table 2. In Table 2, we see that the 10-term MHPM solutions on the time step agree with the RK4 solutions at least 10 decimal places, while the 10-term classical HPM solutions are only valid for . For the chaotic Genesio system we observe that the MHPM has the advantage over the RK4 on achieving a good accuracy with a larger time step.

Table 1: A determination of the accuracy of RK4 for the Genesio system (1.1).

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

Table 2: Differences between 10-term HPM and 10-term MHPM with RK4 solutions on .

1

0.0002355

2
0.01626
0.04539
0.2211

3
0.8815
1.912
10.73

4
14.83
28.01
154.9

5
126.8
236.3
1193

6
703.5
1373
6293

7
2910
6063

8
9790

9

10

11

12

13

14

15

16

17

18

19

20

In Figure 1 we reproduce the well-known , , , and phase portraits of the chaotic Genesio system using the 10-term MHPM solutions on . Having determined the benchmark timestep, we can now investigate the accuracy of the HPM and MHPM as compared to RK4. In Figure 2 we plot the 10-term HPM solutions and 10-term MHPM solutions (on ) against the RK4 solutions on . Both the 10-term MHPM solutions on and RK4 solutions on seem to overlap on the scale used in figure, but the numerical results from the standard HPM start to stay away at about for , , and .

Figure 1: Phase portraits using 10-term MHPM.

Figure 2: Graphical comparisons between HPM, MHPM (), and RK4 () for , , and .

5. Conclusions

In this work, the MHPM was applied to the solutions of the well-known Genesio system. The MHPM is only a simple modification of the standard HPM. Comparisons between the HPM and MHPM solutions and the fourth-order Runge-Kutta (RK4) numerical solutions were made. For the chaotic Genesio system studied we found that the 10-term MHPM solutions on a larger time step achieved comparable accuracy compared with the RK4 solutions on a much smaller time step. We note that the MHPM solutions were computed via a simple algorithm with less amount of computations and without any need for perturbation techniques, special transformations, linearization, or discretization.

Acknowledgments

The authors would like to acknowledge the financial supports received from the Ministry of Higher Education via the FRGS Grant no. FRGS0409-106 and Universiti Putra Malaysia.

References

R. Genesio and A. Tesi, “Harmonic balance methods for the analysis of chaotic dynamics in nonlinear systems,” Automatica, vol. 28, no. 3, pp. 531–548, 1992.

J.-H. He, “A coupling method of a homotopy technique and a perturbation technique for non-linear problems,” International Journal of Non-Linear Mechanics, vol. 35, no. 1, pp. 37–43, 2000.

M. J. Ablowitz, B. M. Herbst, and C. Schober, “Homotopy perturbation method and axisymmetric flow over a stretching sheet,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 7, no. 4, pp. 399–406, 2006.

L. Cveticanin, “The homotopy-perturbation method applied for solving complex-valued differential equations with strong cubic nonlinearity,” Journal of Sound and Vibration, vol. 285, no. 4-5, pp. 1171–1179, 2005.

J. H. He, Non-Perturbative Methods for Strongly Nonlinear Problems, Die Deutsche Bibliothek, Berlin, Germany, 2006.

J.-H. He, “Some asymptotic methods for strongly nonlinear equations,” International Journal of Modern Physics B, vol. 20, no. 10, pp. 1141–1199, 2006.

M. A. Noor and S. T. Mohyud-Din, “An efficient algorithm for solving fifth-order boundary value problems,” Mathematical and Computer Modelling, vol. 45, no. 7-8, pp. 954–964, 2007.

M. El-Shahed, “Application of He's homotopy perturbation method to Volterra's integro-differential equation,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 6, no. 2, pp. 163–168, 2005.

J. H. He, “Analytical methods for thermal science-an elementary introduction,” Thermal Science, vol. 15, supplement 1, pp. S1–S3, 2011.

H. E. Ji-Huan, “A note on the homotopy perturbation method,” Thermal Science, vol. 14, no. 2, pp. 565–568, 2010.

M. S. H. Chowdhury, I. Hashim, and S. Momani, “The multistage homotopy-perturbation method: a powerful scheme for handling the Lorenz system,” Chaos, Solitons and Fractals, vol. 40, no. 4, pp. 1929–1937, 2009.

M. S. H. Chowdhury and I. Hashim, “Application of multistage homotopy-perturbation method for the solutions of the Chen system,” Nonlinear Analysis. Real World Applications, vol. 10, no. 1, pp. 381–391, 2009.

I. Hashim and M. S. H. Chowdhury, “Adaptation of homotopy-perturbation method for numeric-analytic solution of system of ODEs,” Physics Letters A, vol. 372, no. 4, pp. 470–481, 2008.

I. Hashim, M. S. H. Chowdhury, and S. Mawa, “On multistage homotopy-perturbation method applied to nonlinear biochemical reaction model,” Chaos, Solitons and Fractals, vol. 36, no. 4, pp. 823–827, 2008.