Multi-term time-fractional advection diffusion equations are vital for simulating a wide range of physical phenomena, including fluid dynamics and environmental transport processes. However, due to their natural complexity, these equations pose challenges for conventional numerical approaches. In this article, we develop a high order accurate method to solve the multi-term time-fractional advection diffusion equations. We combine the Laplace transform (LT) to integrate the considered equations in time, with Chebyshev spectral method (CSM) for spatial terms. The proposed method produces highly accurate solutions with remarkably low computational
cost as compared to finite difference method. The propose numerical scheme first employs the LT which reduces the considered problem into a finite set of elliptic equations which may be solved in parallel. Then, the CSM is employed for the disctrezation of spatial operators, which makes it possibly to accurately represent the solution chebyshev grid. Finally, numerical inversion of LT is used to convert the obtain solution from the Laplace domain into the real domain. This work utilizes the modified Talbot's method and Stehfest's method for numerical inversion of the LT. To measure the performance, efficiency, and accuracy of the suggested approach, numerical approximations of three models are acquired and verified against the exact solution. The outcomes presented in tables and figures demonstrate that the modified Talbot's method performed better as compared to Stehfest's method.