More Like this

1. Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes

   https://doi.org/10.1051/0004-6361/202142557  

   Andrassy, R. ; Higl, J. ; Mao, H. ; Mocák, M. ; Vlaykov, D. G. ; Arnett, W. D. ; Baraffe, I. ; Campbell, S. W. ; Constantino, T. ; Edelmann, P. V. ; et al ( March 2022 , Astronomy & Astrophysics)

   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 ≲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 simulationsmore »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
2. Multidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star

   https://doi.org/10.1051/0004-6361/202140825  

   Horst, L. ; Hirschi, R. ; Edelmann, P. V. ; Andrássy, R. ; Röpke, F. K. ( September 2021 , Astronomy & Astrophysics)

   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 tomore »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.« less
3. 3D1D hydro-nucleosynthesis simulations – I. Advective–reactive post-processing method and its application to H ingestion into He-shell flash convection in rapidly accreting white dwarfs

   https://doi.org/10.1093/mnras/stab500  

   Stephens, David ; Herwig, Falk ; Woodward, Paul ; Denissenkov, Pavel ; Andrassy, Robert ; Mao, Huaqing ( April 2021 , Monthly Notices of the Royal Astronomical Society)

   ABSTRACT We present two mixing models for post-processing of 3D hydrodynamic simulations applied to convective–reactive i-process nucleosynthesis in a rapidly accreting white dwarf (RAWD) with [Fe/H] = −2.6, in which H is ingested into a convective He shell. A 1D advective two-stream model adopts physically motivated radial and horizontal mixing coefficients constrained by 3D hydrodynamic simulations. A simpler approach uses diffusion coefficients calculated from the same simulations. All 3D simulations include the energy feedback of the 12C(p, γ)13N reaction from the H entrainment. Global oscillations of shell H ingestion in two of the RAWD simulations cause bursts of entrainment of H and non-radial hydrodynamic feedback. With the same nuclear network as in the 3D simulations, the 1D advective two-stream model reproduces the rate and location of the H burning within the He shell closely matching the 3D simulation predictions, as well as qualitatively displaying the asymmetry of the XH profiles between the upstream and downstream. With a full i-process network the advective mixing model captures the difference in the n-capture nucleosynthesis in the upstream and downstream. For example, 89Kr and 90Kr with half-lives of $3.18\,\,\mathrm{\mathrm{min}}$ and $32.3\,\,\mathrm{\mathrm{s}}$ differ by a factor 2–10 in the two streams. In this particular applicationmore »the diffusion approach provides globally the same abundance distribution as the advective two-stream mixing model. The resulting i-process yields are in excellent agreement with observations of the exemplary CEMP-r/s star CS31062-050.« less
4. A Parameterization with a Constrained Potential Energy Conversion Rate of Vertical Mixing Due to Langmuir Turbulence

   https://doi.org/10.1175/JPO-D-18-0258.1  

   Reichl, Brandon G. ; Li, Qing ( November 2019 , Journal of Physical Oceanography)

   In this study we develop a new parameterization for turbulent mixing in the ocean surface boundary layer (OSBL), including the effect of Langmuir turbulence. This new parameterization builds on a recent study ( Reichl and Hallberg 2018 , hereafter RH18 ), which predicts the available energy for turbulent mixing against stable stratification driven by shear and convective turbulence. To investigate the role of Langmuir turbulence in the framework of RH18 , we utilize data from a suite of previously published large-eddy simulation (LES) experiments ( Li and Fox-Kemper 2017 , hereafter LF17 ) with and without Langmuir turbulence under different idealized forcing conditions. We find that the parameterization of RH18 is able to reproduce the mixing simulated by the LES in the non-Langmuir cases, but not the Langmuir cases. We therefore investigate the enhancement of the integrated vertical buoyancy flux within the entrainment layer in the presence of Langmuir turbulence using the LES data. An additional factor is introduced in the RH18 framework to capture the enhanced mixing due to Langmuir turbulence. This additional factor depends on the surface-layer averaged Langmuir number with a reduction in the presence of destabilizing surface buoyancy fluxes. It is demonstrated that including this factormore »within the RH18 OSBL turbulent mixing parameterization framework captures the simulated effect of Langmuir turbulence in the LES, which can be used for simulating the effect of Langmuir turbulence in climate simulations. This new parameterization is compared to the KPP-based Langmuir entrainment parameterization introduced by LF17 , and differences are explored in detail.

   « less
5. An Atlas of Convection in Main-sequence Stars

   https://doi.org/10.3847/1538-4365/ac7cee  

   Jermyn, Adam S. ; Anders, Evan H. ; Lecoanet, Daniel ; Cantiello, Matteo ( August 2022 , The Astrophysical Journal Supplement Series)

   Convection is ubiquitous in stars and occurs under many different conditions. Here we explore convection in main-sequence stars through two lenses: dimensionless parameters arising from stellar structure and parameters that emerge from the application of mixing length theory. We first define each quantity in terms familiar to both the 1D stellar evolution community and the hydrodynamics community. We then explore the variation of these quantities across different convection zones, different masses, and different stages of main-sequence evolution. We find immense diversity across stellar convection zones. Convection occurs in thin shells, deep envelopes, and nearly spherical cores; it can be efficient or inefficient, rotationally constrained or not, transsonic or deeply subsonic. This atlas serves as a guide for future theoretical and observational investigations by indicating which regimes of convection are active in a given star, and by describing appropriate model assumptions for numerical simulations.