We performed 3-D simulations of proton-rich material entrainment into czw-rich He-shell flash convection and the subsequent H-ingestion flash that took place in the post-AGB star Sakurais object. Observations of the transient nature and anomalous abu
ndance features are available to validate our method and assumptions, with the aim to apply them to very low metallicity stars in the future. We include nuclear energy feedback from H burning and cover the full $4pi$ geometry of the shell. Runs on $768^3$ and $1536^3$ grids agree well with each other and have been followed for $1500mathrm{min}$ and $1200mathrm{min}$. After a $850mathrm{min}$ long quiescent entrainment phase the simulations enter into a global non-spherical oscillation that is launched and sustained by individual ignition events of H-rich fluid pockets. Fast circumferential flows collide at the antipode and cause the formation and localized ignition of the next H-overabundant pocket. The cycle repeats for more than a dozen times while its amplitude decreases. During the global oscillation the entrainment rate increases temporarily by a factor $approx 100$. Entrained entropy quenches convective motions in the upper layer until the burning of entrained H establishes a separate convection zone. The lower-resolution run hints at the possibility that another global oscillation, perhaps even more violent will follow. The location of the H-burning convection zone agrees with a 1-D model in which the mixing efficiency is calibrated to reproduce the light curve. The simulations have been performed at the NSF Blue Waters supercomputer at NCSA.
We present the first 3-dimensional, fully compressible gas-dynamics simulations in $4pi$ geometry of He-shell flash convection with proton-rich fuel entrainment at the upper boundary. This work is motivated by the insufficiently understood observed c
onsequences of the H-ingestion flash in post-AGB stars (Sakurais object) and metal-poor AGB stars. Our investigation is focused on the entrainment process at the top convection boundary and on the subsequent advection of H-rich material into deeper layers, and we therefore ignore the burning of the proton-rich fuel in this study. We find that, for our deep convection zone, coherent convective motions of near global scale appear to dominate the flow. At the top boundary convective shear flows are stable against Kelvin-Helmholtz instabilities. However, such shear instabilities are induced by the boundary-layer separation in large-scale, opposing flows. This links the global nature of thick shell convection with the entrainment process. We establish the quantitative dependence of the entrainment rate on grid resolution. With our numerical technique simulations with $1024^3$ cells or more are required to reach a numerical fidelity appropriate for this problem. However, only the result from the $1536^3$ simulation provides a clear indication that we approach convergence with regard to the entrainment rate. Our results demonstrate that our method, which is described in detail, can provide quantitative results related to entrainment and convective boundary mixing in deep stellar interior environments with veryvstiff convective boundaries. For the representative case we study in detail, we find an entrainment rate of $4.38 pm 1.48 times 10^{-13}M_odot mathrm{/s}$.
Depending on mass and metallicity as well as evolutionary phase, stars occasionally experience convective-reactive nucleosynthesis episodes. We specifically investigate the situation when nucleosynthetically unprocessed, H-rich material is convective
ly mixed with a He-burning zone, for example in convectively unstable shell on top of electron-degenerate cores in AGB stars, young white dwarfs or X-ray bursting neutron stars. Such episodes are frequently encountered in stellar evolution models of stars of extremely low or zero metal content [...] We focus on the convective-reactive episode in the very-late thermal pulse star Sakurais object (V4334 Sagittarii). Asplund etal. (1999) determined the abundances of 28 elements, many of which are highly non-solar, ranging from H, He and Li all the way to Ba and La, plus the C isotopic ratio. Our simulations show that the mixing evolution according to standard, one-dimensional stellar evolution models implies neutron densities in the He that are too low to obtain a significant neutron capture nucleosynthesis on the heavy elements. We have carried out 3D hydrodynamic He-shell flash convection [...] we assume that the ingestion process of H into the He-shell convection zone leads only after some delay time to a sufficient entropy barrier that splits the convection zone [...] we obtain significantly higher neutron densities (~few 10^15 1/cm^3) and reproduce the key observed abundance trends found in Sakurais object. These include an overproduction of Rb, Sr and Y by about 2 orders of magnitude higher than the overproduction of Ba and La. Such a peculiar nucleosynthesis signature is impossible to obtain with the mixing predictions in our one-dimensional stellar evolution models. [...] We determine how our results depend on uncertainties of nuclear reaction rates, for example for the C13(alpha, n)O16 reaction.
The He-shell flash convection in AGB stars is the site for the high-temperature component of the s-process in low- and intermediate mass giants, driven by the Ne22 neutron source. [...] The upper convection boundary plays a critical role during the H
-ingestion episode that may lead to neutron-bursts in the most metal-poor AGB stars. We address these problems through global 3-dimensional hydrodynamic simulations including the entire spherical He-shell flash convection zone (as oposed to the 3D box-in-a-star simulations). An important aspect of our current effort is to establish the feasibility of our appoach. We explain why we favour the explicit treatment over the anelastic approximation for this problem. The simulations presented in this paper use a Cartesian grid of 512^3 cells and have been run on four 8-core workstations for four days to simulate ~5000s, which corresponds to almost ten convective turn-over times. The convection layer extends radially at the simulated point in the flash evolution over 7 H_p pressure scale-heights and exceeds the size of the underlying core. Convection is dominated by large convective cells that fill more than an entire octant. [...]
