Does simultaneous solution matter for stellar evolution codes?
arXiv:astro-ph/0605583 · doi:10.1111/j.1365-2966.2006.10595.x
Abstract
A version of the STARS stellar evolution code has been developed that uses a non-simultaneous solution of the equations of stellar structure and evolution. In all other respects it is identical to the normal, fully simultaneous version. It is therefore possible to test the dependence of the solution on how the equations are solved. Two cases are investigated: a 5Msol and a 3Msol star, both of metallicity Z=0.02. Prior to the asymptotic giant branch, the models are almost identical. However once thermal pulses start, the two methods of solution yield diverging results with the non-simultaneous technique predicting longer interpulse periods. This is traced to difficulties associated with hydrogen burning caused by the use of a moving mesh. It is shown that, with careful control of the temporal resolution, the results of the simultaneous technique can be recovered.
7 pages, 7 figures. Accepted for publication in MNRAS
Cited by in corpus (5)
- The death of massive stars - I. Observational constraints on the progenitors of type II-P supernovae
- CLES, Code Liegeois d'Evolution Stellaire
- Mass loss and yield uncertainty in low-mass asymptotic giant branch stars
- Light element abundances in carbon-enhanced metal-poor stars
- A New, Efficient Stellar Evolution Code for Calculating Complete Evolutionary Tracks