THE ATHENA ASTROPHYSICAL MAGNETOHYDRODYNAMICS CODE IN CYLINDRICAL GEOMETRY
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 (CT), 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 CT 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.
- OSTI ID:
- 21454857
- Journal Information:
- Astrophysical Journal, Supplement Series, Vol. 188, Issue 1; Other Information: DOI: 10.1088/0067-0049/188/1/290; ISSN 0067-0049
- Country of Publication:
- United States
- Language:
- English
Similar Records
DISCO: A 3D MOVING-MESH MAGNETOHYDRODYNAMICS CODE DESIGNED FOR THE STUDY OF ASTROPHYSICAL DISKS
AN EXTENSION OF THE ATHENA++ CODE FRAMEWORK FOR GRMHD BASED ON ADVANCED RIEMANN SOLVERS AND STAGGERED-MESH CONSTRAINED TRANSPORT