7.2
CiteScore
3.7
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
ABUNDANCE ESTIMATION IN AN ARID ENVIRONMENT
Case Study
Editorial
Invited review
Letter to the Editor
Original Article
REVIEW
Review Article
SHORT COMMUNICATION
7.2
CiteScore
3.7
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
ABUNDANCE ESTIMATION IN AN ARID ENVIRONMENT
Case Study
Editorial
Invited review
Letter to the Editor
Original Article
REVIEW
Review Article
SHORT COMMUNICATION
View/Download PDF

Translate this page into:

ORIGINAL ARTICLE
24 (
2
); 187-191
doi:
10.1016/j.jksus.2010.08.002

Coupling of homotopy perturbation and modified Lindstedt–Poincaré methods for traveling wave solutions of the nonlinear Klein–Gordon equation

HITEC University Taxila Cantt, Pakistan
Ege University, Department of Mathematics, 35100 Bornova, İzmir, Turkey
Department of Mathematics, Fırat University, Turkey

*Corresponding author. Tel.: +92 333 5151290 syedtauseefs@hotmail.com (Syed Tauseef Mohyud-Din)

Disclaimer:
This article was originally published by Elsevier and was migrated to Scientific Scholar after the change of Publisher.

Available online 10 August 2010

Abstract

oupled method of He’s homotopy perturbation method and modified Lindstedt–Poincaré method is applied to the search for traveling wave solutions of a variety of Klein–Gordon equations. The results obtained provide confirmation for the validity of the coupled method.

Keywords

Klein–Gordon equation
Homotopy perturbation method
modified Lindstedt–Poincaré method
1

1 Introduction

Periodic wave-trains are found in many physical systems and, they are fundamental solutions related to the elementary solutions in the form of sinusoidal wave-trains in linear theory. In nonlinear theory, the solutions are no longer sinusoidal, but periodic solutions may still exist. For this class of periodic solutions, the main effect of non-linearity becomes visible on amplitude dependence in the dispersion relation. This, off course, leads to a new qualitative behavior, not merely to the correction of the dispersion relation in linear case. In this article, we will consider after Whitham (1965) nonlinear, Klein–Gordon, equation governed by

(1)
u tt - u xx + V ( u ) = 0 , where V(u) is any nonlinear potential function which yields oscillatory solutions and V′(u) is chosen as the derivative of a potential energy. Klein–Gordon equation is not only a useful model but also arises in variety of physical situations (Knobel, 2000; Shen, 1994). It is worth to note that the dispersion relation of a periodic wave-train in (1) is dependent on its amplitude. For this class of problems with waves of small amplitude, perturbation methods based on small amplitude expansions infer the existence of the periodic wave-trains (Nayfeh, 1973). However, the perturbation methods, in principle, works only for the nonlinear problems with small parameters. In addition, analytical solution obtained by the perturbation method has a little range of validity in most cases. Consequently, new approximate-analytical techniques are considered necessary. For example, the sine-Gordon and the double sine-Gordon equations, the sinh-Gordon equation and the double sinh-Gordon equations were investigated by using the standard tanh method (Malfliet, 1992, 1996a,b; Wazwaz, 2005). Sirendaoreji and Sun Jiong (2002) applied a direct method for solving Sinh-Gordon type equation. Fu et al. (2004) have obtained exact solutions for double and triple sinh-Gordon equations. Wazwaz (2005) obtained traveling wave solutions for combined and double combined sine-cosine-Gordon equations by the variable separated ODE method.

More recently, Lim et al. (2001), by coupling linearization of Klein–Gordon equation with the method of harmonic balance, established two general analytical approximate formulas for the dispersion relation which depends on the amplitude of the wave train. In this paper, we will be coupling He’s homotopy perturbation method (He, 2000; Öziş and Yıldırım, 2006) and modified Lindstedt–Poincaré method (He, 2002a,b, 2006; Öziş and Yıldırım, 2007, 2006; Liu, 2005) to obtain the periodic wave-trains in (1) and will use three examples to illustrate the applicability and the effectiveness of the proposed method. For latest developments in this field, the reader is referred to see (He et al., 2006; Noor and Mohyud-Din, 2008; Mohyud-Din et al., 2009a,b,c,d, 2010 Mohyud-Din, 2009; Mohyud-Din and Noor, 2007, 2009; Abdou, 2010a, 2009 Abdou et al., in press; El-Wakil and Abdou, 2008, 2010; El-Wakil et al., in press; He, 1999) and the references therein.

