Complexity
Volume 2018 (2018), Article ID 3249720, 7 pages
https://doi.org/10.1155/2018/3249720

Research Article

New Iterative Method for the Solution of Fractional Damped Burger and Fractional Sharma-Tasso-Olver Equations

Mohammad Jibran Khan,¹ Rashid Nawaz,² Samreen Farid,² and Javed Iqbal²

¹Department of Computer and IT, Sarhad University of Science and Information Technology, Peshawar, Pakistan
²Department of Mathematics, Abdul Wali Khan University Mardan, Khyber Pakhtunkhwa, Pakistan

Correspondence should be addressed to Rashid Nawaz; rashid_uop@yahoo.com

Received 14 May 2018; Revised 19 July 2018; Accepted 30 July 2018; Published 5 September 2018

Academic Editor: Jesus M. Muñoz-Pacheco

Copyright © 2018 Mohammad Jibran Khan et al. 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.

The new iterative method has been used to obtain the approximate solutions of time fractional damped Burger and time fractional Sharma-Tasso-Olver equations. Results obtained by the proposed method for different fractional-order derivatives are compared with those obtained by the fractional reduced differential transform method (FRDTM). The 2nd-order approximate solutions by the new iterative method are in good agreement with the exact solution as compared to the 5th-order solution by the FRDTM.

1. Introduction

Most of the problems arising in the physical and biological area of science are nonlinear in nature, and it is not always possible to find the exact solution of such problems. These problems become more complicated when they involve fractional derivatives and are modelled through mathematical tools from fractional calculus. Fractional partial differential equations (FPDEs) are tremendous instrument and are widely used to describe many significant phenomena and dynamic processes such as engineering, rheology, acoustic, electrical networks, and viscoelasticity [1–6]. Generally, partial differential equations (PDEs) are hard to tackle, and their fractional-order types are more complicated [7, 8]. Therefore, several analytical and approximate methods can be used for finding their approximate solutions such as Adomian decomposition [9], homotopy analysis [10], tau method [11], residual power series method [12], and optimal homotopy asymptotic method [13]. Though the study of FPDEs has been obstructed due to the absence of proficient and accurate techniques, the derivation of approximate solution of FPDEs remains a hotspot and demands to attempt some dexterous and solid plans which are of interest. Daftardar-Gejji and Jafari proposed an iterative method called the new iterative method (NIM) for finding the approximate solution of differential equations [14]. NIM does not require the need for calculation of tedious Adomian polynomials in nonlinear terms like ADM, the need for determination of a Lagrange multiplier in its algorithm like VIM, and the need for discretization like numerical methods. The proposed method handles linear and nonlinear equations in an easy and straightforward way. Recently, the method has been extended for differential equations of the fractional order [15–17].

In the present study, we have implemented NIM for finding the approximate solution of the following fractional-order damped Burger equation.

Second, consider the fractional-order Sharma-Tasso-Olver equation of the following form.

where is the parameter describing the order of fractional derivatives, is the function of and , and are constants. The fractional derivatives are described in the Caputo sense. NIM converges rapidly to the exact solution compared to FRDTM, and only at the 2nd iteration does the proposed method yield very encouraging results. The accuracy of the proposed method can further be increased by taking higher-order approximations.

2. Definitions

In this section, we have stated some definitions which are relevant to our work.

Definition 1. A function , is said to be in space , if there exists a real number , such that , where . The function , is said to be in space if only if .

Definition 2. The R-L fractional integral operator of order of a function , is as follows:

Because of certain disadvantages of R-L fractional derivative operator, Caputo proposed modified fractional differential operator as follow.

Definition 3. Caputo fractional derivative of takes the following form.

where

Definition 4. If , and , then and .

The properties of the operator are shown as follows:

(i) exists for almost every .(ii)(iii)(iv)

In the equations above, , , , and

3. New Iterative Method

The basic mathematical theory of NIM is described as follows.

Let us consider the following nonlinear equation:

where , , is the known function and and are the linear and nonlinear functions of , respectively. According to the basic idea of NIM, the solution of the above equation has the series form.

The linear operator can be decomposed as

The decomposition of the nonlinear operator is as follows:

Hence, the general equation of (6) takes the following form:

From this, we have

The -term series solution of the general equation (6) takes the following form:

