Well-balanced treatment of gravity in astrophysical fluid dynamics simulations at low Mach numbers
- 11 August 2021
- journal article
- research article
- Published by EDP Sciences in Astronomy & Astrophysics
- Vol. 652, A53
Abstract
Context.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.Aims.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.Methods.Three such schemes were applied in the implicit, finite-volume SEVEN-LEAGUEHYDROcode 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.Results.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.Conclusions.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.Keywords
All Related Versions
Funding Information
- Klaus Tschira Foundation
- Deutsche Forschungsgemeinschaft (GRK 1147, KL 566/22-1, RO 3676/3-1)
- U.S. Department of Energy (89233218CNA000001)
This publication has 71 references indexed in Scilit:
- FISH: A THREE-DIMENSIONAL PARALLEL MAGNETOHYDRODYNAMICS CODE FOR ASTROPHYSICAL APPLICATIONSThe Astrophysical Journal Supplement Series, 2011
- Turbulent Convection in Stellar Interiors. I. Hydrodynamic SimulationThe Astrophysical Journal, 2007
- Low Mach Number Modeling of Type Ia Supernovae. I. HydrodynamicsThe Astrophysical Journal, 2006
- Active Carbon and Oxygen Shell Burning HydrodynamicsThe Astrophysical Journal, 2006
- Simulations of Core Convection in Rotating A‐Type Stars: Differential Rotation and OvershootingThe Astrophysical Journal, 2004
- A numerical model for the flooding and drying of irregular domainsInternational Journal for Numerical Methods in Fluids, 2002
- The Evolution of Rotating StarsAnnual Review of Astronomy and Astrophysics, 2000
- Riemann Solvers for some Hyperbolic Problems with a Source TermESAIM: Proceedings, 1999
- Low-Diffusion Flux-Splitting Methods for Flows at All SpeedsAIAA Journal, 1998
- Analysis and implementation of TR-BDF2Applied Numerical Mathematics, 1996