CESAM: Un código para cálculos de evolución estelar

El código CESAM es un conjunto consistente de programas y rutinas que realiza cálculos de evolución estelar cuasi-estática 1D, incluyendo difusión y rotación. La principal novedad es la solución del equilibrio cuasi-estático por el método de colocalización basado en aproximaciones polinómicas por tramos proyectadas en su base de B-splines; lo que permite cálculos estables y robustos y la restitución exacta de la solución no solo en los puntos de la malla. Otra ventaja es el monitoreo por un solo parámetro de la precisión y su mejora por superconvergencia. Se ha diseñado un refinamiento automático de la malla para ajustar las localizaciones de los puntos de malla de acuerdo con los cambios de incógnitas, siendo cada límite entre zonas radiativas y convectivas desplazado al punto de malla más cercano. Para modelos estándar, la evolución de la composición química se resuelve mediante esquemas rígidamente estables de órdenes de hasta cuatro; para modelos no estándar, la solución de la ecuación de difusión emplea el esquema de Petrov-Galerkin, con la mezcla de productos químicos en zonas convectivas realizada por difusión turbulenta fuerte. Se permite una restauración precisa de la atmósfera. CESAM calcula la evolución de estrellas desde la pre-secundaria hasta el comienzo del ciclo de quema de 4He. En este documento se presenta una descripción detallada de los algoritmos.