-
Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes
Authors:
R. Andrassy,
J. Higl,
H. Mao,
M. Mocák,
D. G. Vlaykov,
W. D. Arnett,
I. Baraffe,
S. W. Campbell,
T. Constantino,
P. V. F. Edelmann,
T. Goffrey,
T. Guillet,
F. Herwig,
R. Hirschi,
L. Horst,
G. Leidi,
C. Meakin,
J. Pratt,
F. Rizzuti,
F. K. Roepke,
P. Woodward
Abstract:
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 o…
▽ More
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 $\lesssim 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.
△ Less
Submitted 26 January, 2022; v1 submitted 1 November, 2021;
originally announced November 2021.
-
Multidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star
Authors:
L. Horst,
R. Hirschi,
P. V. F. Edelmann,
R. Andrassy,
F. K. Roepke
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 si…
▽ More
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.
△ Less
Submitted 5 July, 2021;
originally announced July 2021.
-
Well-balanced treatment of gravity in astrophysical fluid dynamics simulations at low Mach numbers
Authors:
P. V. F. Edelmann,
L. Horst,
J. P. Berberich,
R. Andrassy,
J. Higl,
G. Leidi,
C. Klingenberg,
F. K. Roepke
Abstract:
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. We demonstrate how discretization errors on grids of reasonable s…
▽ More
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. 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. 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. 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. 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.
△ Less
Submitted 8 July, 2021; v1 submitted 25 February, 2021;
originally announced February 2021.
-
Photometric detection of internal gravity waves in upper main-sequence stars. II. Combined TESS photometry and high-resolution spectroscopy
Authors:
D. M. Bowman,
S. Burssens,
S. Simón-Díaz,
P. V. F. Edelmann,
T. M. Rogers,
L. Horst,
F. K. Roepke,
C. Aerts
Abstract:
Context. Massive stars are predicted to excite internal gravity waves (IGWs) by turbulent core convection and from turbulent pressure fluctuations in their near-surface layers. These IGWs are extremely efficient at transporting angular momentum and chemical species within stellar interiors, but they remain largely unconstrained observationally. Aims. We aim to characterise the photometric detectio…
▽ More
Context. Massive stars are predicted to excite internal gravity waves (IGWs) by turbulent core convection and from turbulent pressure fluctuations in their near-surface layers. These IGWs are extremely efficient at transporting angular momentum and chemical species within stellar interiors, but they remain largely unconstrained observationally. Aims. We aim to characterise the photometric detection of IGWs across a large number of O and early-B stars in the Hertzsprung-Russell diagram, and explain the ubiquitous detection of stochastic variability in the photospheres of massive stars. Methods. We combined high-precision time-series photometry from the NASA Transiting Exoplanet Survey Satellite with high-resolution ground-based spectroscopy of 70 stars with spectral types O and B to probe the relationship between the photometric signatures of IGWs and parameters such as spectroscopic mass, luminosity, and macroturbulence. Results. A relationship is found between the location of a star in the spectroscopic Hertzsprung-Russell diagram and the amplitudes and frequencies of stochastic photometric variability in the light curves of massive stars. Furthermore, the properties of the stochastic variability are statistically correlated with macroturbulent velocity broadening in the spectral lines of massive stars. Conclusions. The common ensemble morphology for the stochastic low-frequency variability detected in space photometry and its relationship to macroturbulence is strong evidence for IGWs in massive stars, since these types of waves are unique in providing the dominant tangential velocity field required to explain the observed spectroscopy.
△ Less
Submitted 8 June, 2020; v1 submitted 4 June, 2020;
originally announced June 2020.
-
Fully compressible simulations of waves and core convection in main-sequence stars
Authors:
L. Horst,
P. V. F. Edelmann,
R. Andrassy,
F. K. Roepke,
D. M. Bowman,
C. Aerts,
R. P. Ratnasingam
Abstract:
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 rather high artificial viscosity and thermal diffusivity are needed for numerical stability. Direct comparison with obse…
▽ More
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 rather high artificial viscosity and thermal diffusivity are needed for numerical stability. Direct comparison with observations is thus limited. Aims. We explore the capabilities of our compressible multidimensional hydrodynamics code SLH 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\mathrm{M}_\odot$ 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 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 the coupling of gravity and pressure waves to be studied too.
△ Less
Submitted 14 September, 2020; v1 submitted 4 June, 2020;
originally announced June 2020.