Search for: All records

Abstract We present 3D hydrodynamical simulations of core convection with a stably stratified envelope of a 25M_⊙star in the early phase of the main sequence. We use the explicit gas-dynamics codePPMstar, which tracks two fluids and includes radiation pressure and radiative diffusion. Multiple series of simulations with different luminosities and radiative thermal conductivities are presented. The entrainment rate at the convective boundary, internal gravity waves in and above the boundary region, and the approach to dynamical equilibrium shortly after a few convective turnovers are investigated. We perform very long simulations on 896³grids accelerated by luminosity boost factors of 1000, 3162 and 10,000. In these simulations, the growing penetrative convection reduces the initially unrealistically large entrainment. This reduction is enabled by a spatial separation that develops between the entropy gradient and the composition gradient. The convective boundary moves outward much more slowly at the end of these simulations. Finally, we present a 1D method to predict the extent and character of penetrative convection beyond the Schwarzschild boundary. The 1D model is based on a spherically averaged reduced entropy equation that takes the turbulent dissipation as input from the 3D hydrodynamic simulation and takes buoyancy and all other energy sources and sinks into account. This 1D method is intended to be ultimately deployed in 1D stellar evolution calculations and is based on the properties of penetrative convection in our simulations carried forward through the local thermal timescale. 

ABSTRACT A few per cent of red giants are enriched in lithium with $$A(\mathrm{Li}) \gt 1.5$$. Their evolutionary status has remained uncertain because these Li-rich giants can be placed both on the red giant branch (RGB) near the bump luminosity and in the red clump (RC) region. However, thanks to asteroseismology, it has been found that most of them are actually RC stars. Starting at the bump luminosity, RGB progenitors of the RC stars experience extra mixing in the radiative zone separating the H-burning shell from the convective envelope followed by a series of convective He-shell flashes at the RGB tip, known as the He-core flash. The He-core flash was proposed to cause fast extra mixing in the stars at the RGB tip that is needed for the Cameron–Fowler mechanism to produce Li. We propose that the RGB stars are getting enriched in Li by the RGB extra mixing that is getting enhanced and begins to produce Li, instead of destroying it, when the stars are approaching the RGB tip. After a discussion of several mechanisms of the RGB extra mixing, including the joint operation of rotation-driven meridional circulation and turbulent diffusion, the azimuthal magnetorotational instability (AMRI), thermohaline convection, buoyancy of magnetic flux tubes, and internal gravity waves, and based on results of (magneto-) hydrodynamics simulations and asteroseismology observations, we are inclined to conclude that it is the mechanism of the AMRI or magnetically enhanced thermohaline convection, that is most likely to support our hypothesis. 

ABSTRACT Recent photometric observations of massive stars have identified a low-frequency power excess which appears as stochastic low-frequency variability in light-curve observations. We present the oscillation properties of high-resolution hydrodynamic simulations of a $$25\,\,{\rm{M}_\odot }$$ star performed with the PPMstar code. The model star has a convective core mass of $$\approx 12\,\,{\rm{M}_\odot }$$ and approximately half of the envelope simulated. From this simulation, we extract light curves from several directions, average them over each hemisphere, and process them as if they were real photometric observations. We show how core convection excites waves with a similar frequency as the convective time-scale in addition to significant power across a forest of low and high angular degree l modes. We find that the coherence of these modes is relatively low as a result of their stochastic excitation by core convection, with lifetimes of the order of 10s of days. Thanks to the still significant power at higher l and this relatively low coherence, we find that integrating over a hemisphere produces a power spectrum that still contains measurable power up to the Brunt–Väisälä frequency. These power spectra extracted from the stable envelope are qualitatively similar to observations, with the same order of magnitude yet lower characteristic frequency. This work further shows the potential of long-duration, high-resolution hydrodynamic simulations for connecting asteroseismic observations to the structure and dynamics of core convection and the convective boundary. 

ABSTRACT The inner structure of core helium burning (CHeB) stars remains uncertain due to the yet unknown nature of mixing at the boundary of their cores. Large convective cores beyond a bare Schwarzschild model are favoured both from theoretical arguments and from asteroseismological constraints. However, the exact nature of this extra mixing, and in particular the possible presence of semiconvective layers, is still debated. In this work, we approach this problem through a new avenue by performing the first full-sphere 3D hydrodynamics simulations of the interiors of CHeB stars. We use the ppmstar explicit gas dynamics code to simulate the inner 0.45$$\, {\rm M}_{\odot }$$ of a 3 M⊙ CHeB star. Simulations are performed using different Cartesian grid resolutions (7683, 11523, and 17283) and heating rates. We use two different initial states, one based on mesas's predictive mixing scheme (which significantly extends the core beyond the Schwarzschild boundary) and one based on the convective premixing approach (which exhibits a semiconvective interface). The general behaviour of the flow in the convective core and in the stable envelope (where internal gravity waves are observed) is consistent with our recent simulations of core convection in massive main-sequence stars, and so are the various luminosity scaling relations. The semiconvective layers are dominated by strong internal gravity waves that do not produce measurable species mixing, but overshooting motions from the convective core gradually homogenize the semiconvective interface. This process can possibly completely erase the semiconvective layers, which would imply that CHeB stars do not harbour a semiconvection zone. 

