We study the properties of oscillatory double-diffusive convection (ODDC) in the presence of a uniform vertical background magnetic field. ODDC takes place in stellar regions that are unstable according to the Schwarzschild criterion and stable according to the Ledoux criterion (sometimes called semiconvective regions), which are often predicted to reside just outside the core of intermediate-mass main-sequence stars. Previous hydrodynamic studies of ODDC have shown that the basic instability saturates into a state of weak wave-like convection, but that a secondary instability can sometimes transform it into a state of layered convection, where layers then rapidly merge and grow until the entire region is fully convective. We find that magnetized ODDC has very similar properties overall, with some important quantitative differences. A linear stability analysis reveals that the fastest-growing modes are unaffected by the field, but that other modes are. Numerically, the magnetic field is seen to influence the saturation of the basic instability, overall reducing the turbulent fluxes of temperature and composition. This in turn affects layer formation, usually delaying it, and occasionally suppressing it entirely for sufficiently strong fields. Further work will be needed, however, to determine the field strength above which layer formation is actually suppressed in stars. Potential observational implications are briefly discussed.
Stellar evolution models calculate convective boundaries using either the Schwarzschild or Ledoux criterion, but confusion remains regarding which criterion to use. Here we present a 3D hydrodynamical simulation of a convection zone and adjacent radiative zone, including both thermal and compositional buoyancy forces. As expected, regions that are unstable according to the Ledoux criterion are convective. Initially, the radiative zone adjacent to the convection zone is Schwarzschild unstable but Ledoux stable due to a composition gradient. Over many convective overturn timescales, the convection zone grows via entrainment. The convection zone saturates at the size originally predicted by the Schwarzschild criterion, although in this final state the Schwarzschild and Ledoux criteria agree. Therefore, the Schwarzschild criterion should be used to determine the size of stellar convection zones, except possibly during short-lived evolutionary stages in which entrainment persists.
more » « less- Award ID(s):
- 1908338
- NSF-PAR ID:
- 10366803
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal Letters
- Volume:
- 928
- Issue:
- 1
- ISSN:
- 2041-8205
- Format(s):
- Medium: X Size: Article No. L10
- Size(s):
- ["Article No. L10"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
ABSTRACT In this work, we investigate the impact of uncertainties due to convective boundary mixing (CBM), commonly called ‘overshoot’, namely the boundary location and the amount of mixing at the convective boundary, on stellar structure and evolution. For this we calculated two grids of stellar evolution models with the MESA code, each with the Ledoux and the Schwarzschild boundary criterion, and vary the amount of CBM. We calculate each grid with the initial masses of 15, 20, and $25\, \rm {M}_\odot$. We present the stellar structure of the models during the hydrogen and helium burning phases. In the latter, we examine the impact on the nucleosynthesis. We find a broadening of the main sequence with more CBM, which is more in agreement with observations. Furthermore, during the core hydrogen burning phase there is a convergence of the convective boundary location due to CBM. The uncertainties of the intermediate convective zone remove this convergence. The behaviour of this convective zone strongly affects the surface evolution of the model, i.e. how fast it evolves redwards. The amount of CBM impacts the size of the convective cores and the nucleosynthesis, e.g. the 12C to 16O ratio and the weak s-process. Lastly, we determine the uncertainty that the range of parameter values investigated introduces and we find differences of up to $70{{\ \rm per\ cent}}$ for the core masses and the total mass of the star.more » « less
-
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
-
more » « less
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
Abstract One-dimensional stellar evolution calculations produce uncertain predictions for quantities like the age, core mass, core compactness, and nucleosynthetic yields; a key source of uncertainty is the modeling of interfaces between regions that are convectively stable and those that are not. Theoretical and numerical work has demonstrated that there should be numerous processes adjacent to the convective boundary that induce chemical and angular momentum transport, as well as modify the thermal structure of the star. One such process is called convective penetration, wherein vigorous convection extends beyond the nominal convective boundary and alters both the composition and thermal structure. In this work, we incorporate the process of convective penetration in stellar evolution calculations using the stellar evolution software instrument
mesa . We implement convective penetration according to the description presented by Anders et al. to to calculate a grid of models from the pre-main sequence to helium core depletion. The extent of the convective penetration zone is self-consistently calculated at each time step without introducing new free parameters. We find both a substantial penetration zone in all models with a convective core and observable differences to global stellar properties such as the luminosity and radius. We present how the predicted radial extent of the penetration zone scales with the total stellar mass, age, and metallicity of the star. We discuss our results in the context of existing numerical and observational studies.
