اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدThree-Dimensional Hydrodynamic Simulations of Convective Nuclear Burning In Massive Stars Near Iron Core Collapse
115
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
Non-spherical structure in massive stars at the point of iron core collapse can have a qualitative impact on the properties of the ensuing core-collapse supernova explosions and the multi-messenger signals they produce. Strong perturbations can aid successful explosions by strengthening turbulence in the post-shock region. Here, we report on a set of $4pi$ 3D hydrodynamic simulations of O- and Si-shell burning in massive star models of varied initial masses using MESA and the FLASH simulation framework. We evolve four separate 3D models for roughly the final ten minutes prior to, and including, iron core collapse. We consider initial 1D MESA models with masses of 14-, 20-, and 25 $M_{odot}$ to survey a range of O/Si shell density and compositional configurations. We characterize the convective shells in our 3D models and compare them to the corresponding 1D models. In general, we find that the angle-average convective speeds in our 3D simulations near collapse are three to four times larger than the convective speeds predicted by MESA at the same epoch for our chosen mixing length parameter of $alpha_{rm{MLT}}=1.5$. In three of our simulations, we observe significant power in the spherical harmonic decomposition of the radial velocity field at harmonic indices of $ell=1-3$ near collapse. Our results suggest that large-scale modes are common in massive stars near collapse and should be considered a key aspect of pre-supernova progenitor models.
قيم البحث
اقرأ أيضاً
We perform two- (2D) and three-dimensional (3D) hydrodynamics simulations of convective oxygen shell-burning that takes place deep inside a massive progenitor star of a core-collapse supernova. Using one dimensional (1D) stellar evolution code, we fi
rst calculate the evolution of massive stars with an initial mass of 9-40 $M_odot$. Four different overshoot parameters are applied, and CO core mass trend similar to previous works is obtained in the 1D models. Selecting eleven 1D models that have a silicon and oxygen coexisting layer, we perform 2D hydrodynamics simulations of the evolution $sim$100 s until the onset of core-collapse. We find that convection with large-scale eddies and the turbulent Mach number $sim$0.1 is obtained in the models having a Si/O layer with a scale of 10$^8$ cm, whereas most models that have an extended O/Si layer up to a few $times 10^9$ cm exhibit lower turbulent velocity. Our results indicate that the supernova progenitors that possess a thick Si/O layer could provide a preferable condition for perturbation-aided explosions. We perform 3D simulation of a 25 $M_odot$ model, which exhibits large-scale convection in the 2D models. The 3D model develops large ($ell = 2$) convection similar to the 2D model, however, the turbulent velocity is lower. By estimating the neutrino emission properties of the 3D model, we point out that a time modulation of the event rates, if observed in KamLAND and Hyper-Kamiokande, would provide an important information about structural changes in the presupernova convective layer.
It has been suggested based on analytic theory that even in non-rotating supernova progenitors stochastic spin-up by internal gravity waves (IGWs) during the late burning stages can impart enough angular momentum to the core to result in neutron star
birth spin periods below 100ms, and a relatively firm upper limit of 500ms for the spin period. We here investigate this process using a 3D simulation of oxygen shell burning in a $3M_odot$ He star. Our model indicates that stochastic spin-up by IGWs is less efficient than previously thought. We find that the stochastic angular momentum flux carried by waves excited at the shell boundary is significantly smaller for a given convective luminosity and turnover time than would be expected from simple dimensional analysis. This can be explained by noting that the waves launched by overshooting convective plumes contain modes of opposite angular wave number with similar amplitudes, so that the net angular momentum of excited wave packets almost cancels. We find that the wave-mediated angular momentum flux from the oxygen shell follows a random walk, but again dimensional analysis overestimates the random walk amplitudes since the correlation time is only a fraction of the convective turnover time. Extrapolating our findings over the entire life time of the last burning stages prior to collapse, we predict that the core angular momentum from stochastic spin-up would translate into long birth spin periods of several seconds for low-mass progenitors and no less than 100ms even for high-mass progenitors.
Context. Multidimensional hydrodynamic simulations of convection in stellar interiors are numerically challenging, especially for flows at low Mach numbers. Methods. We explore the benefits of using a low-Mach hydrodynamic flux solver and demonstra
te its usability for simulations in the astrophysical context. The time-implicit Seven-League Hydro (SLH) code was used to perform multidimensional simulations of convective helium shell burning based on a 25 M$_odot$ star model. The results obtained with the low-Mach AUSM$^{+}$-up solver were compared to results when using its non low-Mach variant AUSM$_mathrm{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 times 90^2$ to $810 times 540^2$ cells and the nuclear energy release was boosted by factors of $3 times 10^3$, $1 times 10^4$, and $3 times 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^{-2}$ to $10^{-3}$. Standard mixing-length theory (MLT) predicts convective velocities of about $1.6 times 10^{-4}$ if no boosting is applied. 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$^{+}_{mathrm{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.
(Abridged) Stars more massive than $20-30M_{odot}$ are so luminous that the radiation force on the cooler, more opaque outer layers can balance or exceed the force of gravity. These near or super-Eddington outer envelopes represent a long standing ch
allenge for calculating the evolution of massive stars in one dimension, a situation that limits our understanding of the stellar progenitors of some of the most exciting and energetic explosions in the universe. In particular, the proximity to the Eddington limit has been the suspected cause for the variability, large mass loss rate and giant eruptions of an enigmatic class of massive stars: the luminous blue variables (LBVs). When in quiescence, LBVs are usually found on the hot ($T_{eff} approx 2 - 4 times 10^4$ K) S Dor instability strip. While in outburst, most LBVs stay on the cold S Dor instability strip with a $T_{eff} approx 9000$ K. Here we show that physically realistic three dimensional global radiation hydrodynamic simulations of radiation dominated massive stars with the largest supercomputers in the world naturally reproduce many observed properties of LBVs, specifically their location in the Hertzsprung-Russell (HR) diagram and their episodic mass loss with rates of $10^{-7}-10^{-5} M_{odot}/yr$. The helium opacity peak is found to play an important role to determine these properties, which is not realized in the traditional one dimensional models of massive stars. The simulations also predict that convection causes irregular envelope oscillations and 10-30% brightness variations on a typical timescale of a few days. The variability is more prominent in our models that are on the cool part of the S Dor instability. These calculations pave the way to a quantitative understanding of the structure, stability and the dominant mode of mass loss of massive stars.
We present 3D hydrodynamics simulations of shell burning in two progenitors with zero-age main-sequence masses of 22 and 27 $M_{odot}$ for $sim$65 and 200 s up to the onset of gravitational collapse, respectively. The 22 and 27 $M_{odot}$ stars are s
elected from a suite of 1D progenitors. The former and the latter have an extended Si- and O-rich layer with a width of $sim$10$^9$ cm and $sim$5$times 10^9$ cm, respectively. Our 3D results show that turbulent mixing occurs in both of the progenitors with the angle-averaged turbulent Mach number exceeding $sim$0.1 at the maximum. We observe that an episodic burning of O and Ne, which takes place underneath the convection bases, enhances the turbulent mixing in the 22 and 27 $M_odot$ models, respectively. The distribution of nucleosynthetic yields is significantly different from that in 1D simulations, namely, in 3D more homogeneous and inhomogeneous in the radial and angular direction, respectively. By performing a spectrum analysis, we investigate the growth of turbulence and its role of material mixing in the convective layers. We also present a scalar spherical harmonics mode analysis of the turbulent Mach number. This analytical formula would be helpful for supernova modelers to implement the precollapse perturbations in core-collapse supernova simulations. Based on the results, we discuss implications for the possible onset of the perturbation-aided neutrino-driven supernova explosion.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
