Abstract
In this paper, an efficient algorithm is presented by adopting the extrapolation technique to improve the accuracy of finite difference schemes for two-dimensional space-fractional diffusion equations with non-smooth solution. The popular fractional centered difference scheme is revisited and the stability and error estimation of numerical solution are given in maximum norm. Based on the analysis of leading singularity of exact solution for the underlying problem, the extrapolation technique and numerical correction method are exploited to enhance the accuracy and convergence rate of the computation. Two numerical examples are provided to validate the theoretical prediction and efficiency of the algorithm. It is shown that, by using the proposed algorithm, both accuracy and convergence rate of numerical solutions can be significantly improved and the second-order accuracy can even be recovered for the equations with large fractional orders.
Similar content being viewed by others
References
Acosta, G., Borthagaray, J. P.: A fractional Laplace equation: regularity of solutions and finite element approximations. SIAM J. Numer. Anal. 55, 472–495 (2017)
Antunes, P. R. S., Ferreira, R. A. C.: An augmented-RBF method for solving fractional Sturm-Liouville eigenvalue problems. SIAM J. Sci. Comput. 37, A515–A535 (2015)
Bu, W., Tang, Y., Yang, J.: Galerkin finite element method for two-dimensional Riesz space fractional diffusion equations. J. Comput. Phys. 276, 26–38 (2014)
Chen, X., Zeng, F., Karniadakis, G. E.: A tunable finite difference method for fractional differential equations with non-smooth solutions. Comput. Methods Appl. Mech. Engrg. 318, 193–214 (2017)
Deng, W., Li, B., Tian, W., Zhang, P.: Boundary problems for the fractional and tempered fractional operators. Multiscale Model Simul. 16, 125–149 (2018)
Ding, H., Li, C., Chen, Y.: High-order algorithms for Riesz derivative and their applications (II). J. Comput. Phys. 293, 218–237 (2015)
Epps, B. P., Cushman-Roisin, B.: Turbulence modeling via the fractional Laplacian. arXiv:1803.05286v1 (2018)
Ervin, V. J., Heuer, N., Roop, J. P.: Regularity of the solution to 1-D fractional order diffusion equations. Math. Comp. 87, 2273–2294 (2018)
Ervin, V. J., Roop, J. P.: Variational formulation for the stationary fractional advection dispersion equation. Numer. Meth. Part. Diff. Eq. 22, 558–576 (2006)
Gatto, P., Hesthaven, J. S.: Numerical approximation of the fractional Laplacian via hp-finite elements, with an application to image denoising. J. Sci. Comput. 65, 249–270 (2015)
Ghanbari, B., Kumar, S., Kumar, R.: A study of behavior for immune and tumor cells in immunogenetic tumor model with non-singular fractional derivative. Chaos Solitons & Fractals. 133, 109619 (2020)
Gunzburger, M., Jiang, N., Xu, F.: Analysis and approximation of a fractional Laplacian-based closure model for turbulent flows and its connection to Richardson pair dispersion. Comput. Math. Appl. 75, 1973–2001 (2018)
Hao, Z., Cao, W.: An improved algorithm based on finite difference schemes for fractional boundary value problems with nonsmooth solution. J. Sci. Comput. 73, 395–415 (2017)
Hao, Z., Lin, G., Zhang, Z.: Error estimates of a spectral Petrov-Galerkin method for two-sided fractional reaction-diffusion equations. Appl. Math. Comput. 374, 125045 (2020)
Hao, Z., Zhang, Z.: Optimal regularity and error estimates of a spectral Galerkin method for fractional advection-diffusion-reaction equations. SIAM J. Numer. Anal. 58, 211–233 (2020)
Hao, Z., Sun, Z. -Z., Cao, W.: A fourth-order approximation of fractional derivatives with its applications. J. Comput. Phys. 281, 787–805 (2015)
Hao, Z., Zhang, Z., Du, R.: Fractional centered difference scheme for high-dimensional integral fractional Laplacian. Researchgate. https://www.researchgate.net/publication/335888811 (2019)
Jin, B., Zhou, Z.: A finite element method with singularity reconstruction for fractional boundary value problems. ESAIM Math. Model. Numer. Anal. 49, 1261–1283 (2015)
Kopteva, N., Stynes, M.: An efficient collocation method for a Caputo two-point boundary value problem. BIT. 55, 1105–1123 (2015)
Kumar, S., Kumar, R., Singh, J., Nisar, K. S., Kumar, D.: An efficient numerical scheme for fractional model of HIV-1 infection of CD4+ T-cells with the effect of antiviral drug therapy. Alexandria Engineering Journal. https://doi.org/10.1016/j.aej.2019.12.046 (2019)
Laskin, N.: Fractional quantum mechanics and lévy path integrals. Phys. Lett. A. 268, 298–305 (2000)
Liu, F., Anh, V., Turner, I.: Numerical solution of the space fractional Fokker-Planck equation. Proceedings of the International Conference on Boundary and Interior Layers—Computational and Asymptotic Methods (BAIL 2002) 166, 209–219 (2004)
Liu, F., Zhuang, P., Turner, I., Anh, V., Burrage, K.: A semi-alternating direction method for a 2-D fractional FitzHugh-Nagumo monodomain model on an approximate irregular domain. J. Comput. Phys. 293, 252–263 (2015)
Magin, R., Abdullah, O., Baleanu, D., Zhou, X. J.: Anomalous diffusion expressed through fractional order differential operators in the Bloch–Torrey equation. J Magn Reson. 190, 255–270 (2008)
Mao, Z., Chen, S., Shen, J.: Efficient and accurate spectral method using generalized Jacobi functions for solving Riesz fractional differential equations. Appl. Numer. Math. 106, 165–181 (2016)
Song, F., Xu, C.: Spectral direction splitting methods for two-dimensional space fractional diffusion equations. J. Comput. Phys. 299, 196–214 (2015)
Stynes, M.: Singularities. In: Karniadakis, G. E. (ed.) Handbook of Fractional Calculus with Applications, vol. 3, pp 287–305. De Gruyter, Berlin (2019)
Sun, H., Sun, Z. Z., Gao, G. H.: Some high order difference schemes for the space and time fractional Bloch-Torrey equations. Appl. Math. Comput. 281, 356–380 (2016)
Tadjeran, C., Meerschaert, M. M., Scheffler, H. -P.: A second-order accurate numerical approximation for the fractional diffusion equation. J. Comput. Phys. 213, 205–213 (2006)
Tian, W., Zhou, H., Deng, W.: A class of second order difference approximations for solving space fractional diffusion equations. Math. Comp. 84, 1703–1727 (2015)
Wang, H., Wang, K., Sircar, T.: A direct \(o(n\log ^{2N)}\) finite difference method for fractional diffusion equations. J. Comput. Phys. 229, 8095–8104 (2010)
Woyczyński, W. A.: Lévy Processes in the Physical Sciences. In: Barndorff-Nielsen, O.E., Resnick, S. I., Mikosch, T. (eds.) Processes, Lévy, pp 241–266. Birkhäuser, Boston (2001)
Xu, C.: Spectral methods for some kinds of fractional differential equations: traditional and Müntz spectral methods. In: Karniadakis, G. E. (ed.) Handbook of Fractional Calculus with Applications, vol. 3, pp 101–126. De Gruyter, Berlin (2019)
Zayernouri, M., Karniadakis, G. E.: Fractional Sturm-Liouville eigen-problems: theory and numerical approximation. J. Comput. Phys. 252, 495–517 (2013)
Zhao, L., Deng, W.: High order finite difference methods on non-uniform meshes for space fractional operators. Adv. Comput. Math. 42, 425–468 (2016)
Zhao, X., Sun, Z. Z., Hao, Z.: A fourth-order compact ADI scheme for two-dimensional nonlinear space fractional Schrödinger equation. SIAM J. Sci. Comput. 36, A2865–A2886 (2014)
Acknowledgments
The authors would like to thank Prof. Zhi-Zhong Sun and Prof. Rui Du for helpful discussion and suggestion.
Funding
This work was partially supported by the NSF of China (No. 11671083) and the Fundamental Research Funds for the Central Universities.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
About this article
Cite this article
Hao, Z., Cao, W. & Li, S. Numerical correction of finite difference solution for two-dimensional space-fractional diffusion equations with boundary singularity. Numer Algor 86, 1071–1087 (2021). https://doi.org/10.1007/s11075-020-00923-8
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11075-020-00923-8
Keywords
- Non-smooth solution
- Riesz derivative
- Extrapolation technique
- Error estimate in maximum norm
- Convergence rate