A New MHD Code with Adaptive Mesh Refinement and Parallelization for Astrophysics
arXiv:1204.5849 · doi:10.1016/j.cpc.2012.02.030
Abstract
A new code, named MAP, is written in Fortran language for magnetohydrodynamics (MHD) calculation with the adaptive mesh refinement (AMR) and Message Passing Interface (MPI) parallelization. There are several optional numerical schemes for computing the MHD part, namely, modified Mac Cormack Scheme (MMC), Lax-Friedrichs scheme (LF) and weighted essentially non-oscillatory (WENO) scheme. All of them are second order, two-step, component-wise schemes for hyperbolic conservative equations. The total variation diminishing (TVD) limiters and approximate Riemann solvers are also equipped. A high resolution can be achieved by the hierarchical block-structured AMR mesh. We use the extended generalized Lagrange multiplier (EGLM) MHD equations to reduce the non-divergence free error produced by the scheme in the magnetic induction equation. The numerical algorithms for the non-ideal terms, e.g., the resistivity and the thermal conduction, are also equipped in the MAP code. The details of the AMR and MPI algorithms are described in the paper.
44 pages, 16 figures
References in corpus (4)
- PLUTO: a Numerical Code for Computational Astrophysics
- A High Order Godunov Scheme with Constrained Transport and Adaptive Mesh Refinement for Astrophysical MHD
- An Unsplit Godunov Method for Ideal MHD via Constrained Transport in Three Dimensions
- Fan-shaped jets in three dimensional reconnection simulation as a model of ubiquitous solar jets