Relative importance of convective uncertainties in massive stars

Edelmann

et al. 2021

A&A

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 to 10 −2 . Standard mixing-length theory (MLT) 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.

Section: Introductionmentioning

confidence: 65%

Section: Construction Of the Initial Modelmentioning

confidence: 99%

See 1 more Smart Citation

Multidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star

Horst

Edelmann

et al. 2021

A&A

“…The step overshoot model starts with a much larger core. This will leave an imprint on the structure of that part of the star, which will affect the behaviour of the intermediate convective zone (Kaiser et al 2020).…”

Section: Impact Of Entrainment On Helium Core Masses and Lifetimesmentioning

confidence: 99%

Convective core entrainment in 1D main-sequence stellar models

Scott

Monthly Notices of the Royal Astronomical Society

Georgy

et al. 2021

Self Cite

3D hydrodynamics models of deep stellar convection exhibit turbulent entrainment at the convective-radiative boundary which follows the entrainment law, varying with boundary penetrability. We implement the entrainment law in the 1D Geneva stellar evolution code. We then calculate models between 1.5 and 60 M⊙ at solar metallicity (Z = 0.014) and compare them to previous generations of models and observations on the main sequence. The boundary penetrability, quantified by the bulk Richardson number, RiB, varies with mass and to a smaller extent with time. The variation of RiB with mass is due to the mass dependence of typical convective velocities in the core and hence the luminosity of the star. The chemical gradient above the convective core dominates the variation of RiB with time. An entrainment law method can therefore explain the apparent mass dependence of convective boundary mixing through RiB. New models including entrainment can better reproduce the mass dependence of the main sequence width using entrainment law parameters A ∼ 2 × 10−4 and n = 1. We compare these empirically constrained values to the results of 3D hydrodynamics simulations and discuss implications.

“…The extent and implementation of convective overshooting in stellar models has a large impact on the mass of the He and CO cores (e.g. Kaiser et al 2020). Three-dimensional models of lower mass stars favour the existence of such mixing at convective boundaries (e.g.…”

Section: Possibility Of H-he Shell Interactionsmentioning

confidence: 99%

Is GW190521 the merger of black holes from the first stellar generations?

Farrell

Groh

Monthly Notices of the Royal Astronomical Society: Letters

et al. 2020

Self Cite

GW190521 challenges our understanding of the late-stage evolution of massive stars and the effects of the pair instability in particular. We discuss the possibility that stars at low or zero metallicity could retain most of their hydrogen envelope until the pre-supernova stage, avoid the pulsational pair-instability regime, and produce a black hole with a mass in the mass gap by fallback. We present a series of new stellar evolution models at zero and low metallicity computed with the geneva and mesa stellar evolution codes and compare to existing grids of models. Models with a metallicity in the range 0–0.0004 have three properties that favour higher black hole (BH) masses. These are (i) lower mass-loss rates during the post main sequence phase, (ii) a more compact star disfavouring binary interaction, and (iii) possible H–He shell interactions which lower the CO core mass. We conclude that it is possible that GW190521 may be the merger of black holes produced directly by massive stars from the first stellar generations. Our models indicate BH masses up to 70–75 M⊙. Uncertainties related to convective mixing, mass loss, H–He shell interactions, and pair-instability pulsations may increase this limit to ∼85 M⊙.