﻿ Approximate Solution of Non-Linear Fractional Klein-Gordon Equation Using Spectral Collocation Method

Applied Mathematics
Vol.06 No.13(2015), Article ID:61584,7 pages
10.4236/am.2015.613190

Approximate Solution of Non-Linear Fractional Klein-Gordon Equation Using Spectral Collocation Method

Rubayyi T. Alqahtani

Department of Mathematics and Statistics, College of Science, Al-Imam Mohammad Ibn Saud Islamic University (IMSIU), Riyadh, Saudi Arabia   Received 1 September 2015; accepted 27 November 2015; published 30 November 2015 ABSTRACT

In this paper, we implement the spectral collocation method with the help of the Legendre polynomials for solving the non-linear Fractional (Caputo sense) Klein-Gordon Equation (FKGE). We present an approximate formula of the fractional derivative. The Legendre collocation method is used to reduce FKGE to the solution of system of ODEs which is solved by using finite difference method. The results of applying the proposed method to the non-linear FKGE show the simplicity and the efficiency of the proposed method.

Keywords:

Fractional Klein-Gordon Equation, Legendre Spectral Method 1. Introduction

The theory of fractional calculus is initiated by Leibniz, Liouville, Riemann, Grunwald and Letnikov and since then has been found many applications in science and engineering. Finding accurate and efficient method for solving fractional differential equations has been an active research subject. Finding the exact solution for most of these equations is not an easy task, thus analytical and numerical methods must be used.

The Klein-Gordon equation plays a significant role in mathematical physics and many scientific applications such as solid-state physics, nonlinear optics, and quantum field theory  . The equation has attracted much attention in studying solitons and condensed matter physics, in investigating the interaction of solitons in a collisionless plasma and the recurrence of initial states, and in examining the nonlinear wave equations  . Wazwaz has obtained the various exact traveling wave solutions such as compactons, solitons and periodic solutions by using the tanh method  . The study of numerical solutions of the Klein-Gordon equation has been investigated considerably in the last few years. In the previous studies, the most papers have carried out different spatial discretization of the equation (  ).

In this work, we apply spectral collocation method (with the help of Legendre polynomials) to obtain the numerical solution of the non-linear FKGE of the form (1)

where denotes the Caputo fractional derivative of order with respect to x, is unknown function, and , and are known constants with , .

We consider the initial conditions and the boundary conditions as follows: (2) (3)

The existence and the uniqueness of the solution of Equations (1)-(3) are given in ( - ).

For more details about the fractional calculus see ( - ) and for more details about the Legendre collocation method see ( - ).

2. An Approximate Formula of the Fractional Derivative

The well-known Legendre polynomials are defined on the interval and can be determined with the aid of the following recurrence formula In order to use these polynomials on the interval we define the so called shifted Legendre polynomials by introducing the change of variable . Let the shifted Legendre polynomials be denoted by. Then can be obtained as follows:

(4)

where and. The analytic form of the shifted Legendre polynomials of degree k is given by

(5)

The function, which is square integrable in, may be expressed in terms of shifted Legendre polynomials as

(6)

where the coefficients are given by

(7)

In practice, only the first -terms of shifted Legendre polynomials are considered. Then we have

(8)

Theorem 1  .

Let be approximated by shifted Legendre polynomials as (8) and also suppose then

(9)

Theorem 2  .

The error in approximating by is bounded by

(10)

where

3. Numerical Implementation

Example 1.

Consider the fractional-order cubically nonlinear Klein-Gordon problem

(11)

with the initial and boundary conditions

where the source term is given by

The exact solution of this problem is.

In order to use the proposed method, we approximate with as

(12)

From Equation (11) and Theorem 1 we have

(13)

We now collocate Equation (13) at points, as

(14)

For suitable collocation points we use roots of shifted Legendre polynomial.

In this case, the roots of shifted Legendre polynomial, i.e.

Also, by substituting Equation (12) in the boundary conditions we can find

(15)

By using Equation (14) and Equation (15) we obtain the following non-linear system of ODEs:

(16)

(17)

(18)

(19)

where

Now, to use FDM  for solving the system (16)-(19), we will use the following notations: to be the integration time, , for. Define, Then the system (16)-(19), is discretized and takes the following form:

(20)

(21)

(22)

(23)

This system presents the numerical scheme of the proposed problem and is non-linear system of algebraic equations, and by solving this system yields the numerical solution of the non-linear FKGE (11).

The obtained numerical results by means of the proposed method are shown in Table 1 and Figure 1 and Figure 2. In Table 1, the absolute errors between the exact solution and the approximate solution, at with the final time are given. But, in Figure 1, we presented comparison between the solution and the approximate solution using the proposed method at for different values of the final time and 0.5 at time step with. Also, in Figure 2, we presented the behavior approximate of solution at for different values of and 2 at time step with.

Example 2.

Consider the fractional cubically non-linear Klein-Gordon problem:

(24)

with the following initial and boundary conditions:

Table 1. The absolute error between the exact and approximate solutions at, and with.

Figure 1. The exact and approximate solutions at for different values of t.

Figure 2. The approximate solution at for different values of.

(25)

(26)

where the source term is given by

(27)

The exact solution of this problem is.

The obtained numerical results by means of the proposed method are shown in Figure 3 and Figure 4. In Figure 3, we presented comparison between the exact solution and the approximate solution using the proposed method at for different values of the final time and 1.25 at time step with. Also, in Figure 4, we presented the behavior of the approximate solution at for different values of at time step with.

Figure 3. The exact and approximate solutions at for different values of time t.

Figure 4. The approximate solution at for different values of.

4. Conclusion and Remarks

