Abstract
In this paper, we present a computational method to solve the fractional Klein-Gordon equation (FKGE). The proposed technique is the grouping of orthogonal polynomial matrices and collocation method. The benefit of the computational method is that it reduces the FKGE into a system of algebraic equations which makes the problem straightforward and easy to solve. The main reason for using this technique is its high accuracy and low computational cost compared to other methods. The main solution behaviors of these equations are due to fractional orders, which are explained graphically. Numerical results obtained by the proposed computational method are also compared with the exact solution. The results obtained by the suggested technique reveals that the method is very useful for solving FKGE.
Introduction
The standard Klein-Gordon equation (KGE) is written as
where v indicates an unknown function in variables x and t, and h(x, t) stands for the source term. Due to the non-local nature and real-life applications of fractional derivatives, the fractional extension of this equation is very useful [1–12]. The fractional extension of this model handles the initial and boundary conditions of the model very accurately. The non-integer derivative helps in understanding the complete memory effect of the system. A broad literature of models with fractional derivatives can be found in [13–17]. Therefore, motivated by our ongoing research work into this special branch of mathematics (namely, fractional calculus), we study non-integer KGE by changing integer order derivative in both time and space using the Liouville-Caputo derivative of fractional order in the following manner:
having the initial conditions:
and boundary conditions:
The KGE is used in science, plasma (especially in quantum field theory), optical fibers, and dispersive wave-phenomena. Due to the great importance of KGE, many authors have studied it using various numerical and analytical schemes [18–27], each with their own limitations and shortcomings. The operational matrix method [28–38] is also applied to solve problems in fractional calculus. There are several other numerical and analytical methods which have been used to solve non-linear problems pertaining to fractional calculus, which can be found in [39, 40]. Some other applications of orthogonal polynomials-based solutions can be found in [41, 42].
In this paper, we present a computational technique which is a combination of the operational matrix and collocation method. We have used Chebyshev polynomials as a basis function for the construction of operational matrices of differentiations and integrations. In our proposed method, first the unknown function and their derivatives are approximated by taking finite dimensional approximations. Then, by using these approximations along with operational matrices of differentiations and integrations in the FKGE, we obtain a system of equations. Finally, by collocating this system, we get an approximate solution for the FKGE. The efficiency and accuracy of the used technique is shown by making a comparison amongst the results derived by our technique, exact solutions, and numerical results by some existing methods.
Some Basic Definitions
In this paper, we use non-integer order integrals and derivatives in the Riemann-Liouville and Caputo sense, respectively, which are given as:
Definition 2.1: The Riemann-Liouville non-integer integral operator of order α is presented as
I0f(x) = f(x).
Definition 2.2: The Liouville-Caputo non-integer derivative of order β are defined as [1–3]
, l − 1 < β < l, x > 0 and l is a natural number.
Chebyshev polynomial of the third kind of degree i on [0, 1] is given as,
The orthogonal property of these polynomials is given as:
where, , is a weight function and n and m are the degrees of polynomials.
A function can be approximated as
where, and
For any approximation taking n1 = n2 = n then Equation (6), can be written as,
The matrix C in Equation (8), is given as:
where, , is called the matrix of dual.
Theorem 1. If , is Chebyshev vector and we consider v > 0, then
where, I(v) = (e(i, j)), is (n + 1) × (n + 1) matrix of integral of non-integer order v and its entries are given by
Proof. Please see [30, 32, 38].
Theorem 2. If , is Chebyshev vector and we consider β > 0, then
where, D(β) = (s(i, j)), is (n + 1) × (n + 1) matrix of differentiation of non-integer order β and its entries are given by
Proof. Please see [38].
Method of Solution
In this section, we apply our proposed algorithm to solve a fractional model of KGE. We use equal number basis elements i.e. n1 = n2 = n, for any approximations of space and time variables. We initially approximate the time derivative of the unknown function as follows:
Taking integral of order β with respect to t on both sides of Equation (12), we have
where I(β)and I(1) are operational matrices of integration of order β and 1, respectively, and are given by Equation (10) and
where A and B are known square matrices and can be calculated using Equation (9).
Taking the differentiation of order γ on both sides of Equation (13), we get
where, D(γ)is the operational matrix of differentiation of order γ and is given by Equation (11). Further, the inhomogeneous term can be approximated as
where E is the known square matrix and can be calculated using Equation (9).
Grouping Equations (12), (13), (16), (17), and (2), we get
Equation (18), can be written as
Equation (19) is a system of equations which is easy to handle using the collocation method to determine the unknown matrix. By making use of the value of C in Equation (13), we can obtain an approximate solution for FLGE.
Numerical Experiments and Discussion
Example 1. Firstly, we take the time fractional KGE [26] given as
, having the ICs:
v(x, 0) = 1+sin(x), , for 0 ≤ x, t ≤ 1, and boundary conditions:
v(0, t) = cosh(t), v(1, t) = sin(1)+ cosh(t).
The exact solution is v(x, t) = sin(x)+ cosh(t).
In Figure 1, we have shown the three-dimensional trajectory of the approximate solution obtained by our used technique for integer KGE. In Figure 2, we have shown absolute errors by our proposed method for integer order KGE at n = 4.
Figure 1
Figure 2
From Figure 2 it is detected that absolute errors are very low, showing good agreement between the exact and approximate solution. In Figure 3, we have plotted fractional order KGE by changing the values of β and t at x = 0.8. In Figure 4, we have plotted fractional order KGE by changing the values of β and t at x = 1.
Figure 3
Figure 4
From Figures 3, 4, it can be seen that the solution changes consistently from fractional order to integer solution, showing the consistency of the proposed algorithm for time fractional order models.
Example 2. Secondly, taking the space fractional KGE [26] given as
, having the ICs:
v(x, 0) = xγ(1 − x), , for 0 ≤ x, t ≤ 1, and boundary conditions:
v(0, t) = 0, v(1, t) = 0, with source function h(x, t) = xγ(1 − x) exp (−t) − [(γ + 1) −(γ + 2) x] exp (−t) and the exact solution v(x, t) = xγ(1 − x) exp (−t).
In Figure 5, we have shown the three-dimensional trajectory of the approximate solution obtained by our proposed method for integer KGE. In Figures 6–8, we have shown absolute errors by our proposed method for integer order KGE at different values of n = 3, 5, and 7, respectively.
Figure 5
Figure 6
Figure 7
Figure 8
From Figures 6–8, it is detected that absolute errors are very low and show good agreement between the exact and approximate solution. It is also observed that absolute errors decrease when increasing the basis elements. In Figure 9, we have plotted fractional order KGE by changing the values of γ and x at t = 0.5. In Figure 10, we have plotted fractional order KGE by changing the values of γ and x at t = 1.
Figure 9
Figure 10
From Figures 9, 10, it can be seen that the solution changes consistently from fractional order to integer solution, showing the consistency of the proposed algorithm for space fractional order models. In Table 1, we have compared absolute errors by our method and the method used in [26] and observed that our used technique is more accurate in comparison to the technique used in [26].
Table 1
| x | n=3 | n=5 | ||
|---|---|---|---|---|
| Present method | Method in [26] | Present method | Method in [26] | |
| 0.0 | 1.6402e-05 | 1.1234e-03 | 2.2794e-08 | 8.1245e-04 |
| 0.1 | 1.0593e-05 | 2.7894e-03 | 1.4125e-08 | 6.8754e-04 |
| 0.2 | 3.4033e-06 | 4.4561e-03 | 2.6052e-09 | 4.9541e-04 |
| 0.3 | 4.5760e-06 | 1.7418e-03 | 1.0543e-08 | 2.4875e-04 |
| 0.4 | 1.2753e-05 | 7.8527e-03 | 2.4098e-08 | 8.5154e-04 |
| 0.5 | 2.0535e-05 | 5.9634e-03 | 3.6840e-08 | 4.0092e-04 |
| 0.6 | 2.7330e-05 | 6.8527e-03 | 4.7547e-08 | 6.1457e-04 |
| 0.7 | 3.2547e-05 | 3.1237e-03 | 5.4998e-08 | 6.9541e-04 |
| 0.8 | 3.5593e-05 | 1.7595e-03 | 5.7972e-08 | 7.1478e-04 |
| 0.9 | 3.5876e-05 | 3.0030e-03 | 5.5247e-08 | 2.0854e-04 |
| 1.0 | 3.2804e-05 | 0.0129e-03 | 4.5603e-08 | 0.0034e-04 |
Comparison of absolute errors by our method and the method in [26] at t = 1, Example 2.
In Table 2, we have compared our solution with the exact solution for different values of x and t at γ = 2.
Table 2
| (x, t) | Exact solution | Present method | Absolute errors |
|---|---|---|---|
| (0.1, 0.1) | 0.00814353 | 0.00814354 | 5.9041e-09 |
| (0.2, 0.2) | 0.02619938 | 0.02619938 | 4.3087e-09 |
| (0.3, 0.3) | 0.04667154 | 0.04667152 | 2.0611e-08 |
| (0.4, 0.4) | 0.06435072 | 0.06435071 | 1.2501e-08 |
| (0.5, 0.5) | 0.07581633 | 0.07581635 | 2.4796e-08 |
| (0.6, 0.6) | 0.07902887 | 0.07902889 | 2.3619e-08 |
| (0.7, 0.7) | 0.07299803 | 0.07299801 | 2.7123e-08 |
| (0.8, 0.8) | 0.05751410 | 0.05751405 | 4.9922e-08 |
| (0.9, 0.9) | 0.03293214 | 0.03293211 | 2.3894e-08 |
Comparison approximate and exact solution at γ = 2 and n = 5, Example 2.
Concluding Remarks
The key benefit of the used algorithm is that it works for both time and space FKGE. Using the proposed algorithm, we can derive an approximate solution for FKGE when the analytical solutions are not possible. It is also easy for computational purposes because FKGE is reduced into algebraic equations. We can apply this method together for time and space fractional, which reduces the time period of computation. Integer and fractional order behavior of KGE is shown. The outcomes of the present study are very helpful for scientists and engineers working in the mathematical modeling of natural phenomena. In a nutshell, we can say that with the aid of this scheme we can examine FKGE for use in quantum field theory, plasma, optical fibers, and dispersive wave-phenomena.
Statements
Data availability statement
All datasets generated for this study are included in the article/supplementary material.
Author contributions
All authors have worked equally on this manuscript and have read and approved the final manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1.
PodlubnyI. Fractional Differential Equations. San Diego, CA: Academic Press. (1999). p. 340.
2.
CaputoM. Elasticita e Dissipazione. Bologna: Zani-Chelli (1969).
3.
KilbasAASrivastavaHMTrujilloJJ. Theory and Applications of Fractional Differential Equations. Amsterdam: Elsevier (2006). p. 540.
4.
HristovJ. Transient heat diffusion with a non-singular fading memory: from the cattaneo constitutive equation with Jeffrey's kernel to the Caputo-Fabrizio time-fractional derivative. Thermal Sci. (2016) 20:765–70. 10.2298/TSCI160112019H
5.
BaskonusHMMekkaouiTHammouchZBulutH. Active control of a Chaotic fractional order economic system. Entropy. (2015) 17:5771–83. 10.3390/e17085771
6.
SinghH. A new stable algorithm for fractional Navier-Stokes equation in polar coordinate. Int J Appl Comput Math. (2017) 3:3705–22. 10.1007/s40819-017-0323-7
7.
YangXJMachadoJATCattaniCGaoF. On a fractal LC-electric circuit modeled by local fractional calculus. Commun Nonlin Sci Numeric Simul. (2017) 47:200–6. 10.1016/j.cnsns.2016.11.017
8.
YangXJ. A new integral transform operator for solving the heat-diffusion problem. Appl Math Lett. (2017) 64:193–7. 10.1016/j.aml.2016.09.011
9.
SrivastavaHMKumarDSinghJ. An efficient analytical technique for fractional model of vibration equation. Appl Math Model. (2017) 45:192–204. 10.1016/j.apm.2016.12.008
10.
KumarDSinghJBaleanuD. A new analysis for fractional model of regularized long-wave equation arising in ion acoustic plasma waves. Math Methods Appl Sci. (2017) 40:5642–53. 10.1002/mma.4414
11.
KumarDAgarwalRPSinghJ. A modified numerical scheme and convergence analysis for fractional model of Lienard's equation. J Comput Appl Math. (2018) 339:405–13. 10.1016/j.cam.2017.03.011
12.
SinghHSrivastavaHM. Numerical simulation for fractional-order Bloch equation arising in nuclear magnetic resonance by using the Jacobi polynomials. Appl Sci. (2020) 10:2850. 10.3390/app10082850
13.
SinghH. A new numerical algorithm for fractional model of Bloch equation in nuclear magnetic resonance. Alex Engl J. (2016) 55:2863–9. 10.1016/j.aej.2016.06.032
14.
SumelkaW. Non-local Kirchhoff–Love plates in terms of fractional calculus. Arch Civil Mech Eng. (2015) 15:231–42. 10.1016/j.acme.2014.03.006
15.
SinghH. Solution of fractional Lienard equation using Chebyshev operational matrix method. Nonlin Sci Lett A. (2017) 8:397–404.
16.
LazopoulosAK. On fractional peridynamic deformations. Arch Appl Mech. (2016) 86:1987–94. 10.1007/s00419-016-1163-3
17.
YingYPLianYPTangSQLiuWK. High-order central difference scheme for Caputo fractional derivative. Comput Methods Appl Mech Eng. (2017) 317:42–54. 10.1016/j.cma.2016.12.008
18.
AbbasbandyS. Numerical solutions of nonlinear Klein-Gordon equation by variational iteration method. Int J Num Meth Engg. (2007) 70:876–81. 10.1002/nme.1924
19.
Mohyud-DinSTNoorMANoorKI. Some relatively new techniques for nonlinear problems. Math Porb Engg. (2009) 2009:234849. 10.1155/2009/234849
20.
WazwazAM. The modified decomposition method for analytic treatment of differential equations. Appl Math Comput. (2006) 173:165–76. 10.1016/j.amc.2005.02.048
21.
KumarDSinghJKumarSSushila. Numerical computation of Klein-Gordon equations arising in quantum field theory by using homotopy analysis transform method. Alexandria Eng J. (2014) 53:469–74. 10.1016/j.aej.2014.02.001
22.
GolmankhanehAKBaleanuD. On nonlinear fractional Klein-Gordon equation. Signal Proc. (2011) 91:446–51. 10.1016/j.sigpro.2010.04.016
23.
KurulayM. Solving the fractional nonlinear Klein-Gordon equation by means of the homotopy analysis method. Adv Diff Eqs. (2012) 13:529–39. 10.1186/1687-1847-2012-187
24.
GepreelKAMohamedMS. Analytical approximate solution for nonlinear space-time fractional Klein-Gordon equation. Chin Phys B. (2013) 22:010201. 10.1088/1674-1056/22/1/010201
25.
KumarDSinghJBaleanuD. A hybrid computational approach for Klein- Gordon equations on Cantor sets. Nonlin Dyn. (2017) 87:511–7. 10.1007/s11071-016-3057-x
26.
KhaderMMKumarS. An accurate numerical method for solving the linear fractional Klein–Gordon equation. Math Meth Appl Sci. (2013) 37:2972–9. 10.1002/mma.3035
27.
YousifMAMahmoodBA. Approximate solutions for solving the Klein-Gordon sine-Gordon equations. J Assoc Arab Univ Bas Appl Sci. (2017) 22:83–90. 10.1016/j.jaubas.2015.10.003
28.
SinghHSrivastavaHMKumarD. A reliable numerical algorithm for the fractional vibration equation. Chaos Solitons Fractals. (2017) 103:131–8. 10.1016/j.chaos.2017.05.042
29.
SinghCSSinghHSinghVKOm SinghP. Fractional order operational matrix methods for fractional singular integro-differential equation. Appl Math Model. (2016) 40:10705–18. 10.1016/j.apm.2016.08.011
30.
SinghHSrivastavaHM. Jacobi collocation method for the approximate solution of some fractional-order Riccati differential equations with variable coefficients. Phys A. (2019) 523:1130–49. 10.1016/j.physa.2019.04.120
31.
SinghCSSinghHSinghSKumarD. An efficient computational method for solving system of nonlinear generalized Abel integral equations arising in astrophysics. Phys A. (2019) 525:1440–8. 10.1016/j.physa.2019.03.085
32.
SinghHPandeyRKBaleanuD. Stable numerical approach for fractional delay differential equations. Few Body Syst. (2017) 58:156. 10.1007/s00601-017-1319-x
33.
SinghHSinghCS. Stable numerical solutions of fractional partial differential equations using Legendre scaling functions operational matrix. Ain Shams Eng J. (2018) 9:717–25. 10.1016/j.asej.2016.03.013
34.
SinghH. Operational matrix approach for approximate solution of fractional model of Bloch equation. J King Saud Univ Sci. (2017) 29:235–40. 10.1016/j.jksus.2016.11.001
35.
WuJL. A wavelet operational method for solving fractional partial differential equations numerically. Appl Math Comput. (2009) 214:31–40. 10.1016/j.amc.2009.03.066
36.
SinghH. An efficient computational method for the approximate solution of nonlinear Lane-Emden type equations arising in astrophysics. Astrophys Space Sci. (2018) 363:71. 10.1007/s10509-018-3286-1
37.
SinghHSahooMRSinghOP. Numerical method based on Galerkin approximation for the fractional advection-dispersion equation. Int J Appl Comput Math. (2016) 3:2171–87. 10.1007/s40819-016-0233-0
38.
SinghH. Approximate solution of fractional vibration equation using Jacobi polynomials. Appl Math Comput. (2018) 317:85–100. 10.1016/j.amc.2017.08.057
39.
BhatterSMathurAKumarDSinghJ. A new analysis of fractional Drinfeld–Sokolov–Wilson model with exponential memory. Physi A. (2020) 537:122578. 10.1016/j.physa.2019.122578
40.
DubeyVPKumarRKumarDKhanISinghJ. An efficient computational scheme for nonlinear time fractional systems of partial differential equations arising in physical sciences. Adv Diff Equat. (2020) 2020:46. 10.1186/s13662-020-2505-6
41.
ShiriBPerfilievaIAlijaniZ. Classical approximation for fuzzy Fredholm integral equation, fuzzy sets and systems. (2020). 10.1016/j.fss.2020.03.023
42.
AlijaniZShiriBBaleanuD. A Chebyshev based numerical methods for solving fractional fuzzy differential equations (2019).
Summary
Keywords
fractional Klein-Gordon equation, fractional derivative, numerical solution, Chebyshev polynomials, operational matrices
Citation
Singh H, Kumar D and Pandey RK (2020) An Efficient Computational Method for the Time-Space Fractional Klein-Gordon Equation. Front. Phys. 8:281. doi: 10.3389/fphy.2020.00281
Received
25 February 2020
Accepted
22 June 2020
Published
23 September 2020
Volume
8 - 2020
Edited by
Jordan Yankov Hristov, University of Chemical Technology and Metallurgy, Bulgaria
Reviewed by
Babak Shiri, Neijiang Normal University, China; Ndolane Sene, Cheikh Anta Diop University, Senegal
Updates
Copyright
© 2020 Singh, Kumar and Pandey.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Devendra Kumar devendra.maths@gmail.com
This article was submitted to Mathematical and Statistical Physics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.