skip to main content


Title: 3D hydrodynamic simulations of C ingestion into a convective O shell
Abstract Interactions between convective shells in evolved massive stars have been linked to supernova impostors, to the production of the odd-Z elements Cl, K, and Sc, and they might also help generate the large-scale asphericities that are known to facilitate shock revival in supernova explosion models. We investigate the process of ingestion of C-shell material into a convective O-burning shell, including the hydrodynamic feedback from the nuclear burning of the ingested material. Our 3D hydrodynamic simulations span almost 3 dex in the total luminosity $L_\rm {tot}$. All but one of the simulations reach a quasi-stationary state with the entrainment rate and convective velocity proportional to $L_\rm {tot}$ and $L_\rm {tot}^{1/3}$, respectively. Carbon burning provides 14 – $33\%$ of the total luminosity, depending on the set of reactions considered. Equivalent simulations done on 7683 and 11523 grids are in excellent quantitative agreement. The flow is dominated by a few large-scale convective cells. An instability leading to large-scale oscillations with Mach numbers in excess of 0.2 develops in an experimental run with the energy yield from C burning increased by a factor of 10. This run represents most closely the conditions expected in a violent O-C shell merger, which is a potential production site for odd-Z elements such as K and Sc and which may seed asymmetries in the supernova progenitor. 1D simulations may underestimate the energy generation from the burning of ingested material by as much as a factor two owing to their missing the effect of clumpiness of entrained material on the nuclear reaction rate.  more » « less
Award ID(s):
1814181
NSF-PAR ID:
10166281
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
491
ISSN:
0035-8711
Page Range / eLocation ID:
972 - 992
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    We present a detailed near-infrared chemical abundance analysis of 10 red giant members of the Galactic open cluster NGC 752. High-resolution (R ≃ 45000) near-infrared spectral data were gathered with the Immersion Grating Infrared Spectrograph, providing simultaneous coverage of the complete H and K bands. We derived the abundances of H-burning (C, N, O), α (Mg, Si, S, Ca), light odd-Z (Na, Al, P, K), Fe-group (Sc, Ti, Cr, Fe, Co, Ni), and neutron-capture (Ce, Nd, Yb) elements. We report the abundances of S, P, K, Ce, and Yb in NGC 752 for the first time. Our analysis yields solar-metallicity and solar abundance ratios for almost all of the elements heavier than the CNO group in NGC 752. O and N abundances were measured from a number of OH and CN features in the H band, and C abundances were determined mainly from CO molecular lines in the K band. High-excitation $\rm{C\,\small {I}}$ lines present in both near-infrared and optical spectra were also included in the C abundance determinations. Carbon isotopic ratios were derived from the R-branch band heads of first overtone (2−0) and (3−1) 12CO and (2−0) 13CO lines near 23 440 Å and (3−1) 13CO lines at about 23 730 Å. The CNO abundances and 12C/13C ratios are all consistent with our giants having completed ‘first dredge-up’ envelope mixing of CN-cyle products. We independently assessed NGC 752 stellar membership from Gaia astrometry, leading to a new colour–magnitude diagram for this cluster. Applications of Victoria isochrones and MESA models to these data yield an updated NGC 752 cluster age (1.52 Gyr) and evolutionary stage indications for the programme stars. The photometric evidence and spectroscopic light element abundances all suggest that the most, perhaps all of the programme stars are members of the helium-burning red clump in this cluster.

     
    more » « less
  2. null (Ed.)
    ABSTRACT We present two mixing models for post-processing of 3D hydrodynamic simulations applied to convective–reactive i-process nucleosynthesis in a rapidly accreting white dwarf (RAWD) with [Fe/H] = −2.6, in which H is ingested into a convective He shell. A 1D advective two-stream model adopts physically motivated radial and horizontal mixing coefficients constrained by 3D hydrodynamic simulations. A simpler approach uses diffusion coefficients calculated from the same simulations. All 3D simulations include the energy feedback of the 12C(p, γ)13N reaction from the H entrainment. Global oscillations of shell H ingestion in two of the RAWD simulations cause bursts of entrainment of H and non-radial hydrodynamic feedback. With the same nuclear network as in the 3D simulations, the 1D advective two-stream model reproduces the rate and location of the H burning within the He shell closely matching the 3D simulation predictions, as well as qualitatively displaying the asymmetry of the XH profiles between the upstream and downstream. With a full i-process network the advective mixing model captures the difference in the n-capture nucleosynthesis in the upstream and downstream. For example, 89Kr and 90Kr with half-lives of $3.18\,\,\mathrm{\mathrm{min}}$ and $32.3\,\,\mathrm{\mathrm{s}}$ differ by a factor 2–10 in the two streams. In this particular application the diffusion approach provides globally the same abundance distribution as the advective two-stream mixing model. The resulting i-process yields are in excellent agreement with observations of the exemplary CEMP-r/s star CS31062-050. 
    more » « less
  3. ABSTRACT

    Our knowledge of stellar evolution is driven by one-dimensional (1D) simulations. 1D models, however, are severely limited by uncertainties on the exact behaviour of many multidimensional phenomena occurring inside stars, affecting their structure and evolution. Recent advances in computing resources have allowed small sections of a star to be reproduced with multi-D hydrodynamic models, with an unprecedented degree of detail and realism. In this work, we present a set of 3D simulations of a convective neon-burning shell in a 20 M⊙ star run for the first time continuously from its early development through to complete fuel exhaustion, using unaltered input conditions from a 321D-guided 1D stellar model. These simulations help answer some open questions in stellar physics. In particular, they show that convective regions do not grow indefinitely due to entrainment of fresh material, but fuel consumption prevails over entrainment, so when fuel is exhausted convection also starts decaying. Our results show convergence between the multi-D simulations and the new 321D-guided 1D model, concerning the amount of convective boundary mixing to include in stellar models. The size of the convective zones in a star strongly affects its structure and evolution; thus, revising their modelling in 1D will have important implications for the life and fate of stars. This will thus affect theoretical predictions related to nucleosynthesis, supernova explosions, and compact remnants.

     
    more » « less
  4. ABSTRACT

    We present a first 3D magnetohydrodynamic (MHD) simulation of oxygen, neon, and carbon shell burning in a rapidly rotating $16\hbox{-}\mathrm{M}_\odot$ core-collapse supernova progenitor. We also run a purely hydrodynamic simulation for comparison. After $\mathord \approx 180\mathrm{s}$ ($\mathord \approx$ 15 and 7 convective turnovers, respectively), the magnetic fields in the oxygen and neon shells achieve saturation at 1011 and 5 × 1010 G. The strong Maxwell stresses become comparable to the radial Reynolds stresses and eventually suppress convection. The suppression of mixing by convection and shear instabilities results in the depletion of fuel at the base of the burning regions, so that the burning shell eventually move outward to cooler regions, thus reducing the energy generation rate. The strong magnetic fields efficiently transport angular momentum outwards, quickly spinning down the rapidly rotating convective oxygen and neon shells and forcing them into rigid rotation. The hydrodynamic model shows complicated redistribution of angular momentum and develops regions of retrograde rotation at the base of the convective shells. We discuss implications of our results for stellar evolution and for the subsequent core-collapse supernova. The rapid redistribution of angular momentum in the MHD model casts some doubt on the possibility of retaining significant core angular momentum for explosions driven by millisecond magnetars. However, findings from multidimensional models remain tentative until stellar evolution calculations can provide more consistent rotation profiles and estimates of magnetic field strengths to initialize multidimensional simulations without substantial numerical transients. We also stress the need for longer simulations, resolution studies, and an investigation of non-ideal effects.

     
    more » « less
  5. null (Ed.)
    Context. A realistic parametrization of convection and convective boundary mixing in conventional stellar evolution codes is still the subject of ongoing research. To improve the current situation, multidimensional hydrodynamic simulations are used to study convection in stellar interiors. Such simulations are numerically challenging, especially for flows at low Mach numbers which are typical for convection during early evolutionary stages. Aims. We explore the benefits of using a low-Mach hydrodynamic flux solver and demonstrate its usability for simulations in the astrophysical context. Simulations of convection for a realistic stellar profile are analyzed regarding the properties of convective boundary mixing. Methods. The time-implicit Seven-League Hydro (SLH) code was used to perform multidimensional simulations of convective helium shell burning based on a 25  M ⊙ star model. The results obtained with the low-Mach AUSM + -up solver were compared to results when using its non low-Mach variant AUSM B + -up. We applied well-balancing of the gravitational source term to maintain the initial hydrostatic background stratification. The computational grids have resolutions ranging from 180 × 90 2 to 810 × 540 2 cells and the nuclear energy release was boosted by factors of 3 × 10 3 , 1 × 10 4 , and 3 × 10 4 to study the dependence of the results on these parameters. Results. The boosted energy input results in convection at Mach numbers in the range of 10 −3 –10 −2 . Standard mixing-length theory predicts convective velocities of about 1.6 × 10 −4 if no boosting is applied. The simulations with AUSM + -up show a Kolmogorov-like inertial range in the kinetic energy spectrum that extends further toward smaller scales compared with its non low-Mach variant. The kinetic energy dissipation of the AUSM + -up solver already converges at a lower resolution compared to AUSM B + -up. The extracted entrainment rates at the boundaries of the convection zone are well represented by the bulk Richardson entrainment law and the corresponding fitting parameters are in agreement with published results for carbon shell burning. However, our study needs to be validated by simulations at higher resolution. Further, we find that a general increase in the entropy in the convection zone may significantly contribute to the measured entrainment of the top boundary. Conclusion. This study demonstrates the successful application of the AUSM + -up solver to a realistic astrophysical setup. Compressible simulations of convection in early phases at nominal stellar luminosity will benefit from its low-Mach capabilities. Similar to other studies, our extrapolated entrainment rate for the helium-burning shell would lead to an unrealistic growth of the convection zone if it is applied over the lifetime of the zone. Studies at nominal stellar luminosities and different phases of the same convection zone are needed to detect a possible evolution of the entrainment rate and the impact of radiation on convective boundary mixing. 
    more » « less