We have implemented Legendre spectral method for solving the non-linear FKGE. The proposed method gives excellent results when it is applied to FKGE. Absolute error by the method decreases while increasing iterations or level of resolution or both, as shown in Figures 1-4. It is evident that the overall errors can be made smaller by adding new terms from the series (8). Comparisons are made between approximate solutions and exact solutions to illustrate the validity and the great potential of the technique.

Acknowledgements

We thank the editor and the referee for their comments.

Cite this paper

Rubayyi T.Alqahtani, (2015) Approximate Solution of Non-Linear Fractional Klein-Gordon Equation Using Spectral Collocation Method. Applied Mathematics,06,2175-2181. doi: 10.4236/am.2015.613190

References

1. 1. Wazwaz, A.M. (2006) Compacton Solitons and Periodic Solutions for Some Forms of Nonlinear Klein-Gordon Equations. Chaos, Solitons and Fractals, 28, 1005-1013.
http://dx.doi.org/10.1016/j.chaos.2005.08.145

2. 2. El-Sayed, S.M. (2003) The Decomposition Method for Studying the Klein-Gordon Equation. Chaos, Solitons and Fractals, 18, 1026-1030.
http://dx.doi.org/10.1016/S0960-0779(02)00647-1

3. 3. Yusufoglu, E. (2008) The Variational Iteration Method for Studying the Klein-Gordon Equation. Applied Mathematics Letters, 21, 669-674.
http://dx.doi.org/10.1016/j.aml.2007.07.023

4. 4. El-Sayed, A.M.A., Elsaid, A. and Hammad, D. (2012) A Reliable Treatment of Homotopy Perturbation Method for Solving the Nonlinear Klein-Gordon Equation of Arbitrary (Fractional) Orders. Journal of Applied Mathematics, 2012, 1-13.
http://dx.doi.org/10.1155/2012/581481

5. 5. Doha, E.H., Bhrawy, A.H. and Ezz-Eldien, S.S. (2011) A Chebyshev Spectral Method Based on Operational Matrix for Initial and Boundary Value Problems of Fractional Order. Computers and Mathematics with Applications, 62, 2364- 2373.
http://dx.doi.org/10.1016/j.camwa.2011.07.024

6. 6. El-Sayed, A.M.A., Elsaid, A., El-Kalla, I.L. and Hammad, D. (2012) A Homotopy Perturbation Technique for Solving Partial Differential Equations of Fractional Order in Finite Domains. Applied Mathematics and Computation, 218, 8329-8340.
http://dx.doi.org/10.1016/j.amc.2012.01.057

7. 7. Oldham, K.B. and Spanier, J. (1974) The Fractional Calculus. Academic Press, New York.

8. 8. Podlubny, I. (1999) Fractional Differential Equations. Academic Press, New York.

9. 9. Khader, M.M. (2011) On the Numerical Solutions for the Fractional Diffusion Equation. Communications in Nonlinear Science and Numerical Simulations, 16, 2535-2542.
http://dx.doi.org/10.1016/j.cnsns.2010.09.007

10. 10. Khader, M.M. (2015) An Efficient Approximate Method for Solving Fractional Variational Problems. Applied Mathematical Modelling, 39, 1643-1649.
http://dx.doi.org/10.1016/j.apm.2014.09.012

11. 11. Khader, M.M. (2015) Fractional Chebyshev Finite Difference Method for Solving the Fractional-Order Delay BVPs. International Journal of Computational Methods, 12, Article ID: 1550033.
http://dx.doi.org/10.1142/s0219876215500334

12. 12. Khader, M.M. (2014) On the Numerical Solution and Convergence Study for System of Non-Linear Fractional Diffusion Equations. Canadian Journal of Physics, 92, 1658-1666.
http://dx.doi.org/10.1139/cjp-2013-0464

13. 13. Khader, M.M., Sweilam, N.H. and Mahdy, A.M.S. (2011) An Efficient Numerical Method for Solving the Fractional Diffusion Equation. Journal of Applied Mathematics and Bioinformatics, 1, 1-12.

14. 14. Khader, M.M. and Hendy, A.S. (2012) The Approximate and Exact Solutions of the Fractional-Order Delay Differential Equations Using Legendre Pseudo-Spectral Method. International Journal of Pure and Applied Mathematics, 74, 287-297.

15. 15. Sweilam, N.H., Khader, M.M. and Al-Bar, R.F. (2007) Numerical Studies for a Multi-Order Fractional Differential Equation. Physics Letters A, 371, 26-33.
http://dx.doi.org/10.1016/j.physleta.2007.06.016

16. 16. Sweilam, N.H. and Khader, M.M. (2010) A Chebyshev Pseudo-Spectral Method for Solving Fractional Integro-Diffe- rential Equations. ANZIAM, 51, 464-475.
http://dx.doi.org/10.1017/S1446181110000830

17. 17. Sweilam, N.H., Khader, M.M. and Nagy, A.M. (2011) Numerical Solution of Two-Sided Space-Fractional Wave Equation Using Finite Difference Method. Journal of Computational and Applied Mathematics, 235, 2832-2841.
http://dx.doi.org/10.1016/j.cam.2010.12.002

18. 18. Sweilam, N.H., Khader, M.M. and Mahdy, A.M.S. (2012) Crank-Nicolson Finite Difference Method for Solving Time-Fractional Diffusion Equation. Journal of Fractional Calculus and Applications, 2, 1-9.

19. 19. Bell, W.W. (1968) Special Functions for Scientists and Engineers. Great Britain, Butler and Tanner Ltd, Frome and London.

20. 20. Smith, G.D. (1965) Numerical Solution of Partial Differential Equations. Oxford University Press, Oxford.