4. Applications

Example 1 (damped Burger equation).
Consider the damped Burger equation

together with IC

where is a constant. The exact solution of (13) is of the following form:

Using the operator on both sides of (13) using the initial condition and Definition 4 yields

where and .
According to (11), we have

The three-term approximate solution of the above equation is

Example 2 (Sharma-Tasso-Olver equation).
One can consider

together with IC

where is a constant. The exact solution of (19) for is of the following form:

Using the operator on both sides of (19) using the initial condition and Definition 4 yields

where and .
According to (11), we have

The three-term approximate solution of the above equation is as follows:

5. Results and Discussion

We have implemented NIM for finding the approximate solutions of the fractional damped Burger equation and fractional Sharma-Tasso-Olver equation. Tables 1 and 2 show the numerical results of the 2nd-order NIM which are compared with those of the 5th-order fractional reduced differential transform method (FRDTM) solution [18] for the fractional-order damped Burger equation. Tables 3 and 4 show the comparison of the proposed scheme with the FRDTM for the fractional-order Sharma-Tasso-Olver equation. Figure 1 shows the comparison of 2D plot of the approximate and exact solution by NIM for the classical damped Burger equation. Figure 2 shows the comparison of the approximate solution for different values of with the exact solution at . In Figures 3 and 4, 3D plots of approximate and exact solutions by NIM for the damped Burger equation are given. In Figure 5, the 2D plots of the approximate and exact solution for the classical Sharma-Tasso-Olver equation are given. Figure 6 shows the comparison of the approximate solution for different values of with the exact solution at . The 3D plots of approximate and exact solutions for the Sharma-Tasso-Olver equation are given in Figures 7 and 8. Throughout computations, we take and .

Table 1: Comparison of numerical results of NIM and FRDTM at and .

5th-order FRDTM [18]
2nd-order NIM
Exact solution
Absolute error

−5
0.002
−4.980060
−4.98006
−4.98006
1.19501 × 10⁻⁷

0.004
−4.960239
−4.96024
−4.96024
9.52046 × 10⁻⁷

0.006
−4.940535
−4.94054
−4.94054
3.19985 × 10⁻⁷

0.008
−4.920949
−4.92096
−4.92095
7.55346 × 10⁻⁶

−3
0.002
−2.988036
−2.98804
−2.98804
7.17009 × 10⁻⁸

0.004
−2.976143
−2.97614
−2.97614
5.71228 × 10⁻⁷

0.006
−2.964321
−2.96432
−2.96432
1.91991 × 10⁻⁶

0.008
−2.952569
−2.95257
−2.95257
4.53208 × 10⁻⁶

3
0.002
2.988036
2.98804
2.98804
7.17009 × 10⁻⁸

0.004
2.976143
2.97614
2.97614
5.71228 × 10⁻⁷

0.006
2.964321
2.96432
2.96432
1.91991 × 10⁻⁶

0.008
2.952569
2.95257
2.95257
4.53208 × 10⁻⁶

5
0.002
4.980060
4.98006
4.98006
1.19501 × 10⁻⁷

0.004
4.960239
4.96024
4.96024
9.52046 × 10⁻⁷

0.006
4.940535
4.94054
4.94054
3.19985 × 10⁻⁶

0.008
4.920949
4.92096
4.92095
7.55346 × 10⁻⁶

Table 2: Comparison of numerical results of the 2nd-order NIM and the 5th-order FRDTM for different values of .

FRDTM [18] ()
FRDTM [18] ()
NIM ()
NIM ()

−5
0.002
−4.211841
−4.876244
−4.55366
−4.89911

0.004
−3.938004
−4.794473
−4.4015
−4.8326

0.006
−3.744961
−4.724456
−4.29706
−4.7758

0.008
−3.592556
−4.661485
−4.21704
−4.72486

−3
0.002
−2.527105
−2.925746
−2.7322
−2.93946

0.004
−2.362803
−2.876684
−2.6409
−2.89956

0.006
−2.246977
−2.834674
−2.57824
−2.86548

0.008
−2.155534
−2.796891
−2.53022
−2.83492

3
0.002
2.527105
2.925746
2.7322
2.93946

0.004
2.362803
2.876684
2.6409
2.89956

0.006
2.246977
2.834674
2.57824
2.86548

