**Applied Mathematics**

Vol.06 No.03(2015), Article ID:54974,11 pages

10.4236/am.2015.63056

Implementation of the Homotopy Perturbation Sumudu Transform Method for Solving Klein-Gordon Equation

Amr M. S. Mahdy^{1}, Adel S. Mohamed^{1}, Ahmad A. H. Mtawa^{2}

^{1}Department of Mathematics, Faculty of Science, Zagazig University, Zagazig, Egypt

^{2}Department of Mathematics, Faculty of Science, Benghize University, Almarj, Libya

Email: amr-mahdy85@yahoo.com, 3adel@live.nl, hussanahmad65@yahoo.com

Copyright © 2015 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

Received 16 February 2015; accepted 23 March 2015; published 24 March 2015

ABSTRACT

This paper extends the homotopy perturbation Sumudu transform method (HPSTM) to solve linear and nonlinear fractional Klein-Gordon equations. To illustrate the reliability of the method, some examples are presented. The convergence of the HPSTM solutions to the exact solutions is shown. As a novel application of homotopy perturbation sumudu transform method, the presented work showed some essential difference with existing similar application four classical examples also highlighted the significance of this work.

**Keywords:**

Mittag-Leffler Functions, Caputo Derivative, Sumudu Transform, Homotopy Perturbation Method, Klein-Gordon Equation

1. Introduction

Nonlinear phenomena that appear in many areas of scientific fields such as solid state physics, plasma physics, fluid dynamics, mathematical biology and chemical kinetics are modeled in terms of nonlinear partial differen- tial equations and in many scientific and engineering applications one of the corner stones of modeling are partial differential equations. For example, the Klein-Gordon equation which is of the form

(1)

with initial conditions

(2)

appears in modeling of problems in quantum field theory, relavistic physics, dispersive wave phenomena, plasma physic, nonlinear optics and applied physical sciences. The complexity of the equations though requires the use of numerical and analytical methods in most cases. A broad class of analytical solution and numerical solution methods were used to handle these problems. The topic of fractional partial differential equations has attracted a great atteation in the recent years. There are several analytical have been presented in the literature to solve fractional partial differential equations (FPDEs), such as the Fourier transform method [1] , the fractional Greens function method [2] , the Mellin transform method and the Laplace transform method [3] [4] , the Su- mudu transform method [5] .

Recently, several numerical methods have been introduced for this purpose, such as: the homotopy pertur- bation method (HPM) has first proposed by He [6] - [8] , the Modified homotopy perturbation method (MHPM) [9] , the differential transform method (DTM) [10] , the variational iteration method (VIM) [11] [12] , the ho- motopy analysis method (HAM) [13] [14] , the Sumudu decomposition method [15] , the Adomian decomposi- tion method [16] [17] .

The homotopy perturbation method (HPM) is extended to drive the exact solutions for linear (nonlinear) ordinary (partial) differential equations of fractional order. The homotopy perturbation method is also combined with the vartional iteration method [18] , to produce ahighly effective technique for handling many nonlinear problems. An also the homotopy perturbation method (HPM) is also combined with the laplace transform me- thod [19] . The advantage of this methods for obtaining exact and approximate solutions for nonlinear equa- tions.

The homotopy perturbation method (HPM) was also investigated by many researchers to handle partial differential equations arising in science and engineering [20] [21] . In addition, some numerical methods use a combination of utilizing specific transformation and obtaining series with converge to the exact solutions. An example of such a method is homotopy analysis Sumudu transform method (HASTD) which is a combination of the homotopy analysis method and the Sumudu transformation method [22] . Another such a combination is the which is the Sumudu decomposition method (SDM), which is constructed by combining two powerful methods, namely, the Sumudu transform method and Adomian decomposition method [23] . An efficent such approach is proposed combining the Sumudu transformation method with the homotopy perturbation method, which gives a new method called the homotopy perturbation Sumudu transform method (HPSTM) [24] . Recently, the ho- motopy perturbation Sumudu transform method (HPSTM) is frequently used for solving linear and nonlinear equations which are PDEs of integer order to obtain the exact solution.

In this paper, we applied homotopy perturbation Sumudu transform method (HPSTM) to obtain the analytical exact and approximate solutions for the fractional Klein-Gordon equation with time-fractional derivatives of the form:

(3)

where is parameters describing the order of the time fractional derivatives of, respectively, and they setisfy, is constants, and is the initial conditions and is is the source term and try to show the convergence of homotopy perturbation Sumudu transform method in solving this equation.

