Advances in Mathematical Physics
Volume 2021 (2021), Article ID 5528928, 8 pages
https://doi.org/10.1155/2021/5528928

Research Article

Solution of Space-Time-Fractional Problem by Shehu Variational Iteration Method

Suleyman Cetinkaya, Ali Demir, and Hulya Kodal Sevindir

Department of Mathematics, Kocaeli University, Kocaeli 41380, Turkey

Correspondence should be addressed to Suleyman Cetinkaya; suleyman.cetinkaya@kocaeli.edu.tr

Received 27 January 2021; Revised 26 March 2021; Accepted 19 April 2021; Published 30 April 2021

Academic Editor: Marianna Ruggieri

Copyright © 2021 Suleyman Cetinkaya 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.

In this study, we deal with the problem of constructing semianalytical solution of mathematical problems including space-time-fractional linear and nonlinear differential equations. The method, called Shehu Variational Iteration Method (SVIM), applied in this study is a combination of Shehu transform (ST) and variational iteration method (VIM). First, ST is utilized to reduce the time-fractional differential equation with fractional derivative in Liouville-Caputo sense into an integer-order differential equation. Later, VIM is implemented to construct the solution of reduced differential equation. The convergence analysis of this method and illustrated examples confirm that the proposed method is one of best procedures to tackle space-time-fractional differential equations.

1. Introduction

Last couple of decades, employing fractional differential equations in modelling of processes such as dynamical systems, biology, fluid flow, signal processing, electrical networks, reaction and diffusion procedure, and advection–diffusion–reaction process [1–4] has gained great importance since these models reflect the behaviour of the processes better than integer-order differential equations.

Consequently, a great deal of methods such as [3, 4] are established to construct analytical and numerical solutions of fractional differential equations. Moreover, their existence, uniqueness, and stability have been studied by many scientists.

One of the significant integral transformations is Shehu transformation proposed by Maitama and Zhao [5]. This linear transformation is a generalization of Laplace transformation. However, the Laplace transformation is obtained by substituting in Shehu transformation. By this transformation, differential equations are reduced into simpler equations.

Various methods such as the homotopy perturbation method (HPM) and VIM are utilized to establish approximate solutions of differential equations of any kind [6, 7]. As a result, it is employed widely to deal with differential equations in various branches of science [8–11]. VIM has been modified by many researchers to improve this method. By modified VIM, the approximate solutions of initial value problems can be established by making use of an initial condition.

2. Preliminaries

In this section, preliminaries, notations, and features of the fractional calculus are given [12, 13]. Riemann-Liouville time-fractional integral of a real valued function is defined as

where denotes the order of the integral.

The -order Liouville-Caputo time-fractional derivative operator of is defined as

The function

is called Mittag-Leffler function depending on two parameters and .

The following set of functions has Shehu transformation:

and it is defined as

which has the following property:

The inverse Shehu inverse transform of is defined as

where [5].

For the -order of Liouville-Caputo time-fractional derivative of , the Shehu transformation has the following form [14]:

3. Main Results

3.1. Fractional Shehu Variational Iteration Method

To reveal the fundamental notions of this method, let us take the following space-time-fractional initial value problem in the Liouville-Caputo fractional derivative:

where , , and denote the nonlinear, linear part of the differential equation, and the source function, respectively.

Utilizing Shehu transformation for Equation (9), we have

Employing the inverse Shehu transformation for Equation (11) leads to

where , and so

The following recurrence relation is established by VIM:

Alternately,

is called the -order of truncated solution.

If

exists, then the analytical solution .

3.2. Convergence Theorem

Now, the convergence of VIM is investigated and required conditions and error estimate [15] are established for Equation (9).

The operator is introduced as

where denote the components of the solution satisfying

Theorem 1 [16]. Let , defined in (5), be an operator from a Banach space BS to BS. The series solution as defined in (6) converges if exists such that , (i.e., ), .

Theorem 1, obtained from the Banach fixed-point theorem, is utilized to establish a sufficient condition for the convergence of fractional VIM.

Theorem 2 [16]. The exact solution of nonlinear problem (9) exists under the condition that the series solution defined in (18) converges.

Theorem 3 [16]. Suppose that the series solution defined in (18) converges to the solution . The maximum error for the approximate solution satisfies the following inequality:

The series solution of problem (9) is convergent to an exact solution , if the conditions

, hold where the parameters for are introduced as

Furthermore, the maximum absolute truncation error satisfies the inequality

where .

4. Illustrative Examples

Example 1. Let us consider following space-time-fractional initial value problem

Step 1. Implementing Shehu transform for (23), we have

