Journal Details Format
Journal
First Published
01 Jan 2016
Publication timeframe
2 times per year
Languages
English
Copyright
© 2020 Sciendo

# A Crank-Nicolson Approximation for the time Fractional Burgers Equation

###### Accepted: 29 Jul 2019
Journal Details Format
Journal
First Published
01 Jan 2016
Publication timeframe
2 times per year
Languages
English
Copyright
© 2020 Sciendo

In the present manuscript, Crank Nicolson finite difference method is going to be applied to get the approximate solutions for the fractional Burgers equation. The fractional derivative used in this equation is going to be taken into consideration in the Caputo sense. The L1 type discretization formula is going to be applied to this equation. For checking the efficiency of proposed methods, the error norms L2 and L have at the same time been calculated. Those newly got solutions using the presented method illustrate the easy usage and efficiency of the approach presented in this manuscript.

#### MSC 2010

Introduction

Fractional order integral and derivative are the generalizations of classical integral and derivative concepts which are examined in detail by Leibniz and Newton. The concepts of fractional integral and derivative are as old as integer order integral and derivative concepts, and the fractional derivative expression is first mentioned by Leibniz's letter to L’Hospital in 1695 . In the letter, Leibniz's question was ’Can the integer order derivatives be generalized to fractional order derivative’. This is known as the first emergence of the concept of fractional differential. In addition to Leibniz, many scientists such as Liouville, Riemann, Weyl, Lagrange, Laplace, Fourier, Euler and Abel have worked on the same subject . Many definitions are given in the literature for fractional derivative. Some of these are Riemann-Liouville, Caputo, Grünwald-Letnikov, Wely, Riesy fractional derivatives . Some studies have shown that these definitions are equivalent under certain conditions. There is more than one derivative definition in the fractional analysis, making it possible to use the most suitable one according to the problem and thus to get the best solution for this problem.

However, if the derivative is described as how the derivative of the fractional order is defined, the expression that when the order is selected being equal to the integer is the same as the integer order of the derivative.

The description of Caputo fractional derivative was first introduced by the Italian mathematician M.Caputo in the 1960s to eliminate the problem of the calculation of Riemann-Liouville definition of initial values in Laplace transform applications. The fundamental advantage of the Caputo approach is the fact that the appropriate initial conditions defined for Caputo fractional differential equations are identical. Therefore, in recent studies in the literature, in the exact and approximate solutions of fractional differential equations, instead of Riemann-Liouville fractional derivative operator, Caputo fractional derivative operator has been more preferred. Recently, studies on the solution of fractional differential equations have increased. There are several studies about fractional problems and their computational accuracy in the literature [YoussefI. K.El DewaikM. H.2017Solving Poisson's Equations with fractional order using Haar wavelet21271284' href="#j_amns.2020.2.00023_ref_004_w2aab3b7e1224b1b6b1ab2b1b4Aa">4, 5, 6]. Since the exact solution for many of the fractional differential equations are not found, various methods have been developed to find approximate or numerical solutions.

Definition 1

f (t) function continuously and can be differentiated, and for n − 1 ≤ γ < n Caputo means for the fractional order derivatives; $aCDtγ f(t)=1Γ(n−γ) ∫atf(n)(τ)(t−τ)γ−n+1 dτ$\,_a^CD_t^\gamma \;f(t) = {1 \over {\Gamma (n - \gamma)}}\;\;\int_a^t {{{f^{(n)}}(\tau)} \over {{{(t - \tau)}^{\gamma - n + 1}}}}\;\;d\tau is defined as above.

Definition 2

Finite difference methods are widely used in the solution of many linear and nonlinear partial differential equations. In general, the following way is followed in applying a finite difference method to a partial differential equation:

The given solution area of the problem is divided into meshes with geometric shapes and approximate solution for the problem is calculated on the nodes of each mesh. Proper finite difference approaches are obtained by using Taylor series instead of derivatives in differential equations. Thus, the present problem of solution of the differential equation has been converted into the problem of the solution of an algebraic system of equation consisting of difference equations. Thus the algebraic equation system obtained now may be solved easily by one of the direct or iterative methods.

Numerical Solution of the Model Problem