The paper is organized as follows: in Section 2, we recall some definitions of fractional calculus theory. In Section 3, we describe the homotopy perturbation Sumudu transform method. In Section 4, contains the main results and an examples to show the efficiency of using HPSTM to solve fractional-time Klein-Gordon equa- tions. Conclusions are given in Section 5.

2. Basic Definitions of Fractional Calculus

In this section, we mention the following basic definitions and properties of the fractional calculus theory and Sumudu transform.

Definition 1 The Riemann-Liouville fractional integral operator of order, of a function, , is defined as:

(4)

Definition 2 The fractional derivative of in the Caputo sense is defined as [13] [18]

(5)

for, , and is the Gamma function.

Definition 3 The Mittag-Leffler function which is ageneralization of exponential function (see [25] ) is de- fined as:

(6)

(7)

Some special cases of the Mittag-Leffler function are as follows:

1.

2.

Definition 4 The Sumudu transform is defined over the set of functions:

(8)

by the following formula:

(9)

Some special properties of the Sumudu transform are as follows:

1.;

2.;;

Other properties of the Sumudu transform can be found in [26] .

Definition 5 The Sumudu transform of the Caputo fractional derivative is defined as follows [5] :

(10)

3. The Homotopy Perturbation Sumudu Transform Method (HPSTM)

To illustrate the basic idea of this method, we consider a general fractional partial differential equation with the initial condition of the form:

(11)

with, and subject to the initial condition

(12)

where is the Caputo fractional derivative of the function, is the source term, is the linear differential operator and is the general nonlinear differential operator.

Applying the Sumudu transform (denoted in this paper by) on both sides of Equation (11), we get

(13)

Using the differentiation property of the Sumudu transform and the initial conditions in Equation (12), we have

(14)

Operating with the Sumudu inverse on both sides of Equation (14) gives

(15)

where represent the prescribed initial conditions. Now we apply the HPM.

(16)

and the nonlinear term can be decomposed as

(17)

for some Adomian’s polynomials that are given by [27]

(18)

Substituting Equation (16) and Equation (17) in Equation (15), we get

(19)

Equating the terms with identical powers of, we can obtain a series of equations as the follows:

(20)

proceeding in the same manner, the rest of the components can be completely found and the series solution is thus entirely determined. We approximate the analytical solution by truncated series as:

(21)

4. Applications

In this section, in order to asses the applicability and the accuracy of the fractional homotopy Sumudu transform method the following four examples.

Example 1 Consider the time-fractional partial differential Klein-Gordon equation

(22)

subject to the initial conditions

(23)

Taking the Sumudu transform on both sides of Equation (22), thus we get

and

Using the property of the Sumudu transform and the initial condition in Equation (23), we have

(24)

Operating with the Sumudu inverse on both sides of Equation (24) we get

(25)

By applying the homotopy perturbation method, and substituting Equation (16) in Equation (25) we have

(26)

Equating the terms with identical powers of, we get

Thus the solution of Equation (22) is given by

(27)

If we put in Equation (27) or solve Equations (22) and (23) with, we obtain the exact solu- tion

Which is in full agreement with the result in Reference [28] .

Example 2 Consider the inhomogeneous linear time-fractional partial differential Klein-Gordon equation

(28)

subject to the initial conditions

(29)

Taking the Sumudu transform on both sides of Equation (28), thus we get

and

Using the property of the Sumudu transform and the initial condition in Equation (29), we have

(30)

Operating with the Sumudu inverse on both sides of Equation (30) we get

(31)

By applying the homotopy perturbation method, and substituting Equation (16) in Equation (31) we have

(32)

Equating the terms with identical powers of, we get

Thus the solution of Equation (36) is given by

(33)

If we put in Equation (33) or solve Equations (28) and (29) with, we obtain the exact so- lution

Which is in full agreement with the result in Reference [28] .

Example 3 Consider the non-linear time-fractional partial differential Klein-Gordon equation

(34)

subject to the initial conditions

(35)

Taking the Sumudu transform on both sides of Equation (34), thus we get

and

Using the property of the Sumudu transform and the initial condition in Equation (35), we have

(36)

Operating with the Sumudu inverse on both sides of Equation (36) we get

(37)

By applying the homotopy perturbation method, and substituting Equations (16) in (37) we have

(38)

Equating the terms with identical powers of, we get

