ON APPROXIMATE SOLUTIONS FOR TIME-FRACTIONAL DIFFUSION EQUATION

Abdulkafi Mohammed Saeed

1Department of Mathematics, College of Science, Qassim University, Saudi Arabia

ABSTRACT

In the last decades differential equations involving fractional derivatives and integrals have been studied by many researchers. Due to their ability to model more adequately some phenomena, fractional partial differential equations have been used in numerous areas such as finance, hydrology, porous media, engineering and control systems, etc. Numerical schemes based on rotated finite difference approximation have been proven to work well in solving standard diffusion equations. However, the formulation of these strategies on time fractional diffusion counterpart is still at its infancy. A well-designed preconditioning for these types of problems reduces the number of iterations to reach convergence. In this research work, we have derived new preconditioned fractional rotated finite difference method for solving 2D time-fractional diffusion equation. Numerical experiments are conducted to examine the effectiveness of the proposed method.

Keywords:Preconditioned rotated method Time-fractional diffusion Equation.

ARTICLE HISTORY: Received:23 August 2018Revised:28 September 2018 Accepted:5 November 2018Published:30 November 2018.

Contribution/ Originality:This study contributes in the existing literature about the foundation of fast iterative schemes from the preconditioned methods for solving the time-fractional diffusion equation. It is one of the few studies which combine a suitable pre-conditioner matrix with the rotated iterative scheme as a way to further improve the convergence rate of the method in solving the 2D time-fractional diffusion equation.

1. INTRODUCTION

