Evolution and Mass Decay of a Gaussian Pulse in Time-Fractional Diffusion: An L1-Implicit Scheme Analysis
Pakhlavon Yovkochev
1,*
,
Pierros Ntelis
2,3,+
1 Institute of Fundamental and Applied Research, National Research University TIIAME , Kori Niyoziy 39, Tashkent 100000, Uzbekistan
2 School of Physics, Harbin Institute of Technology , Harbin 150001, People’s Republic of China
3 Institute of Theoretical Physics, National University of Uzbekistan , Tashkent 100174, Uzbekistan
Abstract
We apply the L1-implicit finite difference scheme to solve the time-fractional diffusion equation with a Gaussian initial condition. The Caputo fractional derivative of order α ∈ (0,1] is discretized using the L1 approximation on a uniform temporal mesh, while the spatial second derivative is approximated by a central difference scheme. Numerical simulations are performed to study the effect of the fractional order α on the solution profile and to investigate mass conservation properties. The results demonstrate the characteristic subdiffusive behavior and the non-conservation of mass typically observed in fractional diffusion processes.
References
[1] Ross, B., Hist. Math. 4, 75 (1977)
[2] Podlubny, I., Math. Sci. Eng. (Academic Press, New York, 1999)
[3] Kilbas, A. A., Srivastava, H. M. and Trujillo, J. J., Theory and Applications of Fractional Differential Equations (Elsevier, 2006)
[4] Mainardi, F., Fractional Calculus and Waves in Linear Viscoelasticity (Imperial College Press, 2010)
[5] Samko, S. G., Kilbas, A. A. and Marichev, O. I., Fractional Integrals and Derivatives: Theory and Applications (Gordon and Breach, 1993)
[6] Diethelm, K., The Analysis of Fractional Differential Equations (Springer, 2010)
[7] Li, C. and Zeng, F., Numerical Methods for Fractional Calculus (CRC Press, 2015)
[8] Metzler, R. and Klafter, J., Phys. Rep. 339, 1 (2000)
[9] Metzler, R., Jeon, J.-H., Cherstvy, A. et al., Phys. Chem. Chem. Phys. 16, 24128 (2014)
[10] Caputo, M., Geophys. J. Int. 13, 529 (1967)
[11] Berkowitz, B., Cortis, A., Dentz, M. and Scher, H., Rev. Geophys. 44, RG2003 (2006)
[12] Scalas, E., Gorenflo, R. and Mainardi, F., Physica A 284, 376 (2000)
[13] Lin, Y. and Xu, C., J. Comput. Phys. 225, 1533 (2007)
[14] Meerschaert, M. and Tadjeran, C., J. Comput. Appl. Math. 172, 65 (2004)
[15] Li, X. and Xu, C., J. Comput. Phys. 229, 3983 (2010)
[16] Jiang, Y. and Ma, J., J. Comput. Appl. Math. 235, 3285 (2011)
[17] Mustapha, K., Abdallah, B., Furati, K. M. et al., Numer. Algorithms 73, 517 (2016)
[18] Guo, Y., Song, J. and Song, X., arXiv:2303.11939 [math.PR] (2023)
[19] Guo, Y., Song, J., Wang, R. et al., J. Differ. Equ. 446, 113604 (2025)
[20] Pandey, P., Das, S., Craciun, E. et al., Meccanica 55, 99 (2020)
[21] Lischke, A., Kelly, J. and Meerschaert, M., Fract. Calc. Appl. Anal. 23, 1159 (2020)
[22] Povstenko, Y. and Kyrylych, T., Entropy 24, 1002 (2022)
[23] Hu, J. and Li, G., J. Comput. Phys. 408, 109231 (2020)
[24] Yang, F., Wang, N. and Li, X.-X., J. Appl. Anal. Comput. 10, 514 (2020)
[25] Alodat, T., Le Gia, Q. T. and Sloan, I., J. Comput. Appl. Math. 446, 115863 (2024)
[26] Alodat, T., Le Gia, Q. T. and Sloan, I., arXiv:2412.05817 (2024)
Yovkochev, P. and Ntelis, P., Evolution and Mass Decay of a Gaussian Pulse in Time-Fractional Diffusion: An L1-Implicit Scheme Analysis, Turanian J. Vol. 2, No. 1 (020103), 2026