The Mezcal AMR/MHD Code

Introduction

Over the last few years I have developed a finite-volume numerical code in one, two and three dimensions in cartesian (1/2/3d), cylindrical (1/2 d) and spherical (1d) coordinates. The code employs a (non-standard) hybrid block/cell-based AMR and it is parallelized using MPI, including an efficient load balancing algorithm. The modular AMR framework can in principle be coupled with any kind of hyperbolic system of equations of arbitrary spatial and temporal order. Currently, the code is used to solve the MHD equations, while an implementation of the relativistic MHD equations is under development.

Code characteristics

  • Geometry: 1/2/3 D cartesian; 1/2 D cylindrical; 1 D spherical.
  • Equations: MHD, RMHD
  • Equation of state: ideal gas, Helmholtz, variable &Gamma (SRMHD)
  • Source terms: external gravity, thermal conduction (isotropic+anisotropic), heating and cooling (from a simple table to a chemistry network)
  • Time integration: II order Runge-Kutta.
  • Space reconstruction: II order (except in shocks), with slope limiters on primitive variables
  • Flux integration: Lax-Friedrichs modified method, Lax-Wendroff, HLL, HLLD, MUSTA method, ``Riemann Solver'' in MHD; HLL in SRMHD
  • Source terms integration: fractional step method; some cooling functions integrated semi-implicitly;
  • Artificial viscosity: ``Lapidus viscosity''.
  • Div B mantained close to zero by constrained trasnport method (with or without staggered mesh)
  • Grid: uniform or AMR

    • AMR Code description

    coming soon!

    Tests

    coming soon!