In this manuscript, we will deal with the nonlinear time fractional Burgers equation with the given initial and boundary conditions as a test problem presented as $Dtγ u+uux−νuxx=f(x,t)u(x,0)=g(x), 0≤x≤1.u(0,t)=h1(t),u(1,t)=h2(t), t≥0$\matrix{{D_t^\gamma \;u + u{u_x} - \nu {u_{xx}} = f(x,t)} \hfill \cr {u\left({x,0} \right) = g\left(x \right),\;\;\;\;\;\;\;\;\;0 \le x \le 1.} \hfill \cr {u(0,t) = {h_1}(t),u\left({1,t} \right) = {h_2}\left(t \right),\;t \ge 0} \hfill \cr} in which ν is the viscosity parameter, u(x, t) represents the speed of fluid media at time-space position (x, t) and $Dtγu(x,t)=1Γ(1−γ)∫at(t−τ)−γ∂u(x,τ)∂τdτ, 0<γ≤1$D_t^\gamma u\left({x,t} \right) = {1 \over {\Gamma (1 - \gamma)}}\int_a^t {(t - \tau)^{- \gamma}}{{\partial u(x,\tau)} \over {\partial \tau}}d\tau,\quad 0 < \gamma \le 1 is the fractional derivative given in the Caputo's form [1, 7]. The Burgers’ equation is the one of the most important differential equations arising in applied sciences and mathematical physics and has interesting applications in physics and astrophysics. The Burgers’ equation examines the modeling of fluid mechanics, diffusive waves in fluid dynamics and also it has many application in the theory of shock waves, mathematical modeling of turbulent fluid, sound waves in a viscous medium and so on. In Mathematical modelling, fractional derivatives provide more accurate and applicable models for real life problems. For the same reasons, recent researches have concentrated on investigating and proposing models on Fractional Burgers Equations. During the years, this attractive modelling ability of the equation have taken attention of most of scientific people, many important and useful papers have brighten science world. In this part of our manuscript, we want to remind some of them; Esen and Tasbozan  have investigated numerical solutions of the equation via quadratic B-spline Galerkin method which is an useful and efficient type of finite element approach, then they have obtained more accurate numerical solutions with using collocation points and raising the degree of basis in . Yaseen and Abbas, in , have used collocation method with the help of cubic trigonometric B-splines basis. Qui et al  have proposed an implicit difference scheme with L1 algorithm which is a specific time discretization of the Caputo fractional derivative. A numerical method focused on a finite difference scheme in terms of time and the Chebyshev spectral collocation method in terms of space is used to get approximate solutions of the Burgers’ equation in . Saad and Al-Sharif  have applied variational iteration method for solving the equation considering several initial conditions. Asgari and Hosseini  have focused on generalized time fractional Burger type equation, they put forward two semi implicit Fourier pseudospectral approximations for seeking solutions of the equation. As a different view to the mentioned equation, Khan et al  are used the generalized version of the differential transform method and homotopy perturbation method. Lombard and Matignon  present an article for better understanding the competition between nonlinear effects and nonlocal relaxation.

Throughout this manuscript, in order to contribute to literature, we will consider numerical solutions of the time fractional Burgers equations using finite difference approach. For the presented numerical solutions, to get a Crank Nicolson finite difference scheme to solve the time fractional Burgers equation as utilized in explicit difference method in Ref. , we will also discretize the derivative respect to time using the widely-known L1 formula  $∂γf∂tγ|tn≅(Δt)−γΓ(2−γ)∑k=0m−1bkγ[f(tn−k)−f(tn−k−1)], 0<γ≤1${\left. {{{{\partial ^\gamma}f} \over {\partial {t^\gamma}}}} \right|_{{t_n}}} \cong {{{{\left({\Delta t} \right)}^{- \gamma}}} \over {\Gamma \left({2 - \gamma} \right)}}\sum\limits_{k = 0}^{m - 1} b_k^\gamma \left[ {f\left({{t_{n - k}}} \right) - f\left({{t_{n - k - 1}}} \right)} \right],\quad 0 < \gamma \le 1 where $bkγ=(k+1)1−γ−k1−γ.$b_k^\gamma = {\left({k + 1} \right)^{1 - \gamma}} - {k^{1 - \gamma}}.

Crank Nicolson Finite Difference Scheme

