Para comprender procesos de transición de fase como la solidificación, los modelos de campo de fase se emplean con frecuencia. En este documento, estudiamos un esquema Crank-Nicolson-ADI de diferencias finitas para una clase de sistemas isotrópicos parabólicos no lineales. Establecemos una estimación de error para este esquema, demostrando su efectividad en la resolución de modelos de campo de fase. Nuestro análisis proporciona una justificación matemática rigurosa para la fiabilidad del método numérico en la simulación de transiciones de fase.