Abstract

Fractional derivative is nonlocal, which is more suitable to simulate physical phenomena and provides more accurate models of physical systems such as earthquake vibration and polymers. The present study suggested a new numerical approach for the fractional diffusion-wave equation (FDWE). The fractional-order derivative is in the Riemann-Liouville (R-L) sense. Discussed the theoretical analysis of stability, consistency, and convergence. The numerical examples demonstrate that the method is more workable and excellently holds the theoretical analysis, showing the scheme’s feasibility.

1. Introduction

Nowadays, fractional calculus gains attention because of many applications in different sciences, such as physical and chemical phenomena [113]. Therefore, many numerical methods are used for the fractional differential equations [1421]. The fractional diffusion-wave equations (FDWE) have many applications, for example, continuous-time random walks, unification of diffusion, and wave propagation phenomenon [22, 23]. To find the solution of FDWE, different researchers used different numerical methods; for example, Hu and Zhanga [24] implemented a higher-order difference approach on the system of the FDWE and proved its stability, solvability, and convergence. Their numerical solutions reported that the high-order difference method is more efficient than the Crank-Nicolson method. Yang et al. [25] first converted FDWE into the integrodifferential equation and then apply the fractional multisteps method to solve the FDWE. They discretized the transform equation’s integration term and Lubich’s fractional multistep approach for space derivative and central difference schemes, respectively. Also, they show that the proposed method is convergent stable and solvable for . Brunner et al. [26] considered the numerical simulation of the 2-D time FDWE. They introduced artificial boundary conditions to transform the problem into the initial boundary value problem (IBVP) and produced a stable and convergence scheme using finite difference approximation. Sweilam et al. [27] solved two-sided space FDWE using a finite difference method. They studied the consistency and stability of both explicit and implicit methods. Also, they concluded that finite difference methods are relatively easy to use and can apply to different fractional differential equations easily. Later, Sweilam et al. [28] solved FDWE with the help of a weighted average finite difference method, which is an extension of the weighted average for the ordinary DWE. The noninteger order derivative is discretized using RL derivative and found that the suggested approach is stable. Liu et al. [29] used the novel finite difference scheme to solve the IBVPs for the FDWE. They reported that the proposed technique is unconditionally stable and convergent, having the order of convergence , where and represent the time step and mesh size, respectively. Dehghan and Abbaszadeh [30] applied A finite difference method to solve one and 2-D FDWE both space and tempered fractional forms. For fractional space derivative, they used Riemann-Liouville’s fractional derivative. They show using error estimate that the suggested scheme is convergence and unconditionally stable.

Sun and Wu [31] researched the fully discrete scheme to solve the solution of the BVP of FDWE. They determined that the proposed numerical approach is solvable, stable, and convergent with convergence order using the energy method. Furthermore, Du et al. [32] constructed a higher-order approximation to find the solution of FDWE with the order , where the fractional derivative is approximated using Caputo fractional derivative. The numerical results show that the Caputo fractional derivative gives convergence , which can be many apply to other fractional differential equations. Zhang et al. [33] proposed two new schemes, namely, the Crank-Nicolson (C-N), alternating direction implicit (ADI), and compact ADI to get a solution of 2-D FDWE. The close ADI scheme used the ADI approach for the time, and the relative difference technique is used for the spatial discretizations. Furthermore, they reported that the compact ADI approach is solvable, unconditionally stable, and convergent. The experiments confirmed that the ADI approach is computationally efficient in terms of CPU timing compared to the C-N ADI scheme for different values of without affecting the results. Huang et al. [34] used two approximations for the initial-boundary value time FDWEs, which is based on its integrodifferential equations. The proposed schemes were convergent with first and second order in temporal and spatial directions, respectively. Wei [35] used the semidiscrete strategy with the help of the finite difference method and fully implicit discrete scheme using the local Galerkin method for space time-fractional order derivatives, respectively, for the FDWE. The proposed approach results in a fully discrete scheme, which is convergent and unconditionally stable. Ali and Abdullah [36] formulated a numerical scheme for the fractional-order wave equation. They provided numerical solutions for the feasibility of the proposed scheme.

The above-cited literature means that fractional calculus is still a novel topic and needs more efficient numerical techniques to investigate the more feasible fractional differential equations (FDEs). This study is aimed at examining the highly accurate and robust numerical approach for fractional-order diffusion-wave equation. So far, the attempt has not been made to discretize the R-L integral and implement in the Riemann-Liouville fractional order derivative to approximate FDWE.

