This work investigates the properties of convection in stars with particular emphasis on entrainment across the upper convective boundary (CB). Idealised simulations of turbulent convection in the O-burning shell of a massive star are performed in 4π geometry on 768 3 and 1536 3 grids, driven by a representative heating rate. A heating series is also performed on the 768 3 grid. The 1536 3 simulation exhibits an entrainment rate at the upper CB of 1.33 × 10 −6 M s −1 . The 768 3 simulation with the same heating rate agrees within 17 per cent. The entrainment rate at the upper convective boundary is found to scale linearly with the driving luminosity and with the cube of the shear velocity at the upper boundary, while the radial RMS fluid velocity scales with the cube root of the driving luminosity, as expected. The mixing is analysed in a 1D diffusion framework, resulting in a simple model for CB mixing. The analysis confirms previous findings that limiting the MLT mixing length to the distance to the CB in 1D simulations better represents the spherically-averaged radial velocity profiles from the 3D simulations and provides an improved determination of the reference diffusion coefficient D 0 for the exponential diffusion CB mixing model in 1D. From the 3D simulation data we adopt as the convective boundary the location of the maximum gradient in the horizontal velocity component which has 2σ spatial fluctuations of ≈ 0.17H P . The exponentially decaying diffusion CB mixing model with f = 0.03 reproduces the spherically-averaged 3D abundance profiles.

We have modelled the multicycle evolution of rapidly-accreting CO white dwarfs (RAWDs) with stable H burning intermittent with strong He-shell flashes on their surfaces for 0.7 ≤ M RAWD /M ≤ 0.75 and [Fe/H] ranging from 0 to −2.6. We have also computed the i-process nucleosynthesis yields for these models. The i process occurs when convection driven by the He-shell flash ingests protons from the accreted H-rich surface layer, which results in maximum neutron densities N n,max ≈ 10 13 -10 15 cm −3 . The H-ingestion rate and the convective boundary mixing (CBM) parameter f top adopted in the one-dimensional nucleosynthesis and stellar evolution models are constrained through 3D hydrodynamic simulations. The mass ingestion rate and, for the first time, the scaling laws for the CBM parameter f top have been determined from 3D hydrodynamic simulations. We confirm our previous result that the highmetallicity RAWDs have a low mass retention efficiency (η < ∼ 10%). A new result is that RAWDs with [Fe/H] < ∼ − 2 have η > ∼ 20%, therefore their masses may reach the Chandrasekhar limit and they may eventually explode as SNeIa. This result and the good fits of the i-process yields from the metal-poor RAWDs to the observed chemical composition of the CEMP-r/s stars suggest that some of the present-day CEMP-r/s stars could be former distant members of triple systems, orbiting close binary systems with RAWDs that may have later exploded as SNeIa.

Context. Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and requires rather high artificial viscosity and thermal diffusivity for numerical stability. A direct comparison with observations is thus limited. Aims. We explore the capabilities of our compressible multidimensional Seven-League Hydro (SLH) code to simulate stellar oscillations. Methods. We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic 3 M⊙ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results. Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with the observed spectra of massive stars. Conclusion. The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows for the coupling of gravity and pressure waves in the outer convective envelopes of evolved stars to be studied in the future.

Interactions between convective shells in evolved massive stars have been linked to supernova impostors, to the production of the odd-Z elements Cl, K, and Sc, and they might also help generate the large-scale asphericities that are known to facilitate shock revival in supernova explosion models. We investigate the process of ingestion of Cshell material into a convective O-burning shell, including the hydrodynamic feedback from the nuclear burning of the ingested material. Our 3D hydrodynamic simulations span almost 3 dex in the total luminosity L tot . All but one of the simulations reach a quasi-stationary state with the entrainment rate and convective velocity proportional to L tot and L 1/3 tot , respectively. Carbon burning provides 14 -33% of the total luminosity, depending on the set of reactions considered. Equivalent simulations done on 768 3 and 1152 3 grids are in excellent quantitative agreement. The flow is dominated by a few large-scale convective cells. An instability leading to large-scale oscillations with Mach numbers in excess of 0.2 develops in an experimental run with the energy yield from C burning increased by a factor of 10. This run represents most closely the conditions expected in a violent O-C shell merger, which is a potential production site for odd-Z elements such as K and Sc and which may seed asymmetries in the supernova progenitor. 1D simulations may underestimate the energy generation from the burning of ingested material by as much as a factor two owing to their missing the effect of clumpiness of entrained material on the nuclear reaction rate.

We address the deficiency of odd-Z elements P, Cl, K and Sc in galactic chemical evolution models through an investigation of the nucleosynthesis of interacting convective O-and C shells in massive stars. 3D hydrodynamic simulations of O-shell convection with moderate C-ingestion rates show no dramatic deviation from spherical symmetry. We derive a spherically averaged diffusion coefficient for 1D nucleosynthesis simulations which show that such convective-reactive ingestion events can be a production site for P, Cl, K and Sc. An entrainment rate of 10 −3 M /s features overproduction factors OP s ≈ 7. Full O-C shell mergers in our 1D stellar evolution massive star models have overproduction factors OP m > 1dex but for such cases 3D hydrodynamic simulations suggest deviations from spherical symmetry. γ-process species can be produced with overproduction factors of OP m > 1dex, e.g. for 130,132 Ba. Using the uncertain prediction of the 15 M , Z = 0.02 massive star model (OP m ≈ 15) as representative for merger or entrainment convective-reactive events involving O-and C-burning shells, and assume that such events occur in more than 50% of all stars, our chemical evolution models reproduce the observed Galactic trends of the odd-Z elements.

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 to 10 −2 . Standard mixing-length theory (MLT) 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.

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 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.

Context. Accurate simulations of flows in stellar interiors are crucial to improving our understanding of stellar structure and evolution. Because the typically slow flows are merely tiny perturbations on top of a close balance between gravity and the pressure gradient, such simulations place heavy demands on numerical hydrodynamics schemes. Aims. We demonstrate how discretization errors on grids of reasonable size can lead to spurious flows orders of magnitude faster than the physical flow. Well-balanced numerical schemes can deal with this problem. Methods. Three such schemes were applied in the implicit, finite-volume Seven-League Hydro (SLH) code in combination with a low-Mach-number numerical flux function. We compare how the schemes perform in four numerical experiments addressing some of the challenges imposed by typical problems in stellar hydrodynamics. Results. We find that the α-β and deviation well-balancing methods can accurately maintain hydrostatic solutions provided that gravitational potential energy is included in the total energy balance. They accurately conserve minuscule entropy fluctuations advected in an isentropic stratification, which enables the methods to reproduce the expected scaling of convective flow speed with the heating rate. The deviation method also substantially increases accuracy of maintaining stationary orbital motions in a Keplerian disk on long timescales. The Cargo-LeRoux method fares substantially worse in our tests, although its simplicity may still offer some merits in certain situations. Conclusions. Overall, we find the well-balanced treatment of gravity in combination with low Mach number flux functions essential to reproducing correct physical solutions to challenging stellar slow-flow problems on affordable collocated grids.