Athena: A New Code for Astrophysical MHD

Abstract

A new code for astrophysical magnetohydrodynamics (MHD) is described. The code has been designed to be easily extensible for use with static and adaptive mesh refinement. It combines higher-order Godunov methods with the constrained transport (CT) technique to enforce the divergence-free constraint on the magnetic field. Discretization is based on cell-centered volume-averages for mass, momentum, and energy, and face-centered area-averages for the magnetic field. Novel features of the algorithm include (1) a consistent framework for computing the time- and edge-averaged electric fields used by CT to evolve the magnetic field from the time- and area-averaged Godunov fluxes, (2) the extension to MHD of spatial reconstruction schemes that involve a dimensionally-split time advance, and (3) the extension to MHD of two different dimensionally-unsplit integration methods. Implementation of the algorithm in both C and Fortran95 is detailed, including strategies for parallelization using domain decomposition. Results from a test suite which includes problems in one-, two-, and three-dimensions for both hydrodynamics and MHD are given, not only to demonstrate the fidelity of the algorithms, but also to enable comparisons to other methods. The source code is freely available for download on the web.

61 pages, 36 figures. accepted by ApJS

References in corpus (7)

• PLUTO: a Numerical Code for Computational Astrophysics
• An Unsplit Godunov Method for Ideal MHD via Constrained Transport in Three Dimensions
• The Magnetic Rayleigh-Taylor Instability in Three Dimensions
• WHAM: A WENO-based general relativistic numerical scheme I: Hydrodynamics
• The Magnetohydrodynamics of Shock-Cloud Interaction in Three Dimensions
• Nonlinear Evolution of the Magnetohydrodynamic Rayleigh-Taylor Instability
• Saturation of the Magnetothermal Instability in Three Dimensions

Cited by in corpus (10)

• Simulating Magnetohydrodynamical Flow with Constrained Transport and Adaptive Mesh Refinement; Algorithms & Tests of the AstroBEAR Code
• Anisotropic Thermal Conduction and the Cooling Flow Problem in Galaxy Clusters
• Density Probability Distribution Functions in Supersonic Hydrodynamic and MHD Turbulence
• Simulations of Magnetorotational Turbulence with a Higher-Order Godunov Scheme
• The Magnetothermal Instability in the Intracluster Medium
• Simulations of MHD Instabilities in Intracluster Medium Including Anisotropic Thermal Conduction
• Local Star formation triggered by SN shocks in magnetized diffuse neutral clouds
• Simulating Supersonic Turbulence in Magnetized Molecular Clouds
• Extremum-Preserving Limiters for MUSCL and PPM
• Spectral Methods for Time-Dependent Studies of Accretion Flows. III. Three-Dimensional, Self-gravitating, Magnetohydrodynamic Disks