3 research outputs found
A M\"untz-Collocation spectral method for weakly singular volterra integral equations
In this paper we propose and analyze a fractional Jacobi-collocation spectral
method for the second kind Volterra integral equations (VIEs) with weakly
singular kernel . First we develop a family of fractional
Jacobi polynomials, along with basic approximation results for some weighted
projection and interpolation operators defined in suitable weighted Sobolev
spaces. Then we construct an efficient fractional Jacobi-collocation spectral
method for the VIEs using the zeros of the new developed fractional Jacobi
polynomial. A detailed convergence analysis is carried out to derive error
estimates of the numerical solution in both - and weighted
-norms. The main novelty of the paper is that the proposed method is
highly efficient for typical solutions that VIEs usually possess. Precisely, it
is proved that the exponential convergence rate can be achieved for solutions
which are smooth after the variable change for a
suitable real number . Finally a series of numerical examples are
presented to demonstrate the efficiency of the method
An approach to construct higher order time discretisation schemes for time fractional partial differential equations with nonsmooth data
Invited review article for Anniversary Edition of Journal.In this paper, we shall review an approach by which we can seek higher order time discretisation schemes for solving time fractional partial differential equations with nonsmooth data. The low regularity of the solutions of time fractional partial differential equations implies standard time discretisation schemes only yield first order accuracy. To obtain higher order time discretisation schemes when the solutions of time fractional partial differential equations have low regularities, one may correct the starting steps of the standard time discretisation schemes to capture the singularities of the solutions. We will consider these corrections of some higher order time discretisation schemes obtained by using Lubich's fractional multistep methods, L1 scheme and its modification, discontinuous Galerkin methods, etc. Numerical examples are given to show that the theoretical results are consistent with the numerical results