In this research article, we develop a powerful algorithm for numerical solutions to variable-order partial differential equations (PDEs). For the said method, we utilize properties of shifted Legendre polynomials to establish some operational matrices of variable-order differentiation and integration. With the help of the aforementioned operational matrices, we reduce the considered problem to a matrix type equation (equations). The resultant matrix equation is then solved by using computational software like Matlab to get the required numerical solution. Here it should be kept in mind that the proposed algorithm omits discretization and collocation which save much of time and memory. Further the numerical scheme based on operational matrices is one of the important procedure of spectral methods. The mentioned scheme is increasingly used for numerical analysis of various problems of differential as well as integral equations in previous many years. Pertinent examples are given to demonstrate the validity and efficiency of the method. Also some error analysis and comparison with traditional Haar wavelet collocations (HWCs) method is also provided to check the accuracy of the proposed scheme.
Citation: Kamal Shah, Hafsa Naz, Muhammad Sarwar, Thabet Abdeljawad. On spectral numerical method for variable-order partial differential equations[J]. AIMS Mathematics, 2022, 7(6): 10422-10438. doi: 10.3934/math.2022581
In this research article, we develop a powerful algorithm for numerical solutions to variable-order partial differential equations (PDEs). For the said method, we utilize properties of shifted Legendre polynomials to establish some operational matrices of variable-order differentiation and integration. With the help of the aforementioned operational matrices, we reduce the considered problem to a matrix type equation (equations). The resultant matrix equation is then solved by using computational software like Matlab to get the required numerical solution. Here it should be kept in mind that the proposed algorithm omits discretization and collocation which save much of time and memory. Further the numerical scheme based on operational matrices is one of the important procedure of spectral methods. The mentioned scheme is increasingly used for numerical analysis of various problems of differential as well as integral equations in previous many years. Pertinent examples are given to demonstrate the validity and efficiency of the method. Also some error analysis and comparison with traditional Haar wavelet collocations (HWCs) method is also provided to check the accuracy of the proposed scheme.
[1] | Y. A. Rossikhin, M. V. Shitikova, Application of fractional calculus for dynamic problems of solid mechanics: novel trends and recent results, ASME. Appl. Mech. Rev., 63 (2010), 010801. https://doi.org/10.1115/1.4000563 doi: 10.1115/1.4000563 |
[2] | R. Hilfer, Applications of fractional calculus in physics, World Scientific, Singapor, 2000. https://doi.org/10.1142/3779 |
[3] | V. E. Tarasov, On history of mathematical economics: Application of fractional calculus, Mathematics, 7 (2019), 509. https://doi.org/10.3390/math7060509 doi: 10.3390/math7060509 |
[4] | J. V. D. C. Sousa, E.C. de Oliveira, On two new operators in fractional calculus and application, arXiv preprint arXiv, 2017. |
[5] | B. J. West, Fractal physiology and the fractional calculus: a perspective, Front. physiol, 1 (2010), 12. |
[6] | M. P. Aghababa, Chaos in a fractional-order micro-electro-mechanical resonator and its suppression, Chinese Phys. B, 21 (2012), 100505. https://doi.org/10.1088/1674-1056/21/10/100505 doi: 10.1088/1674-1056/21/10/100505 |
[7] | A. Carpinteri, F. Mainardi, Fractals and Fractional Calculus in Continuum Mechanics, Springer, Verlage, 2014. |
[8] | R. Gorenflo, F. Mainardi, E. Scalas, M. Raberto, Fractional calculus and continuous time finance III, The diffusion limit, Mathe. Finance, Birkhäuser, Basel, 2001 (2001), 171–180. |
[9] | F. Liu, V. Anh, I. Turner, Numerical solution of the space fractional Fokker–Planck equation, J. Comp. Appl. Math., 66 (2005), 209–219. |
[10] | R. Metzler, J. Klafter, The random walks guide to anomalous diffusion: A fractional dynamics approach, Phys. Reports, 339 (2000), 1–77. https://doi.org/10.1016/S0370-1573(00)00070-3 doi: 10.1016/S0370-1573(00)00070-3 |
[11] | R. Metzler, J. Klafter, Boundary value problems for fractional diffusion equations, Phys. A: Stat. Mech. Appl., 278 (2005), 107–125. |
[12] | A. Saadatmandi, M. Deghan, Numerical solution of the a mathematical formation in tumer angionesis via tau method, Commun. Numer. Methods Eng., 24 (2008), 1467–1474. https://doi.org/10.1002/cnm.1045 doi: 10.1002/cnm.1045 |
[13] | R.A. Khan, M. Rehman, Existence of multiple positive solutions for a general system of fractional differential equations, Commun. Appl. Nonlinear Anal., 18 (2011), 25–35. |
[14] | K. I. Isife, Positive solutions of a class of nonlinear boundary value fractional differential equations, J. Frac. Calc. Nonl. Sys., 2 (2021), 12–30. |
[15] | M. Rehman, R.A. Khan, Existence and uniqueness of solutions for multi-point boundary value problems for fractional differential equations, Appl. Math. Lett. 23 (2010), 1038–1044. https://doi.org/10.1016/j.aml.2010.04.033 |
[16] | E. Ziada, Numerical solution for multi-term fractional delay differential equations, J. Frac. Calc. Nonl. Sys., 2 (2021), 1–12. |
[17] | M. Rehman, R. A. Khan, J. Henderson, Existence and uniqueness of solutions for nonlinear fractional differential equations with integral boundary conditions, Fract. Dyn. Syst., 1 (2011), 29–43. https://doi.org/10.7153/fdc-01-02 doi: 10.7153/fdc-01-02 |
[18] | S. G. Samko, B. Ross, Integration and differentiation to a variable fractional order, Integral Transforms Spec. Funct., 1 (1993), 277–300. |
[19] | H. Dehestani, Y. Ordokhani, M. Razzaghi, Numerical solution of variable-order time fractional weakly singular partial integro-differential equations with error estimation, Mathematical Modelling and Analysis, 25 (2020), 680–701. https://doi.org/10.3846/mma.2020.11692 doi: 10.3846/mma.2020.11692 |
[20] | H. Dehestani, Y. Ordokhani, M. Razzaghi, A novel direct method based on the Lucas multiwavelet functions for variable-order fractional reaction-diffusion and subdiffusion equations, Numer. Linear Algebra Appl., 28 (2021), e2346. https://doi.org/10.1002/nla.2346 doi: 10.1002/nla.2346 |
[21] | H. Dehestani, Y. Ordokhani, M. Razzaghi, Modified wavelet method for solving multitype variable-order fractional partial differential equations generated from the modeling of phenomena, Math. Sci., (2021), 1–17. |
[22] | J. Jiang, J. L. Guirao, T. Saeed, The existence of the extremal solution for the boundary value problems of variable fractional order differential equation with causal operator, Fractals, 28 (2020), 2040025. https://doi.org/10.1142/S0218348X20400253 doi: 10.1142/S0218348X20400253 |
[23] | S. Sarwar, On the existence and stability of variable order caputo type fractional differential equations, Fractal Fract., 6 (2022), 51. https://doi.org/10.3390/fractalfract6020051 doi: 10.3390/fractalfract6020051 |
[24] | C. F. M. Coimbra, Mechanics with variable-order differential operators, Ann. Phys., 12 (2003), 692–703. https://doi.org/10.1002/andp.200310032 doi: 10.1002/andp.200310032 |
[25] | P. Zhuang, F. Liu, V. Anh, I. Turner, Numerical methods for the variable-order fractional advection-diffusion equation with a nonlinear source term, SIAM J. Numer. Anal., 47 (2009), 1760–1781. https://doi.org/10.1137/080730597 doi: 10.1137/080730597 |
[26] | C. Canuto, M. Y. Hussaini, A. Quarteroni, T. A. Zang, Spectral Methods: Fundamentals in Single Domains, Springer Science and Business Media, Berlin, 2007. |
[27] | D. Gottlieb, L. Lustman, E.Tadmor, Stability analysis of spectral methods for hyperbolic initialboundary value systems, SIAM J. Numer. Anal., 24 (1987), 241–256. https://doi.org/10.1137/0724020 doi: 10.1137/0724020 |
[28] | K. Shah, M. Akram, Numerical treatment of non-integer order partial differential equations by omitting discretization of data, Comput. Appl. Math., 37 (2018), 6700–6718. https://doi.org/10.1007/s40314-018-0706-3 doi: 10.1007/s40314-018-0706-3 |
[29] | Y. Chen, L. Liu, B. Li, Y. Sun, Numerical solution for the variable order linear cable equation with Bernstein polynomials, Appl. Math. Comput., 238 (2014), 329–341. https://doi.org/10.1016/j.amc.2014.03.066 doi: 10.1016/j.amc.2014.03.066 |
[30] | D. Valério, J. S. Da Costa, Variable-order fractional derivatives and their numerical approximations, Signal process., 91 (2011), 470–483. https://doi.org/10.1016/j.sigpro.2010.04.006 doi: 10.1016/j.sigpro.2010.04.006 |
[31] | Y. Chen, Y. Wei, D. Liu, H. Yu, Numerical solution for a class of nonlinear variable order fractional differential equations with Legendre wavelets, Appl. Math. Lett., 46 (2015), 83–88. |
[32] | M. A. Zaky, D. Baleanu, J. F. Alzaidy, E. Hashemizadeh, Operational matrix approach for solving the variable-order nonlinear Galilei invariant advection–diffusion equation, Advances in Difference Equations, 2018 (1), 1–11. |
[33] | Y. Yang, Y. Ma, L. Wang, Legendre polynomials operational matrix method for solving fractional partial differential equations with variable coefficients, Math. Probl. Eng., 2015 (2015), 915195. https://doi.org/10.1155/2015/915195 doi: 10.1155/2015/915195 |
[34] | K. Shah, F. Jarad, T. Abdeljawad, Stable numerical results to a class of time-space fractional partial differential equations via spectral method, J. Adv. Res., 25 (2020), 39–48. https://doi.org/10.1016/j.jare.2020.05.022 doi: 10.1016/j.jare.2020.05.022 |
[35] | I. Singh, S. Arora, S. Kumar, Numerical solution of wave equation using Haar wavelet, Int. J.Pure and Appl. Math., 98 (2015), 457–469. |
[36] | Y. Wang, Y. Chen, Shifted Legendre Polynomials algorithm used for the dynamic analysis of viscoelastic pipes conveying fluid with variable fractional order model, Appl. Math. Model., 81 (2020), 159–176. https://doi.org/10.1016/j.apm.2019.12.011 doi: 10.1016/j.apm.2019.12.011 |