0.008
2.155534
2.796891
2.53022
2.83492

5
0.002
4.211841
4.876244
4.55366
4.89911

0.004
3.938004
4.794473
4.4015
4.8326

0.006
3.744961
4.724456
4.29706
4.7758

0.008
3.592556
4.661485
4.21704
4.72486

Table 3: Comparison of numerical results of NIM and FRDTM at and .

5th-order FRDTM [18]
2nd-order NIM
Exact solution
Absolute error

−5
0.002
−0.493320
−0.49332
−0.49332
7.13918 × 10⁻¹²

0.004
−0.493334
−0.493334
−0.493334
5.70817 × 10⁻¹¹

0.006
−0.493347
−0.493347
−0.493347
1.92542 × 10⁻¹⁰

0.008
−0.493360
−0.49336
−0.49336
4.56141 × 10⁻¹⁰

−3
0.002
−0.452664
−0.452664
−0.452664
6.65773 × 10⁻¹²

0.004
−0.452755
−0.452755
−0.452755
5.333 × 10⁻¹¹

0.006
−0.452844
−0.452844
−0.452844
1.80218 × 10⁻¹⁰

0.008
−0.452934
−0.452934
−0.452934
4.27726 × 10⁻¹⁰

3
0.002
0.452484
0.452484
0.452484
6.64085 × 10⁻¹²

0.004
0.452393
0.452393
0.452393
5.30587 × 10⁻¹¹

0.006
0.452302
0.452302
0.452302
1.78845 × 10⁻¹⁰

0.008
0.452211
0.452211
0.452211
4.23386 × 10⁻¹⁰

5
0.002
0.493294
0.493294
0.493294
7.14723 × 10⁻¹²

0.004
0.493281
0.493281
0.493281
5.721 × 10⁻¹¹

0.006
0.493267
0.493267
0.493267
1.93192 × 10⁻¹⁰

0.008
0.493254
0.493254
0.493254
4.58194 × 10⁻¹⁰

Table 4: Comparison of numerical results of the 2nd-order NIM and the 5th-order FRDTM for different values of .

FRDTM [18] ()
FRDTM [18] ()
NIM ()
NIM ()

−5
0.002
−0.493876
−0.49339
−0.49363
−0.493375

0.004
−0.494098
−0.493447
−0.493755
−0.493421

0.006
−0.494262
−0.493496
−0.493849
−0.493461

0.008
−0.494397
−0.49354
−0.493926
−0.493497

−3
0.002
−0.456455
−0.453141
−0.454774
−0.453036

0.004
−0.457972
−0.453523
−0.455639
−0.453348

0.006
−0.459102
−0.453856
−0.456286
−0.45362

0.008
−0.460032
−0.45416
−0.456819
−0.453867

3
0.002
0.448366
0.452001
0.450211
0.452106

0.004
0.446521
0.451607
0.449182
0.451784

0.006
0.445064
0.451259
0.448372
0.4515

0.008
0.443806
0.450937
0.447674
0.451237

5
0.002
0.492686
0.493223
0.492959
0.493238

0.004
0.492412
0.493165
0.492806
0.493191

0.006
0.492194
0.493113
0.492687
0.493149

0.008
0.492007
0.493066
0.492583
0.49311

Figure 1: Numerical solution of the classical damped Burger equation with the exact solution at .

Figure 2: Numerical solution of the fractional damped Burger equation with the exact solution for different values of at .

Figure 3: 3D plot of for (DB) equation at .

Figure 4: 3D plot of the exact solution for the DB equation.

Figure 5: Numerical solution of the classical Sharma-Tasso-Olver equation with the exact solution at .

Figure 6: Numerical solution of the fractional Sharma-Tasso-Olver equation with the exact solution for different values of at .

Figure 7: 3D plot of for the STO equation at .

Figure 8: 3D plot of the exact solution of the STO equation.

By forming the numerical values and graphs, it is clear that NIM is a very powerful tool for the solution of fractional partial differential equations. The accuracy of the NIM can further be increased by taking higher-order approximations.

6. Conclusion

We have successfully applied NIM to time fractional (DB) and (STO) equations. Results reveal that NIM converges to the desired solution in lesser iteration compared to FRDTM. We can conclude that NIM computationally handles many physical and engineering problems in a simple and straightforward way. The accuracy of this method is also better than that of many methods which are computationally difficult to use.

