Well-balanced treatment of gravity in astrophysical fluid dynamics simulations at low Mach numbers

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 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 $\alpha$-$\beta$ 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.

arXiv commentAccepted for publication in A&A

References in corpus (30)

• The evolution of rotating stars
• Turbulent Convection in Stellar Interiors. I. Hydrodynamic Simulation
• Internal Gravity Waves in Massive Stars: Angular Momentum Transport
• Simulations of core convection in rotating A-type stars: Differential rotation and overshooting
• Low Mach Number Modeling of Type Ia Supernovae
• The Last Minutes of Oxygen Shell Burning in a Massive Star
• Active Carbon and Oxygen Shell Burning Hydrodynamics
• Mapping Initial Hydrostatic Models in Godunov Codes
• Turbulent convection in stellar interiors. III. Mean-field analysis and stratification effects
• Three-Dimensional Simulations of Massive Stars: I. Wave Generation and Propagation
• 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
• An advanced leakage scheme for neutrino treatment in astrophysical simulations
• 3D hydrodynamic simulations of C ingestion into a convective O shell
• Fully compressible simulations of waves and core convection in main-sequence stars
• Well balanced Arbitrary-Lagrangian-Eulerian finite volume schemes on moving nonconforming meshes for the Euler equations of gasdynamics with gravity
• Dependence of Convective Boundary Mixing on Boundary Properties and Turbulence Strength
• Third and fourth order well-balanced schemes for the shallow water equations based on the CWENO reconstruction
• Calibrating Core Overshooting Parameters With Two-dimensional Hydrodynamical Simulations
• Benchmarking the Multi-dimensional Stellar Implicit Code MUSIC
• Spherical-shell boundaries for two-dimensional compressible convection in a star
• FISH: A 3D parallel MHD code for astrophysical applications
• 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
• A high performance and portable all-Mach regime flow solver code with well-balanced gravity. Application to compressible convection
• Multidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star
• An Extension of the Athena++ Framework for Fully Conservative Self-Gravitating Hydrodynamics
• Second order finite volume scheme for Euler equations with gravity which is well-balanced for general equations of state and grid systems
• A well-balanced scheme for the simulation tool-kit A-MaZe: implementation, tests, and first applications to stellar structure

Cited by in corpus (10)

• Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes
• Multidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star
• A finite-volume scheme for modeling compressible magnetohydrodynamic flows at low Mach numbers in stellar interiors
• A well-balanced and exactly divergence-free staggered semi-implicit hybrid finite volume/finite element scheme for the incompressible MHD equations
• Multidimensional Simulations of Core Convection
• Turbulent dynamo action and its effects on the mixing at the convective boundary of an idealized oxygen-burning shell
• Toward fully compressible numerical simulations of stellar magneto-convection with the RAMSES code
• 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
• DISPATCH methods: an approximate, entropy-based Riemann solver for ideal magnetohydrodynamics