2

2 Brief concept of coupling of homotopy perturbation and modified Lindstedt–Poincaré methods

Recognizing He’s homotopy perturbation method; the homotopy with and imbedding parameter p ∈ [0, 1] is constructed, and the imbedding parameter is considered as a “small parameter”, so the method is called homotopy perturbation method and is proceed as the standard perturbation method but taking the full advantage of the traditional perturbation methods and the homotopy techniques. The main merit of the homotopy perturbation method is that the perturbation equation can be easily constructed (therefore is problem dependent) by homotopy in topology and the initial approximation can also be freely selected. On the other hand, in He’s modified Lindstedt–Poincaré method the coefficient of the second term, i.e., u is also expanded into a series besides the assumed solution. For further reading, refer to the comprehensive book by He et al. (2006) and the references therein. To our view, if He’s modified Lindstedt–Poincaré method is applicable to the perturbed equation constructed by He’s homotopy perturbation method then the obtained solution would be exceedingly accurate and may possibly be applicable to wide range of physical systems. For better illustration of coupled method of He’s homotopy perturbation method and modified Lindstedt–Poincaré method, and making the underlying idea clear, we demonstrate three examples. By doing so, we will try to validate the applicability, accuracy and effectiveness of the proposed method.

Example 3.1

Consider the Klein–Gordon equation governed by

(2)
u tt - α 2 u xx + γ 2 u = β u 3 We first connect the independent variables x and t (Malfliet, 1992, 1996a,b) into one wave variable by
(3)
ξ = x - ct .
Substituting the transformation (3) into Eq. (2) gives
(4)
u + γ 2 ( c 2 - α 2 ) u - β ( c 2 - α 2 ) u 3 = 0 , ( c 2 α 2 )
where prime denotes differentiation with respect to ξ. Physical considerations show that Eq. (4) has a periodic solution. In order to look for the periodic solution, we construct the following homotopy:
(5)
u + γ 2 ( c 2 - α 2 ) u - β ( c 2 - α 2 ) pu 3 = 0 ,
It is apparent that when p = 0, (5) reduces to linear equation and when p = 1, it becomes the original nonlinear one. Therefore, the embedding parameter p monotonically increases from zero to unit as the linear operator is continuously deformed to the nonlinear problem given. We assume that the periodic solution to Eq. (5) may be written as a power series in p:
(6)
u = u 0 + pu 1 + p 2 u 2 +
If higher order approximate solution is required, the He’s modified Lindstedt–Poincaré method can be applied. Hence, we expand the coefficient of the linear term into a series of p:
(7)
γ 2 ( c 2 - α 2 ) = ω 2 + p ω 1 + p 2 ω 2 +
Substituting (6) and (7) into Eq. (5), and processing as the standard perturbation method, we have:
(8)
u 0 + ω 2 u 0 = 0 ,
(9)
u 1 + ω 2 u 1 + ω 1 u 0 - β ( c 2 - α 2 ) u 0 3 = 0
The initial approximation u0 can be freely chosen, hereby we set
(10)
u 0 = A cos ω ξ ,
which satisfies the Eq. (8). The substitution of (10) in (9) yields
(11)
u 1 + ω 2 u 1 + ω 1 A cos ω ξ - β ( c 2 - α 2 ) A 3 cos 3 ω ξ = 0
or
(12)
u 1 + ω 2 u 1 + A cos ω ξ ω 1 - 3 β A 2 4 ( c 2 - α 2 ) - β A 3 4 ( c 2 - α 2 ) cos 3 ω ξ = 0
No secular term in u1 requires that
(13)
ω 1 = 3 β A 2 4 ( c 2 - α 2 ) .
Eq. (12) becomes
(14)
u 1 + ω 2 u 1 - β A 3 4 ( c 2 - α 2 ) cos 3 ω ξ = 0
We write down its special solution:
(15)
u 1 = - β A 3 32 ω 2 ( c 2 - α 2 ) ( cos 3 ω ξ - cos ω ξ )
Setting p = 1, in (7) gives
(16)
γ 2 ( c 2 - α 2 ) = ω 2 + 3 β A 2 4 ( c 2 - α 2 ) .
Hence, the frequency is:
(17)
ω = 4 γ 2 - 3 β A 2 4 ( c 2 - α 2 ) .
If the first-order approximate solution is adequate, setting p = 1, we have
(18)
u = u 0 + u 1 = A cos ω ξ - β A 3 32 ω 2 ( c 2 - α 2 ) ( cos 3 ω ξ - cos ω ξ )
which is the solution of Klein–Gordon equation in (2).

