In this study, we present a numerical method for solving two--dimensional space--time fractional partial differential equations (FPDEs), where the solutions of the FPDEs are expanded in terms of the shifted Chebyshev polynomials. The numerical approximations are evaluated at the Chebyshev-Gauss-Lobatto points. In the case when the FPDE is nonlinear, we employ a Newton-Raphson approach to linearize the equation. Both the linear and nonlinear cases lead to a consistent system of linear algebraic equations. The scheme is tested on selected FPDEs and the numerical results show that the proposed numerical scheme is accurate and computationally efficient in terms of CPU times. To establish the accuracy of the method, we also present an error analysis which shows the convergence of the numerical method. These positive attributes make the proposed method a good approach for solving two--dimensional fractional partial differential equations with both space and time fractional orders.
This is an open access article distributed under the Creative Commons Attribution License which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.