ABSTRACT We performed 3D hydrodynamic simulations of the inner $$\approx 50{{\ \rm per\ cent}}$$ radial extent of a $$25\,\,\mathrm{\mathrm{M}_\odot }$$ star in the early phase of the main sequence and investigate core convection and internal gravity waves in the core-envelope boundary region. Simulations for different grid resolutions and driving luminosities establish scaling relations to constrain models of mixing for 1D applications. As in previous works, the turbulent mass entrainment rate extrapolated to nominal heating is unrealistically high ($$1.58\times 10^{-4}\,\,\mathrm{\mathrm{M}_\odot \, {\mathrm{yr}}^{-1}}$$), which is discussed in terms of the non-equilibrium response of the simulations to the initial stratification. We measure quantitatively the effect of mixing due to internal gravity waves excited by core convection interacting with the boundary in our simulations. The wave power spectral density as a function of frequency and wavelength agrees well with the GYRE eigenmode predictions based on the 1D spherically averaged radial profile. A diffusion coefficient profile that reproduces the spherically averaged abundance distribution evolution is determined for each simulation. Through a combination of eigenmode analysis and scaling relations it is shown that in the N2-peak region, mixing is due to internal gravity waves and follows the scaling relation DIGW-hydro ∝ L4/3 over a $$\gtrapprox 2\,\,\mathrm{\mathrm{dex}}$$ range of heating factors. Different extrapolations of the mixing efficiency down to nominal heating are discussed. If internal gravity wave mixing is due to thermally enhanced shear mixing, an upper limit is $$D_\mathrm{IGW}\lessapprox 2$$ to $$3\times 10^{4}\,\,\mathrm{cm^2\, s^{-1}}$$ at nominal heating in the N2-peak region above the convective core. 

ABSTRACT Supermassive stars are Population III stars with masses exceeding $$10^4\, {\rm M}_{\odot }$$ that could be the progenitors of the first supermassive black holes. Their interiors are in a regime where radiation pressure dominates the equation of state. In this work, we use the explicit gas dynamics code ppmstar to simulate the hydrogen-burning core of a $$10^4\, {\rm M}_{\odot }$$ supermassive main-sequence star. These are the first three-dimensional hydrodynamics simulations of core convection in supermassive stars. We perform a series of 10 simulations at different heating rates and on Cartesian grids with resolutions of 7683, 11523, and 17283. We examine different properties of the convective flow, including its large-scale morphology, its velocity spectrum, and its mixing properties. We conclude that the radiation pressure-dominated nature of the interior does not noticeably affect the behaviour of convection compared to the case of core convection in a massive main-sequence star where gas pressure dominates. Our simulations also offer support for the use of mixing-length theory in one-dimensional models of supermassive stars. 

ABSTRACT We present the first 3D hydrodynamics simulations of the excitation and propagation of internal gravity waves (IGWs) in the radiative interiors of low-mass stars on the red giant branch (RGB). We use the ppmstar explicit gas dynamics code to simulate a portion of the convective envelope and all the radiative zone down to the hydrogen-burning shell of a $$1.2\, {\rm M}_{\odot }$$ upper RGB star. We perform simulations for different grid resolutions (7683, 15363, and 28803), a range of driving luminosities, and two different stratifications (corresponding to the bump luminosity and the tip of the RGB). Our RGB tip simulations can be directly performed at the nominal luminosity, circumventing the need for extrapolations to lower luminosities. A rich, continuous spectrum of IGWs is observed, with a significant amount of total power contained at high wavenumbers. By following the time evolution of a passive dye in the stable layers, we find that IGW mixing in our simulations is weaker than predicted by a simple analytical prescription based on shear mixing and not efficient enough to explain the missing RGB extra mixing. However, we may be underestimating the efficiency of IGW mixing given that our simulations include a limited portion of the convective envelope. Quadrupling its radial extent compared to our fiducial set-up increases convective velocities by up to a factor 2 and IGW velocities by up to a factor 4. We also report the formation of a $$\sim 0.2\, H_P$$ penetration zone and evidence that IGWs are excited by plumes that overshoot into the stable layers. 

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. 

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

Abstract 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 C-shell 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_\rm {tot}$$. All but one of the simulations reach a quasi-stationary state with the entrainment rate and convective velocity proportional to $$L_\rm {tot}$$ and $$L_\rm {tot}^{1/3}$$, respectively. Carbon burning provides 14 – $$33\%$$ of the total luminosity, depending on the set of reactions considered. Equivalent simulations done on 7683 and 11523 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.