Step 2. Taking the inverse Shehu transform of (25), we get

and so

Step 3. Employing the variational iteration method, we obtain

Based on the iteration formula (28), we have

By using the recurrence relation, we obtained the approximate solution of (23) as follows:

As a result, the analytical solution of (23) is reached by taking the limit of (30):

where is the two-parameter Mittag-Leffler function.

Notice from Table 1 and Figure 1 that the values of the solution for and exact solution are the same which implies that the method implemented in this study is one of the best one for the solution of space-time-fractional differential equations of any order. Moreover, it is clear from Figure 1 that as and tend to 1, the corresponding solutions tend to exact solution. Three-dimensional graphs of exact solution and a truncated solution are given in Figure 2.

Example 2. Let us consider the space-time-fractional equation

with the condition at .

Step 1. Carrying out Shehu transform of (32), we have

Step 2. Enforcing inverse Shehu transform of (34), we obtain

and so

Step 3. Utilizing the variational iteration method, we have

Table 1: Comparison of the exact solution with the truncated solutions by SVIM for various and for Example 1.

0,2
-5
-0,0728794142905183
-0,0556900475433218
0,00822974704902003
0,00822974704902003

0
0,464049675672513
0,865372126175867
1,22140275816017
1,22140275816017

5
272,995384456584
235,378704685654
181,272241875151
181,272241875151

0,4
-5
-0,102568289434788
-0,0717844261439092
0,0100518357446336
0,0100518357446336

0
0,653089516564494
1,11546396921542
1,49182469764127
1,49182469764127

5
384,205469814638
303,402959554171
221,406416204187
221,406416204187

0,6
-5
-0,140869106433447
-0,0919274211102744
0,0122773399030684
0,0122773399030684

0
0,896964716156101
1,42846758746576
1,82211880039051
1,82211880039051

5
527,674601164535
388,539034541079
270,426407426153
270,426407426153

0,8
-5
-0,190052579642384
-0,117094221086887
0,0149955768204777
0,0149955768204777

0
1,21013373669836
1,81953651567707
2,22554092849247
2,22554092849247

5
711,908534824601
494,908647082167
330,299559909649
330,299559909649

1
-5
-0,252961550186961
-0,148489952829492
0,0183156388887342
0,0183156388887342

0
1,61069797918433
2,30739731539736
2,71828182845905
2,71828182845905

5
947,556128411521
627,605367523717
403,428793492735
403,428793492735

Figure 1: 6th order of truncated solutions for various values of and and exact solution at for Example 1.

Figure 2: 6th order of truncated solutions for and exact solution for Example 1.

Based on the iteration formula (37), we have

By using the recurrence relation, the approximate solution of (32) is obtained as follows:

Hence, the analytical solution of (32) and (33) is reached by taking the limit of (39):

which is the same as obtained in [17].

As in Example 1, it is obvious from Table 2 and Figure 3 that the values of the solution for and exact solution are the same and as and tend to 1, the corresponding solutions tend to exact solution which indicates that the method employed in this research is a good choice for the solution of space-time-fractional differential equations of any order in Figure 4, 3-dimensional graphs of exact solution and a truncated solution are presented.

Table 2: Comparison of the exact solution with the truncated solutions by SVIM for various and for Example 2.

0,2
-5
2,97723189881641
2,95468850942983
3,00924465929121
3,00924465929121

0
3,14497276182554
3,70409889494893
4,37202909023507
4,37202909023507

5
88,2859013271367
194,512854240003
206,627171662420
206,627171662420

0,4
-5
3,01929874786381
2,98503015564532
3,00507357659074
3,00507357659074

0
2,87711786964444
3,23261761498140
3,75298552978051
3,75298552978051

5
-69,2902227453216
66,2713155938424
114,752961233252
114,752961233252

0,6
-5
3,03346340884022
2,99928895133517
3,00278443786961
3,00278443786961

0
2,78692633344582
3,01104904237009
3,41324722055397
3,41324722055397

5
-122,348923979388
6,00530743066654
64,3313254927734
64,3313254927734

0,8
-5
3,03543594562253
3,00548196320081
3,00152813190282
3,00152813190282

0
2,77436647600195
2,91481533308062
3,22679488322353
3,22679488322353

5
-129,737751709177
-20,1699819673740
36,6593450875042
36,6593450875042

1
-5
3,03265269180591
3,00773966021300
3,00083865656976
3,00083865656976

0
2,79208846297849
2,87973279769261
3,12446767091966
3,12446767091966

5
-119,312099237844
-29,7123296891665
21,4726402473266
21,4726402473266

