This work was devoted to the construction of a numerical algorithm for solving the initial boundary value problem for the subdiffusion equation with nonlocal boundary conditions. For the case of not strongly regular boundary conditions, the well-known methods cannot be used. We applied an algorithm that consists of reducing the nonlocal problem to a sequential solution of two subproblems with local boundary conditions. The solution to the original problem was summed up from the solutions of the subproblems. To solve the subproblems, we constructed implicit difference schemes on the basis of the L1 formula for approximating the Caputo fractional derivative and central difference for approximating the space derivatives. Stability and convergence of the schemes were established. The Thomas algorithm was used to solve systems of linear algebraic equations. Numerical experiments were conducted to study the constructed algorithm. In terms of accuracy and stability, the algorithm performs well. The results of experiments confirmed that the convergence order of the method coincides with the theoretical one, $ O(\tau^{2-\alpha}+h^2) $.
Citation: Murat A. Sultanov, Vladimir E. Misilov, Makhmud A. Sadybekov. Numerical method for solving the subdiffusion differential equation with nonlocal boundary conditions[J]. AIMS Mathematics, 2024, 9(12): 36385-36404. doi: 10.3934/math.20241726
This work was devoted to the construction of a numerical algorithm for solving the initial boundary value problem for the subdiffusion equation with nonlocal boundary conditions. For the case of not strongly regular boundary conditions, the well-known methods cannot be used. We applied an algorithm that consists of reducing the nonlocal problem to a sequential solution of two subproblems with local boundary conditions. The solution to the original problem was summed up from the solutions of the subproblems. To solve the subproblems, we constructed implicit difference schemes on the basis of the L1 formula for approximating the Caputo fractional derivative and central difference for approximating the space derivatives. Stability and convergence of the schemes were established. The Thomas algorithm was used to solve systems of linear algebraic equations. Numerical experiments were conducted to study the constructed algorithm. In terms of accuracy and stability, the algorithm performs well. The results of experiments confirmed that the convergence order of the method coincides with the theoretical one, $ O(\tau^{2-\alpha}+h^2) $.
[1] | H. Sun, Y. Zhang, D. Baleanu, W. Chen, Y. Chen, A new collection of real world applications of fractional calculus in science and engineering, Commun. Nonlinear Sci., 64 (2018), 213–231. https://doi.org/10.1016/j.cnsns.2018.04.019 doi: 10.1016/j.cnsns.2018.04.019 |
[2] | A. A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and applications of fractional differential equations, Amsterdam: Elsevier, 2006. |
[3] | J. T. Machado, V. Kiryakova, F. Mainardi, Recent history of fractional calculus, Commun. Nonlinear Sci., 16 (2011), 1140–1153. https://doi.org/10.1016/j.cnsns.2010.05.027 doi: 10.1016/j.cnsns.2010.05.027 |
[4] | R. Metzler, J. H. Jeon, A. G. Cherstvy, E. Barkai, Anomalous diffusion models and their properties: non-stationarity, non-ergodicity, and ageing at the centenary of single particle tracking, Phys. Chem. Chem. Phys., 16 (2014), 24128–24164. https://doi.org/10.1039/C4CP03465A doi: 10.1039/C4CP03465A |
[5] | M. J. Buckingham, On pore-fluid viscosity and the wave properties of saturated granular materials including marine sediments, J. Acoust. Soc. Am., 122 (2007), 1486–1501. https://doi.org/10.1121/1.2759167 doi: 10.1121/1.2759167 |
[6] | W. Chen, S. Hu, W. Cai, A causal fractional derivative model for acoustic wave propagation in lossy media, Arch. Appl. Mech., 86 (2016), 529–539. https://doi.org/10.1007/s00419-015-1043-2 doi: 10.1007/s00419-015-1043-2 |
[7] | M. Du, Z. Wang, H. Hu, Measuring memory with the order of fractional derivative, Sci. Rep., 3 (2013), 3431. https://doi.org/10.1038/srep03431 doi: 10.1038/srep03431 |
[8] | J. Zhang, Z. Wei, L. Xiao, Adaptive fractional-order multi-scale method for image denoising, J. Math. Imaging Vis., 43 (2012), 39–49. https://doi.org/10.1007/s10851-011-0285-z doi: 10.1007/s10851-011-0285-z |
[9] | Q. Yang, D. Chen, T. Zhao, Y. Chen, Fractional calculus in image processing: a review, Fract. Calc. Appl. Anal., 19 (2016), 1222–1249. https://doi.org/10.1515/fca-2016-0063 doi: 10.1515/fca-2016-0063 |
[10] | R. Magin, Fractional calculus in bioengineering, part 1, Crit. Rev. Biomed. Eng., 32 (2004), 1–104. https://doi.org/10.1615/CritRevBiomedEng.v32.i1.10 doi: 10.1615/CritRevBiomedEng.v32.i1.10 |
[11] | C. Pinto, A. Carvalho, Fractional complex-order model for HIV infection with drug resistance during therapy, J. Vib. Control, 22 (2016), 2222–2239. https://doi.org/10.1177/1077546315574964 doi: 10.1177/1077546315574964 |
[12] | Y. Zhang, B. Baeumer, L. Chen, D. Reeves, H. Sun, A fully subordinated linear flow model for hillslope subsurface stormflow, Water Resour. Res., 53 (2017), 3491–3504. https://doi.org/10.1002/2016WR020192 doi: 10.1002/2016WR020192 |
[13] | V. E. Tarasov, V. V. Tarasova, Long and short memory in economics: fractional-order difference and differentiation, IRA-International Journal of Management and Social Sciences, 5 (2016), 327–334. https://dx.doi.org/10.21013/jmss.v5.n2.p10 doi: 10.21013/jmss.v5.n2.p10 |
[14] | A. V. Bitsadze, A. A. Samarskii, Some elementary generalizations of linear elliptic boundary value problems, Dokl. Akad. Nauk SSSR, 185 (1969), 739–740. |
[15] | V. A. Il'in, Existence of a reduced system of eigen- and associated functions for a nonselfadjoint ordinary differential operator, Trudy Mat. Inst. Steklov., 142 (1976), 148–155. |
[16] | N. I. Ionkin, The solution of a certain boundary value problem of the theory of heat conduction with a nonclassical boundary condition, Differ. Uravn., 13 (1977), 294–304. |
[17] | N. I. Ionkin, E. I. Moiseev, A problem for a heat equation with two-point boundary conditions, Differ. Uravn., 15 (1979), 1284–1295. |
[18] | V. P. Mikhailov, On Riesz bases in ${L}_2(0, 1)$, Dokl. Akad. Nauk SSSR, 144 (1962), 981–984. |
[19] | M. A. Naimark, Linear differential operators, London: Harrap, 1967. |
[20] | N. Dunford, J. Schwartz, Linear operators, part III: spectral operators, New York: John Wiley & Sons Inc, 1957. |
[21] | A. Mokin, On a family of initial-boundary value problems for the heat equation, Diff. Equat., 45 (2009), 126–141. https://doi.org/10.1134/S0012266109010133 doi: 10.1134/S0012266109010133 |
[22] | M. A. Sadybekov, Initial-boundary value problem for a heat equation with not strongly regular boundary conditions, In: Functional analysis in interdisciplinary applications, Cham: Springer, 2017,330–348. https://doi.org/10.1007/978-3-319-67053-9_32 |
[23] | A. S. Makin, On spectral decompositions corresponding to non-self-adjoint Sturm-Liouville operators, Dokl. Math., 73 (2006), 15–18. https://doi.org/10.1134/S1064562406010042 doi: 10.1134/S1064562406010042 |
[24] | A. V. Gulin, On the spectral stability in subspaces for difference schemes with nonlocal boundary conditions, Diff. Equat., 49 (2013), 815–823. https://doi.org/10.1134/S0012266113070045 doi: 10.1134/S0012266113070045 |
[25] | R. Ashurov, Y. Fayziev, On the nonlocal problems in time for time-fractional subdiffusion equations, Fractal Fract., 6 (2022), 41. https://doi.org/10.3390/fractalfract6010041 doi: 10.3390/fractalfract6010041 |
[26] | E. Ozbilge, F. Kanca, E. Özbilge, Inverse problem for a time fractional parabolic equation with nonlocal boundary conditions, Mathematics, 10 (2022), 1479. https://doi.org/10.3390/math10091479 doi: 10.3390/math10091479 |
[27] | M. A. Almalahi, M. S. Abdo, S. K. Panchal, Periodic boundary value problems for fractional implicit differential equations involving Hilfer fractional derivative, Probl. Anal. Issues Anal., 9 (2020), 16–44. https://doi.org/10.15393/j3.art.2020.7410 doi: 10.15393/j3.art.2020.7410 |
[28] | M. A. Almalahi, A. B. Ibrahim, A. Almutairi, O. Bazighifan, T. A. Aljaaidi, J. A. Awrejcewicz, A qualitative study on second-order nonlinear fractional differential evolution equations with generalized ABC operator, Symmetry, 14 (2022), 207. https://doi.org/10.3390/sym14020207 doi: 10.3390/sym14020207 |
[29] | A. A. Alikhanov, Stability and convergence of difference schemes approximating a two-parameter nonlocal boundary value problem for time-fractional diffusion equation, Comput. Math. Model., 26 (2015), 252–272. https://doi.org/10.1007/s10598-015-9271-4 doi: 10.1007/s10598-015-9271-4 |
[30] | N. B. Kerimov, M. I. Ismailov, Direct and inverse problems for the heat equation with a dynamic-type boundary condition, IMA J. Appl. Math., 80 (2015), 1519–1533. https://doi.org/10.1093/imamat/hxv005 doi: 10.1093/imamat/hxv005 |
[31] | N. H. Tuan, N. A. Triet, N. H. Luc, N. D. Phuong, On a time fractional diffusion with nonlocal in time conditions, Adv. Differ. Equ., 2021 (2021), 204. https://doi.org/10.1186/s13662-021-03365-1 doi: 10.1186/s13662-021-03365-1 |
[32] | Z. Sun, G. Gao, Fractional differential equations: finite difference methods, Berlin: De Gruyter, 2020. https://doi.org/10.1515/9783110616064 |
[33] | M. A. Sadybekov, I. N. Pankratova, Correct and stable algorithm for numerical solving nonlocal heat conduction problems with not strongly regular boundary conditions, Mathematics, 10 (2022), 3780. https://doi.org/10.3390/math10203780 doi: 10.3390/math10203780 |
[34] | M. A. Sultanov, E. N. Akimova, V. E. Misilov, Y. Nurlanuly, Parallel direct and iterative methods for solving the time-fractional diffusion equation on multicore processors, Mathematics, 10 (2022), 323. https://doi.org/10.3390/math10030323 doi: 10.3390/math10030323 |
[35] | E. N. Akimova, M. A. Sultanov, V. E. Misilov, Y. Nurlanuly, Parallel algorithm for solving the inverse two-dimensional fractional diffusion problem of identifying the source term, Fractal Fract., 7 (2023), 801. https://doi.org/10.3390/fractalfract7110801 doi: 10.3390/fractalfract7110801 |
[36] | A. A. Alikhanov, A time-fractional diffusion equation with generalized memory kernel in differential and difference settings with smooth solutions, Comput. Meth. Appl. Math., 17 (2017), 647–660. https://doi.org/10.1515/cmam-2017-0035 doi: 10.1515/cmam-2017-0035 |
[37] | M. A. Sultanov, M. I. Akylbaev, R. Ibragimov, Conditional stability of a solution of a difference scheme for an ill-posed Cauchy problem, Electron. J. Differ. Eq., 2018 (2018), 33. |
[38] | K. Cao, D. Lesnic, M. I. Ismailov, Determination of the time-dependent thermal grooving coefficient, J. Appl. Math. Comput., 65 (2021), 199–221. https://doi.org/10.1007/s12190-020-01388-7 doi: 10.1007/s12190-020-01388-7 |
[39] | N. H. Luc, D. Baleanu, R. P. Agarwal, L. D. Long, Identifying the source function for time fractional diffusion with non-local in time conditions, Comp. Appl. Math., 40 (2021), 159. https://doi.org/10.1007/10.1007/s40314-021-01538-y doi: 10.1007/10.1007/s40314-021-01538-y |
[40] | M. Slodička, Uniqueness for an inverse source problem of determining a space-dependent source in a non-autonomous time-fractional diffusion equation, Fract. Calc. Appl. Anal., 23 (2020), 1702–1711. https://doi.org/10.1515/fca-2020-0084 doi: 10.1515/fca-2020-0084 |
[41] | E. N. Akimova, V. E. Misilov, M. A. Sultanov, Regularized gradient algorithms for solving the nonlinear gravimetry problem for the multilayered medium, Math. Method. Appl. Sci., 45 (2022), 8760–8768. https://doi.org/10.1002/mma.7012 doi: 10.1002/mma.7012 |
[42] | Y. Zhang, A finite difference method for fractional partial differential equation, Appl. Math. Comput., 215 (2009), 524–529. https://doi.org/10.1016/j.amc.2009.05.018 doi: 10.1016/j.amc.2009.05.018 |
[43] | I. Podlubny, Fractional differential equations, London: Academic Press, 1999. |
[44] | A. Samarskii, E. Nikolaev, Numerical methods for grid equations, volume I: direct methods, Basel: Birkhäuser, 1989. https://doi.org/10.1007/978-3-0348-9272-8 |
[45] | Y. Dimitrov, R. Miryanov, V. Todorov, Asymptotic expansions and approximations for the Caputo derivative, Comp. Appl. Math., 37 (2018), 5476–5499. https://doi.org/10.1007/s40314-018-0641-3 doi: 10.1007/s40314-018-0641-3 |
[46] | G. Gao, Z. Sun, H. Zhang, A new fractional numerical differentiation formula to approximate the Caputo fractional derivative and its applications, J. Comput. Phys., 259 (2014), 33–50. https://doi.org/10.1016/j.jcp.2013.11.017 doi: 10.1016/j.jcp.2013.11.017 |
[47] | A. A. Alikhanov, A new difference scheme for the time fractional diffusion equation, J. Comput. Phys., 280 (2015), 424–438. https://doi.org/10.1016/j.jcp.2014.09.031 doi: 10.1016/j.jcp.2014.09.031 |