pith. sign in

arxiv: 2606.13398 · v1 · pith:267XQZKZnew · submitted 2026-06-11 · 🌌 astro-ph.IM · physics.plasm-ph

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

Giancarlo Mattia , Mario Flock , David Melon Fuksman , Anastasia Tzouvanou , Vittoria Berta , Daniele Crocco This is my paper

Pith reviewed 2026-06-27 05:41 UTC · model grok-4.3

classification 🌌 astro-ph.IM physics.plasm-ph
keywords super-time-steppingOhmic diffusionambipolar diffusionnon-ideal MHDGegenbauer polynomialsRunge-Kutta methodsprotoplanetary disksmagnetic reconnection
0
0 comments X

The pith

A Runge-Kutta-Gegenbauer super-time-stepping scheme achieves stable super-timestepping for Ohmic and ambipolar diffusion even with strong anisotropy.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces a super-time-stepping method for the diffusion terms that appear in non-ideal magnetohydrodynamic simulations of astrophysical systems. Conventional explicit integrators are restricted to tiny timesteps set by the fastest diffusion rate, while many existing substepping techniques become unstable when resistivity varies sharply or near domain boundaries. The new scheme draws on the stability region of Gegenbauer polynomials to advance the diffusion operators over multiple internal steps inside a single Runge-Kutta stage. Tests in the PLUTO code show that the method keeps the efficiency gain of super-time-stepping while remaining stable under the anisotropic resistivities typical of protoplanetary disks and collapsing cores. A sympathetic reader would therefore see a route to longer, more accurate simulations of magnetic flux transport without the usual stability penalties.

Core claim

The 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.

What carries the argument

The Runge-Kutta-Gegenbauer scheme, which uses the stability properties of Gegenbauer polynomials to advance diffusion terms over multiple substeps while preserving overall stability.

If this is right

  • Diffusion terms can be integrated with effective timesteps larger than the explicit CFL limit without sacrificing stability.
  • Strongly anisotropic resistivity no longer forces a reduction to the smallest local diffusion timescale.
  • Truncation errors near boundaries produce less instability than in conventional super-time-stepping methods.
  • The scheme reproduces known solutions for magnetic reconnection and the magnetorotational instability at the same accuracy level as explicit methods.
  • Large-scale simulations of protoplanetary disks and dense cores become feasible with non-ideal MHD physics included throughout the domain.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The method could be combined with existing Godunov-type MHD solvers to treat the full induction equation without operator splitting.
  • Extension to Hall diffusion would require only a change in the diffusion operator while retaining the same substepping framework.
  • Adaptive selection of the number of Gegenbauer stages based on local resistivity eigenvalues could further improve efficiency.
  • The same polynomial construction might apply to other parabolic operators such as thermal conduction in stratified atmospheres.

Load-bearing premise

The stability properties of Gegenbauer polynomials will continue to hold when resistivity is strongly anisotropic and when truncation errors appear near boundaries.

What would settle it

A dedicated test run with extreme anisotropy in the resistivity tensor that produces growing oscillations or unphysical field reversals would show the scheme is not robust.

Figures

Figures reproduced from arXiv: 2606.13398 by Anastasia Tzouvanou, Daniele Crocco, David Melon Fuksman, Giancarlo Mattia, Mario Flock, Vittoria Berta.