Let's assume the fact that the solution domain for the present problem 0 ≤ x ≤ 1 is discretized into regular grids with equal length Δx in the x-direction and also with equal time intervals Δt over time t such that xj = jΔx, j = 1(1)M − 1 and tn = nΔt, n = 0(1)N and the numerical solution of u at the grid point (jΔx, nΔt) will denote by $Ujn$U_j^n throughout the study. Using L1 formula in Eq. (2) instead of Caputo derivative in Eq.(1) and utilizing the following discretization in place of the terms uux and uxx respectively: $∂γu(x,t)∂tγ|tn≅(Δt)−γΓ(2−γ)∑k=0m−1bkγ[Ujn−k−Ujn−k−1], 0<γ≤1uux≅Ujn+12(Uj+1n−Uj−1n2Δx)+Ujn2(Uj+1n+1−Uj−1n+12Δx)$\matrix{{{{\left. {{{{\partial ^\gamma}u(x,t)} \over {\partial {t^\gamma}}}} \right|}_{{t_n}}} \cong {{{{(\Delta t)}^{- \gamma}}} \over {\Gamma (2 - \gamma)}}\sum\limits_{k = 0}^{m - 1} b_k^\gamma \left[ {U_j^{n - k} - U_j^{n - k - 1}} \right],\quad 0 < \gamma \le 1} \cr {u{u_x} \cong {{U_j^{n + 1}} \over 2}\left({{{U_{j + 1}^n - U_{j - 1}^n} \over {2\Delta x}}} \right) + {{U_j^n} \over 2}\left({{{U_{j + 1}^{n + 1} - U_{j - 1}^{n + 1}} \over {2\Delta x}}} \right)} \cr} and $uxx≅Uj+1n+1−2Ujn+1+Uj−1n+12(Δx)2+Uj+1n−2Ujn+Uj−1n2(Δx)2.${u_{xx}} \cong {{U_{j + 1}^{n + 1} - 2U_j^{n + 1} + U_{j - 1}^{n + 1}} \over {2(\Delta x{)^2}}} + {{U_{j + 1}^n - 2U_j^n + U_{j - 1}^n} \over {2(\Delta x{)^2}}}.

We are able to get the following system of algebraic equations $Uj−1n+1(−Sν2(Δx)2−S4ΔxUjn)+Ujn+1(1+Sν(Δx)2+S4Δx(Uj+1n−Uj−1n))+Uj+1n+1(−Sν2(Δx)2+S4ΔxUjn)=Sν2(Δx)2Uj−1n+Ujn(1−Sν(Δx)2)+Sν2(Δx)2Uj+1n+Sf(x,t)−∑k=1m−1bkγ[Ujn+1−k−Ujn−k]$\matrix{{U_{j - 1}^{n + 1}\left({- {{S\nu} \over {2{{\left({\Delta x} \right)}^2}}} - {S \over {4\Delta x}}U_j^n} \right) + U_j^{n + 1}\left({1 + {{S\nu} \over {{{\left({\Delta x} \right)}^2}}} + {S \over {4\Delta x}}\left({U_{j + 1}^n - U_{j - 1}^n} \right)} \right) + U_{j + 1}^{n + 1}\left({- {{S\nu} \over {2{{\left({\Delta x} \right)}^2}}} + {S \over {4\Delta x}}U_j^n} \right)} \hfill \cr {= {{S\nu} \over {2{{\left({\Delta x} \right)}^2}}}U_{j - 1}^n + U_j^n\left({1 - {{S\nu} \over {{{\left({\Delta x} \right)}^2}}}} \right) + {{S\nu} \over {2{{\left({\Delta x} \right)}^2}}}U_{j + 1}^n + Sf(x,t) - \sum\limits_{k = 1}^{m - 1} b_k^\gamma \left[ {U_j^{n + 1 - k} - U_j^{n - k}} \right]} \hfill \cr} where S = ΔtγΓ(2 − γ).

Numerical Results

Numerical results obtained for the Eq. (1) are got using the Crank Nicolson finite difference methods. The efficiency of the present methods are tested using the error norm L2$L2=‖uexact−(UN)j‖2=Δx∑j=0M(ujexact−(UN)j)2${L_2} = {\left\| {{u^{exact}} - {{\left({{U_N}} \right)}_j}} \right\|_2} = \sqrt {\Delta x\sum\limits_{j = 0}^M {{\left({u_j^{exact} - {{\left({{U_N}} \right)}_j}} \right)}^2}} and also the maximum error norm L$L∞=‖uexact−(UN)j‖∞=maxj|ujexact−(UN)j|.${L_\infty} = {\left\| {{u^{exact}} - {{\left({{U_N}} \right)}_j}} \right\|_\infty} = \mathop {max}\limits_j \left| {u_j^{exact} - {{\left({{U_N}} \right)}_j}} \right|.