Example 3.2

As a second example, consider the sine-Gordon equation governed by

(19)
u tt - ku xx + 2 α sin u = 0 Using the transformation (3), Eq. (19) gives:
(20)
( c 2 - k ) u + 2 α sin u = 0 .
Approximating sin u u - u 3 6 , Eq. (20) reduces to:
(21)
u + 2 α ( c 2 - k ) u - u 3 6 = 0 , ( c 2 k )
We, in this case, construct the homotopy in the form:
(22)
u + 2 α ( c 2 - k ) u - p α 3 ( c 2 - k ) u 3 = 0 .
Supposing the solution in the form of (6) and if the coefficient of second term of (22) can be expressed as
(23)
2 α ( c 2 - k ) = ω 2 + p ω 1 + p 2 ω 2 +
respectively, and substituting (6) and (23) into (22), collecting terms of the same powers of p, we have
(24)
u 0 + ω 2 u 0 = 0 ,
(25)
u 1 + ω 2 u 1 + ω 1 u 0 - α 3 ( c 2 - k ) u 0 3 = 0
We start initial approximation with u0 = Acosωξ, as in Example 1, hence Eq. (25) yields
(26)
u 1 + ω 2 u 1 + ω 1 A cos ω ξ - α 3 ( c 2 - k ) ( A cos ω ξ ) 3 = 0 .
Making simplification, we obtain:
(27)
u 1 + ω 2 u 1 + A cos ω ξ w 1 - α A 2 4 ( c 2 - k ) - α A 3 12 ( c 2 - k ) cos 3 ω ξ = 0 .
Again, elimination of the secular terms requires
(28)
ω 1 = α A 2 4 ( c 2 - k )
and Eq. (27) reduces to
(29)
u 1 + ω 2 u 1 - α A 3 12 ( c 2 - k ) cos 3 ω ξ = 0 .
We obtain a particular solution of Eq. (29), which reads
(30)
u 1 = - α A 3 96 ω 2 ( c 2 - k ) ( cos 3 ω ξ - cos ω ξ )
Setting p = 1, in (23) gives
(31)
2 α ( c 2 - k ) = ω 2 + α A 2 4 ( c 2 - k )
and therefore, the frequency is
(32)
ω = 8 α - α A 2 4 ( c 2 - k )
If, for example, first-order approximate solution is adequate, setting p = 1, we have
(33)
u = u 0 + u 1 = A cos ω ξ - α A 3 96 ω 2 ( c 2 - k ) ( cos 3 ω ξ - cos ω ξ )
which is the solution of the sine-Gordon equation in (19).

Example 3.3

Consider, now, combined sine-cosine-Gordon equation:

(34)
u tt - ku xx + α sin u + β cos u = 0 Using the transformation (3) again, Eq. (34) gives
(35)
( c 2 - k ) u + α sin u + β cos u = 0 .
Approximating the above equation by sin u u - u 3 6 and cos u 1 - u 2 2 , we obtain following equation:
(36)
u + α ( c 2 - k ) u - u 3 6 + β ( c 2 - k ) 1 - u 2 2 = 0 , ( c 2 k )
We, in a similar manner, construct the homotopy in the form:
(37)
u + α ( c 2 - k ) u + p β ( c 2 - k ) - β 2 ( c 2 - k ) u 2 - α 6 ( c 2 - k ) u 3 = 0 .
Supposing the solution in the form of (6) and if the coefficient of second term of (37) can be expressed as
(38)
α ( c 2 - k ) = ω 2 + p ω 1 + p 2 ω 2 +
respectively, and substituting (6) and (38) into (37), collecting terms of the same powers of p, we have
(39)
u 0 + ω 2 u 0 = 0 ,
(40)
u 1 + ω 2 u 1 + ω 1 u 0 + β ( c 2 - k ) - β 2 ( c 2 - k ) u 0 2
(41)
- α 6 ( c 2 - k ) u 0 3 = 0 .
We start initial approximation, again, with u0 = Acosωξ, as in Examples 2 and 3, hence (41) yields
(42)
u 1 + ω 2 u 1 + ω 1 A cos ω ξ + β ( c 2 - k ) - β 2 ( c 2 - k ) ( A cos ω ξ ) 2 - α 6 ( c 2 - k ) ( A cos ω ξ ) 3 = 0 .
Making simplification, we obtain:
(43)
u 1 + ω 2 u 1 + A cos ω ξ ω 1 - α A 2 8 ( c 2 - k ) + 4 β - β A 2 4 ( c 2 - k ) - β A 2 4 ( c 2 - k ) cos 2 ω ξ - α A 3 24 ( c 2 - k ) cos 3 ω ξ = 0
No secular term in u1 requires that ω 1 = α A 2 4 ( c 2 - k ) and A = 2 and (43) reduces to
(44)
u 1 + ω 2 u 1 - β A 2 4 ( c 2 - k ) cos 2 ω ξ - α A 3 24 ( c 2 - k ) cos 3 ω ξ = 0
We obtain a particular solution of Eq. (44), which reads
(45)
u 1 = - 16 β A 2 + α A 3 192 ω 2 ( c 2 - k ) cos ω ξ + β A 2 12 ω 2 ( c 2 - k ) cos 2 ω ξ + α A 3 192 ω 2 ( c 2 - k ) cos 3 ω ξ
Setting p = 1, in (38) gives
(46)
α ( c 2 - k ) = ω 2 + α A 2 8 ( c 2 - k )
and the frequency is
(47)
ω = 8 α - α A 2 8 ( c 2 - k )
If, for example, first-order approximate solution is adequate, setting p = 1, we have
(48)
u = u 0 + u 1 = A cos ω ξ - 16 β A 2 + α A 3 192 ω 2 ( c 2 - k ) cos ω ξ + β A 2 12 ω 2 ( c 2 - k ) cos 2 ω ξ + α A 3 192 ω 2 ( c 2 - k ) cos 3 ω ξ

