Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userMultidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star
119
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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 demonstrate 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.
rate research
Read More
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 present the first 3-dimensional, fully compressible gas-dynamics simulations in $4pi$ geometry of He-shell flash convection with proton-rich fuel entrainment at the upper boundary. This work is motivated by the insufficiently understood observed consequences of the H-ingestion flash in post-AGB stars (Sakurais object) and metal-poor AGB stars. Our investigation is focused on the entrainment process at the top convection boundary and on the subsequent advection of H-rich material into deeper layers, and we therefore ignore the burning of the proton-rich fuel in this study. We find that, for our deep convection zone, coherent convective motions of near global scale appear to dominate the flow. At the top boundary convective shear flows are stable against Kelvin-Helmholtz instabilities. However, such shear instabilities are induced by the boundary-layer separation in large-scale, opposing flows. This links the global nature of thick shell convection with the entrainment process. We establish the quantitative dependence of the entrainment rate on grid resolution. With our numerical technique simulations with $1024^3$ cells or more are required to reach a numerical fidelity appropriate for this problem. However, only the result from the $1536^3$ simulation provides a clear indication that we approach convergence with regard to the entrainment rate. Our results demonstrate that our method, which is described in detail, can provide quantitative results related to entrainment and convective boundary mixing in deep stellar interior environments with veryvstiff convective boundaries. For the representative case we study in detail, we find an entrainment rate of $4.38 pm 1.48 times 10^{-13}M_odot mathrm{/s}$.
We find significant fluctuations of angular momentum within the convective helium shell of a pre-collapse massive star - a core-collapse supernova progenitor - which may facilitate the formation of accretion disks and jets that can explode the star. The convective flow in our model of an evolved M_ZAMS=15Msun star, computed with the sub-sonic hydrodynamic solver MAESTRO, contains entire shells with net angular momentum in different directions. This phenomenon may have important implications for the late evolutionary stages of massive stars, and for the dynamics of core-collapse.
Stars of sufficiently low mass are convective throughout their interiors, and so do not possess an internal boundary layer akin to the solar tachocline. Because that interface figures so prominently in many theories of the solar magnetic dynamo, a widespread expectation had been that fully convective stars would exhibit surface magnetic behavior very different from that realized in more massive stars. Here I describe how recent observations and theoretical models of dynamo action in low-mass stars are partly confirming, and partly confounding, this basic expectation. In particular, I present the results of 3--D MHD simulations of dynamo action by convection in rotating spherical shells that approximate the interiors of 0.3 solar-mass stars at a range of rotation rates. The simulated stars can establish latitudinal differential rotation at their surfaces which is solar-like at ``rapid rotation rates (defined within) and anti-solar at slower rotation rates; the differential rotation is greatly reduced by feedback from strong dynamo-generated magnetic fields in some parameter regimes. I argue that this ``flip in the sense of differential rotation may be observable in the near future. I also briefly describe how the strength and morphology of the magnetic fields varies with the rotation rate of the simulated star, and show that the maximum magnetic energies attained are compatible with simple scaling arguments.
(abridged) Context: The ratio of kinematic viscosity to thermal diffusivity, the Prandtl number, is much smaller than unity in stellar convection zones. Aims: To study the statistics of convective flows and energy transport as functions of the Prandtl number. Methods: Three-dimensional numerical simulations convection in Cartesian geometry are used. The convection zone (CZ) is embedded between two stably stratified layers. Statistics and transport properties of up- and downflows are studied separately. Results: The rms velocity increases with decreasing Prandtl number. At the same time the filling factor of downflows decreases and leads to stronger downflows at lower Prandtl numbers, and to a strong dependence of overshooting on the Prandtl number. Velocity power spectra do not show marked changes as a function of Prandtl number. At the highest Reynolds numbers the velocity power spectra are compatible with the Bolgiano-Obukhov $k^{-11/5}$ scaling. The horizontally averaged convected energy flux ($overline{F}_{rm conv}$) is independent of the Prandtl number within the CZ. However, the upflows (downflows) are the dominant contribution to the convected flux at low (high) Prandtl number. These results are similar to those from Rayleigh-Benard convection in the low Prandtl number regime where convection is vigorously turbulent but inefficient at transporting energy. Conclusions: The current results indicate a strong dependence of convective overshooting and energy flux on the Prandtl number. Numerical simulations of astrophysical convection often use Prandtl number of unity. The current results suggest that this can lead to misleading results and that the astrophysically relevant low Prandtl number regime is qualitatively different from the parameters regimes explored in typical simulations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
