Solving time fractional Burgers’ and Fisher’s equations using cubic B-spline approximation method

This article presents a numerical algorithm for solving time fractional Burgers’ and Fisher’s equations using cubic B-spline finite element method. The L1 formula with Caputo derivative is used to discretized the time fractional derivative, whereas the Crank–Nicolson scheme based on cubic B-spline functions is used to interpolate the solution curve along the spatial grid. The numerical scheme has been implemented on three test problems. The obtained results indicate that the proposed method is a good option for solving nonlinear fractional Burgers’ and Fisher’s equations. The error norms L2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$L_{2}$\end{document} and L∞\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$L_{\infty }$\end{document} have been calculated to validate the efficiency and accuracy of the presented algorithm.


Introduction
In recent years, most of the practical problems arising in different fields of science like biology, chemistry, physics, engineering, and mathematics lead to nonlinear fractional partial differential equations [1][2][3]. To solve these nonlinear mathematical models, several techniques have been developed. However, still it is noted that a large number of nonlinear fractional differential equations that are encountered in these areas are very difficult to solve numerically. The numerical solution of such kind of models has been the key interest of researchers due to their wide range of applications in real life.
In open literature, a lot of work is available on solving linear/nonlinear partial/fractional partial differential equations numerically by means of different techniques. The authors in [4] used the Galerkin finite element method to solve fractional diffusion and fractional diffusion wave equations. A B-spline collocation method was used by [5] for numerical treatment of fractional diffusion and fractional diffusion wave equations. Mehbobi et al. [6] employed radial basis functions to solve time fractional nonlinear Schrodinger equation. The authors in [7,8] investigated the approximate solutions of fractional telegraph equation by means of radial basis functions and a fully discrete local discontinuous Galerkin method respectively. Nonlinear Fisher's equation has been studied in [9,10] by means of variational iteration method (VIM) and homotopy perturbation method (HPM) respectively. Yokus and Kaya [11] employed the Cole-Hopf transformation and the expansion methods to solve time fractional Burgers' equation. Jin et al. [12] presented the numerical analysis of nonlinear sub-diffusion equations. Li et al. [13] proposed Galerkin's finite element approach for numerical solution of nonlinear time fractional parabolic equations. A Galerkin scheme based on Crank-Nicolson method has been used in [14] for solving nonlocal coupled parabolic problems. A cubic B-spline based Galerkin approach has been applied by Iqbal et al. [15] to solve Schrodinger equation. In general, the finite element method based on third degree basis spline functions has not been widely used for solving time fractional problems and this will be the focus of our paper. In this work, we consider the following generalized time fraction partial differential equation: subject to the term of initial and boundary conditions where η, ν, β are parameters, H(z, t) is the source term, q is a positive integer, and g(z), g i (t) are known functions. There are several interpretations of fractional order derivatives; however, we use Caputo's definition as follows [16]: where M is the smallest integer exceeding α. If we take η = 0 and q = 1, equation (1) takes the following form: The above equation, known as time fractional Fisher's equation [17], is a mathematical model for the temporal spatial propagation of a virile gene in an infinite domain. This model is also used in chemical kinetics [18], auto catalytic chemical reactions [19], neurophysiology [20], nuclear reactor theory [21], branching Brownian motion process [22], and flame propagation [23]. If we set β = 0 and q = 1 in (1), it reduces to Burgers' equation This is the simplest nonlinear model in fluid dynamics for diffusive waves. It appears in various physical issues such as viscous sound waves, waves in fluid-filled viscous elastic pipes, magneto hydrodynamic waves in a medium with finite electrical conductivity [24] and one-dimensional turbulence.
where D j (t n+1 )s are unknowns to be identified and B j (z) are third degree basis spline functions defined as follows [27]: In this way, equation (10) can be written as follows: For the sake of simplicity, we denote Y (z r , t n+1 ) by Y n+1 r and discretize equation (13) along the spatial grid as follows: From (11) and (12), we obtain By plugging (15) into (14), we obtain the following set of equations: where a n r = φ n r . Two more relations are obtained from end conditions (3): Equations (16)- (17) can be expressed in a matrix form as follows: The above system is solved for (D n+1 Using this vector in (11), we get the approximate solution at (n + 1)th time level for n = 0, 1, 2, . . . , N . However, to start the iterative process, we need to evaluate the initial vector. The initial conditions of the problem under consideration will help to compute the initial vector ) T as follows: Discretization of (2) results in the following set of equations [28]: Using (15) in the above system, we get the following matrix form: This above equation, (N + 3) × (N + 3) matrix system, is solved using Thomas algorithm. For numerical computations, MATLAB software has been used.

