-
3D simulations of convective shell Neon-burning in a massive star
Authors:
C. Georgy,
F. Rizzuti,
R. Hirschi,
V. Varma,
W. D. Arnett,
C. Meakin,
M. Mocak,
A. StJ. Murphy,
T. Rauscher
Abstract:
The treatment of convection remains a major weakness in the modelling of stellar evolution with one-dimensional (1D) codes. The ever increasing computing power makes now possible to simulate in 3D part of a star for a fraction of its life, allowing us to study the full complexity of convective zones with hydrodynamics codes. Here, we performed state-of-the-art hydrodynamics simulations of turbulen…
▽ More
The treatment of convection remains a major weakness in the modelling of stellar evolution with one-dimensional (1D) codes. The ever increasing computing power makes now possible to simulate in 3D part of a star for a fraction of its life, allowing us to study the full complexity of convective zones with hydrodynamics codes. Here, we performed state-of-the-art hydrodynamics simulations of turbulence in a neon-burning convective zone, during the late stage of the life of a massive star. We produced a set of simulations varying the resolution of the computing domain (from 1283 to 10243 cells) and the efficiency of the nuclear reactions (by boosting the energy generation rate from nominal to a factor of 1000). We analysed our results by the mean of Fourier transform of the velocity field, and mean-field decomposition of the various transport equations. Our results are in line with previous studies, showing that the behaviour of the bulk of the convective zone is already well captured at a relatively low resolution (2563), while the details of the convective boundaries require higher resolutions. The different boosting factors used show how various quantities (velocity, buoyancy, abundances, abundance variances) depend on the energy generation rate. We found that for low boosting factors, convective zones are well mixed, validating the approach usually used in 1D stellar evolution codes. However, when nuclear burning and turbulent transport occur on the same timescale, a more sophisticated treatment would be needed. This is typically the case when shell mergers occur.
△ Less
Submitted 31 May, 2024;
originally announced May 2024.
-
3D stellar evolution: hydrodynamic simulations of a complete burning phase in a massive star
Authors:
F. Rizzuti,
R. Hirschi,
W. D. Arnett,
C. Georgy,
C. Meakin,
A. StJ. Murphy,
T. Rauscher,
V. Varma
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 multi-dimensional 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…
▽ More
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 multi-dimensional 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$_\odot$ 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.
△ Less
Submitted 23 May, 2023;
originally announced May 2023.
-
Realistic 3D hydrodynamics simulations find significant turbulent entrainment in massive stars
Authors:
F. Rizzuti,
R. Hirschi,
C. Georgy,
W. D. Arnett,
C. Meakin,
A. StJ. Murphy
Abstract:
Our understanding of stellar structure and evolution coming from one-dimensional (1D) stellar models is limited by uncertainties related to multi-dimensional processes taking place in stellar interiors. 1D models, however, can now be tested and improved with the help of detailed three-dimensional (3D) hydrodynamics models, which can reproduce complex multi-dimensional processes over short timescal…
▽ More
Our understanding of stellar structure and evolution coming from one-dimensional (1D) stellar models is limited by uncertainties related to multi-dimensional processes taking place in stellar interiors. 1D models, however, can now be tested and improved with the help of detailed three-dimensional (3D) hydrodynamics models, which can reproduce complex multi-dimensional processes over short timescales, thanks to the recent advances in computing resources. Among these processes, turbulent entrainment leading to mixing across convective boundaries is one of the least understood and most impactful. Here we present the results from a set of hydrodynamics simulations of the neon-burning shell in a massive star, and interpret them in the framework of the turbulent entrainment law from geophysics. Our simulations differ from previous studies in their unprecedented degree of realism in reproducing the stellar environment. Importantly, the strong entrainment found in the simulations highlights the major flaws of the current implementation of convective boundary mixing in 1D stellar models. This study therefore calls for major revisions of how convective boundaries are modelled in 1D, and in particular the implementation of entrainment in these models. This will have important implications for supernova theory, nucleosynthesis, neutron stars and black holes physics.
△ Less
Submitted 16 August, 2022; v1 submitted 7 July, 2022;
originally announced July 2022.
-
Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes
Authors:
R. Andrassy,
J. Higl,
H. Mao,
M. Mocák,
D. G. Vlaykov,
W. D. Arnett,
I. Baraffe,
S. W. Campbell,
T. Constantino,
P. V. F. Edelmann,
T. Goffrey,
T. Guillet,
F. Herwig,
R. Hirschi,
L. Horst,
G. Leidi,
C. Meakin,
J. Pratt,
F. Rizzuti,
F. K. Roepke,
P. Woodward
Abstract:
Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility o…
▽ More
Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within $\lesssim 3σ$ of the mean of all simulations on a given grid ($128^3$ and $256^3$ grid cells), where $σ$ describes the statistical variation due to the flow's time dependence. They also agree well with a $512^3$ reference run. The $128^3$ and $256^3$ simulations agree within $9\%$ and $4\%$, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study.
△ Less
Submitted 26 January, 2022; v1 submitted 1 November, 2021;
originally announced November 2021.
-
Constraints on stellar rotation from the evolution of Sr and Ba in the Galactic halo
Authors:
F. Rizzuti,
G. Cescutti,
F. Matteucci,
A. Chieffi,
R. Hirschi,
M. Limongi,
A. Saro
Abstract:
Recent studies show that the chemical evolution of Sr and Ba in the Galaxy can be explained if different production sites, hosting r- and s-processes, are taken into account. However, the question of unambiguously identifying these sites is still unsolved. Massive stars are shown to play an important role in the production of s-material if rotation is considered. In this work, we study in detail t…
▽ More
Recent studies show that the chemical evolution of Sr and Ba in the Galaxy can be explained if different production sites, hosting r- and s-processes, are taken into account. However, the question of unambiguously identifying these sites is still unsolved. Massive stars are shown to play an important role in the production of s-material if rotation is considered. In this work, we study in detail the contribution of rotating massive stars to the production of Sr and Ba, in order to explain their chemical evolution, but also to constrain the rotational behaviour of massive stars. A stochastic chemical evolution model was employed to reproduce the enrichment of the Galactic halo. We developed new methods for model-data comparison which help to objectively compare the stochastic results to the observations. We employed these methods to estimate the value of free parameters which describe the rotation of massive stars, assumed to be dependent on the stellar metallicity. We constrain the parameters using the observations for Sr and Ba. Employing these parameters for rotating massive stars in our stochastic model, we are able to correctly reproduce the chemical evolution of Sr and Ba, but also Y, Zr and La. The data supports a decrease of both the mean rotational velocities and their dispersion with increasing metallicity. Our results show that a metallicity-dependent rotation is a necessary assumption to explain the s-process in massive stars. Our novel methods of model-data comparison represent a promising tool for future galactic chemical evolution studies.
△ Less
Submitted 13 January, 2021;
originally announced January 2021.
-
The contribution from rotating massive stars to the enrichment in Sr and Ba of the Milky Way
Authors:
Federico Rizzuti,
Gabriele Cescutti,
Francesca Matteucci,
Alessandro Chieffi,
Raphael Hirschi,
Marco Limongi
Abstract:
Most neutron capture elements have a double production by r- and s-processes, but the question of production sites is complex and still open. Recent studies show that including stellar rotation can have a deep impact on nucleosynthesis. We studied the evolution of Sr and Ba in the Milky Way. A chemical evolution model was employed to reproduce the Galactic enrichment. We tested two different nucle…
▽ More
Most neutron capture elements have a double production by r- and s-processes, but the question of production sites is complex and still open. Recent studies show that including stellar rotation can have a deep impact on nucleosynthesis. We studied the evolution of Sr and Ba in the Milky Way. A chemical evolution model was employed to reproduce the Galactic enrichment. We tested two different nucleosynthesis prescriptions for s-process in massive stars, adopted from the Geneva group and the Rome group. Rotation was taken into account, studying the effects of stars without rotation or rotating with different velocities. We also tested different production sites for the r-process: magneto rotational driven supernovae and neutron star mergers. The evolution of the abundances of Sr and Ba is well reproduced. The comparison with the most recent observations shows that stellar rotation is a good assumption, but excessive velocities result in overproduction of these elements. In particular, the predicted evolution of the [Sr/Ba] ratio at low metallicity does not explain the data at best if rotation is not included. Adopting different rotational velocities for different stellar mass and metallicity better explains the observed trends. Despite the differences between the two sets of adopted stellar models, both show a better agreement with the data assuming an increase of rotational velocity toward low metallicity. Assuming different r-process sources does not alter this conclusion.
△ Less
Submitted 10 September, 2019;
originally announced September 2019.