-
Effect of Rotation on Wave Mixing in Intermediate Mass Stars
Authors:
A. Varghese,
R. P. Ratnasingam,
R. Vanon,
P. V. F. Edelmann,
S. Mathis,
T. M. Rogers
Abstract:
Internal gravity waves (IGWs) are likely to cause mixing in stellar interiors. Studies show that the mixing by these waves changes drastically across age and mass (Varghese et al. 2023, arXiv:2211.06432). Here, we study the effect of rotation on this wave mixing by considering a 7 M$_{\odot}$ model at ZAMS and midMS. We compare the mixing profiles at a range of rotation rates ($1\times 10^{-5}$,…
▽ More
Internal gravity waves (IGWs) are likely to cause mixing in stellar interiors. Studies show that the mixing by these waves changes drastically across age and mass (Varghese et al. 2023, arXiv:2211.06432). Here, we study the effect of rotation on this wave mixing by considering a 7 M$_{\odot}$ model at ZAMS and midMS. We compare the mixing profiles at a range of rotation rates ($1\times 10^{-5}$, $2\times 10^{-5}$, $3\times 10^{-5}$, $4\times 10^{-5}$ and $1\times 10^{-4}$ rad.s$^{-1}$) and observe that the mixing decreases with decreasing Rossby number. This can be attributed to the effect of rotation on convection which influences the amplitude with which the waves are excited near the convective-radiative interface.
△ Less
Submitted 6 June, 2024;
originally announced June 2024.
-
Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows
Authors:
G. Leidi,
R. Andrassy,
W. Barsukow,
J. Higl,
P. V. F. Edelmann,
F. K. Röpke
Abstract:
High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so t…
▽ More
High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. In this work, we use our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We consider Mach numbers in the range from $10^{-3}$ to $10^{-1}$ in a well-posed, 2D, Kelvin--Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. We find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a "low-dissipation" Riemann solver and a sextic reconstruction scheme, (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same, (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate, (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow, (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.
△ Less
Submitted 26 February, 2024;
originally announced February 2024.
-
Turbulent dynamo action and its effects on the mixing at the convective boundary of an idealized oxygen-burning shell
Authors:
G. Leidi,
R. Andrassy,
J. Higl,
P. V. F. Edelmann,
F. K. Röpke
Abstract:
Convection is one of the most important mixing processes in stellar interiors. Hydrodynamic mass entrainment can bring fresh fuel from neighboring stable layers into a convection zone, modifying the structure and evolution of the star. Under some conditions, strong magnetic fields can be sustained by the action of a turbulent dynamo, adding another layer of complexity and possibly altering the dyn…
▽ More
Convection is one of the most important mixing processes in stellar interiors. Hydrodynamic mass entrainment can bring fresh fuel from neighboring stable layers into a convection zone, modifying the structure and evolution of the star. Under some conditions, strong magnetic fields can be sustained by the action of a turbulent dynamo, adding another layer of complexity and possibly altering the dynamics in the convection zone and at its boundaries. In this study, we used our fully compressible Seven-League Hydro code to run detailed and highly resolved three-dimensional magnetohydrodynamic simulations of turbulent convection, dynamo amplification, and convective boundary mixing in a simplified setup whose stratification is similar to that of an oxygen-burning shell in a star with an initial mass of $25\ M_\odot$. We find that the random stretching of magnetic field lines by fluid motions in the inertial range of the turbulent spectrum (i.e., a small-scale dynamo) naturally amplifies the seed field by several orders of magnitude in a few convective turnover timescales. During the subsequent saturated regime, the magnetic-to-kinetic energy ratio inside the convective shell reaches values as high as $0.33$, and the average magnetic field strength is ${\sim}10^{10}\,\mathrm{G}$. Such strong fields efficiently suppress shear instabilities, which feed the turbulent cascade of kinetic energy, on a wide range of spatial scales. The resulting convective flows are characterized by thread-like structures that extend over a large fraction of the convective shell. The reduced flow speeds and the presence of magnetic fields with strengths up to $60\%$ of the equipartition value at the upper convective boundary diminish the rate of mass entrainment from the stable layer by ${\approx}\,20\%$ as compared to the purely hydrodynamic case.
△ Less
Submitted 29 September, 2023;
originally announced September 2023.
-
Multidimensional Simulations of Core Convection
Authors:
Daniel Lecoanet,
Philipp V. F. Edelmann
Abstract:
The cores of main sequence intermediate- and high-mass stars are convective. Mixing at the radiative-convective boundary, waves excited by the convection, and magnetic fields generated by convective dynamos all influence the main sequence and post-main sequence evolution of these stars. These effects must be understood to accurately model the structure and evolution of intermediate- and high-mass…
▽ More
The cores of main sequence intermediate- and high-mass stars are convective. Mixing at the radiative-convective boundary, waves excited by the convection, and magnetic fields generated by convective dynamos all influence the main sequence and post-main sequence evolution of these stars. These effects must be understood to accurately model the structure and evolution of intermediate- and high-mass stars. Unfortunately, there are many challenges in simulating core convection due to the wide range of temporal and spatial scales, as well as many important physics effects. In this review, we describe the latest numerical strategies to address these challenges. We then describe the latest state-of-the-art simulations of core convection, summarizing their main findings. These simulations have led to important insights into many of the processes associated with core convection. Two outstanding problems with multidimensional simulations are, 1. it is not always straightforward to extrapolate from simulation parameters to the parameters of real stars; and 2. simulations using different methods sometimes appear to arrive at contradictory results. To address these issues, next generation simulations of core convection must address how their results depend on stellar luminosity, dimensionality, and turbulence intensity. Furthermore, code comparison projects will be essential to establish robust parameterizations that will become the new standard in stellar modeling.
△ Less
Submitted 28 July, 2023;
originally announced July 2023.
-
Three-Dimensional Simulations of Massive Stars: II. Age Dependence
Authors:
R. Vanon,
P. V. F. Edelmann,
R. P. Ratnasingam,
A. Varghese,
T. M. Rogers
Abstract:
We present 3D full star simulations, reaching up to 90% of the total stellar radius, for three $7M_\odot$ stars of different ages (ZAMS, midMS and TAMS). A comparison with several theoretical prescriptions shows the generation spectra for all three ages are dominated by convective plumes. Two distinct overshooting layers are observed, with most plumes stopped within the layer situated directly abo…
▽ More
We present 3D full star simulations, reaching up to 90% of the total stellar radius, for three $7M_\odot$ stars of different ages (ZAMS, midMS and TAMS). A comparison with several theoretical prescriptions shows the generation spectra for all three ages are dominated by convective plumes. Two distinct overshooting layers are observed, with most plumes stopped within the layer situated directly above the convective boundary (CB); overshooting to the second, deeper layer becomes increasingly more infrequent with stellar age. Internal gravity wave (IGW) propagation is significantly impacted in the midMS and TAMS models as a result of some IGWs getting trapped within their Brunt-Väisälä frequency spikes. A fundamental change in the wave structure across radius is also observed, driven by the effect of density stratification on IGW propagation causing waves to become evanescent within the radiative zone, with older stars being affected more strongly. We find that the steepness of the frequency spectrum at the surface increases from ZAMS to the older models, with older stars also showing more modes in their spectra.
△ Less
Submitted 27 July, 2023;
originally announced July 2023.
-
Towards a self-consistent model of the convective core boundary in upper-main-sequence stars
Authors:
R. Andrassy,
G. Leidi,
J. Higl,
P. V. F. Edelmann,
F. R. N. Schneider,
F. K. Roepke
Abstract:
There is strong observational evidence that convective cores of intermediate-mass and massive main-sequence stars are substantially larger than standard stellar-evolution models predict. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-time-scale process th…
▽ More
There is strong observational evidence that convective cores of intermediate-mass and massive main-sequence stars are substantially larger than standard stellar-evolution models predict. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-time-scale process that is likely to be particularly relevant during the slow evolution on the main sequence. We use our low-Mach-number Seven-League Hydro (SLH) code to study this process in 2.5D and 3D geometries. Starting with a chemically homogeneous model of a $15$ M$_\odot$ zero-age main-sequence star, we construct a series of simulations with the luminosity increased and opacity decreased by the same factor ranging from $10^3$ to $10^6$. After reaching thermal equilibrium, all of our models show a clear penetration layer. Its thickness becomes statistically constant in time and it is shown to converge upon grid refinement. As the luminosity is decreased, the penetration layer becomes nearly adiabatic with a steep transition to a radiative stratification. This structure corresponds to the adiabatic ,,step overshoot'' model often employed in stellar-evolution calculations. The thickness of the penetration layer slowly decreases with decreasing luminosity. Depending on how we extrapolate our 3D data to the actual luminosity of the initial stellar model, we obtain penetration distances ranging from $0.09$ to $0.44$ pressure scale heights, which are broadly compatible with observations.
△ Less
Submitted 8 July, 2023;
originally announced July 2023.
-
Internal Gravity Waves in Massive Stars II: Frequency Analysis Across Stellar Mass
Authors:
R. P. Ratnasingam,
T. M. Rogers,
S. Chowdhury,
G. Handler,
R. Vanon,
A. Varghese,
P. V. F. Edelmann
Abstract:
Stars with masses above 1.6 solar masses generally possess convective cores and radiative envelopes, which allows the propagation of outward-travelling internal gravity waves. We have studied the generation and propagation of IGWs in such stars using two-dimensional, fully nonlinear hydrodynamical simulations with realistic stellar reference states from the one-dimensional stellar evolution code,…
▽ More
Stars with masses above 1.6 solar masses generally possess convective cores and radiative envelopes, which allows the propagation of outward-travelling internal gravity waves. We have studied the generation and propagation of IGWs in such stars using two-dimensional, fully nonlinear hydrodynamical simulations with realistic stellar reference states from the one-dimensional stellar evolution code, Modules for Stellar Astrophysics. Compared to previous similar works, this study utilises radius-dependent thermal diffusivity profiles for 5 different stellar masses at the middle of main sequence: 3 - 13 solar masses. From the simulations, we find that the surface perturbations are larger for higher masses, but no noticeable trends are observed for the frequency slopes with different stellar masses. The slopes are also similar to the results from previous works. We compare our simulation results with stellar photometric data from a recent survey and find that for frequency intervals above 8 microHz, there is a good agreement between the temperature frequency slopes from the simulations and the surface brightness variations of these observed stars, indicating that the brightness variations are caused by core-generated IGWs.
△ Less
Submitted 10 May, 2023;
originally announced May 2023.
-
Chemical Mixing Induced by Internal Gravity Waves in Intermediate Mass Stars
Authors:
A. Varghese,
R. P. Ratnasingam,
R. Vanon,
P. V. F. Edelmann,
T. M. Rogers
Abstract:
Internal gravity waves (IGWs) can cause mixing in the radiative interiors of stars. We study this mixing by introducing tracer particles into two - dimensional (2D) hydrodynamic simulations. Following the work of Rogers & McElwaine (2017), arXiv:1709.04920, we extend our study to different masses (3 M$_{\odot}$, 7 M$_{\odot}$ and 20 M$_{\odot}$) and ages (ZAMS, midMS and TAMS). The diffusion profi…
▽ More
Internal gravity waves (IGWs) can cause mixing in the radiative interiors of stars. We study this mixing by introducing tracer particles into two - dimensional (2D) hydrodynamic simulations. Following the work of Rogers & McElwaine (2017), arXiv:1709.04920, we extend our study to different masses (3 M$_{\odot}$, 7 M$_{\odot}$ and 20 M$_{\odot}$) and ages (ZAMS, midMS and TAMS). The diffusion profiles of these models are influenced by various parameters such as the Brunt-Väisälä frequency, density, thermal damping, the geometric effect and the frequencies of waves contributing to these mixing profiles. We find that the mixing profile changes dramatically across age. In younger stars, we noted that the diffusion coefficient increases towards the surface, whereas in older stars the initial increase in the diffusion profile is followed by a decreasing trend. We also find that mixing is stronger in more massive stars. Hence, future stellar evolution models should include this variation. In order to aid the inclusion of this mixing in one-dimensional (1D) stellar evolution models, we determine the dominant waves contributing to these mixing profiles and present a prescription that can be included in 1D models.
△ Less
Submitted 11 November, 2022;
originally announced November 2022.
-
A finite-volume scheme for modeling compressible magnetohydrodynamic flows at low Mach numbers in stellar interiors
Authors:
G. Leidi,
C. Birke,
R. Andrassy,
J. Higl,
P. V. F. Edelmann,
G. Wiest,
C. Klingenberg,
F. K. Röpke
Abstract:
Fully compressible magnetohydrodynamic (MHD) simulations are a fundamental tool for investigating the role of dynamo amplification in the generation of magnetic fields in deep convective layers of stars. The flows that arise in such environments are characterized by low (sonic) Mach numbers (M_son < 0.01 ). In these regimes, conventional MHD codes typically show excessive dissipation and tend to b…
▽ More
Fully compressible magnetohydrodynamic (MHD) simulations are a fundamental tool for investigating the role of dynamo amplification in the generation of magnetic fields in deep convective layers of stars. The flows that arise in such environments are characterized by low (sonic) Mach numbers (M_son < 0.01 ). In these regimes, conventional MHD codes typically show excessive dissipation and tend to be inefficient as the Courant-Friedrichs-Lewy (CFL) constraint on the time step becomes too strict. In this work we present a new method for efficiently simulating MHD flows at low Mach numbers in a space-dependent gravitational potential while still retaining all effects of compressibility. The proposed scheme is implemented in the finite-volume Seven-League Hydro (SLH) code, and it makes use of a low-Mach version of the five-wave Harten-Lax-van Leer discontinuities (HLLD) solver to reduce numerical dissipation, an implicit-explicit time discretization technique based on Strang splitting to overcome the overly strict CFL constraint, and a well-balancing method that dramatically reduces the magnitude of spatial discretization errors in strongly stratified setups. The solenoidal constraint on the magnetic field is enforced by using a constrained transport method on a staggered grid. We carry out five verification tests, including the simulation of a small-scale dynamo in a star-like environment at M_son ~ 0.001 . We demonstrate that the proposed scheme can be used to accurately simulate compressible MHD flows in regimes of low Mach numbers and strongly stratified setups even with moderately coarse grids.
△ Less
Submitted 4 October, 2022;
originally announced October 2022.
-
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.
-
Two-Dimensional Simulations of Internal Gravity Waves in The Radiation Zones of Intermediate-Mass Stars
Authors:
R. P. Ratnasingam,
P. V. F. Edelmann,
T. M. Rogers
Abstract:
Intermediate-mass main sequence stars have large radiative envelopes overlying convective cores. This configuration allows internal gravity waves (IGWs) generated at the convective-radiative interface to propagate towards the stellar surface. The signatures of these waves can be observed in the photometric and spectroscopic data from stars. We have studied the propagation of these IGWs using two-d…
▽ More
Intermediate-mass main sequence stars have large radiative envelopes overlying convective cores. This configuration allows internal gravity waves (IGWs) generated at the convective-radiative interface to propagate towards the stellar surface. The signatures of these waves can be observed in the photometric and spectroscopic data from stars. We have studied the propagation of these IGWs using two-dimensional fully-non-linear hydrodynamical simulations with realistic stellar reference states from the one-dimensional stellar evolution code, Modules for Stellar Astrophysics (MESA). When a single wave is forced, we observe wave self-interaction. When two waves are forced, we observe non-linear interaction (i.e. triadic interaction) between these waves forming waves at different wavelengths and frequencies. When a spectrum of waves similar to that found in numerical simulations is forced, we find that the surface IGW frequency slope is consistent with recent observations. This power law is similar to that predicted by linear theory for the wave propagation, with small deviations which can be an effect of nonlinearities. When the same generation spectrum is applied to 3 solar mass models at different stellar rotation and ages, the surface IGW spectrum slope is very similar to the generation spectrum slope.
△ Less
Submitted 7 August, 2020;
originally announced August 2020.
-
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.
-
Low-frequency gravity waves in blue supergiants revealed by high-precision space photometry
Authors:
Dominic M. Bowman,
Siemen Burssens,
May G. Pedersen,
Cole Johnston,
Conny Aerts,
Bram Buysschaert,
Mathias Michielsen,
Andrew Tkachenko,
Tamara M. Rogers,
Philipp V. F. Edelmann,
Rathish P. Ratnasingam,
Sergio Simón-Díaz,
Norberto Castro,
Ehsan Moravveji,
Benjamin J. S. Pope,
Timothy R. White,
Peter De Cat
Abstract:
Almost all massive stars explode as supernovae and form a black hole or neutron star. The remnant mass and the impact of the chemical yield on subsequent star formation and galactic evolution strongly depend on the internal physics of the progenitor star, which is currently not well understood. The theoretical uncertainties of stellar interiors accumulate with stellar age, which is particularly pe…
▽ More
Almost all massive stars explode as supernovae and form a black hole or neutron star. The remnant mass and the impact of the chemical yield on subsequent star formation and galactic evolution strongly depend on the internal physics of the progenitor star, which is currently not well understood. The theoretical uncertainties of stellar interiors accumulate with stellar age, which is particularly pertinent for the blue supergiant phase. Stellar oscillations represent a unique method of probing stellar interiors, yet inference for blue supergiants is hampered by a dearth of observed pulsation modes. Here we report the detection of diverse variability in blue supergiants using the K2 and TESS space missions. The discovery of pulsation modes or an entire spectrum of low-frequency gravity waves in these stars allow us to map the evolution of hot massive stars towards the ends of their lives. Future asteroseismic modelling will provide constraints on ages, core masses, interior mixing, rotation and angular momentum transport. The discovery of variability in blue supergiants is a step towards a data-driven empirical calibration of theoretical evolution models for the most massive stars in the Universe.
△ Less
Submitted 6 May, 2019;
originally announced May 2019.
-
Three-Dimensional Simulations of Massive Stars: I. Wave Generation and Propagation
Authors:
P. V. F. Edelmann,
R. P. Ratnasingam,
M. G. Pedersen,
D. M. Bowman,
V. Prat,
T. M. Rogers
Abstract:
We present the first three-dimensional (3D), hydrodynamic simulations of the core convection zone (CZ) and extended radiative zone spanning from 1% to 90% of the stellar radius of an intermediate mass (3 $\mathrm{M}_\odot$) star. This allows us to self-consistently follow the generation of internal gravity waves (IGWs) at the convective boundary and their propagation to the surface. We find that c…
▽ More
We present the first three-dimensional (3D), hydrodynamic simulations of the core convection zone (CZ) and extended radiative zone spanning from 1% to 90% of the stellar radius of an intermediate mass (3 $\mathrm{M}_\odot$) star. This allows us to self-consistently follow the generation of internal gravity waves (IGWs) at the convective boundary and their propagation to the surface. We find that convection in the core is dominated by plumes. The frequency spectrum in the CZ and that of IGW generation is a double power law as seen in previous two-dimensional (2D) simulations. The spectrum is significantly flatter than theoretical predictions using excitation through Reynolds stresses induced by convective eddies alone. It is compatible with excitation through plume penetration. An empirically determined distribution of plume frequencies generally matches the one necessary to explain a large part of the observed spectrum. We observe waves propagating in the radiation zone and excited standing modes, which can be identified as gravity and fundamental modes. They show similar frequencies and node patterns to those predicted by the stellar oscillation code GYRE. The continuous part of the spectrum fulfills the IGW dispersion relation. A spectrum of tangential velocity and temperature fluctuations close to the surface is extracted, which are directly related to observable brightness variations in stars. Unlike 2D simulations we do not see the high frequencies associated with wave breaking, likely because these 3D simulations are more heavily damped.
△ Less
Submitted 22 March, 2019;
originally announced March 2019.
-
Onset of Nonlinear Internal Gravity Waves in Intermediate-Mass Stars
Authors:
R. P. Ratnasingam,
P. V. F. Edelmann,
T. M. Rogers
Abstract:
Internal gravity waves (IGW) propagate in the radiation zones of all stars. During propagation, their amplitudes are affected by two main features: radiative diffusion and density stratification. We have studied the implications of these two features on waves traveling within the radiative zones of non-rotating stars with stellar parameters obtained from the one dimensional stellar evolution code,…
▽ More
Internal gravity waves (IGW) propagate in the radiation zones of all stars. During propagation, their amplitudes are affected by two main features: radiative diffusion and density stratification. We have studied the implications of these two features on waves traveling within the radiative zones of non-rotating stars with stellar parameters obtained from the one dimensional stellar evolution code, MESA. As a simple measure of induced wave dynamics, we define a criterion to see if waves can become nonlinear and if so, under what conditions. This was done to understand the role IGW may play in angular momentum transport and mixing within stellar interiors. We find that the IGW generation spectrum, convective velocities and the strength of density stratification all play major roles in whether waves become nonlinear. With increasing stellar mass, there is an increasing trend in nonlinear wave energies. The trends with different metallicities and ages depend on the generation spectrum.
△ Less
Submitted 3 December, 2018;
originally announced December 2018.
-
Photometric detection of internal gravity waves in early-type stars observed by CoRoT
Authors:
D. M. Bowman,
C. Aerts,
C. Johnston,
M. G. Pedersen,
T. M. Rogers,
P. V. F. Edelmann,
S. Simón-Díaz,
T. Van Reeth,
B. Buysschaert,
A. Tkachenko,
S. A. Triana
Abstract:
Early-type stars are predicted to excite an entire spectrum of internal gravity waves (IGWs) at the interface of their convective cores and radiative envelopes. Numerical simulations of IGWs predict stochastic low-frequency variability in photometric observations, yet the detection of IGWs in early-type stars has been limited by a dearth of high-quality photometric time series. We present observat…
▽ More
Early-type stars are predicted to excite an entire spectrum of internal gravity waves (IGWs) at the interface of their convective cores and radiative envelopes. Numerical simulations of IGWs predict stochastic low-frequency variability in photometric observations, yet the detection of IGWs in early-type stars has been limited by a dearth of high-quality photometric time series. We present observational evidence of stochastic low-frequency variability in the CoRoT photometry of a sample of O, B, A and F stars. The presence of this stochastic low-frequency variability in stars across the upper main-sequence cannot be universally explained as granulation or stellar winds, but its morphology is found to be consistent with predictions from IGW simulations.
△ Less
Submitted 30 November, 2018;
originally announced November 2018.
-
On the photometric detection of Internal Gravity Waves in upper main-sequence stars I. Methodology and application to CoRoT targets
Authors:
D. M. Bowman,
C. Aerts,
C. Johnston,
M. G. Pedersen,
T. M. Rogers,
P. V. F. Edelmann,
S. Simón-Díaz,
T. Van Reeth,
B. Buysschaert,
A. Tkachenko,
S. A. Triana
Abstract:
Context. Main sequence stars with a convective core are predicted to stochastically excite Internal Gravity Waves (IGWs), which effectively transport angular momentum throughout the stellar interior and explain the observed near-uniform interior rotation rates of intermediate-mass stars. However, there are few detections of IGWs, and fewer still made using photometry, with more detections needed t…
▽ More
Context. Main sequence stars with a convective core are predicted to stochastically excite Internal Gravity Waves (IGWs), which effectively transport angular momentum throughout the stellar interior and explain the observed near-uniform interior rotation rates of intermediate-mass stars. However, there are few detections of IGWs, and fewer still made using photometry, with more detections needed to constrain numerical simulations. Aims. We aim to formalise the detection and characterisation of IGWs in photometric observations of stars born with convective cores (M > 1.5 M$_{\odot}$) and parameterise the low-frequency power excess caused by IGWs. Methods. Using the most recent CoRoT light curves for a sample of O, B, A and F stars, we parameterise the morphology of the flux contribution of IGWs in Fourier space using an MCMC numerical scheme within a Bayesian framework. We compare this to predictions from IGW numerical simulations and investigate how the observed morphology changes as a function of stellar parameters. Results. We demonstrate that a common morphology for the low-frequency power excess is observed in early-type stars observed by CoRoT. Our study shows that a background frequency-dependent source of astrophysical signal is common, which we interpret as IGWs. We provide constraints on the amplitudes of IGWs and the shape of their detected frequency spectrum across a range of mass, which is the first ensemble study of stochastic variability in such a diverse sample of stars. Conclusions. The evidence of a low-frequency power excess across a wide mass range supports the interpretation of IGWs in photometry of O, B, A and F stars. We also discuss the prospects of observing hundreds of massive stars with the Transiting Exoplanet Survey Satellite (TESS) in the near future.
△ Less
Submitted 19 November, 2018;
originally announced November 2018.
-
Testing a one-dimensional prescription of dynamical shear mixing with a two-dimensional hydrodynamic simulation
Authors:
Philipp V. F. Edelmann,
Friedrich K. Roepke,
Raphael Hirschi,
Cyril Georgy,
Samuel Jones
Abstract:
The treatment of mixing processes is still one of the major uncertainties in 1D stellar evolution models. This is mostly due to the need to parametrize and approximate aspects of hydrodynamics in hydrostatic codes. In particular, the effect of hydrodynamic instabilities in rotating stars, for example, dynamical shear instability, evades consistent description. We intend to study the accuracy of th…
▽ More
The treatment of mixing processes is still one of the major uncertainties in 1D stellar evolution models. This is mostly due to the need to parametrize and approximate aspects of hydrodynamics in hydrostatic codes. In particular, the effect of hydrodynamic instabilities in rotating stars, for example, dynamical shear instability, evades consistent description. We intend to study the accuracy of the diffusion approximation to dynamical shear in hydrostatic stellar evolution models by comparing 1D models to a first-principle hydrodynamics simulation starting from the same initial conditions. We chose an initial model calculated with the stellar evolution code GENEC that is just at the onset of a dynamical shear instability but does not show any other instabilities (e.g., convection). This was mapped to the hydrodynamics code SLH to perform a 2D simulation in the equatorial plane. We compare the resulting profiles in the two codes and compute an effective diffusion coefficient for the hydro simulation. Shear instabilities develop in the 2D simulation in the regions predicted by linear theory to become unstable in the 1D model. Angular velocity and chemical composition is redistributed in the unstable region, thereby creating new unstable regions. Eventually the 2D simulation settles in a symmetric, steady state, which is Richardson stable everywhere, whereas the instability remains for longer in the 1D model due to current limitations in the 1D code. A spatially resolved diffusion coefficient is extracted by comparing the initial and final profiles of mean atomic mass. The presented simulation gives a first insight on hydrodynamics of shear instabilities in a real stellar environment and even allows us to directly extract an effective diffusion coefficient. We see evidence for a critical Richardson number of 0.25 as regions above this value remain stable for the course of the simulation.
△ Less
Submitted 20 April, 2017;
originally announced April 2017.
-
A numerical scheme for the compressible low-Mach number regime of ideal fluid dynamics
Authors:
Wasilij Barsukow,
Philipp V. F. Edelmann,
Christian Klingenberg,
Fabian Miczek,
Friedrich K. Roepke
Abstract:
Based on the Roe solver a new technique that allows to correctly represent low Mach number flows with a discretization of the compressible Euler equations was proposed in Miczek et al.: New numerical solver for flows at various mach numbers. A&A 576, A50 (2015). We analyze properties of this scheme and demonstrate that its limit yields a discretization of the continuous limit system. Furthermore w…
▽ More
Based on the Roe solver a new technique that allows to correctly represent low Mach number flows with a discretization of the compressible Euler equations was proposed in Miczek et al.: New numerical solver for flows at various mach numbers. A&A 576, A50 (2015). We analyze properties of this scheme and demonstrate that its limit yields a discretization of the continuous limit system. Furthermore we perform a linear stability analysis for the case of explicit time integration and study the performance of the scheme under implicit time integration via the evolution of its condition number. A numerical implementation demonstrates the capabilities of the scheme on the example of the Gresho vortex which can be accurately followed down to Mach numbers of ~1e-10 .
△ Less
Submitted 12 December, 2016;
originally announced December 2016.
-
Extreme Scale-out SuperMUC Phase 2 - lessons learned
Authors:
Nicolay Hammer,
Ferdinand Jamitzky,
Helmut Satzger,
Momme Allalen,
Alexander Block,
Anupam Karmakar,
Matthias Brehm,
Reinhold Bader,
Luigi Iapichino,
Antonio Ragagnin,
Vasilios Karakasis,
Dieter Kranzlmüller,
Arndt Bode,
Herbert Huber,
Martin Kühn,
Rui Machado,
Daniel Grünewald,
Philipp V. F. Edelmann,
Friedrich K. Röpke,
Markus Wittmann,
Thomas Zeiser,
Gerhard Wellein,
Gerald Mathias,
Magnus Schwörer,
Konstantin Lorenzen
, et al. (14 additional authors not shown)
Abstract:
In spring 2015, the Leibniz Supercomputing Centre (Leibniz-Rechenzentrum, LRZ), installed their new Peta-Scale System SuperMUC Phase2. Selected users were invited for a 28 day extreme scale-out block operation during which they were allowed to use the full system for their applications. The following projects participated in the extreme scale-out workshop: BQCD (Quantum Physics), SeisSol (Geophysi…
▽ More
In spring 2015, the Leibniz Supercomputing Centre (Leibniz-Rechenzentrum, LRZ), installed their new Peta-Scale System SuperMUC Phase2. Selected users were invited for a 28 day extreme scale-out block operation during which they were allowed to use the full system for their applications. The following projects participated in the extreme scale-out workshop: BQCD (Quantum Physics), SeisSol (Geophysics, Seismics), GPI-2/GASPI (Toolkit for HPC), Seven-League Hydro (Astrophysics), ILBDC (Lattice Boltzmann CFD), Iphigenie (Molecular Dynamic), FLASH (Astrophysics), GADGET (Cosmological Dynamics), PSC (Plasma Physics), waLBerla (Lattice Boltzmann CFD), Musubi (Lattice Boltzmann CFD), Vertex3D (Stellar Astrophysics), CIAO (Combustion CFD), and LS1-Mardyn (Material Science). The projects were allowed to use the machine exclusively during the 28 day period, which corresponds to a total of 63.4 million core-hours, of which 43.8 million core-hours were used by the applications, resulting in a utilization of 69%. The top 3 users were using 15.2, 6.4, and 4.7 million core-hours, respectively.
△ Less
Submitted 6 September, 2016;
originally announced September 2016.
-
Do electron-capture supernovae make neutron stars? First multidimensional hydrodynamic simulations of the oxygen deflagration
Authors:
Samuel Jones,
Friedrich K. Roepke,
Ruediger Pakmor,
Ivo R. Seitenzahl,
Sebastian T. Ohlmann,
Philipp V. F. Edelmann
Abstract:
In the classical picture, electron-capture supernovae and the accretion-induced collapse of oxygen-neon white dwarfs (ONeWDs) undergo an oxygen deflagration phase before gravitational collapse produces a neutron star (NS). These types of core collapse events are postulated to explain several astronomical phenomena. In this work, the deflagration phase is simulated for the first time using multidim…
▽ More
In the classical picture, electron-capture supernovae and the accretion-induced collapse of oxygen-neon white dwarfs (ONeWDs) undergo an oxygen deflagration phase before gravitational collapse produces a neutron star (NS). These types of core collapse events are postulated to explain several astronomical phenomena. In this work, the deflagration phase is simulated for the first time using multidimensional hydrodynamics, with the aim of gaining new insight into the explosive deaths of $8-10~M_\odot$ stars and ONeWDs that accrete material from a binary companion star. The main aim is to determine whether these events are thermonuclear or core-collapse supernova explosions, and hence whether NSs are formed by such phenomena. The deflagration is simulated in ONe cores with three different central ignition densities. The intermediate density case is perhaps the most realistic, being based on recent nuclear physics calculations and 1D stellar models. The 3D hydrodynamic simulations presented in this work begin from a centrally confined flame structure using a level-set-based flame approach and are performed in $256^3$ and $512^3$ numerical resolutions. In the simulations with intermediate and low ignition density, the cores do not appear to collapse into NSs. Instead, almost a solar mass of material becomes unbound from the cores, leaving bound remnants. These simulations represent the case in which semiconvective mixing during the electron-capture phase preceding the deflagration is inefficient. The masses of the bound remnants double when Coulomb corrections are included in the EoS, however they still do not exceed the effective Chandrasekhar mass and, hence, would not collapse into NSs. The simulations with the highest ignition density ($\log_{10}ρ_{\rm c}=10.3$), representing the case where semiconvective mixing is very efficient, show clear signs that the core will collapse into a NS.
△ Less
Submitted 13 July, 2016; v1 submitted 18 February, 2016;
originally announced February 2016.
-
New numerical solver for flows at various Mach numbers
Authors:
F. Miczek,
F. K. Roepke,
P. V. F. Edelmann
Abstract:
Many problems in stellar astrophysics feature flows at low Mach numbers. Conventional compressible hydrodynamics schemes frequently used in the field have been developed for the transonic regime and exhibit excessive numerical dissipation for these flows. While schemes were proposed that solve hydrodynamics strictly in the low Mach regime and thus restrict their applicability, we aim at developing…
▽ More
Many problems in stellar astrophysics feature flows at low Mach numbers. Conventional compressible hydrodynamics schemes frequently used in the field have been developed for the transonic regime and exhibit excessive numerical dissipation for these flows. While schemes were proposed that solve hydrodynamics strictly in the low Mach regime and thus restrict their applicability, we aim at developing a scheme that correctly operates in a wide range of Mach numbers. Based on an analysis of the asymptotic behavior of the Euler equations in the low Mach limit we propose a novel scheme that is able to maintain a low Mach number flow setup while retaining all effects of compressibility. This is achieved by a suitable modification of the well-known Roe solver. Numerical tests demonstrate the capability of this new scheme to reproduce slow flow structures even in moderate numerical resolution. Our scheme provides a promising approach to a consistent multidimensional hydrodynamical treatment of astrophysical low Mach number problems such as convection, instabilities, and mixing in stellar evolution.
△ Less
Submitted 6 February, 2015; v1 submitted 29 September, 2014;
originally announced September 2014.