Diffusion in stellar interiors: critical tests of three numerical methods

Abstract

We describe and discuss the properties of three numerical methods for solving the diffusion equation for the transport of the chemical species and of the angular momentum in stellar interiors. We study through numerical experiments both their accuracy and their ability to provide physical solutions. On the basis of new tests and analyses applied to the stellar astrophysical context, we show that the most robust method to follow the secular evolution is the implicit finite differences method. The importance of correctly estimating the diffusion coefficient between mesh points is emphasized and a procedure for estimating the average diffusion coefficient between a convective and a radiative zone is described.

15 pages, 8 figures. accepted for publication in Astronomy and Astrophysics

References in corpus (1)

• Stellar evolution with rotation X: Wolf-Rayet star populations at solar metallicity

Cited by in corpus (11)

• Evolution of Thermally Pulsing Asymptotic Giant Branch Stars I. The COLIBRI Code
• Stellar evolution with rotation XII: Pre-supernova models
• Evolution and Nucleosynthesis of Extremely Metal Poor & Metal-Free Low- and Intermediate-Mass Stars I: Stellar Yield Tables and the CEMPs
• Can Extra Mixing in RGB and AGB Stars Be Attributed to Magnetic Mechanisms?
• Evolution and nucleosynthesis of extremely metal-poor and metal-free low- and intermediate-mass stars II. s-process nucleosynthesis during the core He flash
• Grids of stellar models with rotation VI: Models from 0.8 to 120 $M_{\odot }$ at a metallicity Z = 0.006
• The Gaia-ESO survey: 3D NLTE abundances in the open cluster NGC 2420 suggest atomic diffusion and turbulent mixing at the origin of chemical abundance variations
• Constraints on rotational mixing from surface evolution of light elements in massive stars
• Confronting uncertainties in stellar physics II. exploring differences in main-sequence stellar evolution tracks
• The effects of surface fossil magnetic fields on massive star evolution: III. The case of $\tau$ Sco
• Numerical Tests of Rotational Mixing in Massive Stars with the new Population Synthesis Code BONNFIRES