En este trabajo, se presentan resultados de la simulación numérica de la descomposición de fases para sistemas binarios y ternarios de aleaciones metálicas. Estos resultados se obtuvieron a partir de la ecuación diferencial parcial no lineal de Cahn y Hilliard. Se propone el uso del Modelo de Solución Regular para la Energía libre. La ecuación se resuelve por el método explícito de Diferencias finitas, dando como resultado la evolución microestructural y aproximaciones en la cinética de descomposición de fases de aleaciones binarias y ternarias. Los resultados de las simulaciones permiten predecir y analizar la cinética y morfologia de la descomposición de fases en las primeras etapas de envejecido de aleaciones metálicas, para analizar el efecto de los diferentes parámetros de estas ecuaciones sobre la descomposición de fases durante un tratamiento térmico.