We performed 3-D simulations of proton-rich material entrainment into 12 C-rich He-shell flash convection and the subsequent H-ingestion flash that took place in the post-AGB star Sakurai's object. Observations of the transient nature and anomalous abundance 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 4π geometry of the shell. Runs on 768 3 and 1536 3 grids agree well with each other and have been followed for 1500min and 1200min. After a 850min 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 ≈ 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 4π geometry of Heshell flash convection with proton-rich fuel entrainment at the upper boundary. This work is motivated by the insufficiently understood observed consequences of the H-ingestion flash in post-AGB stars (Sakurai's 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 boundarylayer 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 very stiff convective boundaries. For the representative case we study in detail, we find an entrainment rate of 4.38 ± 1.48 × 10 −13 M /s.

FES will improve the swallowing function of NPC patients with dysphagia and bring about better quality of life.

We report on numerical simulations of the detailed evolution of the single mode Rayleigh-Taylor [Lord Rayleigh, Scientific Papers II (Cambridge University Press, Cambridge, 1900), p. 200; G. I. Taylor, “The instability of liquid surfaces when accelerated in a direction perpendicular to their plane,” Proc. R. Soc. London, Ser. A 201, 192 (1950)10.1098/rspa.1950.0052; S. Chandrasekhar, Hydrodynamic and Hydromagnetic Stability (Oxford University Press, Oxford, 1961)] instability to late times and high aspect ratios. In contrast to established potential flow models that predict a terminal velocity and a constant Froude number at low Atwood numbers, we observe a complex sequence of events that can be summarized in four stages: I. Exponential growth of imposed perturbations, II. Saturation to terminal velocity, III. Reacceleration to a higher Froude number, and IV. Chaotic mixing. The observed reacceleration away from the Froude number predicted by potential flow theory is attributed to the appearance of secondary Kelvin–Helmholtz structures, and described with a modification to the potential flow model proposed by Betti and Sanz [R. Betti and J. Sanz, “Bubble acceleration in the ablative Rayleigh-Taylor instability,” Phys. Rev. Lett. 97, 205002 (2006)10.1103/PhysRevLett.97.205002]. The secondary KH instability is in turn sensitive to several parameters, and can be suppressed at large Atwood numbers, as well as viscosity (physical or numerical), with the bubble/spike velocity in each case reverting to the potential flow value. Our simulations delineate the change in dynamics of the primary and secondary instabilities due to changes in these flow parameters. When the flow is allowed to evolve to late times, further instability is observed, resulting in chaotic mixing which is quantified here. The increased atomic mixing due to small-scale structures results in a dramatic drop in the late-time Froude number. Spike behavior resembles bubbles at low A, while for large A, spikes approach free-fall – thus, the notion of a terminal velocity appears not to be applicable to spikes at any density difference. We expect the results to be relevant to turbulent mix models that are based on bubble growth and interaction.

No abstract