Data Availability

All the data and the metadata regarding the finding of the manuscript have been given in the research manuscript.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

Sarhad University of Science and Information Technology, Peshawar, Pakistan, will pay the publication charges of the manuscript.

References

I. Podlubny, Fractional Differential Equations, Mathematics in Science and Engineering, Academic Press, San Diego, California, USA, 1999.

X.-J. Yang, J. A. Tenreiro Machado, and H. M. Srivastava, “A new numerical technique for solving the local fractional diffusion equation: two-dimensional extended differential transform approach,” Applied Mathematics and Computation, vol. 274, pp. 143–151, 2016.

R. S. Sengar, M. Sharma, and A. Trivedi, “Fractional calculus applied in solving instability phenomenon in fluid dynamics,” Journal Impact Factor, vol. 6, no. 5, pp. 34–44, 2015.

X. J. Yang, F. Gao, and H. M. Srivastava, “New rheological models within local fractional derivatives,” Romanian Reports in Physics, vol. 69, no. 113, 2017.

Y. Yu, P. Perdikaris, and G. E. Karniadakis, “Fractional modeling of viscoelasticity in 3D cerebral arteries and aneurysms,” Journal of Computational Physics, vol. 323, pp. 219–242, 2016.

V. Gill, K. Modi, and Y. Singh, “Analytic solutions of fractional differential equation associated with RLC electrical circuit,” Journal of Statistics and Management Systems, vol. 21, no. 4, pp. 575–582, 2018.

S. Das, “Analytical solution of a fractional diffusion equation by variational iteration method,” Computers & Mathematics with Applications, vol. 57, no. 3, pp. 483–487, 2009.

A. Elsaid, “The variational iteration method for solving Riesz fractional partial differential equations,” Computers & Mathematics with Applications, vol. 60, no. 7, pp. 1940–1947, 2010.

M. Ahmad, A. A. Bhat, and R. Jain, “Space time fractional telegraph equation and its application by using Adomian decomposition method,” Journal of New Theory, vol. 22, pp. 73–81, 2018.

F. Mabood, W. A. Khan, and A. I. M. Ismail, “MHD flow over exponential radiating stretching sheet using homotopy analysis method,” Journal of King Saud University - Engineering Sciences, vol. 29, no. 1, pp. 68–74, 2017.

S. Salahshour, T. Allahviranloo, and S. Abbasbandy, “Solving fuzzy fractional differential equations by fuzzy Laplace transforms,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 3, pp. 1372–1381, 2012.

S. Momani, O. Abu Arqub, M.’m. Abu Hammad, and Z. S. Abo-Hammour, “A residual power series technique for solving systems of initial value problems,” Applied Mathematics & Information Sciences, vol. 10, no. 2, pp. 765–775, 2016.

S. Sarwar and M. M. Rashidi, “Approximate solution of two-term fractional-order diffusion, wave-diffusion and telegraph models arising in mathematical physics using optimal homotopy asymptotic method,” Waves in Random and Complex Media, vol. 26, no. 3, pp. 365–382, 2016.

V. Daftardar-Gejji and H. Jafari, “An iterative method for solving nonlinear functional equations,” Journal of Mathematical Analysis and Applications, vol. 316, no. 2, pp. 753–763, 2006.

A. M. S. Mahdy and N. A. H. Mukhtar, “New iterative method for solving nonlinear partial differential equations,” Journal of Progressive Research in Mathematics, vol. 11, no. 3, pp. 1701–1711, 2017.

M. S. Al-Luhaibi, “New iterative method for fractional gas dynamics and coupled Burger’s equations,” The Scientific World Journal, vol. 2015, Article ID 153124, 8 pages, 2015.

B. R. Sontakke and A. Shaikh, “Approximate solutions of time fractional Kawahara and modified Kawahara equations by fractional complex transform,” Communications in Numerical Analysis, vol. 2016, no. 2, pp. 218–229, 2016.

M. S. Rawashdeh, “An efficient approach for time-fractional damped Burger & time-Sharma-Tasso-Olver equation using the FRDTM,” Applied Mathematics & Information Sciences, vol. 9, no. 3, pp. 1239–1246, 2015.