Advances in Mathematical Physics
Volume 2015 (2015), Article ID 932905, 5 pages
http://dx.doi.org/10.1155/2015/932905

Research Article

Comparative Solution of Nonlinear Quintic Cubic Oscillator Using Modified Homotopy Perturbation Method

Muhammad Suleman1,2 and Qingbiao Wu¹

¹Department of Mathematics, Zhejiang University, Hangzhou, China
²Department of Mathematics, Comsat Institute of Information Technology, Islamabad, Pakistan

Received 1 April 2015; Revised 14 May 2015; Accepted 17 May 2015

Academic Editor: Nikos Mastorakis

Copyright © 2015 Muhammad Suleman and Qingbiao Wu. 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.

We use modified homotopy perturbation method to find out the solution of nonlinear cubic quintic equation. Besides this method solution of the problem with the following methods is discussed, Energy Balance Method and He’s Frequency Formulation method and then compare the results with each other and Global Error Method. The results show that these three methods are effective as global error method for nonlinear cubic quintic oscillator equation with multiple nonlinear terms and have different effects on the solution. In particular, the homotopy perturbation solution is quite surprising. A cubic quintic nonlinear oscillator is used as an example to compare the results.

1. Introduction

In the present era large numbers of engineering problems are nonlinear and it is almost difficult to solve them analytically. Perturbation methods are among the known methods to solve the nonlinear problems, which are based on the existence of small/large parameters so-called perturbation parameters. These so-called perturbation parameters are not present in many nonlinear problems so we cannot use perturbation methods. To overcome this difficulty researchers introduced such methods which do not require small parameters such as homotopy perturbation method (HPM) [1, 2], book keeping parameter perturbation method [3], He’s energy balance method (EBM) [4–7], global error minimization method [8], Variational Iterational method [4], and many others. In this paper, solution of nonlinear cubic quintic oscillatory equation is solved with iterative methods, the basic idea of the homotopy perturbation is used, and its applications in some nonlinear oscillator equations are studied. In most engineering problems, it is easy to find angular frequency of a solution to a nonlinear equation. He presented the method for the first time. These results are valid for weakly nonlinear systems as well as strongly nonlinear systems. In Section 2, solution with energy balance method is explained. Solution with He’s frequency formulation method is explained in Section 3. In Section 4, HPM solution is explained and some conclusions are given in Section 5.

2. Energy Balance Method (EBM)

Consider the nonlinear quintic cubic equation

Its variational principle can be easily obtained:

Its Hamiltonian, therefore, can be written in the form

In (3) and (4), the kinetic energy () and potential energy () can be, respectively, expressed as throughout the oscillation; it holds that . We use the following trial function to determine the angular frequency . Let

Being the trial function, substituting (5) into (4), we obtain the following residual equation:

If, by chance, the exact solution had been chosen as the trial function, then it would be possible to make residual zero for all values of “” by appropriate choice of . Since (5) is only an approximation to the exact solution, cannot be made zero everywhere.

Collocation at gives

3. He’s Frequency Formulation Method

Equation (1) can be written in the general form

We use two trial functions and , which are, respectively, the solutions of the following linear oscillator equations:

where is assumed to be the frequency of the nonlinear oscillator. The residuals for (8) are

The original frequency-amplitude formulation reads [9]

We just used the following formulation:

Geng and Cai improved the formulation by choosing another location point [10]:

The accuracy depends upon the chosen location point, and we have no general rule for the choice.

To illustrate this shortcoming, we consider the cubic quintic nonlinear equation [8]:

The residual of the equations is

where frequency can be calculated at 0 or :

Now locating (where , , for (15) and for (16)), we obtain

It shows that the accuracy depends on the point of location.

3.1. The Method of Weighted Residual

We use the method of weighted residuals to overcome the shortcoming. To this end, we introduce two new residual variables and defined as

We can approximately calculate in the form

Now, by using (20) for frequency formulation, we find the value of for cubic quintic nonlinear oscillator:

Its approximate frequency reads

4. Modified Homotopy Perturbation Method with Two Expanding Parameters

Consider the following nonlinear cubic quintic oscillator [8]:

where and are constant. Equation (23) is exactly solvable and its solution is given in [11]. According to homotopy perturbation method [1, 2], homotopy equation can be constructed in the form

where is unknown frequency of oscillator to be determined. When , we have

This can be approximately describing the solution properties. Embedding parameter increases monotonically from zero to unit as the trivial problem. Equation (25) is continuously deformed to the original one. According to standard homotopy perturbation method [1, 2], solution is expanded in a series of :

Generally one iteration is enough to obtain an approximate solution with relatively high accuracy. If higher accuracy of solution is needed, homotopy equation is suggested in the following [1]:

Solution is also expanded in the same form of (26) and coefficient 1 of the linear term in (27) is also expanded in series of :

This method is also used in parameter-expansion method. In [6] a homotopy equation with an auxiliary term is suggested which vanishes for both cases with homotopy parameters and in constructed homotopy equation. A detailed solution process is given and a relatively comprehensive survey on concepts, theory, and applications of homotopy perturbation method is reported [1–3, 12]. In this paper, we suggest a homotopy perturbation method with two expanding parameters.

