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.

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

References in corpus (18)

• 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
• Three-Dimensional Simulations of Massive Stars: I. Wave Generation and Propagation
• 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
• 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
• Fully compressible simulations of waves and core convection in main-sequence stars
• Convective core entrainment in 1D main sequence stellar models
• 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
• Calibrating Core Overshooting Parameters With Two-dimensional Hydrodynamical Simulations
• Benchmarking the Multi-dimensional Stellar Implicit Code MUSIC
• Multidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star

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