Development and Analysis of Optimized Numerical Scheme for Three-Dimensional Linear Time-Fractional Diffusion Equations

Document Type : Research article

Author

Basic Science Group, Golpayegan College of Engineering, Isfahan University of Technology, Golpayegan, 87717-67498, Iran.

Abstract

In this paper, an optimized numerical method for solving time-fractional diffusion equations in three-dimensional space is proposed. Fractional-order differential equations can model complex physical processes with memory effects more accurately than differential equations with integer order derivatives. On the other hand, numerical solutions of this type of problem are usually challenging due to the existence of a singularity near the initial time. In the proposed method, in order to overcome this issue, the $L1$ formula has been used on a graded mesh to maintain the accuracy of calculations near the initial time. Also, a fourth-order compact operator has been used to discretize the Laplace operator in three-dimensional space. Stability and convergence analyses show that this method has a convergence order proportional to the grading parameter in time and a fourth-order convergence order in space. The obtained results from numerical simulations confirm the high accuracy of the proposed method and the efficiency of the optimal selection of the scaling parameter in achieving the suitable convergence order, even in the presence of initial singularities.

Keywords

Main Subjects

Full Text

 

Article PDF

References

[1] F. Safari, M. A. Zaky, M. Fardi, and Ch. Hong, Fast and accurate high-order discretization of nonlinear fractional integro-differential equations in multi-dimensional settings, International Communications in Heat and Mass Transfer, 175, 111204, (2026).
[2] K. Diethelm and N. J. Ford, The analysis of fractional differential equations, Lecture Notes in Mathematics, 2004, (2010).
[3] B. Azarnavid, M. Fardi, and H. Emami, Efficient polynomial kernel-based method for multi-term time-fractional diffusion systems in regular and irregular domains, Chaos, Solitons & Fractals, 199, 116666, (2025).
[4] M. Fardi, M. A. Zaky, and B. Azarnavid, Fractional convection-diffusion systems in complex 2D and 3D geometries: A Bernoulli polynomial-based kernel method, Applied Numerical Mathematics, (2025).
[5] M. Fardi and B. Azarnavid, A stable numerical method for distributed order time-fractional diffusion equations, International Journal of Numerical Modelling: Electronic Networks, Devices and Fields, 38(2), e70012, (2025).
[6] H. Gündodu and H. Joshi, Numerical analysis of time-fractional cancer models with different types of net killing rate, Mathematics, 13(3), 536, (2025).
[7] M. Fardi, F. Safari, and X. Zhang, Compact radial basis function partition of unity method for 2D neutral reactiondiffusion equations with temporal delays, International Communications in Heat and Mass Transfer, 172, 110330, (2026).
[8] B. Azarnavid, The Bernoulli polynomials reproducing kernel method for nonlinear Volterra integro-differential equations of fractional order with convergence analysis, Computational & Applied Mathematics, 42(1), (2023).
[9] B. Azarnavid, A kernel-based method for fractional integro-differential equations with a weakly singular kernel in multi-dimensional complex domains, Engineering Analysis with Boundary Elements, 159, 1–10, (2024).
[10] M. Stynes, E. O’Riordan, and J. L. Gracia, Error analysis of a finite difference method on graded meshes for a time-fractional diffusion equation, SIAM Journal on Numerical Analysis, 55(2), 1057–1079, (2017).
[11] Z. Z. Sun and G. H. Gao, Finite difference methods for fractional differential equations, Science Press, Beijing, (2015).
[12] I. Podlubny, Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some in Science and Engineering, Academic Press, (1998).
[13] K. B. Oldham and J. Spanier, The Fractional Calculus: Theory and Applications of Differentiation and Integration to Arbitrary Order, Dover Publications, (2006).
[14] R. Hilfer, Applications of Fractional Calculus in Physics, World Scientific, (2000).
[15] M. A. Zaky and M. Fardi, An accurate kernel functions-based method for coupled advection-diffusion-reaction models under Neumann, periodic, and no flux conditions, International Journal of Computer Mathematics, 102(8), 1228–1249, (2025).
[16] F. Razavi Sosan, M. Ghasemi, and M. Fardi, Fourier spectral approximation for the time-fractional reaction-diffusion Brusselator system, Journal of Advanced Mathematical Modeling, 15(2), 87–112, (2025).
[17] S. Irandoust-Pakchin, M. H. Derakhshan, S. Rezapour, and M. Adel, An efficient numerical method for the distributed-order time-fractional diffusion equation with error analysis and stability properties, Mathematical Methods in the Applied Sciences, 48(3), 2743–2765, (2025).
[18] M. Cai and W. Tong, Nonuniform L1/spectral element algorithm for the time fractional diffusion equation, Mathematics and Computers in Simulation, 239, 361–375, (2026).
[19] M. A. Zaky, A linearized two-dimensional Galerkin-L1 spectral method with diagonalization for time-fractional diffusion equations with delay, Applied Numerical Mathematics, (2025).
[20] W. M. Abd-Elhameed, Y. H. Youssri, and A. G. Atta, Adopted spectral tau approach for the time-fractional diffusion equation via seventh-kind Chebyshev polynomials, Boundary Value Problems, 2024(1), 102, (2024).
[21] M. Moustafa, Y. H. Youssri, and A. G. Atta, Explicit ChebyshevGalerkin scheme for the time-fractional diffusion equation, International Journal of Modern Physics C, 35(01), 2450002, (2024).
[22] H. Liao, D. Li, and J. Zhang, Sharp error estimate of the nonuniform L1 formula for linear reaction-subdiffusion equations, SIAM Journal on Numerical Analysis, 56(2), 1112–1133, (2018).
[23] H. Liu, Y. Ma, H. Li, and W. Zhang, Combination of discrete technique on graded meshes with barycentric rational interpolation for solving a class of time-dependent partial integro-differential equations with weakly singular kernels, Computers & Mathematics with Applications, 141, 159–169, (2023).
[24] J. Ma and Y. Jiang, On a graded mesh method for a class of weakly singular Volterra integral equations, Journal of Computational and Applied Mathematics, 231(2), 807–814, (2009).
[25] H. Majidian, Modified Eulers method with a graded mesh for a class of Volterra integral equations with weakly singular kernel, Numerical Algorithms, 67(2), 405–422, (2014).
[26] M. Fardi, S. Mohammadi, and B. Azarnavid, A fast method based on variable time steps for 2D nonlinear time-fractional generalized BenjaminBonaMahonyBurgers equation: error and stability analysis, Journal of Computational Science, 85, 102517, (2025).
[27] H. Brunner, Collocation Methods for Volterra Integral and Related Functional Differential Equations, Cambridge University Press, (2004). 
Analytical and Numerical Solutions for Nonlinear Equations
Volume 11, Issue 2
September 2026
Pages 238-247

Files

History

  • Receive Date: 14 March 2026
  • Revise Date: 23 April 2026
  • Accept Date: 02 June 2026
  • Publish Date: 04 September 2026

Share

How to cite

Statistics