A robust super-time-stepping scheme for Ohmic and ambipolar diffusion

@misc{mattia2026robust,
  title = {{A robust super-time-stepping scheme for Ohmic and ambipolar diffusion}},
  author = {Giancarlo Mattia and Mario Flock and David Melon Fuksman and Anastasia Tzouvanou and Vittoria Berta and Daniele Crocco},
  howpublished = {arXiv},
  year = {2026},
  month = jun,
  doi = {10.48550/arxiv.2606.13398},
  eprint = {2606.13398},
  archivePrefix = {arXiv},
  url = {https://arxiv.org/abs/2606.13398},
}

TY  - UNPB
TI  - A robust super-time-stepping scheme for Ohmic and ambipolar diffusion
AU  - Giancarlo Mattia
AU  - Mario Flock
AU  - David Melon Fuksman
AU  - Anastasia Tzouvanou
AU  - Vittoria Berta
AU  - Daniele Crocco
PY  - 2026
DA  - 2026/06/11
PB  - arXiv
DO  - 10.48550/arxiv.2606.13398
UR  - https://arxiv.org/abs/2606.13398
AB  - Context. Non-ideal magnetohydrodynamics (MHD) is a key tool for modeling magnetic flux transport in astrophysical systems such as molecular clouds, protostellar cores, and protoplanetary disks. Conventional explicit methods for non-ideal MHD diffusion are severely limited by timestep constraints, while substepping approaches can be unstable due to truncation errors near boundaries and strong magnetic-field gradients. Aims. Our main goal is to address these limitations by developing robust super-time-stepping methods for Ohmic and ambipolar diffusion. Methods. We present a super-time-stepping method based on the stability of the Gegenbauer polynomials. The method is designed to enhance robustness in the presence of strongly anisotropic resistivity and to reduce sensitivity to truncation errors near boundaries. We implement the scheme in the PLUTO code and assess its performance through dedicated Ohmic and ambipolar diffusion tests. We also compare this novel numerical scheme against two common astrophysical problems, namely magnetic reconnection and the magnetorotational instability. Results. The novel Runge-Kutta-Gegenbauer scheme retains computational efficiency beyond purely explicit schemes while providing excellent stability compared with other traditional substepping methods. It remains stable in the presence of strongly anisotropic diffusion, enabling accurate magnetic-field evolution in regimes characteristic of protoplanetary disks and collapsing dense cores. Benchmark tests, including magnetic reconnection and magnetorotational-instability setups, confirm the method's accuracy, efficiency, and suitability for large-scale non-ideal MHD simulations.
ER  -