Numerical simulation method of lines together with a pseudospectral method for solving space-time partial differential equations with space left- and right-sided fractional derivative

Purpose –This study aims to use new formula derived based on the shifted Jacobi functions have been defined and some theorems of the leftand right-sided fractional derivative for them have been presented. Design/methodology/approach – In this article, the authors apply the method of lines (MOL) together with the pseudospectral method for solving space-time partial differential equations with space leftand right-sided fractional derivative (SFPDEs). Then, using the collocation nodes to reduce the SFPDEs to the system of ordinary differential equations, which can be solved by the ode45 MATLAB toolbox. Findings – Applying the MOL method together with the pseudospectral discretization method converts the space-dependent on fractional partial differential equations to the system of ordinary differential equations. Originality/value – This paper contributes to gain choosing the shifted Jacobi functions basis with special parameters a, b and give the authors this opportunity to obtain the leftand right-sided fractional differentiation matrices for this basis exactly. The results of the examples are presented in this article. The authors found that the method is efficient and provides accurate results, and the authors found significant implications for success in the science, technology, engineering and mathematics domain.


Introduction
The recent development in the last few decades has shown that most of the complex system in engineering and other several phenomena can be accurately modeled using partial differential equations with fractional order. This contributed to a great development in several areas such as biotechnology, chemistry, signal and image processing, finance and many others [1][2][3][4][5]. The main aim of this paper is to introduce an efficient numerical method to approximate the fractional partial differential equation (FPDE) of the form vuðx; tÞ vt ¼ sðx; tÞ þ c þ ðx; tÞ R 0 D α x uðx; tÞ þ c − ðx; tÞ R x D α ' uðx; tÞ; (1) we also assume the initial and Dirichlet boundary conditions: uðx; 0Þ ¼ FðxÞ; (2) uð0; tÞ ¼ uð'; tÞ ¼ 0; on the domain of space 0 < x < l and time 0 ≤ t ≤ T, and we consider that the parameter α is the fractional order where 1 < α ≤ 2. And a source or sink term is the function s(x, t). The functions c þ (x, t) ≥ 0 and c − ðx; tÞÞ ≥ 0 represented the interpreted as transport-related coefficients or the advection, and the diffusion coefficients. Many approximation methods in the numerical analysis have been a survey to solve space-dependent on fractional partial differential equations (SFPDEs), and the target of the main subject of these methods in terms of convergence to real solutions, the stability of methods, and order of accuracy and value of error.
One of the best types of methods for solving SFPDEs numerically is by discretization of the space variable without the time variable. These kinds of method are referred to as method of lines (MOL). In this method, the spatial dimensions of the space variable can be discretized by using diverse techniques such as mesh methods, or meshless methods [6][7][8][9][10][11][12]. In general, these methods convert the SFPDEs to a system of ordinary differential equations, or differential-algebraic equations based on the type of boundary conditions [13]. There are many types of partial differential problems that have been solved by the MOL, we refer to [7,10,[14][15][16][17][18][19], and therein. Also, many different methods have been discussed for SFPDEs, we refer to [20][21][22][23][24], and therein.
Our target of this work is to use the advantage of the pseudospectral method based on the shifted Jacobi functions together with MOL and employed the collocation method to approximate the solution to the SFPDE (1)- (3).
The structure of this paper was arranged in the following way: In Section 2, preliminaries, the definitions of fractional derivatives and some notations of Jacobi-Gauss (JG) nodes. In Section 3, the new numerical technique for solving SFPDE (1)-(3) is presented. In Section 4, the illustrative examples were included to demonstrate the validity and applicability of the proposed method. In Section 5, a brief conclusion and some remarks.

Preliminaries and definitions
We introduce several important basic definitions and properties of fractional Riemann-Liouville integrals and derivatives and JG nodes. where r is any real number. For more information see [25][26][27][28].
Definition 2.3. The fractional Riesz operator of the order α is defined as [29]: where x belongs to the finite interval 0 ≤ x ≤ l and n À 1 < α ≤ n such that α ≠ 1.
According to the definition of Jacobi polynomials, various types of functions with fractional order have been constructed based on them. For more information see [26].
We insert some properties of the classical Jacobi polynomials in the following: And the following theorem plays a cornerstone to establish our proposed method.
Then for all x belong to the interval [0, 1] we have: It is clear from combined properties and theorems that have been inserted as it hold the proof of theorem in References [30][31][32]. And there are many applications of the above theorem that can be found in References [33][34][35][36][37][38][39][40].

Th presented pseudospectral discretization method
The important step of our method includes to discretize the space variable of the unknown function u(x, t) that appears in SFPDE (1)-(3). To get started, we use the pseudospectral method based on the shifted Jacobi function of the parameters a 5 1, b 5 1, together with the collocation JG nodes. Let the domain space of the function u(x, t) belongs to the interval [0, 1], then the JG nodes will be correspondent to the interval [0, 1] as: Let us begin to approximate the unknown function u(x, t) byũ n ðx; tÞ as follows: where such that ' j (t), j 5 1, 2, 3, . . ., n, are the Lagrange basis polynomials based on the JG nodes From the Kronecker properties in the JG nodes, the ' j (x) and w j (t) are satisfied as follows: Let us denote the shifted Jacobi polynomials b P ða;bÞ k ðtÞ, k 5 0, 1, . . ., n, which can be written as: b P ða;bÞ k ðtÞ ¼ P ða;bÞ k ð2t À 1Þ; k ¼ 0; 1; 2; 3; . . . ; n: Since the degree of polynomial ' j (x) is n À 1. Then, the expansion of ' j (x) in terms of the shifted Jacobi polynomials b P ð1;1Þ k−1 ðxÞ, k 5 1, . . ., n, can be written as follows: Now, by multiplying both sides of Eqn (15) by xð1 − xÞ b P ð1;1Þ k−1 ðxÞ, and using the orthogonality property (4) to the shifted Jacobi polynomials in the interval (0, 1), get that By using the JG quadrature rule (6), and then using the Kronecker property (14), we get By using above Eqn (17), then we can define w j (t) as: where Now dependent on the fact of Eqns (19), (18), (15) and (11) Numerical simulation method Similarly, we approximate the initial and boundary conditions as follows: uðx; 0Þ 'ũ n ðx; 0Þ ¼ X n k¼0 a k ð0Þw k ðxÞ ¼ FðxÞ; (21) uð0; tÞ 'ũ n ð0; tÞ ¼ 0; uð'; tÞ 'ũ n ð'; tÞ ¼ 0: By substituting the approximations functions (20), (21) and (22) into the problems (1)-(3), we get: . . . ; n be the zeros of b P ð1;1Þ nþ1 ðxÞ. By using collocating nodes at . . . ; n the above Eqn (23), will change to the system of algebra equations dependent on the time variable: In general, Eqn (24) can be rewritten in the following matrix form:  ; . . . ; a n ð0Þ T . In the final result, the main problems (1)-(3) will be reduced to the following system of ODEs with initial conditions:

On the derivation of the left and right fractional differentiation matrices
The coefficients d þ ij and d − ij , appear in the system of Eqn (25), we need to compute them accurately and efficiently. In the following, we present an efficient method to compute the coefficients d þ ij and d − ij . It is interesting to point out that, we choose w k (x) as the basis functions in the pseudospectral method. The best feature of these bases is that, we can obtain a closed form of the left and right fractional derivatives, by the use of the next two theorems 2 and 3, which may simplify the discretization stage.
Theorem 2. The left fractional differentiation matrix D α þ ¼ ½d þ ij can be obtained as the element of it in the explicit form as follows: and ζ k ðxÞd Γð1 þ kÞ Γðk þ 1 À αÞ where b τ i are the shifted JG nodes defined in Eqn. (10).

ðxÞ:
Proof: Depending on the fact we get it in Eqn (19), and by using property (8) for the special parameter a 5 2 and b 5 1, we obtain: k−2 ðxÞ: Now, combining Theorem 1 with the above equation, we obtain:

AJMS
The inverse of matrix P in the explicit form can be found by using the next theorem.
Theorem 3. The element of the matrix P can be written in the explicit form as follows: where w ð1;1Þ s ; s ¼ 0; 1; 2; . . . ; n are the JG weights with respect to the weight function w (1,1) Proof: Start the proof by noting the fact that: On the other hand, we can approximate the above integral by using the JG quadrature rule with respect to the weight function t(1 À t) on [0, 1] as follows: j, k 5 0, 1, . . ., n to reach the mass matrix P, we can rewrite the above equation as: Thus by rewriting the above equation in the matrix form, we have: which concludes the proof. " The above theorem shows that the matrix P is nonsingular, and we can rewrite the implicit initial value problem (26) to the following explicit initial value problem: It should be noted that the IVP (35) can be solved by various well-known software. We solve this IVP by the ode45 MATLAB toolbox. In the next section, we provide some numerical examples to check the efficiency of the proposed method.
Numerical simulation method

Numerical results and comparisons
Our goal to solve problems (1)-(3) and find the unknown function u(x, t), by approximate the space variables via pseudospectral method based on JG nodes, and then solving the new IVP by the ode45 MATLAB toolbox. We will begin in the first step by transforming interval (0, ') to interval [0, 1], where the function u(x, t) defined on the interval x ∈ (0, '), t ∈ (0, T), and the collocation points {τ i , i 5 0, 1, 2, . . ., m} belong to the interval [À1, 1], then x ¼ 1 2 ½ð'Þτ i þ ' are the corresponding collocation points on [0, ']. Also if we transform this interval to [0, 1], then, In order to confirm the utility of the presented method, we apply the method to solve some IVPs. The proposed method has been implemented with MathWorks MATLAB 2017a in a personal computer 3.5 GHz Core i7 PC with 8 GB of RAM. By applying the presented method with n 5 40, the obtained functionũðx; tÞ is plotted in Figure 1. Moreover, the error function E i;j ðuÞduðb τ i ; t j Þ −ũðb τ i ; t j Þ is plotted in this figure too.
From that, this figure can be seen that our numerical solution is in good agreement with the analytic solution. Consumed CPU time for an accurate solution is obtained in just 21.357 seconds.
By applying the presented method with n 5 20 and α 5 0.2, the obtained functionũðx; tÞ is plotted in Figure 2. Moreover, the error of the obtained functionũðx; tÞ is plotted in this figure too. Consumed CPU time for an accurate solution is obtained in just 23.427 seconds. In Table 1, E 2 n;m ðuÞ is the two-norms of the error for the obtained functionũðx; tÞ, which is defined as, From Table 1, it can be seen that norms of errors for various values of n and α appear that the small numbers of n, with that an accurate solution.
Example 3. We consider the following SFPDEs [29]: vuðx; tÞ vt with initial and boundary conditions: uðx; 0Þ ¼ sinð4xÞ; uð0; tÞ ¼ uðπ; tÞ ¼ 0; Numerical simulation method Figure 3 shows the approximation function ofũðx; tÞ, for order of fractional derivative α 5 1.25 and ℵ 5 1.5 and the time for 0 < t < 1. Observed that if value of α is decreased from 2 to 1 the amplitude of the sinusoidal solution behavior is increased. In Table 2, the values of uðx; 1Þ in various values of x and n with the value of α 5 1.5 Example 4. We consider the following FPDEs [42,43] vuðx; tÞ vt þ uðx; tÞ ¼ v α v x j j α uðx; tÞ þ sðx; tÞ; 1 < α ≤ 2; with initial and boundary conditions: The exact solution is of the form u(x, t) 5 x 2 (L À x) 2 e Àt The problem (42) has been solved in [43] by proposed a class of difference scheme based on the parameter spline function and improved matrix transform method. Also, the authors in [42] have solved this problem (42) by

AJMS
proposing a space-time spectral algorithm based on the shifted Jacobi tau technique. To demonstrate the accuracy of our proposed method, in Table 3, compare the absolute errors ju(x i , 0.1) À u(x i , 0.1)j, i 5 1, . . ., 9 with the numerical method proposed in [42,43] and our results for various choices of α.

Conclusion
In the present paper, we developed an efficient and accurate method for solving SFPDE.
Applying the MOL together with the pseudospectral discretization method converts the SFPDE to the system of ordinary differential equations. Choosing the shifted Jacobi functions as a test basis with special parameters a, b gives us this opportunity to obtain the left-and right-sided fractional differentiation matrices for this basis exactly. Four examples have been solved, and the results are reported. These results show that our method is efficient and provides accurate results, whereas a small number of JG nodes are used based on the collocation method. Obtaining some theoretical estimates for the approximation errors would be desirable. Table 3. Example 4: Comparison of absolute errors of our scheme with scheme in [42,43] AJMS