Control and Stability Analysis of Complex Dynamical Systems involved in Finance, Ecology and EngineeringView this Special Issue
A Numerical Method for Time-Fractional Reaction-Diffusion and Integro Reaction-Diffusion Equation Based on Quasi-Wavelet
In this research work, we focused on finding the numerical solution of time-fractional reaction-diffusion and another class of integro-differential equation known as the integro reaction-diffusion equation. For this, we developed a numerical scheme with the help of quasi-wavelets. The fractional term in the time direction is approximated by using the Crank–Nicolson scheme. The spatial term and the integral term present in integro reaction-diffusion are discretized and approximated with the help of quasi-wavelets. We study this model with Dirichlet boundary conditions. The discretization of these initial and boundary conditions is done with a different approach by the quasi-wavelet-based numerical method. The validity of this proposed method is tested by taking some numerical examples having an exact analytical solution. The accuracy of this method can be seen by error tables which we have drawn between the exact solution and the approximate solution. The effectiveness and validity can be seen by the graphs of the exact and numerical solutions. We conclude that this method has the desired accuracy and has a distinctive local property.
Differential and integral calculus is a vital and broad branch of mathematics [1, 2]. First of all, the integer-order differential and integer calculus came into existence. About a few decades ago, the concept of fractional calculus came into existence. Day by day the new concept, new definitions, and new terminology are emerging. The concepts and ideas of fractional calculus were evolved by Liouville and Abel. The more basic fundamentals and concepts are available in the literature [3, 4]. This branch of mathematics allows us to generalize the order of derivatives from the integer order to the real and arbitrary order. There are so many types available in the literature of fractional calculus. Some of them are Caputo, Hadamard, Riemann–Liouville, and Grunwald–Letnikov. As time passes, in the development of this branch, the variable-order derivative concepts are introduced in which orders of derivatives are functions of variables. The most prevailing derivatives are power law kernel derivatives. But, in recent few years, the derivatives with nonsingular kernel are introduced. Derivatives with the Mittag–Leffler kernel are known as Atangana–Baleanu derivatives, and derivatives with the exponential kernel are known as the Caputo–Fabrizio derivative. The fractional differential equations are in trends to study. There are so many physical phenomena in which fractional differential equations depict and represent the behavior and dynamics of the phenomena in a more accurate way than integer-order differential equations. To find out the solution of the fractional differential equation and fractional partial differential equations is a tough task as it is not possible always to find out the exact solution. The researchers started to develop the methods and schemes for finding the numerical solution of fractional differential equations. Many numerical methods are developed and available in the literature. Some of them are the fractional differential transform method , predictor-corrector method , Adomian decomposition method , homotopy perturbation method , and generalized block pulse operational matrix method . The operational matrix method is a spectral method which is easy to be applicable for many types of fractional differential equations whether they are of integer, fractional, and variable type. Its results are so accurate, and it is easy to apply. Some of them are Gegenbauer wavelet , Haar wavelets , Chebyshev wavelets , and Legendre wavelets . The stability of discrete-time delayed switched systems and time delayed systems is given in the literature [13, 14].
FDEs, FPDEs, and fractional integro-differential equations are used to model many physical phenomena [15–17], which have arisen in physics, chemistry, biology, electrical engineering, nuclear reactor, viscoelastic, fluid dynamics, population dynamics, and medical science. To show the memory effect of the system, in the branch of nuclear reactor and thermo elasticity, we use the fractional differential equations. But, while modeling, the effect of history is being ignored. This can be taken into account by adding a term of integration in these fractional PDEs. The integro-differential equations have so many applications in control theory, financial market, fluid dynamics, aerospace systems, and viscoelastic mechanism. The most occurring physical process in nature is the diffusion process. In this process of diffusion, the molecules or energy or any physical quantity is transferred to the low concentration area from the high concentration area. If another process known as reaction is added to the process of diffusion, then it is called the reaction-diffusion process. In this process, molecules react with the fluids. This phenomenon of diffusion can be depicted by the following PDE known as the reaction-diffusion partial differential equation:where represents the reaction term and the first term on the right-hand side of the equation represents the diffusion process with as the diffusion coefficient.
We divided our article into many sections as follows: Section 1 contains the preliminary definitions which will be used throughout the article. This section contains the theory and methods of the quasi-wavelet to find out an approximation of an unknown function. The taken model of integro reaction-diffusion and reaction-diffusion is given in Section 3. How to deal with discretization of spatial and time-fractional terms with Dirichlet boundary conditions is also given in Section 3. The numerical results which depict the effectiveness and accuracy of our proposed method are present in Section 4. The conclusion of this paper is contained in Section 5.
2. Preliminary Definition about Fractional Derivative and Quasi-Wavelets
In this section, we described the definitions of fractional derivatives which will be used throughout the article. There are two types of definitions available in the literature of the fractional derivative, i.e., Riemann–Liouville and Caputo definitions. We have taken Caputo’s one because it is more reliable in modeling.
2.1. Integration and Differentiation of Riemann–Liouville Type
The R-L-type definition of fractional integration is given as follows :
The definition of fractional differentiation in the R-L sense is given by the following definition:
2.2. Definition of Fractional Derivative in Caputo Sense
The following equation represents the definition of fractional differentiation in Caputo sense:where the parameter represents the integer. A resemblance of the Caputo definition with the integer-order derivative is as follows:with being a constant.where represents the floor function. All these types of fractional operators follow the linear property:with and being constants. A useful relation between the R-L definition and Caputo is given as follows:
2.3. Approximation of Functions with the Help of Quasi-Wavelets
The method based upon quasi-wavelets nowadays is developing in the form of the local spectral collocation method. It is used for seeking the numerical solution of the fractional-order partial differential equation and integral equations. The singular convolution used in this method is a significant transformation. We define this transformation as follows:where is known as the test function and is a singular kernel. We know that the construction of the family of wavelets is done with the help of a function known as the mother wavelet . We use the operations of dilation and transformation to form the family of wavelets:
We have used two parameters here. The parameter is the dilation parameter, and is the translation parameter. The kernel known as Shannon’s delta sequence is used in the construction of quasi-wavelets. It is defined aswhere and is a well-known Dirac delta function. Shannon’s delta sequence kernel produces the basis for the Paley–Wiener reproducing kernel Hilbert space  for a given which is a subspace of . We can produce uniquely a function as follows:
The sampling scaling function can be written in another form as follows:
We have denoted the sampling points around the space by which are sampling points centered around . The discrete form of function which belongs to the space is presented by the following infinite series expansion:
The above approximation can be rewritten with the help of Shannon sampling theorem at the Nyquist frequency in the band-limited signal in :where the value of parameter is and the grid size in spatial direction is denoted by . We improved the localized asymptotic behavior of Dirichlet’s delta sequence kernel by using the regularizer . In the literature, many regularizers are introduced, and we use here the following Gaussian type regularizer:where represents the width parameter of the Gaussian envelope. The parameters and are related by the expression , with as the computation parameter. The sampling scaling function is modified as follows by using the Gaussian regularizer:wherewhere the Gaussian-regularized sampling scaling function follows the property of the quasi-scaling function. So, we can approximate an unknown function belonging to the space with the help of the quasi-scaling function:
Here, the infinite sampling points have been used, but while commutating in practical manner, we cannot use these infinite points. By seeing this problem, the domain of computation is restricted to the finite sampling points. For this purpose, we select . The simplified form of the above equation is as follows:
Now, for the approximation of integer-order derivatives of order of the function ,
The sampling points vary from to , i.e., there are sampling points. We give some of the formula of , and  for the calculation purpose:
3. Mathematical Model and Used Proposed Method
In this section, we propose the integrated mathematical model of the reaction-diffusion and integro reaction-diffusion equation. In this section, we develop the Crank–Nicolson scheme with the combination of the quasi-wavelet-based numerical method to find the numerical solution of the taken model:
We have taken the following Dirichlet boundary conditions and initial conditions:where the term denotes the forced function and is known as a kernel.
3.1. Using Crank–Nicolson Scheme for the Time Discretization
To discretize the time derivatives present in model, we will use the Crank–Nicolson scheme. We divided the time interval into discrete points as where and time length is noted by . The notations and represent the values of and at time . The term from the model is approximated as at the grid point :
In aberration form, we can rewrite as
Hence, we obtainwhere
The bound of is as follows after simplifying:
The trapezoidal technique is implemented to deal with the integral term :
Taking the exponential form of the kernel as , then
3.2. Use of Quasi-Wavelet-Based Approximation for the Discretization of Spatial Derivatives
By using the approximation by quasi-wavelets given in Section 2, we discretize the unknown function and its spatial derivatives. The spatial step is chosen as . We denote the value of the unknown function at spatial point and time by with and . By using the approximation from equation (21), we have
Considering , we have the following:
This is the full discrete form of the taken model in case of the kernel . For other kernels and , the discrete form is obtained as follows:
The discretization of initial and boundary conditions is performed by adopting a technique because the value of is not defined outside the domain :
3.3. Theorem (see )
Considering the function which is a member of the space and band limited to , , , , . Then, we have the following required result:where
4. Numerical Results and Validation
In this section, we take some numerical examples of the reaction-diffusion equation and the integro reaction-diffusion equation. We show the validity and feasibility by showing the absolute error tables. The formula of finding errors between exact and numerical solutions is as follows:where denotes the numerical solution at point . We do all numerical computation in Wolfram Mathematica version-11.3.
Example 1. We take the following particular case of our model with parameters , , and :The initial and boundary conditions are taken as follows:We determine the force function by putting the value of the exact solution in the taken particular case of the model. The graph is plotted for the exact and approximate solutions with parameters’ value which can be seen in Figure 1. The variation of absolute errors between exact and numerical solutions is depicted in Table 1.
Example 2. Taking the following integro diffusion equation as a particular case of our model with parametric value , , and :The boundary and initial conditions are taken as follows:The exact solution of the above problem is taken as while force function is chosen with the help of exact solutions and the given problem. The graph is plotted for the numerical and exact solutions with parameters’ value which can be seen in Figure 2. The variation of absolute errors between exact and numerical solutions is depicted in Table 2.
Example 3. If we take , , , and kernel , then we getWith the following initial and boundary conditions:The function is an exact solution of the problem. The graph is plotted for the numerical and exact solutions with parameters’ value which can be seen in Figure 3. The variation of absolute errors between exact and numerical solutions is depicted in Table 3.
Example 4. If we take , , , and , then our model (42) is reduced toThe following equations are chosen as initial and boundary conditions:with the exact solution .
The graph is plotted for the exact and approximate solutions with parameters’ value which can be seen in Figure 4. The variation of absolute errors between exact and numerical solutions is depicted in Table 4. The graphs and error tables clearly show that our numerical results are in feasible and desirable accuracy.
Example 5. Let us take another example of the Burger–Fisher equation:with following initial and boundary conditions:where the exact solution of the above problem is . Now, our obtained numerical results are compared with the results given in .
In this research work, we developed a numerical scheme for time-fractional partial differential equations and integro fractional partial differential equations. The discretization of the time-fractional derivative is done by using the Crank–Nicolson scheme, and the discretization of the unknown function and its spatial derivatives is performed with the help of the quasi-wavelet-based numerical method. For different kernels, we show the validity and effectiveness of our method. Our method has a good accuracy for both fractional PDEs and fractional integro PDEs as seen from the error tables. The same can be seen by the plotted graph for the existing exact solution and numerical solution.
No data were used to support this study.
Conflicts of Interest
The authors declare that there are no conflicts of interest.
This work was supported by the National Natural Science Foundation of China (61673247) and the Research Fund for Distinguished Young Scholars and Excellent Young Scholars of Shandong Province (JQ201719).
A. Kilbas, H. Srivastava, and J. J. Trujillo, Theory and Applications of the Fractional Differential Equations, vol. 204, Elsevier, Amsterdam, Netherlands, 2006.
I. Podlubny, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications, Fractional Differential Equations: An Introduction to Fractional Derivatives, Academic Press, San Diego, CA, USA, 1998.
S. Kumar, P. Pandey, and S. Das, “Operational matrix method for solving nonlinear space-time fractional order reaction-diffusion equation based on genocchi polynomial,” Special Topics & Reviews in Porous Media: An International Journal, vol. 11, no. 1, pp. 33–47, 2020.View at: Publisher Site | Google Scholar