Resumen de Aproximación de la temperatura en modelos no lineales de masas de hielo politérmico

Natividad Calvo Ruibal , José Durany Castrillo , C. Vázquez

• En este trabajo se propone un método numérico para resolver un subproblema que aparece en la modelización de placas de hielo politérmico. Dicho problema presenta tres términos no lineales: un término de reacción debido a la disipación viscosa, una condición de frontera de tipo desigualdad asociada al flujo de calor geotérmico y un término de tipo entalpía originado por la formulación de tipo Stefan para el cambio de fase. La temperatura estacionaria se obtiene como límite de la solución de problemas evolutivos, discretizados en tiempo mediante un algoritmo de tipo transporte-difusión. Los aspectos no lineales se resuelven mediante técnicas de dualidad y el método de Newton. Tras la descripción del algoritmo, se presentan algunos resultados obtenidos para problemas con datos reales.