The importance of this study lies in the various applications of fractional partial differential equations (FPDE's) in finance, physics, image processing and engineering [1 , 2 ].  It is well known that FPDE's is a generalized of the classical partial differential equations (PDE's). As a result of that there is no general method that can be used in solving FPDE's same as classical PDE's. Approximation methods such as finite difference methods have played important role for solving FPDE's in the last few years [3 , 4 ]. It is noteworthy to observe that the finite difference schemes derived from skewed (rotated) difference operators have been extensively investigated over the years for solving FPDE's. These iterative methods have been shown to be much faster than the methods based on the standard five-point formula which is due to the formers’ overall lower computational complexities (Saeed and Ali [5 ]; Ali and Saeed [6 ]; Saeed and Ali [7 ]. In Saeed [8 ] the preconditioned rotated finite difference method applied successfully for solving fractional elliptic partial differential equations and the reveal results was very encouraging.

This work involves an investigation on the utilization of the new preconditioned fractional rotated finite difference method for solving 2D Time-Fractional Diffusion Equations. An outline of this paper is as follows. In Section 2, the proposed accelerated version of fractional rotated five point’s approximation method will be formulated. The numerical results will be presented to show the efficiency of the new proposed methods in Section 3. Finally, Conclusion and Future work are given in Section 4.

2.  FORMULATION OF THE PROPOSED PRECONDITIONED ITERATIVE METHOD

Consider the following time fractional diffusion equation

3. NUMERICAL RESULTS

In this section, we present numerical results for the proposed method applied to two particular examples. The first problem as the following [16]:

4   CONCLUSION AND FUTURE WORK

In this study, we have introduced new preconditioned iterative methods based on fractional rotated finite difference method for solving 2D time-fractional diffusion equation. From observation of all experimental results, it can be conclude that the proposed P2FRFD method requires less time and iterations number when compared to FRFD and P1FRFD methods with same levels of precision. Therefore, the proposed scheme P2FRFD may be a good alternative to solve this type of equations and many other numerical problems.  Numerical results strongly suggest that the efficiency of the proposed preconditioning methods. The convergence analysis of the present iterative method regarding solutions for 2D time-fractional diffusion equation is currently under study. Furthermore, the idea of this proposed method can be extended to group iterative solver which will be reported separately in the future.

Funding: The author gratefully acknowledges Qassim University, represented by the Deanship of Scientific Research, on the material support for this research under the number (1123-cos-2016-12-s) during the academic year 1438 AH/2017AD).
Competing Interests: The author declares that there are no conflicts of interests regarding the publication of this paper.

REFERENCES

[1]          S. I. Muslih and O. P. Agrawal, "Riesz fractional derivatives and fractional dimensional space," International Journal of Theoretical Physics, vol. 49, pp. 270-275, 2010.Available at: https://doi.org/10.1007/s10773-009-0200-1.

[2]          J. Hristov, "Approximate solutions to fractional subdiffusion equations," The European Physical Journal Special Topics, vol. 193, pp. 229-243, 2011.Available at: https://doi.org/10.1140/epjst/e2011-01394-2.

[3]          V. D. Beibalaev and R. P. Meilanov, "The Dirihlet problem for the fractional Poisson’s equation with Caputo derivatives: A finite difference approximation and a numerical solution," Thermal Science, vol. 16, pp. 385-394, 2012.Available at: https://doi.org/10.2298/tsci110421076b.

[4]          Z.-B. Li and J.-H. He, "Fractional complex transform for fractional differential equations," Mathematical and Computational Applications, vol. 15, pp. 970-973, 2010.Available at: https://doi.org/10.3390/mca15050970.

[5]          A. M. Saeed and N. H. M. Ali, "Preconditioned modified explicit decoupled group method in the solution of elliptic PDEs," Applied Mathematical Sciences, vol. 4, pp. 1165-1181, 2010.

[6]          N. H. M. Ali and A. M. Saeed, "Preconditioned modified explicit decoupled group for the solution of steady state navier-stokes equation," Applied Mathematics & Information Sciences, vol. 7, pp. 1837-1844, 2013.Available at: https://doi.org/10.12785/amis/070522.

[7]          A. Saeed and N. H. M. Ali, "On the convergence of the preconditioned group rotated iterative methods in the solution of elliptic PDEs," Applied Mathematics & Information Sciences, vol. 5, pp. 65-73, 2011.

[8]          A. M. Saeed, "Improved rotated finite difference method for solving fractional elliptic partial differential equations," American Scientific Research Journal for Engineering, Technology, and Sciences, vol. 26, pp. 261-270, 2016.

[9]          Y.-n. Zhang and Z.-z. Sun, "Alternating direction implicit schemes for the two-dimensional fractional sub-diffusion equation," Journal of Computational Physics, vol. 230, pp. 8713-8728, 2011.Available at: https://doi.org/10.1016/j.jcp.2011.08.020.

[10]        I. Karatay, N. Kale, and S. Bayramoglu, "A new difference scheme for time fractional heat equations based on the Crank-Nicholson method," Fractional Calculus and Applied Analysis, vol. 16, pp. 892-910, 2013.Available at: https://doi.org/10.2478/s13540-013-0055-2.

[11]        N. H. M. Ali and A. M. Saeed, "Convergence analysis of the preconditioned group splitting methods in boundary value ‎problems," Abstract and Applied Analysis, pp. 1-14, 201‎‏2.Available at: https://doi.org/10.1155/2012/867598.

[12]        A. M. Saeed, "Fast iterative solver for the 2-D convection-diffusion equations," Journal of Advances In Mathematics, vol. 9, pp. 2773-2782, 2014.

[13]        A. M. Saeed, "A numerical study on MEDG iterative method for solving nonlinear steady burgers' equation," Global Journal of EngineeringScience and Research Management, vol. 4, pp. 32-40, 2017.

[14]        A. M. Saeed and A. M. Badiea, "A new approach on numerical solutions of burgers equation using PMEDG iterative method," Journal of Asian Scientific Research, vol. 7, pp. 263-270, 2017.Available at: https://doi.org/10.18488/journal.2.2017.77.263.270.

[15]        A. M. Saeed, "A new version of the accelerated FRFD method for solving time-fractional diffusion equations," Journal of Computer Science and Computational Mathematics, In Press, 2018.

[16]        M. Cui, "Convergence analysis of high-order compact alternating direction implicit schemes for the two-dimensional time fractional diffusion equation," Numerical Algorithms, vol. 62, pp. 383-409, 2013.Available at: https://doi.org/10.1007/s11075-012-9589-3.