3

3 Conclusion

In this paper, coupling of He’s homotopy perturbation and modified Lindstedt–Poincaré methods is applied to solve a variety of Klein–Gordon equations. The advantage of the approach is that it does not need a small parameter in the physical system, leading to wide application in nonlinear wave equations. Moreover, the method is capable of significantly minimizing the size of computational labor compared to other existing techniques. The obtained results are entirely new.

References

  1. , . Approximate solution of nonlinear difference differential equations. Int. J. Comput. Methods. 2009;6(4):1-15.
    [Google Scholar]
  2. Abdou, M.A., in press. Explicit series solutions of nonlinear physical problems via the homotopy analysis method. Inter. J. Modern Phys. B
  3. Abdou, M.A., in press. On the use of homotopy analysis method for solving heat transfer in MHD viscoelastic boundary layer flow. Inter. J. Modern Phys. B.
  4. , . New applications of He’s homotopy perturbation method for nonlinear difference differential equations. Phys. Scripta. 2010;81:015003. 1–8
    [Google Scholar]
  5. Abdou, M.A., in press. Analytical approximate solutions for Zakharov Kuznetsov equation in ion acoustic waves and shallow water equation with source term. Int. J. Modern Phys. B.
  6. , , . New applications of the homotopy analysis method. Z. Naturforsch.. 2008;63A:1-8.
    [Google Scholar]
  7. , , . Solution of coupled system of nonlinear physical problems using homotopy analysis method. Phys. Scripta. 2010;81:015001. 1–8
    [Google Scholar]
  8. El-Wakil, S.A., Degedy, A.R., Abulwafa, E.M., Madkour, M.A. Abdou, M.A., in press. Maximum entropy homotopy perturbation method for solving the moment equations of Fokker–Planck equations. Int. J. Modern Phys. B.
  9. , , , . Exact solutions to double and triple sinh-Gordon equations. Z. Naturforsch.. 2004;59a:933-937.
    [Google Scholar]
  10. , . Homotopy perturbation technique. Comput. Meth. Appl. Mech. Eng.. 1999;178:257-262.
    [Google Scholar]
  11. , . A coupling method of homotopy technique and perturbation technique for nonlinear problems. Int. J. Nonlinear Mech.. 2000;35:37-43.
    [Google Scholar]
  12. , . Modified Lindstedt–Poincaré methods for some strongly non-linear oscillations Part I: expansion of a constant. Int. J. Non-linear Mech.. 2002;37:309-314.
    [Google Scholar]
  13. , . Recent developments in asymptotic methods for nonlinear ordinary equations. Int. J. Comput. Num. Anal. Appl.. 2002;2(2):127-190.
    [Google Scholar]
  14. , . Some asymptotic methods for strongly nonlinear equations. Int. J. Modern Phys. B. 2006;20(10):1141-1199.
    [Google Scholar]
  15. He, J.H., 2006. Non-Perturbative Methods for Strongly nonlinear Problems. Dissertation.de-verlag im Internet GmbH, Berlin.
  16. , . An Introduction to the Mathematical Theory of Waves. Providence, RI: American Mathematical Society; .
  17. , , , . new approximate analytical approach for dispersion relation of the nonlinear Klein–Gordon equation. Chaos. 2001;11(4):843-848.
    [Google Scholar]
  18. , . Approximate period of nonlinear oscillators with discontinuities by modified Lindstedt–Poincare method. Choas, Solitons and Fractals. 2005;23(2):577-579.
    [Google Scholar]
  19. , . Solitary wave solutions of nonlinear wave equations. Am. J. Phys.. 1992;60(7):650-654.
    [Google Scholar]
  20. , . The tanh method: I. Exact solutions of nonlinear evolution and wave equations. Phys. Scr.. 1996;54:63-568.
    [Google Scholar]
  21. , . The tanh method: II. Perturbation technique for conservative systems. Phys. Scr.. 1996;54:569-575.
    [Google Scholar]
  22. , . On the coupling of He’s polynomials and Laplace transformation. Int. J. Nonlinear Sci. Num. Simulat. 2009
    [Google Scholar]
  23. Mohyud-Din, S.T., Noor, M.A., 2007. Homotopy perturbation method for solving fourth-order boundary value problems. Mathematical Problems in Engineering, Hindawi 2007, pp. 1–15. Article ID 98602, doi:10.1155/2007/98602.
  24. , , . Homotopy perturbation method for solving partial differential equations. Z. Naturforsch. A – A J. Phys. Sci.. 2009;64a:157-170.
    [Google Scholar]
  25. , , , . Travelling wave solutions of seventh-order generalized KdV equations using He’s polynomials. Int. J. Nonlinear Sci. Numer. Simulat.. 2009;10(2):223-229.
    [Google Scholar]
  26. Mohyud-Din, S.T., Noor, M.A., Noor, K.I., 2009. Some relatively new techniques for nonlinear problems. Mathematical Problems in Engineering, Hindawi, 2009, Article ID 234849, 25 p., doi:10.1155/2009/234849.
  27. , , , . Solving second-order singular problems using He’s polynomials. World Appl. Sci. J.. 2009;6(6):769-775.
    [Google Scholar]
  28. , , , . Parameter-expansion techniques for strongly nonlinear oscillators. Int. J. Nonlinear Sci. Numer. Simulat.. 2009;10(5):581-583.
    [Google Scholar]
  29. , , , . Traveling wave solutions of Whitham–Broer–Kaup equations by homotopy perturbation method. J. King Saud University (Science) Elsevier. 2010;22:173-176.
    [Google Scholar]
  30. , . Perturbation Methods. New York: Wiley-Interscience; .
  31. , , . Variational iteration method for solving higher-order nonlinear boundary value problems using He’s polynomials. Int. J. Nonlinear Sci. Numer. Simulat.. 2008;9(2):141-157.
    [Google Scholar]
  32. Öziş, T., Yıldırım, A., 2006. A note on He’s homotopy perturbation method for van der Pol oscillator with very strong nonlinearity. Chaos Solitons and Fractals, doi:10.1016/j.chaos.2006.04.013.
  33. Öziş, T., Yıldırım, A., 2006. Determination of periodic solutions for a u1/3 Force by He’s modified Lindstedt–Poincaré method. J. Sound Vib., doi:10.1016/j.jsv.2006.10.001.
  34. , , . Determination of limit cycles by a modified straightforward expansion for nonlinear oscillators. Chaos Solitons and Fractals. 2007;32:445-448.
    [Google Scholar]
  35. , . A Course on Nonlinear Waves. Dordrecht: Kluwer; .
  36. , , . A direct method for solving sinh-Gordon type equation. Phys. Lett. A. 2002;298:133-139.
    [Google Scholar]
  37. Wazwaz, A.M., 2005. Traveling wave solutions for combined and double combined sine-cosine-Gordon equations by the variable separated ODE method. Appl. Math. Comput. doi:10.1016/j.amc.2005.09.104.
  38. , . The tanh method: exact solutions of the sine-Gordon and the sinh-Gordon equations. Appl. Math. Comput.. 2005;167(2):1196-1210.
    [Google Scholar]
  39. , . A general approach to linear and nonlinear waves using a Lagrangian. J. Fluid Mech.. 1965;22:273-283.
    [Google Scholar]
Show Sections