A Multidimensional Jacobi-Based Spectral Framework for 3D Time-Fractional Diffusion and Transport Equations

This work presents a new and efficient numerical framework for solving three-dimensional time-fractional diffusion and mobile–immobile equations in the Caputo sense. The method is formulated using four-variable Jacobi polynomials, constructed systematically via the Kronecker product of one-dimensional Jacobi bases to accurately represent the multidimensional nature of the governing equations. Within a pseudo-operational collocation formulation, these polynomials enable a highly accurate and computationally efficient approximation of the fractional operators in both temporal and spatial directions. From the theoretical standpoint, the existence and uniqueness of the approximate solution are rigorously established through Schauder’s fixed-point theorem. Furthermore, the Ulam–Hyers stability of the numerical solution is verified, demonstrating the robustness of the method with respect to perturbations in the input data. To reinforce the reliability of the approach, an explicit error bound for the residual function is derived in a Jacobi-weighted Sobolev space, offering a firm analytical basis for assessing convergence. Numerical experiments confirm that the proposed approach achieves superior accuracy and efficiency, highlighting its potential as a powerful tool for high-dimensional fractional partial differential equations.