En este artículo, se resuelven las ecuaciones de reacción-difusión fraccionarias en el tiempo utilizando un método de elemento de volumen finito mixto (MFVE) y la fórmula de aproximación de la derivada fraccionaria de Caputo. Se proporciona el análisis de existencia, unicidad y estabilidad incondicional para el esquema MFVE completamente discreto. Se derivan estimaciones de error para la variable desconocida escalar (en norma) y la variable auxiliar vectorial (en norma y en norma). Finalmente, se presentan dos ejemplos numéricos en regiones espaciales unidimensionales y bidimensionales para examinar la viabilidad y efectividad.