We present 3D implicit large eddy simulations (ILES) of the turbulent convection in the envelope of a 5 Msun red giant star and in the oxygen-burning shell of a 23 Msun supernova progenitor. The numerical models are analyzed in the framework of 1D Reynolds-Averaged Navier-Stokes (RANS) equations. The effects of pressure fluctuations are more important in the red giant model, owing to larger stratification of the convective zone. We show how this impacts different terms in the mean-field equations. We clarify the driving sources of kinetic energy, and show that the rate of turbulent dissipation is comparable to the convective luminosity. Although our flows have low Mach number and are nearly adiabatic, our analysis is general and can be applied to photospheric convection as well. The robustness of our analysis of turbulent convection is supported by the insensitivity of the mean-field balances to linear mesh resolution. We find robust results for the turbulent convection zone and the stable layers in the oxygen-burning shell model, and robust results everywhere in the red giant model, but the mean fields are not well converged in the narrow boundary regions (which contain steep gradients) in the oxygen-burning shell model. This last result illustrates the importance of unresolved physics at the convective boundary, which governs the mixing there.

The injection of hydrogen into the convection shell powered by helium burning during the core helium flash is commonly encountered during the evolution of metal-free and extremely metal-poor low-mass stars. With specifically designed multidimensional hydrodynamic simulations, we aim to prove that an entropy barrier is no obstacle for the growth of the helium-burning shell convection zone in the helium core of a metal-rich Pop I star, i.e. convection can penetrate into the hydrogen-rich layers for these stars, too. We further study whether this is also possible in one-dimensional stellar evolutionary calculations. Our hydrodynamical simulations show that the helium-burning shell convection zone in the helium core moves across the entropy barrier and reaches the hydrogen-rich layers. This leads to mixing of protons into the hotter layers of the core and to a rapid increase of the nuclear energy production at the upper edge of the helium-burning convection shell - the hydrogen injection flash. As a result a second convection zone appears in the hydrogen-rich layers. Contrary to 1D models, the entropy barrier separating the two convective shells from each other is largely permeable to chemical transport when allowing for multidimensional flow, and consequently, hydrogen is continuously mixed deep into the helium core. We find it difficult to achieve such a behavior in one-dimensional stellar evolutionary calculations.

Shell convection driven by nuclear burning in a stellar core is a common hydrodynamic event in the evolution of many types of stars. We encounter and simulate this convection (i) in the helium core of a low-mass red giant during core helium flash leading to a dredge-down of protons across an entropy barrier, (ii) in a carbon-oxygen core of an intermediate-mass star during core carbon flash, and (iii) in the oxygen and carbon burning shell above the silicon-sulfur rich core of a massive star prior to supernova explosion. Our results, which were obtained with the hydrodynamics code HERAKLES, suggest that both entropy gradients and entropy barriers are less important for stellar structure than commonly assumed. Our simulations further reveal a new dynamic mixing process operating below the base of shell convection zones.

Mixing processes in stars driven by composition gradients as a result of the Rayleigh-Taylor instability are not anticipated. They are supported only by hydrodynamic studies of stellar convection. We find that such mixing occurs below the bottom edge of convection zones in our multidimensional hydrodynamic shell convection models. It operates at interfaces created by off-center nuclear burning, where less dense gas with higher mean molecular weight is located above denser gas with a lower mean molecular weight. We discuss the mixing under various conditions with hydrodynamic convection models based on stellar evolutionary calculations of the core helium flash in a 1.25 Msun star, the core carbon flash in a 9.3 Msun star, and of oxygen burning shell in a star with a mass of 23 Msun. We simulate the hydrodynamic behavior of shell convection during various phases of stellar evolution with the Eulerian hydrodynamics code HERAKLES in two and three spatial dimensions. Initial models for this purpose are obtained by state-of-the-art stellar evolutionary codes, namely GARSTEC, STAREVOL, and TYCHO for the core helium flash, core carbon flash, and oxygen shell burning, respectively. Most of our analysis is performed for two-dimensional hydrodynamic models of shell convection during the core helium flash at its peak covering approximately 250 convective turnover timescales. The mixing manifests itself in the form of overdense and cold fingers enriched with matter of higher mean molecular weight, originating from density fluctuations at the lower boundary of the convective shell, and "shooting" down into the core. They result from the Rayleigh-Taylor instability at the lower convection zone boundary due to a negative mean molecular weight gradient. They do not appear when the mean molecular weight gradient is positive.

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.

We describe and discuss hydrodynamic simulations of the core helium flash using an initial model of a 1.25 M_sol star with a metallicity of 0.02 near at its peak. Past research concerned with the dynamics of the core helium flash is inconclusive. Its results range from a confirmation of the standard picture, where the star remains in hydrostatic equilibrium during the flash (Deupree 1996), to a disruption or a significant mass loss of the star (Edwards 1969; Cole & Deupree 1980). However, the most recent multidimensional hydrodynamic study (Dearborn 2006) suggests a quiescent behavior of the core helium flash and seems to rule out an explosive scenario. Here we present partial results of a new comprehensive study of the core helium flash, which seem to confirm this qualitative behavior and give a better insight into operation of the convection zone powered by helium burning during the flash. The hydrodynamic evolution is followed on a computational grid in spherical coordinates using our new version of the multi-dimensional hydrodynamic code HERAKLES, which is based on a direct Eulerian implementation of the piecewise parabolic method.

We study turbulent convection during the core helium flash close to its peak by comparing the results of two and three-dimensional hydrodynamic simulations. We use a multidimensional Eulerian hydrodynamics code based on state-of-the-art numerical techniques to simulate the evolution of the helium core of a $1.25 M_{\odot}$ Pop I star. Our three-dimensional hydrodynamic simulations of the evolution of a star during the peak of the core helium flash do not show any explosive behavior. The convective flow patterns developing in the three-dimensional models are structurally different from those of the corresponding two-dimensional models, and the typical convective velocities are smaller than those found in their two-dimensional counterparts. Three-dimensional models also tend to agree better with the predictions of mixing length theory. Our hydrodynamic simulations show the presence of turbulent entrainment that results in a growth of the convection zone on a dynamic time scale. Contrary to mixing length theory, the outer part of the convection zone is characterized by a sub-adiabatic temperature gradient.

We investigate the hydrodynamics of the core helium flash near its peak. Past research concerned with the dynamics of this event is inconclusive. However, the most recent multidimensional hydrodynamic studies suggest a quiescent behavior and seem to rule out an explosive scenario. Previous work indicated, that depending on initial conditions, employed turbulence models, grid resolution, and dimensionality of the simulation, the core helium flash leads either to the disruption of a low-mass star or to a quiescent quasi-hydrostatic evolution. We try to clarify this issue by simulating the evolution with advanced numerical methods and detailed microphysics. Assuming spherical or axial symmetry, we simulate the evolution of the helium core of a $1.25 M_{\odot}$ star with a metallicity Z=0.02 during the core helium flash at its peak with a grid-based hydrodynamics code. We find that the core helium flash neither rips the star apart, nor that it significantly alters its structure, as convection plays a crucial role in keeping the star in hydrostatic equilibrium. In addition, our simulations show the presence of overshooting, which implies new predictions concerning mixing of chemical species in red giants.