Example 1.Firstly, we are going to take time fractional Burgers equation (1) given with the following appropriate boundary conditions$U(0,t)=0,U(1,t)=0,t≥0$U(0,t) = 0,U(1,t) = 0,t \ge 0and the initial condition as$U(x,0)=0,0≤x≤1.$U(x,0) = 0,0 \le x \le 1.

The term f (x, t) is of the form $f(x,t)=2t2−γsin(2πx)Γ(3−γ)+2πt4sin(2πx)cos(2πx)+4νπ2t2sin(2πx).$f(x,t) = {{2{t^{2 - \gamma}}\sin \left({2\pi x} \right)} \over {\Gamma \left({3 - \gamma} \right)}} + 2\pi {t^4}\sin \left({2\pi x} \right)\cos \left({2\pi x} \right) + 4\nu {\pi ^2}{t^2}\sin \left({2\pi x} \right).

The analytic solution for the present problem is described as$U(x,t)=t2sin(2πx).$U\left({x,t} \right) = {t^2}\sin \left({2\pi x} \right).

A comparison for the exact and approximate solutions got by the Crank-Nicolson methods for fractional Burgers equation for the values of γ = 0.5, Δt = 0.00025 and tf = 1 for various values of M is presented in Table 1. As one could clearly see from this table, the decreasing values of the error norms L2and Lvalidate this comment. In Table 2, the approximate solutions obtained for values of M have been compared with those in Ref. . Those error norms for each values of M got by this method are smaller than those presented in Ref. .

The error norms L2 and L of the time fractional Burgers equation problem using the Crank Nicolson finite difference method with ν = 1.0, Δt = 0.00025 and tf =1 for various values of M

N = 10N = 20N = 40M = 80
L2 × 10322.640873265.444834241.220073330.16846258
L × 10330.494450827.700511771.725529150.23826679

A comparison of the errors for Example 1 at tf =1

N = 40 N = 80 N = 100 Present  Present L2 × 103 1.22007333 1.224329 0.16846258 0.177703 0.04239382 0.052299 L∞ × 103 1.72552915 1.730469 0.23826679 0.253053 0.05996900 0.076541

The solutions obtained at different times in Fig. 1 were given graphically with exact solutions. The results obtained at different times were found to be too close to each other on the graphics. A comparison of the analytic and approximate solutions using the Crank-Nicolson method for ν=1.0, N=40, γ=0.5 and Δt = 0.0025 at various values of final time

In Table 3, solutions for different viscosity values were given for values of γ = 0.5, Δt = 0.00025, tf =1.0, N = 40. It was observed that the errors increased gradually at small values of viscosity.

The error norms L2 and L of Example 1 for γ=0.5, Δt=0.00025, tf =1.0, N=40 and various values of ν

ν = 1ν = 0.5ν = 0.1ν = 0.01ν = 0.001
L2 × 1030.417611570.512591611.039281122.138803146.46231244
L × 1030.590407800.723427311.512801854.6570727522.95302718

Table 4 presents the error norms obtained at various selected values of γ.

The error norms L2 and L of Example 1 for N=120, tf = 1.0, Δt = 0.00025 for different values of γ

γ = 0.1γ = 0.25γ = 0.75γ = 0.9
L2 × 1030.024119760.024905180.026695470.02579288
L × 1030.034099050.035211150.037745920.03646791

Example 2.Secondly, we will consider the model problem with the appropriate boundary conditions$U(0,t)=t2,U(1,t)=−t2,t≥0$U(0,t) = {t^2},U(1,t) = - {t^2},t \ge 0and the initial condition as$U(x,0)=0,0≤x≤1.$U(x,0) = 0,0 \le x \le 1.

The term$f(x,t)=2t2−γcos(πx)Γ(3−γ)−πt4cos(πx)sin(πx)+νπ2t2cos(πx)$f(x,t) = {{2{t^{2 - \gamma}}\cos \left({\pi x} \right)} \over {\Gamma \left({3 - \gamma} \right)}} - \pi {t^4}\cos \left({\pi x} \right)\sin \left({\pi x} \right) + \nu {\pi ^2}{t^2}\cos \left({\pi x} \right)and the analytic solution for the problem$U(x,t)=t2cos(πx).$U(x,t) = {t^2}\cos (\pi x).

