This paper addresses the numerical solution of the multi-dimensional variable-order fractional integro-partial differential equations. The upwind scheme and a piecewise linear interpolation, are proposed to approximate the Coimbra variable-order fractional derivatives and integral term with kernel, respectively. Two new approaches via the Sinc collocation method based on single and double exponential transformations are adopted for the temporal and spatial discretizations, respectively. The convergence behaviour of the methods is analysed and the error bounds are provided. In addition, four test problems illustrate the validity and effectiveness of the proposed algorithms.