Stability analysis
For the sake of simplicity, we discuss the stability of the presented scheme for the following linear force free case of problem (1): Employing the method described in Sect. 3 for θ = 1, the fully discretized scheme for equation (21) is given by where Now, we bring the Fourier method into play for stability analysis. Let Ω n denote the growth factor andΩ n be its estimated value. The error Υ n at nth time stage is given by where . (24) The end conditions (2)-(3) are also satisfied by error term as follows: Now, we introduce a mesh function Also, Υ n (z) can be written in terms of Fourier series as follows: where (using Parseval's equality).
Now, we plug Υ n r = ρ n e rhι into (24) and simplify to get the following expression: where ι = √ -1 and = 2mπ b-a . After some simplification, (30) takes the form where . It can be noted that |χ| ≥ 1. Now, we make use of mathematical induction to show that |ρ n | ≤ |ρ 0 | for all n.
Using (29) and (33), we have which concludes that the proposed scheme is unconditionally stable.

Convergence analysis
In this section, we derive the uniform convergence of the proposed scheme.
Let y, Y be the exact and approximate solutions for the given problem andỸ = R+1 l=-1 l B l be the spline approximation for Y with the interpolation conditions Ly(z r , t) = LY (z r , t) = g(z r , t), r = 0 : 1 : R, then LỸ (z, t) = g((z r , t).
At t = t n , the problem can be expressed in terms of difference equation L(Ỹ (z r , t) -Y (z r , t)) as follows: where ℘ n r = D n r -n r , r = 0, 1, . . . , R. From the end conditions, we can write If n r = h 2 [ n rr n ] for r = 0, 1, . . . , R, then using Lemma (4.1), we can write n r = h 2 n rr n ≤ ‫ג‬h 4 .
Let n = max{| n r |; 0 ≤ r ≤ R},Ẽ n r = |℘ n j |, andẼ n = max{|E n r |; 0 ≤ r ≤ R}. At n = 0, equation (36) takes the shape Making use of the initial condition, we have E 0 = 0 ⇒ Using small grid spacing with end conditions, we can write where ‫ג‬ 1 does not contain h. Using mathematical induction, we can easily prove that Using Lemma (4.2) , we get where ‫ג‬ and ς are constants.

Applications and discussion
In this section, we present some numerical experiments. L 2 and L ∞ norms are used to analyze the precision of the suggested technique. and Example 1 Consider the following time fractional Burgers' equation:    The three-dimensional plot of solution surface is shown in Fig. 2. Table 1 presents the comparison of approximate and exact solutions. We can see that the numerical results are reasonably precise. The error standards for L 2 and L ∞ were measured and displayed in Table 2. The numerical order of convergence has been reported in Table 3.
Example 2 Consider the following time fractional Burgers' equation: IC: y(z, 0) = 2z 2 , 0<z < 1, BCs: y(0, t) = 0, y(1, t) = 2(t + 1), t ≥ 0.   The source term is given by The closed form solution to this problem is y(z, t) = 2z 2 (i + t). The exact and approximate solutions of equation (42) Fig. 3(a); from there, we can see that both solutions agree with each other well. Figure 3(b) depicts the symmetry in both solutions for different time levels t = 0.5, 0.75, and 1. Figure 4 displays three-dimensional representation of physical behavior of the solution curve. The exact and numerical outcomes are compared in Table 4. One can see a very close agreement of the exact and approximate solutions in Table 5. From the tabular information, it can be concluded that error reduces as the value of t varies. Table 6 shows the calculations of numerical order of convergence.

Concluding remarks
In this work, the authors have presented a B-spline collocation method based on third degree basis spline functions and finite difference formulation for numerical investigation of time fractional Burgers' and Fisher's equations. Caputo's interpretation for time fractional derivative has been considered. The fractional time derivative has been discretized by L1 formula, and B-spline functions are used to interpolate the solution along the spatial grid. The stability and convergence analysis of the proposed method is carried out. The numerical results of three different test problems reflect the applicability of the proposed scheme. The numerical outcomes validate the accuracy and efficiency of the presented method.