Control and Stability Analysis of Complex Dynamical Systems involved in Finance, Ecology and Engineering
View this Special IssueResearch Article  Open Access
Sachin Kumar, Jinde Cao, Xiaodi Li, "A Numerical Method for TimeFractional ReactionDiffusion and Integro ReactionDiffusion Equation Based on QuasiWavelet", Complexity, vol. 2020, Article ID 3291723, 11 pages, 2020. https://doi.org/10.1155/2020/3291723
A Numerical Method for TimeFractional ReactionDiffusion and Integro ReactionDiffusion Equation Based on QuasiWavelet
Abstract
In this research work, we focused on finding the numerical solution of timefractional reactiondiffusion and another class of integrodifferential equation known as the integro reactiondiffusion equation. For this, we developed a numerical scheme with the help of quasiwavelets. 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 reactiondiffusion are discretized and approximated with the help of quasiwavelets. We study this model with Dirichlet boundary conditions. The discretization of these initial and boundary conditions is done with a different approach by the quasiwaveletbased 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.
1. Introduction
Differential and integral calculus is a vital and broad branch of mathematics [1, 2]. First of all, the integerorder 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 variableorder 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 integerorder 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 [5], predictorcorrector method [6], Adomian decomposition method [7], homotopy perturbation method [8], and generalized block pulse operational matrix method [9]. 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 [9], Haar wavelets [10], Chebyshev wavelets [11], and Legendre wavelets [12]. The stability of discretetime delayed switched systems and time delayed systems is given in the literature [13, 14].
FDEs, FPDEs, and fractional integrodifferential 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 integrodifferential 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 reactiondiffusion process. In this process, molecules react with the fluids. This phenomenon of diffusion can be depicted by the following PDE known as the reactiondiffusion partial differential equation:where represents the reaction term and the first term on the righthand 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 quasiwavelet to find out an approximation of an unknown function. The taken model of integro reactiondiffusion and reactiondiffusion is given in Section 3. How to deal with discretization of spatial and timefractional 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 QuasiWavelets
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 RLtype definition of fractional integration is given as follows [18]:
The definition of fractional differentiation in the RL 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 integerorder 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 RL definition and Caputo is given as follows:
2.3. Approximation of Functions with the Help of QuasiWavelets
The method based upon quasiwavelets nowadays is developing in the form of the local spectral collocation method. It is used for seeking the numerical solution of the fractionalorder 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 quasiwavelets. It is defined aswhere and is a wellknown Dirac delta function. Shannon’s delta sequence kernel produces the basis for the Paley–Wiener reproducing kernel Hilbert space [19] 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 bandlimited 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 Gaussianregularized sampling scaling function follows the property of the quasiscaling function. So, we can approximate an unknown function belonging to the space with the help of the quasiscaling 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 integerorder 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 [20] for the calculation purpose:
3. Mathematical Model and Used Proposed Method
In this section, we propose the integrated mathematical model of the reactiondiffusion and integro reactiondiffusion equation. In this section, we develop the Crank–Nicolson scheme with the combination of the quasiwaveletbased 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
Similarly,
The temporal semidiscrete form is obtained with the help of equations (32)–(34):
3.2. Use of QuasiWaveletBased Approximation for the Discretization of Spatial Derivatives
By using the approximation by quasiwavelets 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
With the help of equation (33), we can rewrite equation (34) as follows:
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 [21])
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 reactiondiffusion equation and the integro reactiondiffusion 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 version11.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 [22].
5. Conclusion
In this research work, we developed a numerical scheme for timefractional partial differential equations and integro fractional partial differential equations. The discretization of the timefractional 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 quasiwaveletbased 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.
Data Availability
No data were used to support this study.
Conflicts of Interest
The authors declare that there are no conflicts of interest.
Acknowledgments
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).
References
 X. Lv, R. Rakkiyappan, R. Rakkiyappan, and X. Li, “mstability criteria for nonlinear differential systems with additive leakage and transmission timevarying delays,” Nonlinear Analysis: Modelling and Control, vol. 23, no. 3, pp. 380–400, 2018. View at: Publisher Site  Google Scholar
 A. Vinodkumar, T. Senthilkumar, and X. Li, “Robust exponential stability results for uncertain infinite delay differential systems with random impulsive moments,” Advances in Difference Equations, no. 1, p. 39, 2018. View at: Publisher Site  Google Scholar
 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.
 P. Darania and A. Ebadian, “A method for the numerical solution of the integrodifferential equations,” Applied Mathematics and Computation, vol. 188, no. 1, pp. 657–668, 2007. View at: Publisher Site  Google Scholar
 K. Diethelm, N. J. Ford, and A. D. Freed, “A predictorcorrector approach for the numerical solution of fractional differential equations,” Nonlinear Dynamics, vol. 29, no. 1–4, pp. 3–22, 2002. View at: Publisher Site  Google Scholar
 L. Suarez and A. Shokooh, “An eigenvector expansion method for the solution of motion containing fractional derivatives,” Journal of Applied Mechanics, vol. 64, p. 9, 1997. View at: Publisher Site  Google Scholar
 S. Kumar and D. Baleanu, “Numerical solution of twodimensional time fractional cable equation with mittagleffler kernel,” Mathematical Methods in the Applied Sciences, vol. 43, no. 15, pp. 8348–8362, 2020. View at: Publisher Site  Google Scholar
 S. Kumar, P. Pandey, and S. Das, “Gegenbauer wavelet operational matrix method for solving variableorder nonlinear reaction–diffusion and galilei invariant advection–diffusion equations,” Computational and Applied Mathematics, vol. 38, no. 4, p. 162, 2019. View at: Publisher Site  Google Scholar
 Y. Li and W. Zhao, “Haar wavelet operational matrix of fractional order integration and its applications in solving the fractional order differential equations,” Applied Mathematics and Computation, vol. 216, no. 8, pp. 2276–2285, 2010. View at: Publisher Site  Google Scholar
 L. Yuanlu, “Solving a nonlinear fractional differential equation using Chebyshev wavelets,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 9, pp. 2284–2292, 2010. View at: Publisher Site  Google Scholar
 H. Jafari, S. A. Yousefi, M. A. Firoozjaee, S. Momani, and C. M. Khalique, “Application of legendre wavelets for solving fractional differential equations,” Computers & Mathematics with Applications, vol. 62, no. 3, pp. 1038–1045, 2011. View at: Publisher Site  Google Scholar
 B. Liu and D. J. Hill, “Stability via hybrideventtime Lyapunov function and impulsive stabilization for discretetime delayed switched systems,” SIAM Journal on Control and Optimization, vol. 52, no. 2, pp. 1338–1365, 2014. View at: Publisher Site  Google Scholar
 B. Liu, D. J. Hill, Z. Sun, and J. Huang, “Eventtriggered control via impulses for exponential stabilization of discrete‐time delayed systems and networks,” International Journal of Robust and Nonlinear Control, vol. 29, no. 6, pp. 1613–1638, 2019. View at: Publisher Site  Google Scholar
 X. Li, X. Yang, and T. Huang, “Persistence of delayed cooperative models: impulsive control method,” Applied Mathematics and Computation, vol. 342, pp. 130–146, 2019. View at: Publisher Site  Google Scholar
 D. Yang, X. Li, J. Shen, and Z. Zhou, “Statedependent switching control of delayed switched systems with stable and unstable modes,” Mathematical Methods in the Applied Sciences, vol. 41, no. 16, pp. 6968–6983, 2018. View at: Publisher Site  Google Scholar
 X. Zhang, X. Li, and X. Han, “Design of hybrid controller for synchronization control of chen chaotic system,” The Journal of Nonlinear Sciences and Applications, vol. 10, no. 6, pp. 3320–3327, 2017. View at: Publisher Site  Google Scholar
 D. Tavares, R. Almeida, and D. F. M. Torres, “Caputo derivatives of fractional variable order: numerical approximations,” Communications in Nonlinear Science and Numerical Simulation, vol. 35, pp. 69–87, 2016. View at: Publisher Site  Google Scholar
 G. W. Wei, “Discrete singular convolution for the solution of the FokkerPlanck equation,” The Journal of Chemical Physics, vol. 110, no. 18, pp. 8930–8942, 1999. View at: Publisher Site  Google Scholar
 X. Yang, D. Xu, and H. Zhang, “Quasiwavelet based numerical method for fourthorder partial integrodifferential equations with a weakly singular kernel,” International Journal of Computer Mathematics, vol. 88, no. 15, pp. 3236–3254, 2011. View at: Publisher Site  Google Scholar
 L. Qian, “On the regularized WhittakerKotel’nikovShannon sampling formula,” Proceedings of the American Mathematical Society, vol. 131, no. 4, pp. 1169–1176, 2003. View at: Publisher Site  Google Scholar
 S. Kumar, P. Pandey, and S. Das, “Operational matrix method for solving nonlinear spacetime fractional order reactiondiffusion 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
Copyright
Copyright © 2020 Sachin Kumar et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.