Chebyshev Pseudospectral Method for Fractional Differential Equations in Non-Overlapping Partitioned Domains

Fractional differential operators are inherently non-local, so global methods, such as spectral methods, are well suited for handling these non-local operators. Long-time integration of differential models such as chaotic dynamical systems poses specific challenges and considerations that make multi-domain numerical methods advantageous when dealing with such problems. This study proposes a novel multi-domain pseudospectral method based on the first kind of Chebyshev polynomials and the Gauss–Lobatto quadrature for fractional initial value problems.The proposed technique involves partitioning the problem’s domain into non-overlapping sub-domains, calculating the fractional differential operator in each sub-domain as the sum of the ‘local’ and ‘memory’ parts and deriving the corresponding differentiation matrices to develop the numerical schemes. The linear stability analysis indicates that the numerical scheme is absolutely stable for certain values of arbitrary non-integer order and conditionally stable for others. Numerical examples, ranging from single linear equations to systems of non-linear equations, demonstrate that the multi-domain approach is more appropriate, efficient and accurate than the single-domain scheme, particularly for problems with long-term dynamics.