Abstract
We construct in this paper two efficient spectral algorithms in the frequency space for solving unsteady advection-reaction-diffusion equations with constant and variable coefficients. We first consider a Jacobi-Galerkin method for solving linear equations with constant coefficients. We then develop the direct solution algorithm for the linear advection-reaction-diffusion equations with variable coefficients using the Jacobi-Galerkin method with numerical integration. The proposed Jacobi-Galerkin methods, both in temporal and spatial discretizations, are successfully developed to handle the two-dimensional unsteady advection-reaction-diffusion equations with constant and variable coefficients and with fractional orders. In these methods, the model solution is expanded in both space and time in terms of polynomials bases built upon a linear combination of Jacobi polynomials. The homogeneous initial and Dirichlet boundary conditions are satisfied exactly by expanding the model solution in terms of these polynomials. The proposed Jacobi-Galerkin methods yield an exponential rate of convergence when the solution is smooth and allow a great flexibility to handle multi-dimensional time fractional advection-reaction-diffusion equations. Finally, a series of numerical examples are presented to demonstrate the efficiency and flexibility of the methods.