Sobre la solución numérica del problema de stefan unidimensional

Abstract

Se presenta un estudio de los métodos de tratamiento numérico de una interfase móvil a través del análisis de un problema de Stefan unidimensional. Se consideran los métodos de malla fija y variable en diferencias finitas. Se discuten sus ventajas e inconvenientes y se propone un nuevo método consistente en la utilización de malla variable, de una aproximación de orden superior en el calculo de la interfase, de una iteración compuesta en el calculo del paso de tiempo y de la utilización de una aproximación lineal basada en los métodos de malla fija como valor inicial en el calculo iterativo del paso de tiempo. Se presentan resultados numéricos que se comparan con los obtenidos por otros autores. De esta comparación surge que el método propuesto presenta la ventaja de ser mas preciso, de converger en un número menor de iteraciones y de utilizar menor tiempo de cómputo .