Abstract:Degenerate ignition of helium in low-mass stars at the end of the red giant branch phase leads to dynamic convection in their helium cores. One-dimensional (1D) stellar modeling of this intrinsically multi-dimensional dynamic event is likely to be inadequate. Previous hydrodynamic simulations imply that the single convection zone in the helium core of metal-rich Pop I stars grows during the flash on a dynamic timescale. This may lead to hydrogen injection into the core, and a double convection zone structure as known from one-dimensional core helium flash simulations of low-mass Pop III stars. We perform hydrodynamic simulations of the core helium flash in two and three dimensions to better constrain the nature of these events. To this end we study the hydrodynamics of convection within the helium cores of a 1.25 \Msun metal-rich Pop I star (Z=0.02), and a 0.85 \Msun metal-free Pop III star (Z=0) near the peak of the flash. These models possess single and double convection zones, respectively. We use 1D stellar models of the core helium flash computed with state-of-the-art stellar evolution codes as initial models for our multidimensional hydrodynamic study, and simulate the evolution of these models with the Riemann solver based hydrodynamics code Herakles which integrates the Euler equations coupled with source terms corresponding to gravity and nuclear burning. The hydrodynamic simulation of the Pop I model involving a single convection zone covers 27 hours of stellar evolution, while the first hydrodynamic simulations of a double convection zone, in the Pop III model, span 1.8 hours of stellar life. We find differences between the predictions of mixing length theory and our hydrodynamic simulations. The simulation of the single convection zone in the Pop I model shows a strong growth of the size of the convection zone due to turbulent entrainment. Hence we predict that for the Pop I model a hydrogen injection phase (i.e. hydrogen injection into the helium core) will commence after about 23 days, which should eventually lead to a double convection zone structure known from 1D stellar modeling of low-mass Pop III stars. Our two and three-dimensional hydrodynamic simulations of the double (Pop III) convection zone model show that the velocity field in the convection zones is different from that predicted by stellar evolutionary calculations. The simulations suggest that the double convection zone decays quickly, the flow eventually being dominated by internal gravity waves.

3D hydrodynamics simulations of core convection in supermassive main-sequence stars

3D hydrodynamic simulations of massive main-sequence stars. III. The effect of radiation pressure and diffusion leading to a 1D equilibrium model

3D hydrodynamics simulations of a 3 $M_{\odot}$ core-helium burning star

3D hydrodynamic simulations of massive main-sequence stars. I. Dynamics and mixing of convection and internal gravity waves

Multidimensional Simulations of Core Convection

3D simulations of convective shell Neon-burning in a massive star

3D hydrodynamic simulations of massive main-sequence stars II. Convective excitation and spectra of internal gravity waves

3D hydrodynamic simulations of massive main-sequence stars – II. Convective excitation and spectra of internal gravity waves

3D simulations of a neon burning convective shell in a massive star

3D Hydrodynamic Simulations of Carbon Burning in Massive Stars

Looking Deep Within an A-type Star: Core Convection Under the Influence of Rotation

Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes

Towards a self-consistent model of the convective core boundary in upper main sequence stars. Part I: 2.5D and 3D simulations

3D stellar evolution: hydrodynamic simulations of a complete burning phase in a massive star

Three-Dimensional Simulations of Massive Stars: I. Wave Generation and Propagation

Hydrodynamic Processes in Massive Stars

The core helium flash revisited III. From Pop I to Pop III stars

Signatures of convection in the atmospheres of cool evolved stars

Simulations of core convection in rotating A-type stars: Magnetic dynamo action

3D Simulations of Magnetoconvection in a Rapidly Rotating Supernova Progenitor