A pseudo-spectral method based on reproducing kernel for solving the time-fractional diffusion-wave equation

In this paper, we focus on the development and study of the finite difference/pseudo-spectral method to obtain an approximate solution for the time-fractional diffusion-wave equation in a reproducing kernel Hilbert space. Moreover, we make use of the theory of reproducing kernels to establish certain reproducing kernel functions in the aforementioned reproducing kernel Hilbert space. Furthermore, we give an approximation to the time-fractional derivative term by applying the finite difference scheme by our proposed method. Over and above, we present an appropriate technique to derive the numerical solution of the given equation by utilizing a pseudo-spectral method based on the reproducing kernel. Then, we provide two numerical examples to support the accuracy and efficiency of our proposed method. Finally, we apply numerical experiments to calculate the quality of our approximation by employing discrete error norms.


Introduction
The use of the fractional-order derivative has become popular due to its nonlocality property, which is an intrinsic property of many complex systems. The fractional-order derivative has recently been applied in modeling different phenomena including viscoelasticity, financial modeling, nanotechnology, control theory of dynamical systems, random walk, anomalous transport, biological modeling, and anomalous diffusion as well. For further applications of the fractional order derivative in the fields of engineering, physical sciences, we may refer to [1][2][3][4][5][6][7][8][9][10][11] and [12][13][14][15][16] and the references cited therein.
The fractional partial differential equations (FPDEs) have been deployed in recent years as a powerful tool in nonlocality and spatial heterogeneity modelling. Many applications of the fractional models can be found in [17][18][19][20][21]. The fractional diffusion equation assumed in this paper covers not only the classic state of the heat equation but also the kernel of many other FPDEs. This equation describes the anomalous diffusion of particles. For applications of this equation, we refer to the one presented in [22], which describes the transfer processes with a long memory and the water transport in the soil model [23].
Many mathematical models are usually based on the diffusion model such as the diffusion equation on fractals [24], Fischer information theory [25], and so on. Moreover, the implication of the fractional in time (sub-diffusion) and fractional in space (super-diffusion) has been observed in the solution profiles in many fractional models [26,27]. The superior capabilities of fractional differential equations to accurately model such processes have raised significant interest in assaying numerical methods for obtaining the solutions to such problems [28][29][30]. The analysis we present in this paper depends on the following time-fractional diffusion-wave problem: is the time interval satisfying 0 < T < +∞; u : × J → R is a sufficiently differentiable function, and 1 < α < 2 is the order of the fractional derivative, and the time-fractional derivative ∂ α t is the Caputo fractional derivative of order 1 < α < 2 defined by The reproducing theory was investigated by Mercer in 1909 [31]. He named functions satisfying the reproducing property as "positive definite kernels". Around 1948, Aronszajn [32] systematized the concept of reproducing kernels. From 1980, Cui and co-workers [33,34] have been pioneers and beginners in the numerical analysis of linear and nonlinear problems using the "reproducing kernel Hilbert space method". Recently, a lot of research has been done to solve several linear and nonlinear problems using the theory of reproducing kernel [35][36][37][38][39][40][41][42][43][44][45][46].
The aim of this paper is to introduce a finite difference/pseudo-spectral method based on reproducing kernel (RK) for solving the time-fractional diffusion-wave equation (1.1)-(1.3). This paper spreads over four sections, including the introduction. In Sect. 2, we present a finite difference/pseudo-spectral method based on a reproducing kernel (RK) for solving the time-fractional diffusion-wave equation. In Sect. 2, we solve some test problems and derive several results. In Sect. 4, we present some concluding remarks.

Implementation of the method 2.1 Discretization of Caputo derivative and semi-discrete scheme
First, we obtain the semi-discrete scheme for (1.1)-(1.3). The discretization of J is performed with a constant time step τ = T N , where N ∈ N * . Denote t n = nτ for N = 0 : N . Let u n = u(x, t n ). For a discrete function {u n } N+1 n=0 , we provide some preliminaries concerning the approximation of the time fractional derivative ∂ α t u(x, t) with 1 < α < 2. A Caputo derivative approximation formula (CDAF) for ∂ α t u(x, t n+1 ) with 1 < α < 2 can be defined as a linear combination of the discrete second time derivatives {∂ 2 u j } n+1 j=1 [47] ∂ α t u n+1 = and R n+1 1 (u) is the local truncation error such that The following lemma summarizes some properties of the coefficients b j which will be used in this paper.

Scheme I Given
For the convenience of discussion, define the linear operator L as follows: Therefore, a semi-discrete problem can be converted into the following equivalent:

A pseudo-spectral kernel-based method
Now, we employ a pseudo-spectral method based on RK to discrete the space direction and obtain a full-discrete scheme of (2.3). To obtain this, we need some notations and preliminaries.
We now give background material and preliminaries, which are used in the following sections. Recall that a real reproducing kernel Hilbert space (RKHS) on a nonempty abstract set is a particular type of a real Hilbert space H of functions that satisfies the following additional property (called reproducing kernel property): for each x ∈ , there exists K(x, ·) ∈ H (R : × − → R) such that, for every u ∈ H, one has (2.4)

Definition 2.2 (See [40]) A Hilbert space H of real functions on a set is called an RKHS
if there exists an RK K(x, ·) of H.  (2.5) and the norm u H is defined by The space H p [0, b] is an RKHS and the RK K y (x) can be denoted by [40] K y (x) = (2.7) With the help of pseudo-spectral method based on RK, we will illustrate how to derive the numerical solution. Now, we will give the representation of a numerical solution to the semi-discrete problem (2.3) be a distinct subset of . We consider the finite-dimensional space The semi-discrete problem can be written into following equivalent form An approximant U n+1 M to U n+1 can be obtained by calculating a truncated series based on trial functions as follows: To determine the interpolation coefficients {α n+1 j } M j=1 , the set of collocation conditions is used by applying (2.3) to B M . Thus where the functional λ i , (1 ≤ i ≤ M) is defined by applying the differential operator followed by a point evaluation at x i ∈ B n . In general, a single set M := {λ i } M i=1 of functionals contains several types of differential operators.
The arising collocation matrix K is unsymmetric and has the ij-entries: where the subscript x in λ x j indicates that λ x j applies to the function of x.
Therefore the unknown coefficients α n+1 j , j = 1, 2, . . . , M, can be obtained by solving the following system: We know that The following matrix vector form is achieved by differentiating (2.12) with respect to x and evaluating it at the point girds x i ∈ B M : and therefore Now, by using (2.13), we can write and

Nonsingularity of the collocation matrix
From Lemma 2.5 it is clear that the functions λ x j K(x, ·) for λ j ∈ M are continuous. Furthermore, note that {x i } i≥1 is a dense set. Therefore M j=1 c j λ x j K(x, ·) = 0, which implies that c j = 0 (j = 1, 2, . . . , M). This completes the proof.  (2.17) which implies that c j = 0 (j = 1, 2, . . . , M), and this completes the proof.
From Lemma 2.8 and Theorem 2.7, we can derive the following theorem.

Error analysis
Applying the Gram-Schmidt orthogonalization process to {ψ 1 , ψ 2 , . . . , ψ M }, we can obtain (2.18) Therefore {ψ 1 , ψ 2 , . . . , ψ n } is an orthonormal basis for U M . Therefore, we can write the interpolant U n+1 M (x) to U n+1 at B n in the following form: (2.20) Proof Using the reproducing property, we have Thus, the proof is completed.

24)
where c is a real constant.
Proof According to Lemma 2.11, we have

Illustrative test problems
We have studied some example tests to illustrate the performance of the proposed methods. We show the stability and accuracy of the proposed method for different values of M and N . As the exact solution is known, the root mean square error L rms and the maximum absolute error L ∞ are measured with the following formulas: Example 3.1 In this example, we deal with the following problem: The exact solution u is given by u(x, t) = t 2 sin(πx). The proposed method in the previous section is tested on this problem with In Tables 1 and 2, we present, the root mean square error L rms , the maximum absolute error L ∞ , and the convergence ratio in the computed solutions for Example 3.1 with α = 1.2, 1.4, 1.8, 1.9. In a considerable number of cases, an exciting agreement between the results is observed, which confirms the excellent validity of the proposed method.
In Tables 3 and 4, we present the root mean square error L rms , the maximum absolute error L ∞ , and the convergence ratio in the computed solutions for Example 3.2 with α = 1.3, 1.5, 1.7, 1.95. In a considerable number of cases, an exciting agreement between the results is observed, which confirms the excellent validity of the proposed method.

Conclusion
In this paper, a finite difference/pseudo-spectral method is presented to solve timefractional diffusion-wave equation. The method is based on a finite difference method in a temporal direction to obtain a semi-discrete configuration, whereas a pseudo-spectral method based on RK introduces spatial discretization. The implementation of the proposed method is very simple and has reasonable accuracy. It can be seen from the error norms and numerical results that the proposed method for solving time-fractional diffusion-wave equation is an excellent method.