Thus the solution of Equation (34) is given by

(39)

If we put in Equation (39) or solve Equations (34) and (35) with, and so on, we can find that

we obtain the exact solution

Which is in full agreement with the result in Reference [28] .

Example 4 Consider the one-dimensional linear inhomogeneous fractional Klein-Gordon equation

(40)

subject to the initial conditions

(41)

Taking the Sumudu transform on both sides of Equation (40), thus we get

and

Using the property of the Sumudu transform and the initial condition in Equation (41), we have

(42)

Operating with the Sumudu inverse on both sides of Equation (42) we get

(43)

By applying the homotopy perturbation method, and substituting Equation (16) in Equation (43) we have

(44)

Equating the terms with identical powers of, we get

Thus the solution of Equation (40) is given by

(45)

If we put in Equation (45) or solve Equations (40) and (41) with, we obtain the exact so- lution

Which is in full agreement with the result in Reference [29] .

As it is presented above in Example 4 we obtained homotopy perturbation Sumudu transform solution of Equation (40) for values of, ,. Figures 1-4 show the approximate solutions for Equation (40) obtained for the three different values of using the homotopy perturbation Sumudu transform

method (HPSTM). The values of is the only case for which we know the exact solution and the results of (HPSTM) are in excellent agreement with the exact solution.