4.1. Solution of Problem Using HPM

In (23) there are two nonlinear terms, which might have different effects on the solution. In order to find the best match with different nonlinear terms, we can construct a homotopy equation in the form

where and are homotopy parameters, , . Solution and coefficient 1 of linear term in (29) are expanded in the form

This expansion is similar to modified Lindstedt-Poincare method with double series expansion [3]. Substituting (30) and (31) into (29), collecting the same power of (), and setting coefficients zero, we can readily obtain the following equations:

Initial conditions for () can be freely chosen, but it must hold that and for final closed approximate solution. Solving (32), we have

Substitution of into (33) results in

Requiring no secular term in from (39), we have

A particular solution for (39) is readily obtained as

Substituting , , and into (34) one gets

In order to remove a secular term from (42), consider

Now obtain a particular solution for (42):

Substituting the obtained results into (35), we obtain equation for and , respectively:

Similarly, we can find out the values of and by using (36) and (37):

The solution process continues in this way. Substituting the values of into (31) and using and , we have

According to Nayfeh’s result [13] for comparison

This result is obtained under the assumption of small amplitude. Consequently, the above methods give the same result as that of Nayfeh [13]. In case of being zero it will become the well-known Duffing equation and its nonlinear angular frequency can be obtained. The result is the same as that obtained in [13] and its maximal relative error is less than 7.5%.

5. Conclusions

From Table 1 it is suggested that these methods not only are effective for the given nonlinear cubic quintic oscillatory equation but also can be extended to various nonlinear problems with multiple nonlinear terms. In particular, the solution procedure of modified homotopy perturbation method with two expanding parameters can be used as a paradigm for many other applications.

Table 1: Comparison of HPM, He’s frequency, energy balance method, global error method, and exact solution.

(amplitude)
Modified HPM
Energy balance method
He’s frequency
Global error minimization method
Exact

0.2
1.0724
1.0826
1.0998
1.0726
1.0720

0.5
1.1954
1.4701
1.4764
1.2438
1.1846

1
1.9947
1.9254
1.9635
2.3262
1.3778

5
2.0017
2.5217
2.5321
2.1980
2.1504

10
2.3434
2.7423
2.8984
2.9463
2.8666

100
7.1033
8.1798
8.1975
8.8327
8.5336

Conflict of Interests

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

Acknowledgments

This work is supported by the National Natural Science Foundation of China (Grant no. 11371320), Zhejiang Natural Science Foundation (Grant no. LZ4A010002), Foundation of Science and Technology Department (Grant no. 2013C31084) of Zhejiang province, and Scientific Research Fund of Zhejiang Provincial Education Department (Grants nos. Y201431077 and Y2013329420).

References

J. H. He, “Homotopy perturbation method with two expanding parameters,” Indian Journal of Physics, vol. 88, no. 2, pp. 193–196, 2014.

J.-H. He, “Bookkeeping parameter in perturbation methods,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 2, no. 3, pp. 257–264, 2001.

H.-L. Zhang, “Periodic solutions for some strongly nonlinear oscillations by He's energy balance method,” Computers & Mathematics with Applications, vol. 58, no. 11-12, pp. 2480–2485, 2009.

N. Jamshidi and D. D. Ganji, “Application of energy balance method and variational iteration method to an oscillation of a mass attached to a stretched elastic wire,” Current Applied Physics, vol. 10, no. 2, pp. 484–486, 2010.

S. S. Ganji, D. D. Ganji, Z. Z. Ganji, and S. Karimpour, “Periodic solution for strongly nonlinear vibration systems by He's energy balance method,” Acta Applicandae Mathematicae, vol. 106, no. 1, pp. 79–92, 2009.

J.-H. He, “Homotopy perturbation method with an auxiliary term,” Abstract and Applied Analysis, vol. 2012, Article ID 857612, 7 pages, 2012.

J.-H. He, “Comment on He's frequency formulation for nonlinear oscillators,” European Journal of Physics, vol. 29, no. 4, 2008.

A. Kargar and M. Akbarzade, “Analytical solution of nonlinear cubic-quintic Duffing oscillator using global error minimization method,” Advanced Studies in Theoretical Physics, vol. 6, no. 9–12, pp. 467–471, 2012.

M. Gorji, D. D. Ganji, and S. Soleimani, “New application of He's homotopy perturbation method,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 8, no. 3, pp. 319–328, 2007.

L. Geng and X.-C. Cai, “He's frequency formulation for nonlinear oscillators,” European Journal of Physics, vol. 28, no. 5, pp. 923–931, 2007.

A. Elías-Zúñiga, “Exact solution of the cubic-quintic Duffing oscillator,” Applied Mathematical Modelling, vol. 37, no. 4, pp. 2574–2579, 2013.

M. Febbo, “A finite extensibility nonlinear oscillator,” Applied Mathematics and Computation, vol. 217, no. 14, pp. 6464–6475, 2011.

A. H. Nayfeh, Introduction to Perturbation Technique, John Wiley & Sons, New York, NY, USA, 1985.