One of the ongoing issues with time fractional diffusion models is the design of efficient high-order numerical schemes for the solutions of limited regularity. We construct in this paper two efficient Galerkin spectral algorithms for solving multi-dimensional time fractional advection–diffusion–reaction equations with constant and variable coefficients. The model solution is discretized in time with a spectral expansion of fractional-order Jacobi orthogonal functions. For the space discretization, the proposed schemes accommodate high-order Jacobi Galerkin spectral discretization. The numerical schemes do not require imposition of artificial smoothness assumptions in time direction as is required for most methods based on polynomial interpolation. We illustrate the flexibility of the algorithms by comparing the standard Jacobi and the fractional Jacobi spectral methods for three numerical examples. The numerical results indicate that the global character of the fractional Jacobi functions makes them well-suited to time fractional diffusion equations because they naturally take the irregular behavior of the solution into account and thus preserve the singularity of the solution.
相似文献This paper presents a spectral collocation technique to solve fractional stochastic Volterra integro-differential equations (FSV-IDEs). The algorithm is based on shifted fractional order Legendre orthogonal functions generated by Legendre polynomials. The shifted fractional order Legendre–Gauss–Radau collocation (SFL-GR-C) method is developed for approximating the FSV-IDEs, with the objective of obtaining a system of algebraic equations. For computational purposes, the Brownian motion function W(x) is discretized by Lagrange interpolation, while the integral terms are interpolated by Legendre–Gauss–Lobatto quadrature. Numerical examples demonstrate the accuracy and applicability of the proposed technique, even when dealing with non-smooth solutions.
相似文献