An Extension of the Athena++ Framework for Fully Conservative Self-Gravitating Hydrodynamics

Abstract

Numerical simulations of self-gravitating flows evolve a momentum equation and an energy equation that account for accelerations and gravitational energy releases due to a time-dependent gravitational potential. In this work, we implement a fully conservative numerical algorithm for self-gravitating flows, using source terms, in the astrophysical magnetohydrodynamics framework Athena++. We demonstrate that properly evaluated source terms are conservative when they are equivalent to the divergence of a corresponding "gravity flux" (i.e., a gravitational stress tensor or a gravitational energy flux). We provide test problems that demonstrate several advantages of the source-term-based algorithm, including second order convergence and round-off error total momentum and total energy conservation. The fully conservative scheme suppresses anomalous accelerations that arise when applying a common numerical discretization of the gravitational stress tensor that does not guarantee curl-free gravity.

34 pages, 9 figures

References in corpus (9)

• Theory of Star Formation
• Athena: A New Code for Astrophysical MHD
• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• An Unsplit Godunov Method for Ideal MHD via Constrained Transport in Three Dimensions
• A Direct Multigrid Poisson Solver for Oct-Tree Adaptive Meshes
• Gravito-Turbulent Disks in 3D: Turbulent Velocities vs. Depth
• Characterizing gravito-turbulence in 3D: turbulent properties and stability against fragmentation
• A New Way to Conserve Total Energy for Eulerian Hydrodynamic Simulations with Self-Gravity
• Three Dimensional Magnetohydrodynamical Simulations of Core Collapse Supernova

Cited by in corpus (2)

• The Athena++ Adaptive Mesh Refinement Framework: Multigrid Solvers for Self-Gravity
• Standing solitary waves as transitions to spiral structures in gravitationally unstable accretion disks