The remaining paper is organized as follows: Section 2 discussed the related preliminaries. Section 3 explains the methodology of the proposed scheme. Section 4 discusses the stability analysis. In Section 5, consistency and convergence are discussed. Section 6 provides numerical results, and finally, the conclusion is presented in Section 7.

Here, we consider the FDWE as following [36]: where the conditions are where is denoted the source term and represented the R-L of order lying between 1 and 2.

2. Preliminaries

The R-L fractional order derivative is defined as follows:

Here, the R-L integral can be defined as:

To approximate the Riemann-Liouville integral in equation (5), following the approach in [37]: using the Jumarie property [38] as: applying the Jumarie property , and

Lemma 1. The order R-L fractional integral of in defined as:

Lemma 2. The coefficient satisfies the below properties [37]: (i)(ii)(iii)Here, the constant as (iv)

3. The Proposed Scheme

The IDS for FDWE is developed in this section. The fractional order part is replaced by Lemma 1, and space derivative is approximated by central difference approximation. The step for space is and the step for time is where , and , , respectively. Using equation (4) in (1) at grid point , as follows

Eliminating the second order time derivative by applying backward difference approximation w.r.t time and then integrating from to and using trapezoidal rule for source term, we obtained substituting Lemma 1 in place of RL fractional integral in equation (11), the simplified approximated difference scheme is obtained for fractional diffusion-wave equation, as follows we know that , therefore where .

4. Stability

The Von Neumann technique is using to find the stability of the suggested scheme. Let represent the exact solution for equation (12), we get

The error is defined as . And satisfies equation (15) as follows:

Consider the growth factor in the form of a one Fourier mode as where and are the mode number and step size, respectively. Equation (17) can be the solution of the above error equation (16), as

Dividing both sides by , then replace , we obtained

Proposition 3. Suppose is the solution of equation (19), then we need to prove that

Proof. Let take in equation (19), to prove the proposition by induction method, as

By utilizing equations (13) and (17), we get and putting in equation (21). After simplification, the above equation becomes obtained the following relation

Let holds for

Using equations (19) and (23) and Lemma 2, we have

From Lemma 2, the value , so it is clear that the value

Here, , also , so it can written as . It shows that the suggested approach is unconditionally stable.

5. Consistency

To find the consistency analysis, suppose is the closed form solution, and is the approximated solution, and the function is the approximated scheme for the proposed equation at the mesh point . Then, denoted the local truncation error at .

Theorem 4. The local transcation error for the suggested scheme is

Proof. The obtained difference scheme in equation (12), can also be written as: using Taylor series, we obtained

The above equation can also be written as

We obtained the following truncation error

The IDS of fractional diffusion-wave equation is consistent if space step and time step approaches zero; then, truncation error approaches zero.

Theorem 5. The consistency and stability analysis are both the necessary and the sufficient for convergence. By Lax-equivalence theorem [33], the proposed scheme is convergent.

6. Numerical Results

Here, the FDWE examples reported to find the accuracy and feasibility of the scheme. The numerical examples are coded in Maple 15, i.e., the maximum error, is defined as follows:

The error is:

The convergence order is:

Example 6. In equations (1)–(3), take the source term and the closed solution .

Example 7. In equations (1)–(3) taking the source term and the closed solution .

7. Discussion

The numerical examples 5.1 and 5.2 are solved by the suggested numerical scheme in equation (12). Tables 14 represent the performance of the numerical results. The accuracy level increases as reducing the space and the time step; please refer to Tables 1 and 2 for different and same space and time steps, respectively. In Table 3, fixed and for various values of and , and find the convergence order completely demonstrates the analysis theoretically. Similarly, Table 4 also shows the accuracy for the second test example 5.2, which converges to the exact solution to reduce the error by reducing the step size. From the above discussion, it is clear that the proposed scheme proved our theoretical analysis.

8. Conclusion

A practical and fast numerical scheme has been developed for FDWE. The approximation is based on the discretization of Riemann-Liouville integral in Lemma 1. We have successfully proved the analysis theoretically of stability by mathematical induction, consistency, and convergence. The numerical result confirmed our theoretical analysis and demonstrated that the proposed scheme is fast convergent and more feasible. This approach can also apply to other types of FDWEs.

Data Availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

Acknowledgments

This research was supported by the Taif University Researchers Supporting Project Number (TURSP-2020/48), Taif University, Taif, Saudi Arabia.