Multidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star

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 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.

arXiv commentAccepted for publication in A&A

References in corpus (43)

• Modules for Experiments in Stellar Astrophysics (MESA)
• Modules for Experiments in Stellar Astrophysics (MESA): Giant Planets, Oscillations, Rotation, and Massive Stars
• Modules for Experiments in Stellar Astrophysics (MESA): Binaries, Pulsations, and Explosions
• Modules for Experiments in Stellar Astrophysics (MESA): Pulsating Variable Stars, Rotation, Convective Boundaries, and Energy Conservation
• Turbulent Convection in Stellar Interiors. I. Hydrodynamic Simulation
• Probing the interior physics of stars through asteroseismology
• Hydrodynamics of core-collapse supernovae and their progenitors
• Sub-Inertial Gravity Modes in the B8V Star KIC 7760680 Reveal Moderate Core Overshooting and Low Vertical Diffusive Mixing
• Internal mixing of rotating stars inferred from dipole gravity modes
• The Last Minutes of Oxygen Shell Burning in a Massive Star
• Phenomenology of buoyancy-driven turbulence: recent results
• Turbulent Convection in Stellar Interiors. II. The Velocity Field
• Beyond Mixing-length Theory: a step toward 321D
• The shape of convective core overshooting from gravity-mode period spacings
• MAESTRO: An Adaptive Low Mach Number Hydrodynamics Algorithm for Stellar Flows
• Idealised hydrodynamic simulations of turbulent oxygen-burning shell convection in 4π geometry
• 3D Hydrodynamic Simulations of Carbon Burning in Massive Stars
• New numerical solver for flows at various Mach numbers
• Toward a consistent use of overshooting parametrizations in 1D stellar evolution codes
• Application of a Theory and Simulation based Convective Boundary Mixing model for AGB Star Evolution and Nucleosynthesis
• 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
• Probing the shape of the mixing profile and of the thermal structure at the convective core boundary through asteroseismology
• 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
• Angular momentum fluctuations in the convective helium shell of massive stars
• Properties of Convective Oxygen and Silicon Burning Shells in Supernova Progenitors
• A Jacobian-free Newton-Krylov method for time-implicit multidimensional hydrodynamics
• Convective boundary mixing in a post-He core burning massive star model
• Calibrating Core Overshooting Parameters With Two-dimensional Hydrodynamical Simulations
• Convective boundary mixing in low- and intermediate-mass stars I. Core properties from pressure-mode asteroseismology
• Turbulent Mixing and Nuclear Burning in Stellar Interiors
• Benchmarking the Multi-dimensional Stellar Implicit Code MUSIC
• MAESTROeX: A Massively Parallel Low Mach Number Astrophysical Solver
• 3D Simulations and MLT: I. Renzini's Critique
• Examining the Relationship Between Convective Core Overshoot and Stellar Properties Using Asteroseismology
• Comparison of different nonlinear solvers for 2D time-implicit stellar hydrodynamics
• 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
• Synergies between Asteroseismology and Three-dimensional Simulations of Stellar Turbulence
• A well-balanced scheme for the simulation tool-kit A-MaZe: implementation, tests, and first applications to stellar structure

Cited by in corpus (13)

• Stellar convective penetration: parameterized theory and dynamical simulations
• Numerical experiments to help understand cause and effect in massive star evolution
• 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
• Is Thermohaline Mixing the Full Story? Evidence for Separate Mixing Events near the Red Giant Branch Bump
• Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes
• Well-balanced treatment of gravity in astrophysical fluid dynamics simulations at low Mach numbers
• A finite-volume scheme for modeling compressible magnetohydrodynamic flows at low Mach numbers in stellar interiors
• The impact of convective criteria on the properties of massive stars
• Toward fully compressible numerical simulations of stellar magneto-convection with the RAMSES code
• 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
• Spectral Difference method with a posteriori limiting: II- Application to low Mach number flows