In Table 5, for various values of N, the error norms obtained for Example 2 have been illustrated. In the table, one can obviously see that the improvement in numerical solutions has been observed as expected for the increasing values of N.

Comparison of results at tf = 1.0 for γ = 0.5, Δt = 0.0001, ν = 1.0 and various mesh sizes

N = 10N = 20N = 40N = 80
L2 × 1030.816670900.235948590.092154640.05901177
L × 1031.136780940.320862370.122458710.09790576

In Table 6, the error norms obtained for different values of Δt have been given for values of ν = 1.0, Δx = 0.025 and γ = 0.5. As the Δt values decreased, it was observed that the errors also decreased accordingly. In Figure 2, the solutions of Problem 2 obtained at different times for values of ν = 1.0, Δx = 0.025, γ = 0.5 and Δt = 0.0025 have been presented together with their analytical values.

Comparison of results at tf = 1.0 for ν=1, N = 40 and various time steps

Δt = 0.0001Δt = 0.00025Δt = 0.0005Δt = 0.001Δt = 0.0025Δt = 0.005
L2 × 1030.092154640.161425900.279069060.515832671.227758342.41506007
L × 1030.122458710.245778980.474783710.932808102.306967154.59744900 A comparison of the analytic and approximate solutions using the Crank-Nicolson method for ν = 1.0, Δx = 0.025, γ = 0.5 and Δt = 0.0025 at different values of final time

Example 3.Lastly, we consider the model problem with the boundary conditions$U(0,t)=t2,U(1,t)=et2,t≥0$U(0,t) = {t^2},U(1,t) = e{t^2},t \ge 0and the initial conditions as$U(x,0)=0,0≤x≤1.$U(x,0) = 0,0 \le x \le 1.

The function f (x, t) is of the following form$f(x,t)=2t2−γexΓ(3−γ)+t4e2x−νt2ex.$f\left({x,t} \right) = {{2{t^{2 - \gamma}}{e^x}} \over {\Gamma \left({3 - \gamma} \right)}} + {t^4}{e^{2x}} - \nu {t^2}{e^x}.

The analytical solution of the problem is given by$U(x,t)=t2ex.$U(x,t) = {t^2}{e^x}.

In Table 7, approximate solutions obtained for values of tf = 1.0 and, Δx = 0.025 have been given at different Δt values. It is obvious that the errors get smaller and smaller as Δt values get smaller and smaller. Figure 3 shows the numerical and exact solutions obtained at different times for Problem 3. It is seen that the results are indiscriminately close to each other on the graphics.

Comparison of results at tf = 1.0 for ν = 1, Δx = 0.025 and various time steps

Δt = 0.0001Δt = 0.00025Δt = 0.0005Δt = 0.001
L2 × 1030.083444540.222522580.454979940.92040447
L × 1030.236483230.597225841.198489412.40107923 A comparison of the analytic and approximate solutions using the Crank-Nicolson method for ν = 1.0, N = 40, γ = 0.5 and Δt = 0.0025 at different values of final time
Conclusions

As a conclusion, in the present study the numerical solutions of the time fractional Burgers equation have been obtained by using finite difference method based on Crank-Nicolson discretization. The obtained results are compared with analytic and some of the numerical results available in the literature. This comparison has shown that the presented method is efficient and effective and can also be used for a wide range of physical and scientific applications. Moreover to illustrate the accuracy of the present method the error norms L2 and L are computed and given in tables. Two test problems have been used to show the accuracy of the present scheme for various values of parameters in the problem. Tables and figures show the results of these various tests and also comparisons with some available results together with the error norms L2 and L. Finally the present method has been shown to be applicable for more widely used fractional differential equations. A comparison of the analytic and approximate solutions using the Crank-Nicolson method for ν=1.0, N=40, γ=0.5 and Δt = 0.0025 at various values of final time A comparison of the analytic and approximate solutions using the Crank-Nicolson method for ν = 1.0, Δx = 0.025, γ = 0.5 and Δt = 0.0025 at different values of final time A comparison of the analytic and approximate solutions using the Crank-Nicolson method for ν = 1.0, N = 40, γ = 0.5 and Δt = 0.0025 at different values of final time