Figure 3: 6th order of truncated solutions for various values of and and exact solution at for Example 2.

Figure 4: 6th order of truncated solutions for and exact solution for Example 2.

5. Conclusions

In this research, the targeted goal is to construct truncated solutions of linear/nonlinear space-time-fractional initial value problem by employing SVIM, the combination of the Shehu transform and variational iteration method. The main advantage of this method is that its implementation is straightforward and fruitful. Moreover, the illustrated examples reveal that the obtained approximate solutions with high precision converge swiftly to exact analytical solutions.

In the future study, this method and its improved modifications are applied to initial value problems including space-time-fractional linear and nonlinear differential equations.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Authors’ Contributions

All authors contributed equally to this work.

Acknowledgments

The first author would like to thank Scientific and Technological Research Council of Turkey (TUBITAK) for the financial support of the 2228-B Fellowship Program.

References

M. H. Heydari, Z. Avazzadeh, and A. Atangana, “Orthonormal shifted discrete Legendre polynomials for solving a coupled system of nonlinear variable-order time fractional reaction-advection-diffusion equations,” Applied Numerical Mathematics, vol. 161, pp. 425–436, 2021.

A. Jannelli, M. Ruggieri, and M. P. Speciale, “Analytical and numerical solutions of time and space fractional advection- diffusion-reaction equation,” Communications in Nonlinear Science and Numerical Simulation, vol. 70, pp. 89–101, 2019.

M. Şenol, O. S. Iyiola, H. D. Kasmaei, and L. Akinyemi, “Efficient analytical techniques for solving time-fractional nonlinear coupled Jaulent–Miodek system with energy-dependent Schrödinger potential,” Advances in Difference Equations, vol. 2019, no. 1, 2019.

L. Akinyemi and O. S. Iyiola, “Exact and approximate solutions of time-fractional models arising from physics via Shehu transform,” Mathematical Methods in the Applied Sciences, vol. 43, no. 12, pp. 7442–7464, 2020.

S. Maitama and W. Zhao, “New integral transform: Shehu transform a generalization of Sumudu and Laplace transform for solving differential equations,” International Journal of Analysis and Applications, vol. 17, no. 2, pp. 167–190, 2019.

J. H. He, “Approximate analytical solution for seepage flow with fractional derivatives in porous media,” Computer Methods in Applied Mechanics and Engineering, vol. 167, no. 1-2, pp. 57–68, 1998.

P. K. Gupta and M. Singh, “Homotopy perturbation method for fractional Fornberg-Whitham equation,” Computers & Mathematcs with Applications, vol. 61, no. 2, pp. 250–254, 2011.

Z. Odibat and S. Momani, “The variational iteration method: an efficient scheme for handling fractional partial differential equations in fluid mechanics,” Computers & Mathematcs with Applications, vol. 58, no. 11–12, pp. 2199–2208, 2009.

B. Ibis and M. Bayram, “Numerical comparison of methods for solving fractional differential-algebraic equations (FDAEs),” Computers & Mathematcs with Applications, vol. 62, no. 8, pp. 3270–3278, 2011.

S. Abbasbandy, “An approximation solution of a nonlinear equation with Riemann-Liouville's fractional derivatives by He's variational iteration method,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 53–58, 2007.

M. Inc, “The approximate and exact solutions of the space- and time-fractional Burgers equations with initial conditions by variational iteration method,” Journal of Mathematical Analysis and Applications, vol. 345, no. 1, pp. 476–484, 2008.

A. A. Kilbas, H. M. Srivastava, and J. J. Trujıllo, Theory and Applications of Fractional Differential Equations, Elsevier, New York, NY, USA, 2006.

I. Podlubny, Fractional Differential Equations, Elsevier, San Diego, CA, USA, 1998.

R. Belgacem, D. Baleanu, and A. Bokhari, “Shehu transform and applications to Caputo-fractional differential equations,” International Journal of Analysis and Applications, vol. 17, no. 6, pp. 917–927, 2019.

M. G. Sakar and H. Ergoren, “Alternative variational iteration method for solving the time-fractional Fornberg-Whitham equation,” Applied Mathematical Modelling, vol. 39, no. 14, pp. 3972–3979, 2015.

Z. M. Odibat, “A study on the convergence of variational iteration method,” Mathematical and Computer Modelling, vol. 51, no. 9-10, pp. 1181–1192, 2010.

S. Choudhary and V. Daftardar-Gejji, “Invariant subspace method: a tool for solving fractional partial differential equations,” Fractional Calculus and Applied Analysis, vol. 20, no. 2, 2017.