An Efficient Spectral Method for a Class of Asymmetric Functional-Order Diffusion–Wave Equations Using Generalized Chelyshkov Wavelets

Asymmetric functional-order (variable-order) fractional diffusion–wave equations (FO-FDWEs) introduce considerable computational challenges, as the fractional order of the derivatives can vary spatially or temporally. To overcome these challenges, a novel spectral method employing generalized fractional-order Chelyshkov wavelets (FO-CWs) is developed to efficiently solve such equations. In this approach, the Riemann–Liouville fractional integral operator of variable order is evaluated in closed form via a regularized incomplete Beta function, enabling the transformation of the governing equation into a system of algebraic equations. This wavelet-based spectral scheme attains extremely high accuracy, yielding significantly lower errors than existing numerical techniques. In particular, numerical results show that the proposed method achieves notably improved accuracy compared to existing methods under the same number of basis functions. Its strong convergence properties allow high precision to be achieved with relatively few wavelet basis functions, leading to efficient computations. The method’s accuracy and efficiency are demonstrated on several practical diffusion–wave examples, indicating its suitability for real-world applications. Furthermore, it readily applies to a wide class of fractional partial differential equations (FPDEs) with spatially or temporally varying order, demonstrating versatility for diverse applications.