Figure 1
Figure 1. Figure 1: Ratio between the RKG steps and the RKL steps (i.e., RKG with α = 0.5, represented by the vertical dashed line) for different values of α and τ. In the top (bottom) panel, the ratio is shown for the 1st- order (2nd-order) algorithms. The solid black lines indicate the τ values below which the RKG method would require fewer than 3 steps. derivatives) to be unity at z = 0 (O’Sullivan 2019). Once the coef￾fic… view at source ↗
Figure 2
Figure 2. Figure 2: L1 errors for RKG-1 and RKG-2 on a decaying sine mode of the one-dimensional scalar decay, versus timestep (top row) and grid resolution under adaptive CFL stepping (bottom row), for different values of the resistivity η and the Gegenbauer parameter α. Note that, since the z−component of the magnetic field is set to zero and the problem is two-dimensional, ideally, the x− and y− resistivity components shou… view at source ↗
Figure 3
Figure 3. Figure 3: Electric current for the resistive Orszag-Tang vortex cases. The top, middle, and bottom panels show, respectively, the cases HΩ, LΩ, and TΩ. In the different columns, from left to right, the EXPl, RKG, STS, and RKL cases are reported, respectively. The black dashed lines represent the slice shown in [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 6
Figure 6. Figure 6: Cuts of the plasma-β for the ambipolar Orszag-Tang vortex tests. Dashed and dotted lines represent, respectively, x- and y-boundaries. The top and bottom panels show, respectively, the results obtained from the cases Had and Lad. given by η = v 2 A γadρi . (28) Solving the dispersion relation yields a damped oscillatory solution provided that (vAk) 2 is larger than the square of the damping rate γd. The la… view at source ↗
Figure 5
Figure 5. Figure 5: Plasma magnetization for the case AΩ of the Orszag-Tang vortex test at times t = π/2 (top panels) and t = 2π (bottom panels) for the RKL (left panel) and RKG (right panels) methods [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 7
Figure 7. Figure 7: Temporal evolution of the root-mean-square transverse magnetic perturbation ⟨B 2 z ⟩ 1/2 for resistive Alfvén waves with ηΩ = 0.01 (top) and ηΩ = 0.05 (bottom). Solid lines show the analytical solution, while sym￾bols denote numerical results (EXPL, RKG, STS, and RKL) sampled at successive oscillation maxima [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Same as [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Comparison between the EXPL (left) and RKG (right) schemes for magnetic reconnection at Lundquist numbers S = 103 (first and third rows) and S = 104 (second and fourth rows). The top two rows show the density ρ, while the bottom two rows display the magnetic energy density B 2 . Snapshots are shown at time t = 3×105 [PITH_FULL_IMAGE:figures/full_fig_p010_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Same as [PITH_FULL_IMAGE:figures/full_fig_p010_10.png] view at source ↗
Figure 12
Figure 12. Figure 12: Temporal evolution of max(Br[r = 1]) for ηad = 0.1, 0.5, and 1.0 (top to bottom), comparing the EXPL (blue) and RKG (red) schemes for ambipolar MRI. The dotted line indicates the theoretical growth rate. Black crosses mark the time interval used to fit the linear growth rate; filled circles and squares denote the time corresponding to [PITH_FULL_IMAGE:figures/full_fig_p011_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Spatial distribution of the local Alfvén speed squared, B 2 /ρ, for ambipolar MRI, shown for η0 = 0.1, 0.5, and 1.0. The left and right columns compare the EXPL and RKG schemes, respectively. Top and bottom panels correspond to the times marked by the filled circles and filled squares in [PITH_FULL_IMAGE:figures/full_fig_p012_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Same as [PITH_FULL_IMAGE:figures/full_fig_p012_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Same as [PITH_FULL_IMAGE:figures/full_fig_p013_15.png] view at source ↗
read the original abstract

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.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 2 minor

Summary. The manuscript introduces a Runge-Kutta-Gegenbauer super-time-stepping scheme for Ohmic and ambipolar diffusion in non-ideal MHD. Implemented in the PLUTO code, the method is tested on dedicated diffusion problems and benchmarked against magnetic reconnection and magnetorotational instability setups, with the central claim that it delivers efficiency gains over explicit schemes, superior stability relative to other substepping approaches, and robustness under strong anisotropy.

Significance. If the reported stability and accuracy hold, the scheme would enable larger timesteps in simulations of protoplanetary disks and collapsing cores where anisotropic non-ideal effects dominate. The approach is grounded in polynomial stability properties rather than fitted parameters, and the inclusion of reconnection and MRI benchmarks directly targets the regimes of interest. The stress-test concern regarding translation of Gegenbauer stability to anisotropic resistivity and boundary truncation errors does not land, as the manuscript provides dedicated tests in those regimes.

minor comments (2)
  1. [Abstract] Abstract: the phrase 'retains computational efficiency beyond purely explicit schemes' would benefit from a quantitative statement of the speedup factor achieved in the reported tests.
  2. [Methods] The description of boundary handling in the super-time-stepping implementation should include a brief statement on how truncation errors are controlled, even if the tests demonstrate robustness.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript and the recommendation of minor revision. No specific major comments were raised in the report.

Circularity Check

0 steps flagged

No significant circularity; derivation anchored in external polynomial stability properties

full rationale

The paper presents a Runge-Kutta-Gegenbauer super-time-stepping scheme explicitly derived from the established stability properties of Gegenbauer polynomials (abstract: 'a super-time-stepping method based on the stability of the Gegenbauer polynomials'). This is a standard external mathematical fact, not constructed from the paper's own fitted results or prior self-citations. No load-bearing steps reduce to self-definition, parameter fitting renamed as prediction, or uniqueness theorems imported from the authors' own prior work. Benchmark tests and comparisons to reconnection/MRI problems provide independent verification outside any internal loop. The central efficiency/stability claims therefore remain self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract provides no explicit free parameters, axioms, or invented entities; the method rests on the mathematical stability of Gegenbauer polynomials and standard MHD equations.

pith-pipeline@v0.9.1-grok · 5817 in / 973 out tokens · 19668 ms · 2026-06-27T05:41:06.650955+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

95 extracted references · 32 linked inside Pith

  1. [1]

    Handbook of Mathematical Functions

  2. [2]

    and Amiez, Genevgve and Gremaud, Pierre , year =

    Alexiades, V. and Amiez, Genevgve and Gremaud, Pierre , year =. Super-time-stepping acceleration of explicit schemes for parabolic problems , volume =. Communications in Numerical Methods in Engineering , doi =

  3. [3]

    , year = 1942, month = oct, volume =

    Existence of Electromagnetic-Hydrodynamic Waves. , year = 1942, month = oct, volume =. doi:10.1038/150405d0 , adsurl =

    doi:10.1038/150405d0 1942

  4. [4]

    , keywords =

    Dynamics of Protoplanetary Disks. , keywords =. doi:10.1146/annurev-astro-081710-102521 , archivePrefix =. 1011.1496 , primaryClass =

    Pith/arXiv arXiv doi:10.1146/annurev-astro-081710-102521

  5. [5]

    , keywords =

    Effect of Ambipolar Diffusion on the Nonlinear Evolution of Magnetorotational Instability in Weakly Ionized Disks. , keywords =. doi:10.1088/0004-637X/736/2/144 , archivePrefix =. 1103.1380 , primaryClass =

    Pith/arXiv arXiv doi:10.1088/0004-637x/736/2/144

  6. [6]

    , keywords =

    Wave Propagation in Molecular Clouds. , keywords =. doi:10.1086/177462 , adsurl =

    doi:10.1086/177462

  7. [7]

    Barenblatt, G. I. , title =. Prikl.\ Mat.\ Mekh. , volume =

  8. [8]

    , keywords =

    Magnetic fields in spiral galaxies. , keywords =. doi:10.1007/s00159-015-0084-4 , archivePrefix =. 1509.04522 , primaryClass =

    Pith/arXiv arXiv doi:10.1007/s00159-015-0084-4

  9. [9]

    Journal of Computational Physics , keywords =

    A 4 ^ th -order accurate finite volume method for ideal classical and special relativistic MHD based on pointwise reconstructions. Journal of Computational Physics , keywords =. doi:10.1016/j.jcp.2023.112701 , archivePrefix =. 2310.11831 , primaryClass =

    doi:10.1016/j.jcp.2023.112701 2023

  10. [10]

    , keywords =

    2D or not 2D? Exploring 3D relativistic magnetic reconnection dynamics with highly accurate numerical simulations. , keywords =. doi:10.1093/mnras/stag286 , archivePrefix =. 2512.15954 , primaryClass =

    doi:10.1093/mnras/stag286

  11. [11]

    Journal of Computational Physics , year = 2008, month = mar, volume =

    An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws. Journal of Computational Physics , year = 2008, month = mar, volume =. doi:10.1016/j.jcp.2007.11.038 , adsurl =

    doi:10.1016/j.jcp.2007.11.038 2008

  12. [12]

    , keywords =

    Astrophysical magnetic fields and nonlinear dynamo theory. , keywords =. doi:10.1016/j.physrep.2005.06.005 , archivePrefix =. astro-ph/0405052 , primaryClass =

    Pith/arXiv arXiv doi:10.1016/j.physrep.2005.06.005 2005

  13. [13]

    Butcher, J. C. , biburl =. The Numerical analysis of ordinary differential equations : Runge-Kutta and general linear methods , url =

  14. [14]

    and Johnston, Craig D

    Caplan, Ronald M. and Johnston, Craig D. and Daldoff, Lars K. S. and Linker, Jon A. , title =. 2024 , month =. doi:10.1088/1742-6596/2742/1/012020 , url =

    doi:10.1088/1742-6596/2742/1/012020 2024

  15. [15]

    doi:10.1007/978-88-470-5280-2 , adsurl =

    Basics of Plasma Astrophysics. doi:10.1007/978-88-470-5280-2 , adsurl =

    doi:10.1007/978-88-470-5280-2

  16. [16]

    , keywords =

    An Explicit Scheme for Incorporating Ambipolar Diffusion in a Magnetohydrodynamics Code. , keywords =. doi:10.1088/0067-0049/181/2/413 , archivePrefix =. 0812.3748 , primaryClass =

    Pith/arXiv arXiv doi:10.1088/0067-0049/181/2/413

  17. [17]

    Mathematische Annalen , year = 1928, month = jan, volume =

    \"U ber die partiellen Differenzengleichungen der mathematischen Physik. Mathematische Annalen , year = 1928, month = jan, volume =. doi:10.1007/BF01448839 , adsurl =

    doi:10.1007/bf01448839 1928

  18. [18]

    , year = 2012, month = sep, volume =

    Magnetic Fields in Molecular Clouds. , year = 2012, month = sep, volume =. doi:10.1146/annurev-astro-081811-125514 , adsurl =

    doi:10.1146/annurev-astro-081811-125514 2012

  19. [19]

    Journal of Computational Physics , year = 2002, month = jan, volume = 175, pages =

    Hyperbolic Divergence Cleaning for the MHD Equations. Journal of Computational Physics , year = 2002, month = jan, volume = 175, pages =. doi:10.1006/jcph.2001.6961 , adsurl =

    doi:10.1006/jcph.2001.6961 2002

  20. [20]

    1991 , issn =

    The coefficients of differentiated expansions and derivatives of ultraspherical polynomials , journal =. 1991 , issn =. doi:https://doi.org/10.1016/0898-1221(91)90089-M , url =

    doi:10.1016/0898-1221(91)90089-m 1991

  21. [21]

    , keywords =

    Simulating hydromagnetic processes in star formation: introducing ambipolar diffusion into an adaptive mesh refinement code. , keywords =. doi:10.1111/j.1365-2966.2008.14026.x , archivePrefix =. 0810.0299 , primaryClass =

    doi:10.1111/j.1365-2966.2008.14026.x 2008

  22. [22]

    Higher Transcendental Functions , year = 1953, pages =

    Higher Transcendental Functions. Higher Transcendental Functions , year = 1953, pages =

    1953

  23. [23]

    , keywords =

    Collimation of Astrophysical Jets: The Role of the Accretion Disk Magnetic Field Distribution. , keywords =. doi:10.1086/507976 , archivePrefix =. astro-ph/0511611 , primaryClass =

    Pith/arXiv arXiv doi:10.1086/507976

  24. [24]

    Magnetically disrupted discs and surface accretion

    A global 3D simulation of magnetospheric accretion - I. Magnetically disrupted discs and surface accretion. , keywords =. doi:10.1093/mnras/stad3712 , archivePrefix =. 2309.15318 , primaryClass =

    doi:10.1093/mnras/stad3712

  25. [25]

    MHD simulations of the magnetorotational instability in a shearing box with zero net flux. I. The issue of convergence. , keywords =. doi:10.1051/0004-6361:20077942 , archivePrefix =. 0705.3621 , primaryClass =

    Pith/arXiv arXiv doi:10.1051/0004-6361:20077942

  26. [26]

    , keywords =

    Array programming with NumPy. , keywords =. doi:10.1038/s41586-020-2649-2 , archivePrefix =. 2006.10256 , primaryClass =

    Pith/arXiv arXiv doi:10.1038/s41586-020-2649-2 2006

  27. [27]

    A Powerful Local Shear Instability in Weakly Magnetized Disks. II. Nonlinear Evolution. , keywords =. doi:10.1086/170271 , adsurl =

    doi:10.1086/170271

  28. [28]

    Computing in Science and Engineering , keywords =

    Matplotlib: A 2D Graphics Environment. Computing in Science and Engineering , keywords =. doi:10.1109/MCSE.2007.55 , adsurl =

    doi:10.1109/mcse.2007.55 2007

  29. [29]

    A Powerful Local Shear Instability in Weakly Magnetized Disks. I. Linear Analysis. , keywords =. doi:10.1086/170270 , adsurl =

    doi:10.1086/170270

  30. [30]

    , keywords =

    3D Radiation Nonideal Magnetohydrodynamical Simulations of the Inner Rim in Protoplanetary Disks. , keywords =. doi:10.3847/1538-4357/835/2/230 , archivePrefix =. 1612.02740 , primaryClass =

    Pith/arXiv arXiv doi:10.3847/1538-4357/835/2/230

  31. [31]

    , keywords =

    Radiation magnetohydrodynamics in global simulations of protoplanetary discs. , keywords =. doi:10.1051/0004-6361/201322451 , archivePrefix =. 1310.5865 , primaryClass =

    Pith/arXiv arXiv doi:10.1051/0004-6361/201322451

  32. [32]

    SIAM Review , keywords =

    Strong Stability-Preserving High-Order Time Discretization Methods. SIAM Review , keywords =. doi:10.1137/S003614450036757X , adsurl =

    doi:10.1137/s003614450036757x

  33. [33]

    , keywords =

    Global Simulations of Protoplanetary Disks With Ohmic Resistivity and Ambipolar Diffusion. , keywords =. doi:10.1088/0004-637X/801/2/84 , archivePrefix =. 1501.05431 , primaryClass =

    Pith/arXiv arXiv doi:10.1088/0004-637x/801/2/84

  34. [34]

    Proceedings of the Royal Society of London Series A , year = 1982, month = sep, volume =

    Eigenvalues of the Barenblatt-Pattle Similarity Solution in Nonlinear Diffusion. Proceedings of the Royal Society of London Series A , year = 1982, month = sep, volume =. doi:10.1098/rspa.1982.0122 , adsurl =

    doi:10.1098/rspa.1982.0122 1982

  35. [35]

    Hairer, E. and N. Solving Ordinary Differential Equations

  36. [36]

    and Lax, P

    Harten, A. and Lax, P. and Leer, B. , title =. SIAM Review , volume =. 1983 , doi =

    1983

  37. [37]

    1992 , issn =

    On the coefficients of differentiated expansions of ultraspherical polynomials , journal =. 1992 , issn =. doi:https://doi.org/10.1016/0168-9274(92)90010-B , url =

    doi:10.1016/0168-9274(92)90010-b 1992

  38. [38]

    , keywords =

    Energy distribution and substructure formation in astrophysical MHD simulations. , keywords =. doi:10.1093/mnras/stad3807 , archivePrefix =. 2308.16062 , primaryClass =

    doi:10.1093/mnras/stad3807

  39. [39]

    , keywords =

    Ambipolar diffusion in the magnetorotational instability. , keywords =. doi:10.1111/j.1365-2966.2004.07383.x , archivePrefix =. astro-ph/0309707 , primaryClass =

    doi:10.1111/j.1365-2966.2004.07383.x 2004

  40. [40]

    High-order Godunov schemes for global 3D MHD simulations of accretion disks. I. Testing the linear growth of the magneto-rotational instability. , keywords =. doi:10.1051/0004-6361/200912443 , archivePrefix =. 0906.5516 , primaryClass =

    Pith/arXiv arXiv doi:10.1051/0004-6361/200912443

  41. [41]

    The combined influence of Ohmic, Hall, and ambipolar diffusion on dead zones

    Thanatology in protoplanetary discs. The combined influence of Ohmic, Hall, and ambipolar diffusion on dead zones. , keywords =. doi:10.1051/0004-6361/201423660 , archivePrefix =. 1402.4133 , primaryClass =

    Pith/arXiv arXiv doi:10.1051/0004-6361/201423660

  42. [42]

    , keywords =

    Gap opening in protoplanetary discs: gas dynamics from global axisymmetric non-ideal MHD simulations with consistent thermochemistry. , keywords =. doi:10.1093/mnras/stad1632 , archivePrefix =. 2304.05972 , primaryClass =

    doi:10.1093/mnras/stad1632

  43. [43]

    , keywords =

    Global three-dimensional simulations of outer protoplanetary discs with ambipolar diffusion. , keywords =. doi:10.1093/mnras/stab2220 , archivePrefix =. 2106.10167 , primaryClass =

    doi:10.1093/mnras/stab2220

  44. [44]

    Global simulations of protoplanetary disks with net magnetic flux. I. Non-ideal MHD case. , keywords =. doi:10.1051/0004-6361/201630056 , archivePrefix =. 1612.00883 , primaryClass =

    Pith/arXiv arXiv doi:10.1051/0004-6361/201630056

  45. [45]

    Turbulence in the Outer Regions of Protoplanetary Disks. I. Weak Accretion with No Vertical Magnetic Flux. , keywords =. doi:10.1088/0004-637X/764/1/66 , archivePrefix =. 1210.4164 , primaryClass =

    Pith/arXiv arXiv doi:10.1088/0004-637x/764/1/66

  46. [46]

    , keywords =

    Exact shearing box solutions of magnetohydrodynamic flows with resistivity, viscosity and cooling. , keywords =. doi:10.1111/j.1365-2966.2007.12270.x , archivePrefix =. 0709.1388 , primaryClass =

    doi:10.1111/j.1365-2966.2007.12270.x 2007

  47. [47]

    Protostars and Planets VII , year = 2023, editor =

    Hydro-, Magnetohydro-, and Dust-Gas Dynamics of Protoplanetary Disks. Protostars and Planets VII , year = 2023, editor =. doi:10.48550/arXiv.2203.09821 , archivePrefix =. 2203.09821 , primaryClass =

    doi:10.48550/arxiv.2203.09821 2023

  48. [48]

    , keywords =

    IDEFIX: A versatile performance-portable Godunov code for astrophysical flows. , keywords =. doi:10.1051/0004-6361/202346005 , archivePrefix =. 2304.13746 , primaryClass =

    doi:10.1051/0004-6361/202346005

  49. [49]

    Journal of Computational Physics , keywords =

    On the divergence-free condition in Godunov-type schemes for ideal magnetohydrodynamics: the upwind constrained transport method. Journal of Computational Physics , keywords =. doi:10.1016/j.jcp.2003.09.016 , archivePrefix =. astro-ph/0310183 , primaryClass =

    Pith/arXiv arXiv doi:10.1016/j.jcp.2003.09.016 2003

  50. [50]

    , keywords =

    Chemical solver to compute molecule and grain abundances and non-ideal MHD resistivities in prestellar core-collapse calculations. , keywords =. doi:10.1051/0004-6361/201526780 , archivePrefix =. 1604.05613 , primaryClass =

    Pith/arXiv arXiv doi:10.1051/0004-6361/201526780

  51. [51]

    Plasma Physics and Controlled Fusion , keywords =

    Magnetic reconnection: from the Sweet-Parker model to stochastic plasmoid chains. Plasma Physics and Controlled Fusion , keywords =. doi:10.1088/0741-3335/58/1/014021 , archivePrefix =. 1507.07756 , primaryClass =

    Pith/arXiv arXiv doi:10.1088/0741-3335/58/1/014021

  52. [52]

    , keywords =

    Non-ideal magnetohydrodynamics on a moving mesh. , keywords =. doi:10.1093/mnras/sty397 , archivePrefix =. 1710.10265 , primaryClass =

    Pith/arXiv arXiv doi:10.1093/mnras/sty397

  53. [53]

    , keywords =

    Incorporating Ambipolar and Ohmic Diffusion in the AMR MHD Code RAMSES. , keywords =. doi:10.1088/0067-0049/201/2/24 , archivePrefix =. 1206.2476 , primaryClass =

    Pith/arXiv arXiv doi:10.1088/0067-0049/201/2/24

  54. [54]

    Ambipolar diffusion in low-mass star formation. I. General comparison with the ideal magnetohydrodynamic case. , keywords =. doi:10.1051/0004-6361/201526371 , archivePrefix =. 1509.05630 , primaryClass =

    Pith/arXiv arXiv doi:10.1051/0004-6361/201526371

  55. [55]

    The Journal of Open Source Software , keywords =

    PyPLUTO: a data analysis Python package for the PLUTO code. The Journal of Open Source Software , keywords =. doi:10.21105/joss.08448 , archivePrefix =. 2501.09748 , primaryClass =

    doi:10.21105/joss.08448

  56. [56]

    , keywords =

    Theory of Star Formation. , keywords =. doi:10.1146/annurev.astro.45.051806.110602 , archivePrefix =. 0707.3514 , primaryClass =

    Pith/arXiv arXiv doi:10.1146/annurev.astro.45.051806.110602

  57. [57]

    , keywords =

    A second-order accurate Super TimeStepping formulation for anisotropic thermal conduction. , keywords =. doi:10.1111/j.1365-2966.2012.20744.x , adsurl =

    doi:10.1111/j.1365-2966.2012.20744.x 2012

  58. [58]

    Journal of Computational Physics , year = 2014, month = jan, volume =

    A stabilized Runge-Kutta-Legendre method for explicit super-time-stepping of parabolic and mixed equations. Journal of Computational Physics , year = 2014, month = jan, volume =. doi:10.1016/j.jcp.2013.08.021 , adsurl =

    doi:10.1016/j.jcp.2013.08.021 2014

  59. [59]

    and Bodo, G

    Mignone, A. and Bodo, G. and Massaglia, S. and Matsakos, T. and Tesileanu, O. and Zanni, C. and Ferrari, A. , year =. The Astrophysical Journal Supplement Series , volume =. doi:10.1086/513316 , urldate =

    doi:10.1086/513316

  60. [60]

    , keywords =

    A conservative orbital advection scheme for simulations of magnetized shear flows with the PLUTO code. , keywords =. doi:10.1051/0004-6361/201219557 , archivePrefix =. 1207.2955 , primaryClass =

    Pith/arXiv arXiv doi:10.1051/0004-6361/201219557

  61. [61]

    Journal of Computational Physics , keywords =

    High-order conservative reconstruction schemes for finite volume methods in cylindrical and spherical coordinates. Journal of Computational Physics , keywords =. doi:10.1016/j.jcp.2014.04.001 , archivePrefix =. 1404.0537 , primaryClass =

    Pith/arXiv arXiv doi:10.1016/j.jcp.2014.04.001 2014

  62. [62]

    Journal of Computational Physics , keywords =

    Systematic construction of upwind constrained transport schemes for MHD. Journal of Computational Physics , keywords =. doi:10.1016/j.jcp.2020.109748 , archivePrefix =. 2004.10542 , primaryClass =

    doi:10.1016/j.jcp.2020.109748 2020

  63. [63]

    2005", month =

    A multi-state HLL approximate Riemann solver for ideal magnetohydrodynamics. Journal of Computational Physics , year = "2005", month = "Sep", volume =. doi:10.1016/j.jcp.2005.02.017 , adsurl =

    doi:10.1016/j.jcp.2005.02.017 2005

  64. [64]

    Nonhomologous contraction and equilibria of self-gravitating, magnetic interstellar clouds embedded in an intercloud medium: star formation. II. Results. , keywords =. doi:10.1086/154478 , adsurl =

    doi:10.1086/154478

  65. [65]

    , keywords =

    Ambipolar diffusion in the Bifrost code. , keywords =. doi:10.1051/0004-6361/202037809 , archivePrefix =. 2004.11927 , primaryClass =

    doi:10.1051/0004-6361/202037809 2004

  66. [66]

    Journal of Computational Physics , keywords =

    Runge-Kutta-Gegenbauer explicit methods for advection-diffusion problems. Journal of Computational Physics , keywords =. doi:10.1016/j.jcp.2019.03.001 , adsurl =

    doi:10.1016/j.jcp.2019.03.001 2019

  67. [67]

    Their origin and their activity

    Cosmical magnetic fields. Their origin and their activity

  68. [68]

    Protostars and Planets V , year = 2007, editor =

    Disk Winds, Jets, and Outflows: Theoretical and Computational Foundations. Protostars and Planets V , year = 2007, editor =. doi:10.48550/arXiv.astro-ph/0603592 , archivePrefix =. astro-ph/0603592 , primaryClass =

    Pith/arXiv arXiv doi:10.48550/arxiv.astro-ph/0603592 2007

  69. [69]

    The MHD case

    On the impact of the numerical method on magnetic reconnection and particle acceleration - I. The MHD case. , keywords =. doi:10.1093/mnras/stab2813 , archivePrefix =. 2109.12858 , primaryClass =

    doi:10.1093/mnras/stab2813

  70. [70]

    , keywords =

    The impact of resistive electric fields on particle acceleration in reconnection layers. , keywords =. doi:10.1093/mnras/stac2807 , archivePrefix =. 2210.01113 , primaryClass =

    doi:10.1093/mnras/stac2807

  71. [71]

    Journal of Computational Physics , year = 1981, month = oct, volume =

    Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes. Journal of Computational Physics , year = 1981, month = oct, volume =. doi:10.1016/0021-9991(81)90128-5 , adsurl =

    doi:10.1016/0021-9991(81)90128-5 1981

  72. [72]

    Rossazza and A

    M. Rossazza and A. Mignone and M. Bugli and S. Truzzi and L. Riha and T. Panoc and O. Vysocky and N. Shukla and A. Romeo and V. Berta , keywords =. The PLUTO code on GPUs: A first look at Eulerian MHD methods , journal =. 2026 , issn =. doi:https://doi.org/10.1016/j.ascom.2026.101076 , url =

    doi:10.1016/j.ascom.2026.101076 2026

  73. [73]

    , keywords =

    Star formation in molecular clouds: observation and theory. , keywords =. doi:10.1146/annurev.aa.25.090187.000323 , adsurl =

    doi:10.1146/annurev.aa.25.090187.000323

  74. [74]

    Journal of the British Astronomical Association , year = 1992, month = aug, volume =

    Books-Received - the Physics of Astrophysics - V.2 - Gas Dynamics. Journal of the British Astronomical Association , year = 1992, month = aug, volume =

    1992

  75. [75]

    Journal of Computational Physics , keywords =

    Super-time-stepping schemes for parabolic equations with boundary conditions. Journal of Computational Physics , keywords =. doi:10.1016/j.jcp.2020.109879 , adsurl =

    doi:10.1016/j.jcp.2020.109879 2020

  76. [76]

    doi:10.1002/9783527617722 , adsurl =

    Physical processes in the interstellar medium. doi:10.1002/9783527617722 , adsurl =

    doi:10.1002/9783527617722

  77. [77]

    , keywords =

    The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers. , keywords =. doi:10.3847/1538-4365/ab929b , archivePrefix =. 2005.06651 , primaryClass =

    doi:10.3847/1538-4365/ab929b 2005

  78. [78]

    , keywords =

    Radiation Magnetohydrodynamic Simulations of Protostellar Collapse: Protostellar Core Formation. , keywords =. doi:10.1088/0004-637X/763/1/6 , archivePrefix =. 1206.3567 , primaryClass =

    Pith/arXiv arXiv doi:10.1088/0004-637x/763/1/6

  79. [79]

    , keywords =

    Radiation Magnetohydrodynamic Simulations of Protostellar Collapse: Nonideal Magnetohydrodynamic Effects and Early Formation of Circumstellar Disks. , keywords =. doi:10.1088/0004-637X/801/2/117 , archivePrefix =. 1501.04102 , primaryClass =

    Pith/arXiv arXiv doi:10.1088/0004-637x/801/2/117

  80. [80]

    Riemann Solvers and Numerical Methods for Fluid Dynamics: A Practical Introduction , journal =

    Toro, Eleuterio , year =. Riemann Solvers and Numerical Methods for Fluid Dynamics: A Practical Introduction , journal =

Showing first 80 references.