Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes

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 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\sigma$ of the mean of all simulations on a given grid ($128^3$ and $256^3$ grid cells), where $\sigma$ 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.

arXiv comment20 pages, 19 figures, new Fig. B.2, new clarifications, electronic materials updated, accepted for publication in A&A

References in corpus (52)

• Grids of stellar models with rotation - I. Models from 0.8 to 120 Msun at solar metallicity (Z = 0.014)
• Comparing the statistics of interstellar turbulence in simulations and observations: Solenoidal versus compressive turbulence forcing
• Turbulent Convection in Stellar Interiors. I. Hydrodynamic Simulation
• Probing the interior physics of stars through asteroseismology
• The AGORA High-Resolution Galaxy Simulations Comparison Project
• Type Ia Supernova Explosion: Gravitationally Confined Detonation
• The Launching of Cold Clouds by Galaxy Outflows I: Hydrodynamic Interactions with Radiative Cooling
• The Role of Turbulence in Neutrino-Driven Core-Collapse Supernova Explosions
• The Stagger-grid: A grid of 3D stellar atmosphere models - III. The relation to mixing-length convection theory
• The spectroscopic Hertzsprung-Russell diagram of Galactic massive stars
• The AGORA High-Resolution Galaxy Simulations Comparison Project. II: Isolated Disk Test
• A Well-Posed Kelvin-Helmholtz Instability Test and Comparison
• Bottleneck effect in three-dimensional turbulence simulations
• Turbulent Convection in Stellar Interiors. II. The Velocity Field
• Turbulent convection in stellar interiors. III. Mean-field analysis and stratification effects
• Global non-spherical oscillations in 3-D $4\pi$ simulations of the H-ingestion flash
• Three-Dimensional Simulations of Massive Stars: I. Wave Generation and Propagation
• A Validated Nonlinear Kelvin-Helmholtz Benchmark for Numerical Hydrodynamics
• Simulations of the solar near-surface layers with the CO5BOLD, MURaM, and Stagger codes
• MESA Models of Classical Nova Outbursts: The Multicycle Evolution and Effects of Convective Boundary Mixing
• Idealised hydrodynamic simulations of turbulent oxygen-burning shell convection in 4π geometry
• The i-process yields of rapidly-accreting white dwarfs from multicycle He-shell flash stellar evolution models with mixing parameterizations from 3D hydrodynamics simulations
• The dependence of convective core overshooting on stellar mass: reality check, and additional evidence
• New numerical solver for flows at various Mach numbers
• Theoretical seismology in 3D : nonlinear simulations of internal gravity waves in solar-like stars
• Application of a Theory and Simulation based Convective Boundary Mixing model for AGB Star Evolution and Nucleosynthesis
• Improvements to Stellar Structure Models, Based on a Grid of 3D Convection Simulations. I. $T(\tau )$-Relations
• Lithium depletion in solar-like stars: effect of overshooting based on realistic multi-dimensional simulations
• 3D hydrodynamic simulations of C ingestion into a convective O shell
• Fully compressible simulations of waves and core convection in main-sequence stars
• Convective core entrainment in 1D main sequence stellar models
• The Aarhus Red Giants Challenge I: Stellar structures in the red giant branch phase
• Dependence of Convective Boundary Mixing on Boundary Properties and Turbulence Strength
• Extreme value statistics for two-dimensional convective penetration in a pre-Main Sequence star
• The treatment of mixing in core helium burning models -- III. Suppressing core breathing pulses with a new constraint on overshoot
• A Jacobian-free Newton-Krylov method for time-implicit multidimensional hydrodynamics
• Coupling 1D stellar evolution with 3D-hydrodynamical simulations on-the-fly I: A New Standard Solar Model
• The Aarhus red giants challenge II. Stellar oscillations in the red giant branch phase
• Calibrating Core Overshooting Parameters With Two-dimensional Hydrodynamical Simulations
• Type Ia Supernova Explosions from Hybrid Carbon-Oxygen-Neon White Dwarf Progenitors
• Turbulent Mixing and Nuclear Burning in Stellar Interiors
• Benchmarking the Multi-dimensional Stellar Implicit Code MUSIC
• Calibration of the mixing-length parameter $\alpha$ for the MLT and FST models by matching with CO${}^5$BOLD models
• 3D1D hydro-nucleosynthesis simulations. I. Advective-reactive post-processing method and its application to H ingestion into He-shell flash convection in rapidly accreting white dwarfs
• Multi-dimensional structure of accreting young stars
• The growth of helium burning cores
• Comparison of two- and three-dimensional compressible convection in a pre-main sequence star
• Testing a one-dimensional prescription of dynamical shear mixing with a two-dimensional hydrodynamic simulation
• Well-balanced treatment of gravity in astrophysical fluid dynamics simulations at low Mach numbers
• Surface Manifestation of Stochastically Excited Internal Gravity Waves
• Multidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star
• Acoustic-gravity wave propagation characteristics in 3D radiation hydrodynamic simulations of the solar atmosphere

Cited by in corpus (8)

• 3D hydrodynamic simulations of massive main-sequence stars. I. Dynamics and mixing of convection and internal gravity waves
• 3D stellar evolution: hydrodynamic simulations of a complete burning phase in a massive star
• Realistic 3D hydrodynamics simulations find significant turbulent entrainment in massive stars
• A finite-volume scheme for modeling compressible magnetohydrodynamic flows at low Mach numbers in stellar interiors
• Multidimensional Simulations of Core Convection
• Turbulent dynamo action and its effects on the mixing at the convective boundary of an idealized oxygen-burning shell
• Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows
• 3D Convective Urca Process in a Simmering White Dwarf