Numerical analysis of an operational Jacobi Tau method for fractional weakly singular integro-differential equations

Abstract The main concern of this paper is to develop and analyze an operational Tau method for obtaining the numerical solution of fractional weakly singular integro-differential equations when the Jacobi polynomials are used as natural basis functions. This strategy is an application of the matrix–vector–product approach in Tau formulation of the problem. We first study the regularity of the exact solution and show that some derivatives of the exact solution have a singularity at origin dependence on both order of fractional derivative and weakly singular kernel function which makes poor convergence results for the Tau discretization of the problem. In order to recover high-order of convergence, we propose a new variable transformation to regularize the given functions and then to approximate the solution via a satisfactory order of convergence using an operational Tau method. Convergence analysis of this novel method is presented and the expected spectral rate of convergence for the proposed method is established. Numerical results are given which confirm both the theoretical predictions obtained and efficiency of the proposed method.

[1]  Paulo B. Vasconcelos,et al.  Towards a Lanczos’ $$\tau $$τ-Method Toolkit for Differential Problems , 2016, Math. Comput. Sci..

[2]  P. Mokhtary Discrete Galerkin Method for Fractional Integro-Differential Equations , 2015, 1501.01111.

[3]  Li Zhu,et al.  SCW method for solving the fractional integro-differential equations with a weakly singular kernel , 2016, Appl. Math. Comput..

[4]  Somayeh Nemati,et al.  A fast numerical algorithm based on the second kind Chebyshev polynomials for fractional integro-differential equations with weakly singular kernels , 2016, J. Comput. Appl. Math..

[5]  Eduardo L. Ortiz,et al.  Numerical solution of differential eigenvalue problems with an operational approach to the Tau method , 2005, Computing.

[6]  Mingxu Yi,et al.  Legendre wavelets method for the numerical solution of fractional integro-differential equations with weakly singular kernel , 2016 .

[7]  Mokhtary Payam,et al.  Convergence analysis of spectral Tau method for fractional Riccati differential equations , 2014 .

[8]  I. Podlubny Fractional differential equations , 1998 .

[9]  K. Ito,et al.  On a Legendre Tau method for fractional boundary value problems with a Caputo derivative , 2016 .

[10]  Jie Shen,et al.  Spectral Methods: Algorithms, Analysis and Applications , 2011 .

[11]  P. Mokhtary Numerical treatment of a well-posed Chebyshev Tau method for Bagley-Torvik equation with high-order of accuracy , 2015, Numerical Algorithms.

[12]  P. Mokhtary,et al.  SPECTRAL COLLOCATION METHOD FOR MULTI-ORDER FRACTIONAL DIFFERENTIAL EQUATIONS , 2014 .

[13]  Jan S. Hesthaven,et al.  Spectral Methods for Time-Dependent Problems: Contents , 2007 .

[14]  René Escalante,et al.  Segmented Tau approximation for a parametric nonlinear neutral differential equation , 2007, Appl. Math. Comput..

[15]  H. Srivastava,et al.  Theory and Applications of Fractional Differential Equations , 2006 .

[16]  Jingjun Zhao,et al.  Collocation methods for fractional integro-differential equations with weakly singular kernels , 2013, Numerical Algorithms.

[17]  Tao Tang,et al.  Convergence analysis of Jacobi spectral collocation methods for Abel-Volterra integral equations of second kind , 2012 .

[18]  Bao-Qing Tang,et al.  Solution of a class of Volterra integral equations with singular and weakly singular kernels , 2008, Appl. Math. Comput..

[19]  Arvet Pedas,et al.  Smoothing transformation and spline collocation for nonlinear fractional initial and boundary value problems , 2017, J. Comput. Appl. Math..

[20]  T. A. Zang,et al.  Spectral Methods: Fundamentals in Single Domains , 2010 .

[21]  G. Vainikko,et al.  Smoothing transformation and spline collocation for weakly singular Volterra integro-differential equations , 2017 .

[22]  P. Mokhtary,et al.  Reconstruction of exponentially rate of convergence to Legendre collocation solution of a class of fractional integro-differential equations , 2015, J. Comput. Appl. Math..

[23]  Farideh Ghoreishi,et al.  The L2-convergence of the Legendre spectral Tau matrix formulation for nonlinear fractional integro differential equations , 2011, Numerical Algorithms.

[24]  René Escalante,et al.  Segmented Tau approximation for a forward-backward functional differential equation , 2011, Comput. Math. Appl..

[25]  Eduardo L. Ortiz,et al.  An operational approach to the Tau method for the numerical solution of non-linear differential equations , 1981, Computing.

[26]  N. Ford,et al.  Analysis of Fractional Differential Equations , 2002 .

[27]  Arvet Pedas,et al.  Smoothing transformation and spline collocation for linear fractional boundary value problems , 2016, Appl. Math. Comput..

[28]  V. V. Zozulya,et al.  Weakly singular, singular and hypersingular integrals in 3-D elasticity and fracture mechanics , 1999 .

[29]  P. Kythe,et al.  Computational Methods for Linear Integral Equations , 2002 .