-
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.
-
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.
-
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.
-
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.
-
Calibrating Core Overshooting Parameters With Two-dimensional Hydrodynamical Simulations
Authors:
Johann Higl,
Ewald Mueller,
Achim Weiss
Abstract:
The extent of mixed regions around convective zones is one of the biggest uncertainties in stellar evolution. 1D overshooting descriptions introduce a free parameter ($f_{ov}$) that is in general not well constrained from observations. Especially in small central convective regions the value is highly uncertain due to its tight connection to the pressure scale height. Long-term multi-dimensional h…
▽ More
The extent of mixed regions around convective zones is one of the biggest uncertainties in stellar evolution. 1D overshooting descriptions introduce a free parameter ($f_{ov}$) that is in general not well constrained from observations. Especially in small central convective regions the value is highly uncertain due to its tight connection to the pressure scale height. Long-term multi-dimensional hydrodynamic simulations can be used to study the size of the overshooting region and the involved mixing processes. Here we show how one can calibrate an overshooting parameter by performing 2D Maestro simulations of Zero-Age-Main-Sequence stars ranging from $1.3$ to $3.5 M_\odot$. The simulations cover the convective cores of the stars and a large fraction of the surrounding radiative envelope. We follow the convective flow for at least 20 convective turnover times, while the longest simulation covers 430 turnover time scales. This allows us to study how the mixing as well as the convective boundary evolve with time, and how the resulting entrainment can be interpreted in terms of overshooting parameters. We find that increasing the overshooting parameter $f_{ov}$ beyond a certain value in the initial model of our simulations, changes the mixing behaviour completely. This result can be used to put limits on the overshooting parameter. We find $0.010 < f_{ov} < 0.017$ to be in good agreement with our simulations of a $3.5 M_\odot$ mass star. We also identify a diffusive mixing component due to internal gravity waves (IGW) that is active throughout the convectively stable layer, but likely overestimated in our simulations. Furthermore, applying our calibration method to simulations of less massive stars suggests a need for a mass-dependent overshooting description where the mixing in terms of the pressure scale height is reduced for small convective cores.
△ Less
Submitted 9 December, 2020;
originally announced December 2020.