#### The error norms L2 and L∞ of the time fractional Burgers equation problem using the Crank Nicolson finite difference method with ν = 1.0, Δt = 0.00025 and tf =1 for various values of M

N = 10N = 20N = 40M = 80
L2 × 10322.640873265.444834241.220073330.16846258
L × 10330.494450827.700511771.725529150.23826679

#### Comparison of results at tf = 1.0 for γ = 0.5, Δt = 0.0001, ν = 1.0 and various mesh sizes

N = 10N = 20N = 40N = 80
L2 × 1030.816670900.235948590.092154640.05901177
L × 1031.136780940.320862370.122458710.09790576

#### A comparison of the errors for Example 1 at tf =1

N = 40 N = 80 N = 100 Present  Present L2 × 103 1.22007333 1.224329 0.16846258 0.177703 0.04239382 0.052299 L∞ × 103 1.72552915 1.730469 0.23826679 0.253053 0.05996900 0.076541

#### The error norms L2 and L∞ of Example 1 for N=120, tf = 1.0, Δt = 0.00025 for different values of γ

γ = 0.1γ = 0.25γ = 0.75γ = 0.9
L2 × 1030.024119760.024905180.026695470.02579288
L × 1030.034099050.035211150.037745920.03646791

#### Comparison of results at tf = 1.0 for ν = 1, Δx = 0.025 and various time steps

Δt = 0.0001Δt = 0.00025Δt = 0.0005Δt = 0.001
L2 × 1030.083444540.222522580.454979940.92040447
L × 1030.236483230.597225841.198489412.40107923

#### The error norms L2 and L∞ of Example 1 for γ=0.5, Δt=0.00025, tf =1.0, N=40 and various values of ν

ν = 1ν = 0.5ν = 0.1ν = 0.01ν = 0.001
L2 × 1030.417611570.512591611.039281122.138803146.46231244
L × 1030.590407800.723427311.512801854.6570727522.95302718

#### Comparison of results at tf = 1.0 for ν=1, N = 40 and various time steps

Δt = 0.0001Δt = 0.00025Δt = 0.0005Δt = 0.001Δt = 0.0025Δt = 0.005
L2 × 1030.092154640.161425900.279069060.515832671.227758342.41506007
L × 1030.122458710.245778980.474783710.932808102.306967154.59744900

Podlubny, I., (1998), Fractional Differential Equations, Acasemic Pres, San Diego.PodlubnyI.1998Fractional Differential EquationsAcasemic PresSan DiegoSearch in Google Scholar

Miller, K.S, Ross, B. (1993), An Introduction to The Fractional Calculus and Differential Equations, John Wiley and Sons Ltd.MillerK.SRossB.1993An Introduction to The Fractional Calculus and Differential EquationsJohn Wiley and Sons Ltd.Search in Google Scholar

Smith, G. D. (1986), Numerical Solution of Partial Differential Equations: Finite Difference Methods, Oxford University Press.SmithG. D.1986Numerical Solution of Partial Differential Equations: Finite Difference MethodsOxford University PressSearch in Google Scholar

Youssef, I. K. and El Dewaik, M. H., (2017) Solving Poisson's Equations with fractional order using Haar wavelet, 2,1, 271–284.YoussefI. K.El DewaikM. H.2017Solving Poisson's Equations with fractional order using Haar wavelet21271284Search in Google Scholar

Brzeziński, D. W., (2017), Comparison of Fractional Order Derivatives Computational Accuracy - Right Hand vs Left Hand Definition, 2,1, 237–248.BrzezińskiD. W.2017Comparison of Fractional Order Derivatives Computational Accuracy - Right Hand vs Left Hand Definition21237248Search in Google Scholar

Brzeziński, D. W., (2018), Review of numerical methods for NumILPT with computational accuracy assessment for fractional calculus, 3, 2, 487–502BrzezińskiD. W.2018Review of numerical methods for NumILPT with computational accuracy assessment for fractional calculus32487502Search in Google Scholar

Kilbas, A.A., Srivastava, H.M., Trujillo, J.J., (2006), Theory and Applications of Fractional Differential Equations, Elsevier, AmsterdamKilbasA.A.SrivastavaH.M.TrujilloJ.J.2006Theory and Applications of Fractional Differential EquationsElsevierAmsterdamSearch in Google Scholar

