<abstract><p>In this article, the meshfree radial basis function method based on the Gaussian function is proposed for some time-space dependent fractional order partial differential equation (PDE) models. These PDE models have significant applications in chemical engineering and physical science. Some main advantages of the proposed method are that it is easy to implement, and the output response is quick and highly accurate, especially in the higher dimension. In this method, the time-dependent derivative terms are treated by Caputo fractional derivative while space-dependent derivative terms are treated by Riesz, Riemann-Liouville, and Grünwald-Letnikov derivatives. The proposed method is tested on some numerical examples and the accuracy is analyzed by $ \|L\|_\infty $.</p></abstract>