Este documento trata sobre las soluciones numéricas de una clase de modelos matemáticos fraccionarios que surgen en ciencias de la ingeniería gobernadas por ecuaciones de advección-difusión-reacción fraccionarias en el tiempo (TF-ADR), que involucran la derivada de Caputo. En particular, estamos interesados en los modelos que vinculan procesos químicos e hidrodinámicos. El objetivo de este documento es proponer un método de diferencias finitas implícito, simple y robusto, incondicionalmente estable para resolver las ecuaciones TF-ADR. Los resultados numéricos muestran que el método propuesto es eficiente, confiable y fácil de implementar desde un punto de vista computacional y puede ser empleado para problemas en ciencias de la ingeniería.