Esen, A., and Tasbozan, O., (2015), Numerical solution of time fractional Burgers equation, Acta Universitatis Sapientiae, Mathematica 7,2, 167–185.EsenA.TasbozanO.2015Numerical solution of time fractional Burgers equationActa Universitatis Sapientiae, Mathematica72167185Search in Google Scholar

Esen, A., and Tasbozan O., (2016), Numerical solution of time fractional burgers equation by cubic B-spline finite elements, Mediterranean Journal of Mathematics 13,3, 1325–1337.EsenA.TasbozanO.2016Numerical solution of time fractional burgers equation by cubic B-spline finite elementsMediterranean Journal of Mathematics13313251337Search in Google Scholar

Qiu, W., Chen H. and Zheng X., (2019), An implicit difference scheme and algorithm implementation for the one-dimensional time-fractional Burgers equations, Mathematics and Computers in Simulation, doi.org/10.1016/j.matcom.2019.05.017.QiuW.ChenH.ZhengX.2019An implicit difference scheme and algorithm implementation for the one-dimensional time-fractional Burgers equationsMathematics and Computers in Simulationdoi.org/10.1016/j.matcom.2019.05.017Search in Google Scholar

Mohebbi, A., (2018), Analysis of a Numerical Method for the Solution of Time Fractional Burgers Equation, Bulletin of the Iranian Mathematical Society 44,2 457–480.MohebbiA.2018Analysis of a Numerical Method for the Solution of Time Fractional Burgers EquationBulletin of the Iranian Mathematical Society442457480Search in Google Scholar

Yaseen, M., and Muhammad A., (2019), An efficient computational technique based on cubic trigonometric B-splines for time fractional Burgers’ equation, International Journal of Computer Mathematics, 1–14, doi.org/10.1080/00207160.2019.1612053.YaseenM.MuhammadA.2019An efficient computational technique based on cubic trigonometric B-splines for time fractional Burgers’ equationInternational Journal of Computer Mathematics114doi.org/10.1080/00207160.2019.1612053Search in Google Scholar

Asgari, Z., and Hosseini, S. M., (2018), Efficient numerical schemes for the solution of generalized time fractional Burgers type equations, Numerical Algorithms, 77,3, 763–792.AsgariZ.HosseiniS. M.2018Efficient numerical schemes for the solution of generalized time fractional Burgers type equationsNumerical Algorithms773763792Search in Google Scholar

Alam K., Najeeb, A. A., and Amir M., (2012), Numerical solutions of time-fractional Burgers equations: a comparison between generalized differential transformation technique and homotopy perturbation method, International Journal of Numerical Methods for Heat and Fluid Flow 22,2, 175–193.AlamK.NajeebA. A.AmirM.2012Numerical solutions of time-fractional Burgers equations: a comparison between generalized differential transformation technique and homotopy perturbation methodInternational Journal of Numerical Methods for Heat and Fluid Flow222175193Search in Google Scholar

Lombard, B., and Matignon, D., (2017), Numerical modeling of a time-fractional Burgers equation, WAVES, 219–220.LombardB.MatignonD.2017Numerical modeling of a time-fractional Burgers equationWAVES219220Search in Google Scholar

Saad, K. M., and Al-Sharif H.F., (2017), Analytical study for time and time-space fractional Burgers’ equation, Advances in Difference Equations, 2017,300, 1–15.SaadK. M.Al-SharifH.F.2017Analytical study for time and time-space fractional Burgers’ equationAdvances in Difference Equations2017300,115Search in Google Scholar

Murillo, J.Q., Yuste, S.B., (2011), An Explicit Difference Method for solving Fractional Diffusion and Diffusion-Wave Equations in the Caputo Form, J. Comput. Nonlinear Dynam. 6, 021014.MurilloJ.Q.YusteS.B.2011An Explicit Difference Method for solving Fractional Diffusion and Diffusion-Wave Equations in the Caputo FormJ. Comput. Nonlinear Dynam.6021014.Search in Google Scholar

Oldham, K.B., Spanier, J., (1974), The Fractional Calculus, Academic, New YorkOldhamK.B.SpanierJ.1974The Fractional CalculusAcademicNew YorkSearch in Google Scholar