This article presents an efficient numerical algorithm based on Legendre wavelets operational matrix for solving some fractional models. Fractional integration operational matrix of Legendre wavelets is derived and it is employed to reduce fractional differential equations into a system of algebraic equations. Mixing problems, Newton law of cooling problems and sugar inversion problems are included to elucidate the applicability and the simplicity of the Legendre wavelet-based numerical algorithm.