LAUSR

This paper deals with developing a fast and robust numerical formulation to simulate a system of fractional PDEs. At the first stage, the time variable is approximated by a finite difference method with first-order accuracy. At the second stage, the spectral Galerkin method based upon the fractional Jacobi polynomials is employed to discretize the spatial variables. We apply a reduced-order method based upon the proper orthogonal decomposition technique to decrease the utilized computational time. The unconditional stability property and the order of convergence of the new technique are analyzed in detail. The proposed numerical technique is well known as the reduced-order spectral Galerkin scheme. Furthermore, by employing the Newton–Raphson method and semi-implicit schemes, the proposed method can be used for solving linear and nonlinear ODEs and PDEs. Finally, some examples are provided to confirm the theoretical results.