DOI: 10.5281/zenodo.17303430
title: Reproduction package for the paper: Resistive relativistic magnetohydrodynamics without Ampere's law

description: Resistive magnetohydrodynamics is thought to play a key role in transient astrophysical phenomena such as black hole flares and neutron star magnetospheres. When performing numerical simulations of resistive magnetohydrodynamics, one is faced with the issue that Amperes law becomes stiff in the high conductivity limit which poses challenges to the numerical evolution. We show that using a description of resistive magnetohydrodynamics based on higher form symmetry, one can perform simulations with a generalized dual Faraday tensor without having to use Amperes Law, thereby avoiding the stiffness problem. We also explain the relation of this dual model to a traditional description of resistive magnetohydrodynamics and how causality is guaranteed by introducing second order corrections. ArXiv: https://arxiv.org/abs/2501.04638ADS: https://ui.adsabs.harvard.edu/abs/2025arXiv250104638LAmerican Physical Society (APS), Physical Review D: https://doi.org/10.1103/stvc-tzfh
description: Reproducion package

Simulation results obtained with the open-source code BHAC (GPL v3)

Porth, O., Olivares, H., Mizuno, Y., Younsi, Z., Rezzolla, L., Moscibrodzka, M., Falcke, H., Kramer, M. (2017). The Black Hole Accretion Code. Computational Astrophysics and Cosmology, 4(1), 1. https://doi.org/10.1186/s40668-017-0020-2

www.bhac.science

This reproduction package includes raw data and plotting scripts. Filetypes: ipynb, vtu, vti, csv

scripts

fig1.ipynb  fig2.ipynb  fig3.ipynb  fig4.ipynb  fig5.ipynb  fig6.ipynb  fig8.ipynb  fig9.ipynb

data folders