Figure 1. Profiles of w(x, t) when α = 2: Exact solution of (40(.

Figure 2. Profiles of w(x, t) when α = 2: Approximate solution of (40(.

Figure 3. Profiles of w(x, t) when α = 1.5: Approximate solution of (40).

Figure 4. Profiles of w(x, t) when α = 1.75: Approximate solution of (40).

5. Conclusion

In this paper, we have introduced a combination of the homotopy perturbation method and the Sumudu transform method for time fractional problems. This combination builds a strong method called the HPSTD. This method has been successfully applied to one-dimensional fractional equations and also for problems of linear and nonlinear partial differential equations. The HPSTD is an analytical method and runs by using the initial conditions only. Thus, it can be used to solve equations with fractional and integer order with respect to time. An important advantage of the new approach is its low computational load.

References

- Magin, R.L. and Ovadia, M. (2008) Modeling the Cardiac Tissue Electrode In-Terface Using Fractional Calculus. Journal of Vibration and Control, 14, 1431-1442. http://dx.doi.org/10.1177/1077546307087439
- Mainardi, F. (1995) Fractional Diffusive Waves in Viscoelastic Solids. In: Wegner, J.L. and Norwood, F.R., Eds., Nonlinear Waves in Solids, ASME Book No. AMR 137, Fairfield, 93-97.
- Odibat, Z. and Momani, S. (2007) A Reliable Treatment of Homotopy Perturbation Method for Klein-Gordon Equations. Physics Letters A, 365, 351-357. http://dx.doi.org/10.1016/j.physleta.2007.01.064
- Podlubny, I. (1999) Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Equations, to Methods of Their Solution and Same of Their Applications. Academic Press, New York.
- Gupta, V.G. and Sharma, B. (2010) Application of Sumudu Transform in Reaction-Diffusion Systems and Nonlinear Waves. Applied Mathematical Sciences, 4, 435-446.
- He, J.H. (1998) Approximate Analytical Solution for Seepage Flow with Fractional Derivatives in Porous Media. Computer Methods in Applied Mechanics and Engineering, 167, 57-68. http://dx.doi.org/10.1016/S0045-7825(98)00108-X
- He, J.H. (2005) Limit Cycle and Bifurcation of Nonlinear Problems. Chaos, Solutions and Fractals, 24, 827-833. http://dx.doi.org/10.1016/j.chaos.2005.03.007
- He, J.H. (1997) A New Approach to Nonlinear Partial Differential Equations. Communications in Nonlinear Science and Numerical Simulation, 2, 230-235. http://dx.doi.org/10.1016/S1007-5704(97)90007-1
- Hesameddini, E. and Latifzadeh, H. (2011) An Optimal Choice of Initial Solutions in the Homotopy Perturbation Method. International Journal of Nonlinear Sciences and Numerical Simulation, 10, 1389-1398.
- Miller, K.S. and Ross, B. (1993) An Introduction to the Fractional Calculus and Fractional Differential Equations. John Wiley and Sons, New York.
- Bhalekar, S. and Daftardar-Gejji, V. (2008) New Iterative Method: Application to Partial Differential Equations. Applied Mathematics and Computation, 203, 778-783. http://dx.doi.org/10.1016/j.amc.2008.05.071
- Daftardar-Gejji, V. and Bhalekar, S. (2010) Solving Fractional Boundary Value Problems with Dirichlet Boundary Conditions Using a New Iterative Method. Computers & Mathematics with Applications, 59, 1801-1809. http://dx.doi.org/10.1016/j.camwa.2009.08.018
- Arafa, A.A.M., Rida, S.Z. and Mohamed, H. (2011) Homotopy Analysis Method for Solving Biological Population Model. Communications in Theoretical Physics, 56, 797-800.
- Hilfe, R., Ed. (2000) Applications of Fractional Calculus in Physics. World Scientific, Singapore City.
- Eltayeb, H. and Kilicman, A. (2012) Application of Sumudu Decomposition Method to Solve Nonlinear System of Partial Differential Equations. Hindawi Publishing Corporation, Abstract and Applied Analysis, 2012, Article ID: 412948. http://dx.doi.org/10.1155/2012/412948
- Adomian, G. (1994) Solving Frontier Problems of Physics: The Decomposition Method. Kluwer Academic Publishers, Boston and London. http://dx.doi.org/10.1007/978-94-015-8289-6
- Cheng, J.F. and Chu, Y.M. (2011) Solution to the Linear Fractional Differential Equation Using Adomian Decomposition Method. Mathematical Problems in Engineering, 2011, 1-14. http://dx.doi.org/10.1155/2011/587068
- Noor, M.A. and Mohyud-Din, S.T. (2008) Variational Homotopy Perturbation Method for Solving Higher Dimensional Initial Boundary Value Problems. Mathematical Problems in Engineering, 2008, Article ID: 696734. http://dx.doi.org/10.1155/2008/696734
- Kumar, S., Yildirin, A. and Wei, L. (2012) A Fractional Model of the Diffusion Equation and Its Analytical Solution Using Laplace Transform. Scientia Iranica, 19, 1117-1123. http://dx.doi.org/10.1016/j.scient.2012.06.016
- Moustafa, O.L. (2003) On the Cauchy Problem for Some Fractional Order Partial Differential Equations. Chaos, Solitons & Fractals, 18, 135-140. http://dx.doi.org/10.1016/S0960-0779(02)00586-6
- Rafei, M. and Ganji, D.D. (2006) Explicit Solutions of Helmholtz Equation and Fifth-Order KdV Equation Using Homotopy Perturbation Method. International Journal of Nonlinear Sciences and Numerical Simulation, 7, 321-329.
- Rathore, S., Kumar, D., Singh, J. and Gupta, S. (2012) Homotopy Analysis Sumudu Transform Method for Nonlinear Equations. International Journal of Industrial Mathematics, 4, 301-314.
- Kumar, D., Singh, J. and Rathore, S. (2012) Sumudu Decomposition Method for Nonlinear Equations. International Mathematical Forum, 7, 515-521.
- Karbalaie, A., Montazeri, M.M. and Muhammed, H.H. (2014) Exact Solution of Time-Fractional Partial Dioerential Equations Using Sumudu Transform. WSEAS Transactions on Mathematics, 13, 142-151.
- Kilbas, A.A., Saigo, M. and Saxena, R.K. (2004) Generalized Mittag-Leffler Function and Generalized Fractional Calculus Operators. Integral Transforms and Special Functions, 15, 31-49. http://dx.doi.org/10.1080/10652460310001600717
- Belgacem, F.B.M. and Karaballi, A.A. (2006) Sumudu Transform Fundamental Properties Investigations and Application. Journal of Applied Mathematics and Stochastic Analysis, 2006, 1-23. http://dx.doi.org/10.1155/JAMSA/2006/91083
- Ghorbani, A. (2009) Beyond Adomian Polynomials: He Polynomials. Chaos, Solitons & Fractals, 39, 1486-1492. http://dx.doi.org/10.1016/j.chaos.2007.06.034
- Marasi, H.R. and Karimi, S. (2014) Convergence of the Variational Iteration Method for Solving Fractional Klein- Gordon Equation. Journal of Mathematical and Computational Science, 4, 257-266.
- Turut, V. and Güzel, N. (2013) On Solving Partial Differential Equations of Fractional Order by Using the Variational Iteration Method and Multivariate Padé Approximations. European Journal of Pure and Applied Mathematics, 6, 147-171.