The Athena Astrophysical MHD Code in Cylindrical Geometry

Aaron Skinner,Eve Ostriker
DOI: https://doi.org/10.1088/0067-0049/188/1/290
2010-04-15
Abstract:A method for implementing cylindrical coordinates in the Athena magnetohydrodynamics (MHD) code is described. The extension follows the approach of Athena's original developers and has been designed to alter the existing Cartesian-coordinates code as minimally and transparently as possible. The numerical equations in cylindrical coordinates are formulated to maintain consistency with constrained transport, a central feature of the Athena algorithm, while making use of previously implemented code modules such as the Riemann solvers. Angular-momentum transport, which is critical in astrophysical disk systems dominated by rotation, is treated carefully. We describe modifications for cylindrical coordinates of the higher-order spatial reconstruction and characteristic evolution steps as well as the finite-volume and constrained transport updates. Finally, we present a test suite of standard and novel problems in one-, two-, and three-dimensions designed to validate our algorithms and implementation and to be of use to other code developers. The code is suitable for use in a wide variety of astrophysical applications and is freely available for download on the web.
Instrumentation and Methods for Astrophysics,Astrophysics of Galaxies,High Energy Astrophysical Phenomena
What problem does this paper attempt to address?