fig1-data└── output-128    ├── amrvac.log    ├── data0000.dat    ├── data0000.vti    ├── data0001.dat    ├── data0001.vti    ├── data0002.dat    ├── data0002.vti    ├── data0003.dat    ├── data0003.vti    ├── data0004.dat    ├── data0004.vti    ├── data0005.dat    ├── data0005.vti    ├── data0006.dat    ├── data0006.vti    ├── data0007.dat    ├── data0007.vti    ├── data0008.dat    ├── data0008.vti    ├── data0009.dat    ├── data0009.vti    ├── data0010.dat    ├── data0010.vti    ├── data0011.dat    ├── data0011.vti    ├── data0012.dat    ├── data0012.vti    ├── data0013.dat    ├── data0013.vti    ├── data0014.dat    ├── data0014.vti    ├── data0015.dat    ├── data0015.vti    ├── data0016.dat    ├── data0016.vti    ├── data0017.dat    ├── data0017.vti    ├── data0018.dat    ├── data0018.vti    ├── data0019.dat    ├── data0019.vti    ├── data0020.dat    └── data0020.vtifig2-data├── output-128│   ├── amrvac.log│   ├── data0000.dat│   ├── data0000.vti│   ├── data0001.dat│   ├── data0001.vti│   ├── data0002.dat│   ├── data0002.vti│   ├── data0003.dat│   ├── data0003.vti│   ├── data0004.dat│   ├── data0004.vti│   ├── data0005.dat│   ├── data0005.vti│   ├── data0006.dat│   ├── data0006.vti│   ├── data0007.dat│   ├── data0007.vti│   ├── data0008.dat│   ├── data0008.vti│   ├── data0009.dat│   ├── data0009.vti│   ├── data0010.dat│   ├── data0010.vti│   └── out├── output-256│   ├── amrvac.log│   ├── data0000.dat│   ├── data0000.vti│   ├── data0001.dat│   ├── data0001.vti│   ├── data0002.dat│   ├── data0002.vti│   ├── data0003.dat│   ├── data0003.vti│   ├── data0004.dat│   ├── data0004.vti│   ├── data0005.dat│   ├── data0005.vti│   ├── data0006.dat│   ├── data0006.vti│   ├── data0007.dat│   ├── data0007.vti│   ├── data0008.dat│   ├── data0008.vti│   ├── data0009.dat│   ├── data0009.vti│   ├── data0010.dat│   ├── data0010.vti│   └── out├── output-512│   ├── amrvac.log│   ├── data0000.dat│   ├── data0000.vti│   ├── data0001.dat│   ├── data0001.vti│   ├── data0002.dat│   ├── data0002.vti│   ├── data0003.dat│   ├── data0003.vti│   ├── data0004.dat│   ├── data0004.vti│   ├── data0005.dat│   ├── data0005.vti│   ├── data0006.dat│   ├── data0006.vti│   ├── data0007.dat│   ├── data0007.vti│   ├── data0008.dat│   ├── data0008.vti│   ├── data0009.dat│   ├── data0009.vti│   ├── data0010.dat│   ├── data0010.vti│   └── out└── output-64    ├── amrvac.log    ├── data0000.dat    ├── data0000.vti    ├── data0001.dat    ├── data0001.vti    ├── data0002.dat    ├── data0002.vti    ├── data0003.dat    ├── data0003.vti    ├── data0004.dat    ├── data0004.vti    ├── data0005.dat    ├── data0005.vti    ├── data0006.dat    ├── data0006.vti    ├── data0007.dat    ├── data0007.vti    ├── data0008.dat    ├── data0008.vti    ├── data0009.dat    ├── data0009.vti    ├── data0010.dat    ├── data0010.vti    └── outfig3-data├── bdnk│   ├── output│   │   ├── amrvac.log│   │   ├── data0000.dat│   │   ├── data0000.vti│   │   ├── data0001.dat│   │   ├── data0001.vti│   │   ├── data0002.dat│   │   ├── data0002.vti│   │   ├── data0003.dat│   │   ├── data0003.vti│   │   ├── data0004.dat│   │   └── data0004.vti│   ├── output-1e-2│   │   ├── amrvac.log│   │   ├── data0000.dat│   │   ├── data0000.vti│   │   ├── data0001.dat│   │   ├── data0001.vti│   │   ├── data0002.dat│   │   ├── data0002.vti│   │   ├── data0003.dat│   │   ├── data0003.vti│   │   ├── data0004.dat│   │   └── data0004.vti│   └── output-1e-4│       ├── amrvac.log│       ├── data0000.dat│       ├── data0000.vti│       ├── data0001.dat│       ├── data0001.vti│       ├── data0002.dat│       ├── data0002.vti│       ├── data0003.dat│       ├── data0003.vti│       ├── data0004.dat│       └── data0004.vti├── rmhd│   └── output│       ├── amrvac.log│       ├── data0000.dat│       ├── data0000.vti│       ├── data0001.dat│       ├── data0001.vti│       ├── data0002.dat│       ├── data0002.vti│       ├── data0003.dat│       ├── data0003.vti│       ├── data0004.dat│       └── data0004.vti└── rrmhd    ├── output    │   ├── amrvac.log    │   ├── data0000.dat    │   ├── data0000.vti    │   ├── data0001.dat    │   ├── data0001.vti    │   ├── data0002.dat    │   ├── data0002.vti    │   ├── data0003.dat    │   ├── data0003.vti    │   ├── data0004.dat    │   └── data0004.vti    ├── output-1e-2    │   ├── amrvac.log    │   ├── data0000.dat    │   ├── data0000.vti    │   ├── data0001.dat    │   ├── data0001.vti    │   ├── data0002.dat    │   ├── data0002.vti    │   ├── data0003.dat    │   ├── data0003.vti    │   ├── data0004.dat    │   └── data0004.vti    └── output-1e-4        ├── amrvac.log        ├── data0000.dat        ├── data0000.vti        ├── data0001.dat        ├── data0001.vti        ├── data0002.dat        ├── data0002.vti        ├── data0003.dat        ├── data0003.vti        ├── data0004.dat        └── data0004.vtifig4-data└── output    ├── amrvac.log    ├── data0000.dat    ├── data0000.vtu    ├── data0001.dat    ├── data0001.vtu    ├── data0002.dat    ├── data0002.vtu    ├── data0003.dat    ├── data0003.vtu    ├── data0004.dat    └── data0004.vtufig6-data├── bdnk│   └── output-1024│       ├── amrvac.log│       ├── data0060.dat│       ├── data0060.vtu│       ├── data0100.dat│       ├── data0100.vtu│       └── data_analysis.csv├── rmhd│   └── output-1024│       ├── amrvac.log│       ├── data0060.dat│       ├── data0060.vtu│       ├── data0100.dat│       ├── data0100.vtu│       └── data_analysis.csv└── rrmhd    └── output-1024        ├── amrvac.log        ├── data0060.dat        ├── data0060.vtu        ├── data0100.dat        ├── data0100.vtu        └── data_analysis.csvfig9-data├── bdnk│   ├── cut0060.csv│   ├── cut.csv│   ├── data0030.vtu│   └── data0060.vtu├── cut.pvsm├── rmhd│   ├── cut0060.csv│   ├── cut.csv│   ├── data0030.vtu│   └── data0060.vtu└── rrmhd    ├── cut0060.csv    ├── cut.csv    ├── data0030.vtu    └── data0060.vtu