pith. sign in

arxiv: 2511.15937 · v3 · pith:LPF5L3TQnew · submitted 2025-11-20 · ⚛️ physics.plasm-ph

Asymptotic-preserving semi-implicit finite volume scheme for Extended Magnetohydrodynamics

Yi Han Toh , Joshua Dolence , Karthik Duraisamy This is my paper

Pith reviewed 2026-05-25 07:24 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph
keywords extended magnetohydrodynamicsfinite volume schemesemi-implicit methodasymptotic preservingHall MHDresistive MHDplasma simulationadaptive mesh refinement
0
0 comments X

The pith

A finite volume scheme for extended magnetohydrodynamics yields accurate solutions in ideal, resistive, and Hall limits by rewriting the equations to retain ideal MHD solvers.

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

The paper develops a finite volume method for the extended magnetohydrodynamic equations that include electron inertia and displacement current. It first rewrites the equations so that standard ideal MHD Riemann solvers and constrained transport for divergence-free magnetic fields can still be used. The added stiffness is treated with a semi-implicit approach that casts the system as a relaxation problem, advanced by explicit second-order Runge-Kutta time stepping with operator splitting on the implicit terms and a density-dependent slope limiter for stability at low densities. The resulting scheme is placed in an adaptive mesh refinement framework and tested on reference problems, where it naturally recovers the ideal MHD limit and performs well in the resistive and Hall limits.

Core claim

The authors construct a semi-implicit finite volume scheme for the extended magnetohydrodynamic equations by first rewriting them to preserve the use of ideal MHD Riemann solvers and the constrained transport method. The model is recast as a relaxation system whose implicit source terms are handled by operator splitting inside an explicit second-order Runge-Kutta integrator. A density-dependent slope limiter augments numerical diffusivity at low densities. Implemented inside a scalable adaptive mesh refinement code, the algorithm is shown on standard test problems to produce accurate results in the ideal, resistive, and Hall MHD regimes while asymptoting to the ideal MHD limit without extra

What carries the argument

Reformulation of the XMHD equations as a relaxation system that retains ideal MHD Riemann solvers and constrained transport while treating electron inertia and displacement current through semi-implicit operator-split updates.

If this is right

  • The scheme asymptotes naturally to the ideal MHD limit without additional modifications.
  • It produces accurate results on reference test problems for resistive and Hall MHD.
  • The density-dependent slope limiter improves stability where non-ideal effects are strong.
  • The method supports scalable adaptive mesh refinement simulations of plasmas.

Where Pith is reading between the lines

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

  • A single implementation could simulate plasmas that transition between MHD regimes without manual switching of solvers.
  • The approach may be examined on problems with strong density gradients to test whether the limiter choice remains robust.
  • Extension to fully three-dimensional configurations with complex boundaries could expose additional stability constraints not visible in the presented tests.

Load-bearing premise

Rewriting the extended magnetohydrodynamic equations permits retention of ideal MHD Riemann solvers and constrained transport without introducing errors that invalidate accuracy claims in the resistive and Hall limits.

What would settle it

Comparison of the scheme's output on a standard Hall MHD reference test against an independent high-resolution solution; large deviation in magnetic field or current density would falsify the accuracy claim in that limit.

Figures

Figures reproduced from arXiv: 2511.15937 by Joshua Dolence, Karthik Duraisamy, Yi Han Toh.

Figure 1
Figure 1. Figure 1: A 3D computational cell on a Cartesian grid is illustrated to indicate the positions of various field [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: 2D diagram illustrating the faces F3 of six computational cells, with their normals directed into the paper along the z-axis. In this configuration, the faces F3(i, j, k − 1 2 ) and F3(i, j − 1, k − 1 2 ) correspond to cell centers at (i, j, k) and (i, j − 1, k), respectively. Each computational cell is bounded by edges that define its geometry: the horizontal arrows represent the edges E1, while the verti… view at source ↗
Figure 3
Figure 3. Figure 3: Brio-Wu shock tube solution from Zhao et. al. [16] using LLF Riemann solver at cell interface [PITH_FULL_IMAGE:figures/full_fig_p020_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: For the standard ideal MHD test case on shock tube, the XMHD (Artemis) is quite similar to ideal [PITH_FULL_IMAGE:figures/full_fig_p021_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: At low density, XMHD (Artemis) begins to deviate significantly from ideal MHD (Athena++) due [PITH_FULL_IMAGE:figures/full_fig_p021_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: For 1D problem at grid resolution of [16, 32, 64, 128, 256, 512, 1024]. 23 [PITH_FULL_IMAGE:figures/full_fig_p023_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: For 2D problem, the code is tested by setting both Hall and resistive terms to zero. 2nd-order [PITH_FULL_IMAGE:figures/full_fig_p024_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: To test for order of convergence in 1D field diffusion problem, several grid resolution and their L1 [PITH_FULL_IMAGE:figures/full_fig_p025_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Results for Huba’s Hall drift problem obtained from XMHD implemented on Artemis and compared [PITH_FULL_IMAGE:figures/full_fig_p028_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Default 2D ideal MHD blast wave problem with slanted magnetic field along x-y direction. [PITH_FULL_IMAGE:figures/full_fig_p029_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: 2D MHD blast wave problem at lower density of [PITH_FULL_IMAGE:figures/full_fig_p030_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Default 2D ideal MHD rotor problem with ρ0 = 1. 32 [PITH_FULL_IMAGE:figures/full_fig_p032_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Default 2D ideal MHD rotor problem with ρ0 = 0.1 when non-ideal effects become more significant. 33 [PITH_FULL_IMAGE:figures/full_fig_p033_13.png] view at source ↗
read the original abstract

A Finite Volume (FV) scheme is developed for solving the extended magnetohydrodynamic (XMHD) equations, yielding accurate results in the ideal, resistive, and Hall MHD limits. This is accomplished by first re-writing the XMHD equations such that it allows the algorithm to retain the use of ideal MHD Riemann solvers and the constrained transport method to preserve divergence-free magnetic fields. Incorporation of electron inertia and displacement current introduces additional numerical stiffness which motivates a semi-implicit FV scheme that re-formulates the XMHD model as a relaxation system. The equations are then advanced in time using an explicit 2nd-order Runge-Kutta scheme with operator splitting applied to the implicit source term updates at each sub-stage. For additional numerical stability, a density-dependent slope limiter is implemented to increase flux diffusivity at low density regions where non-ideal effects become significant. The algorithm is subsequently implemented in a scalable adaptive mesh refinement (AMR) framework. As the new algorithm retains many aspects of the ideal MHD formulations, it asymptotes naturally to the ideal MHD limit. Moreover, it shows promising results at the resistive and Hall MHD limits. This is verified against reference test problems for ideal, resistive and Hall MHD.

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

3 major / 1 minor

Summary. The manuscript develops a finite volume scheme for the extended magnetohydrodynamic (XMHD) equations. The approach rewrites the XMHD system to retain ideal MHD Riemann solvers and constrained transport, reformulates electron inertia and displacement current as a relaxation system, applies a semi-implicit time advance with explicit second-order Runge-Kutta and operator splitting on implicit sources, and adds a density-dependent slope limiter for low-density stability. The scheme is implemented in an AMR framework and is asserted to asymptote naturally to the ideal MHD limit while yielding promising results in resistive and Hall limits, as verified on reference test problems.

Significance. If the reformulation and splitting preserve consistency without O(1) errors when resistivity or Hall parameter are finite, the method would allow reuse of mature ideal-MHD infrastructure for multi-regime plasma problems and could be practically useful in AMR codes. The relaxation-system treatment of stiffness and the AMR implementation are positive engineering features, but the absence of quantitative verification leaves the practical significance unclear.

major comments (3)
  1. [Abstract] Abstract: the verification statement supplies no quantitative error measures, convergence rates, or L1/L2 norms against reference solutions in the resistive or Hall regimes; without these, the central claim that the scheme remains accurate when non-ideal terms are active cannot be assessed.
  2. [Abstract] Abstract / reformulation description: the rewriting that permits retention of ideal MHD Riemann solvers is presented without derivation of the modified numerical fluxes or a consistency analysis of the operator-split scheme for finite resistivity and Hall parameter; this step is load-bearing for the assertion that no O(1) errors are introduced in the non-ideal limits.
  3. [Abstract] Abstract: the density-dependent slope limiter is introduced for stability at low density, yet no analysis of its additional diffusivity on the resistive or Hall terms is given, leaving an uncontrolled source of error unquantified in the regimes where the scheme is claimed to be promising.
minor comments (1)
  1. [Abstract] The abstract would be clearer if it named the specific reference test problems and the order of the scheme.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful review and constructive comments. We address each major comment below and will make revisions to strengthen the manuscript's clarity and rigor.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the verification statement supplies no quantitative error measures, convergence rates, or L1/L2 norms against reference solutions in the resistive or Hall regimes; without these, the central claim that the scheme remains accurate when non-ideal terms are active cannot be assessed.

    Authors: We agree the abstract would benefit from quantitative support. The manuscript body (Sections 4.2–4.3) already contains L1/L2 error norms and observed convergence rates for the resistive and Hall test problems. We will revise the abstract to include a concise summary of these metrics. revision: yes

  2. Referee: [Abstract] Abstract / reformulation description: the rewriting that permits retention of ideal MHD Riemann solvers is presented without derivation of the modified numerical fluxes or a consistency analysis of the operator-split scheme for finite resistivity and Hall parameter; this step is load-bearing for the assertion that no O(1) errors are introduced in the non-ideal limits.

    Authors: The reformulation, modified fluxes, and consistency analysis for finite resistivity/Hall parameters (showing no O(1) splitting errors) are derived in Sections 2.1 and 3.3. We will add one sentence to the abstract referencing these sections and the key consistency result. revision: yes

  3. Referee: [Abstract] Abstract: the density-dependent slope limiter is introduced for stability at low density, yet no analysis of its additional diffusivity on the resistive or Hall terms is given, leaving an uncontrolled source of error unquantified in the regimes where the scheme is claimed to be promising.

    Authors: We agree an explicit analysis of the limiter's diffusivity on non-ideal terms is warranted. We will add a short discussion (with supporting estimates and test data) in Section 3.4 quantifying this effect in the resistive and Hall regimes. revision: yes

Circularity Check

0 steps flagged

No circularity: algorithmic construction verified externally on tests

full rationale

The paper constructs a semi-implicit FV scheme by rewriting XMHD equations to reuse ideal MHD Riemann solvers and constrained transport, then applies a relaxation-system treatment of stiffness with operator-split implicit updates and a density-dependent limiter. These are standard numerical techniques assembled into an algorithm. Claims of asymptotic preservation and accuracy in ideal/resistive/Hall limits rest on verification against independent reference test problems, not on any definitional equivalence, fitted-parameter renaming, or self-citation chain. No load-bearing step reduces to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Abstract-only review yields no explicit free parameters, invented entities, or non-standard axioms; the work implicitly rests on standard assumptions of finite-volume methods and operator splitting.

axioms (2)
  • standard math Standard finite-volume discretization combined with constrained transport preserves divergence-free magnetic fields to machine precision when ideal MHD Riemann solvers are retained.
    Invoked when the authors state that the rewriting allows retention of ideal MHD solvers and constrained transport.
  • domain assumption Operator splitting between explicit transport and implicit source terms in a second-order Runge-Kutta scheme maintains overall second-order accuracy for the relaxation system.
    Implicit in the time-advancement description.

pith-pipeline@v0.9.0 · 5746 in / 1446 out tokens · 40266 ms · 2026-05-25T07:24:49.846576+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

38 extracted references · 38 canonical work pages

  1. [1]

    Extended- magnetohydrodynamics in under-dense plasmas,

    C. A. Walsh, J. P. Chittenden, D. W. Hill, and C. Ridgers, “Extended- magnetohydrodynamics in under-dense plasmas,”Physics of Plasmas, vol. 27, p.022103, 02 2020

    work page 2020

  2. [2]

    Extended magnetohydrodynamics simulations of thin-foil z-pinch implosions with comparison to experiments,

    J. M. Woolstrum, D. A. Yager-Elorriaga, P. C. Campbell, N. M. Jordan, C. E. Seyler, and R. D. McBride, “Extended magnetohydrodynamics simulations of thin-foil z-pinch implosions with comparison to experiments,”Physics of Plasmas, vol. 27, p. 092705, 09 2020. 37 Figure C.14: Atn 0 = 10 30 where both ions and electrons temperatures are expected to equilibra...

    work page 2020

  3. [3]

    Using extended mhd to explore lasers as a trigger for x-pinches,

    J. R. Young, M. B. Adams, H. Hasson, I. West-Abdallah, M. Evans, and P.-A. Gourdain, “Using extended mhd to explore lasers as a trigger for x-pinches,”Physics of Plasmas, vol. 28, p. 102703, 10 2021

    work page 2021

  4. [4]

    Extended-MHD simulations of disruption mitigation via massive gas injection in SPARC,

    A. Kleiner, N. M. Ferraro, R. Sweeney, B. C. Lyons, and M. Reinke, “Extended-MHD simulations of disruption mitigation via massive gas injection in SPARC,”Nuclear Fu- sion, vol. 65, no. 2, p. 026015, 2025. Published 30 December 2024.©2024 The Au- thor(s). Open access

    work page 2025

  5. [5]

    Power flow in pulsed-power systems: The influence of hall physics and modeling of the plasma–vacuum interface,

    N. D. Hamlin and C. E. Seyler, “Power flow in pulsed-power systems: The influence of hall physics and modeling of the plasma–vacuum interface,”IEEE Transactions on Plasma Science, vol. 47, no. 5, pp. 2064–2073, 2019

    work page 2064

  6. [6]

    R. D. Hazeltine,The Framework of Plasma Physics. Boca Raton: CRC Press, 1st ed.,

    work page

  7. [7]

    eBook published March 7, 2018

    work page 2018

  8. [8]

    Krall and A

    N. Krall and A. Trivelpiece,Principles of Plasma Physics. International series in pure and applied physics, McGraw-Hill, 1973

    work page 1973

  9. [9]

    Spitzer,Physics of Fully Ionized Gases

    L. Spitzer,Physics of Fully Ionized Gases. 1956

    work page 1956

  10. [10]

    A 2d high-βhall mhd implicit nonlinear solver,

    L. Chacón and D. Knoll, “A 2d high-βhall mhd implicit nonlinear solver,”Journal of Computational Physics, vol. 188, no. 2, pp. 573–592, 2003

    work page 2003

  11. [11]

    A high resolution wave propagation scheme for ideal two-fluid plasma equations,

    A. Hakim, J. Loverich, and U. Shumlak, “A high resolution wave propagation scheme for ideal two-fluid plasma equations,”Journal of Computational Physics, vol. 219, no. 1, pp. 418–442, 2006

    work page 2006

  12. [12]

    Extended mhd modelling with the ten-moment equations,

    A. H. Hakim, “Extended mhd modelling with the ten-moment equations,”Journal of Fusion Energy, vol. 27, no. 1, pp. 36–43, 2008

    work page 2008

  13. [13]

    Numerical methods: Ideal and hall mhd,

    J. D. Huba, “Numerical methods: Ideal and hall mhd,” inProceedings of ISSS, vol. 7, pp. 26–31, March 2005

    work page 2005

  14. [14]

    A discontinuous galerkin method for the full two-fluid plasma model,

    J. Loverich and U. Shumlak, “A discontinuous galerkin method for the full two-fluid plasma model,”Computer Physics Communications, vol. 169, no. 1, pp. 251–255, 2005. Proceedings of the Europhysics Conference on Computational Physics 2004

    work page 2005

  15. [15]

    Approximate riemann solver for the two-fluid plasma model,

    U. Shumlak and J. Loverich, “Approximate riemann solver for the two-fluid plasma model,”Journal of Computational Physics, vol. 187, no. 2, pp. 620–638, 2003

    work page 2003

  16. [16]

    A new 3d mhd algorithm: the distribution function method,

    J. D. Huba and J. G. Lyon, “A new 3d mhd algorithm: the distribution function method,”Journal of Plasma Physics, vol. 61, no. 3, pp. 391–405, 1999. Received 20 May 1998; revised 17 November 1998

    work page 1999

  17. [17]

    A positivity-preserving semi-implicit discontinuous galerkin scheme for solving extended magnetohydrodynamics equations,

    X. Zhao, Y. Yang, and C. E. Seyler, “A positivity-preserving semi-implicit discontinuous galerkin scheme for solving extended magnetohydrodynamics equations,”Journal of Computational Physics, vol. 278, pp. 400–415, 2014. 40

    work page 2014

  18. [18]

    Relaxation model for extended magnetohydrodynamics: Comparisontomagnetohydrodynamicsfordensez-pinches,

    C. E. Seyler and M. R. Martin, “Relaxation model for extended magnetohydrodynamics: Comparisontomagnetohydrodynamicsfordensez-pinches,”Physics of Plasmas, vol.18, p. 012703, 01 2011

    work page 2011

  19. [19]

    Derivation and numerical approximation of two-temperature euler plasma model,

    A. Sangam, Élise Estibals, and H. Guillard, “Derivation and numerical approximation of two-temperature euler plasma model,”Journal of Computational Physics, vol. 444, p. 110565, 2021

    work page 2021

  20. [20]

    Parthenon—a performance portable block-structured adap- tive mesh refinement framework,

    P.Grete, J.C.Dolence, J.M.Miller, J.Brown, B.Ryan, A.Gaspar, F.Glines, S.Swami- narayan, J. Lippuner, C. J. Solomon, G. Shipman, C. Junghans, D. Holladay, J. M. Stone, and L. F. Roberts, “Parthenon—a performance portable block-structured adap- tive mesh refinement framework,”International Journal of High Performance Comput- ing Applications, vol. 37, 12 2022

    work page 2022

  21. [21]

    The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers,

    J. M. Stone, K. Tomida, C. J. White, and K. G. Felker, “The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers,”The Astrophysical Journal Supplement Series, vol. 249, p. 4, July 2020

    work page 2020

  22. [22]

    A staggered mesh algorithm using high order godunov fluxes to ensure solenoidal magnetic fields in magnetohydrodynamic simulations,

    D. S. Balsara and D. S. Spicer, “A staggered mesh algorithm using high order godunov fluxes to ensure solenoidal magnetic fields in magnetohydrodynamic simulations,”Jour- nal of Computational Physics, vol. 149, no. 2, pp. 270–292, 1999

    work page 1999

  23. [23]

    Weaksolutionsofnonlinearhyperbolicequationsandtheirnumericalcompu- tation,

    P.D.Lax, “Weaksolutionsofnonlinearhyperbolicequationsandtheirnumericalcompu- tation,”Communications on Pure and Applied Mathematics, vol. 7, no. 1, pp. 159–193, 1954

    work page 1954

  24. [24]

    A multi-state HLL approximate Riemann solver for ideal magnetohydrodynamics,

    T. Miyoshi and K. Kusano, “A multi-state HLL approximate Riemann solver for ideal magnetohydrodynamics,”Journal of Computational Physics, vol. 208, pp. 315–344, Sept. 2005

    work page 2005

  25. [25]

    L. N. Trefethen and D. Bau,Numerical Linear Algebra, Twenty-fifth Anniversary Edi- tion. Philadelphia, PA: Society for Industrial and Applied Mathematics, 2022

    work page 2022

  26. [26]

    Efficient implementation of essentially non-oscillatory shock- capturing schemes,

    C.-W. Shu and S. Osher, “Efficient implementation of essentially non-oscillatory shock- capturing schemes,”Journal of Computational Physics, vol. 77, no. 2, pp. 439–471, 1988

    work page 1988

  27. [27]

    Towardstheultimateconservativedifferenceschemeiii.upstream-centered finite-difference schemes for ideal compressible flow,

    B.VanLeer, “Towardstheultimateconservativedifferenceschemeiii.upstream-centered finite-difference schemes for ideal compressible flow,”Journal of Computational Physics, vol. 23, no. 3, pp. 263–275, 1977

    work page 1977

  28. [28]

    Uniformly accurate schemes for hyperbolic systems withrelaxation,

    R. E. Caflisch, S. Jin, and G. Russo, “Uniformly accurate schemes for hyperbolic systems withrelaxation,”SIAM Journal on Numerical Analysis, vol.34, no.1, pp.246–281, 1997

    work page 1997

  29. [29]

    An upwind differencing scheme for the equations of ideal mag- netohydrodynamics,

    M. Brio and C. Wu, “An upwind differencing scheme for the equations of ideal mag- netohydrodynamics,”Journal of Computational Physics, vol. 75, no. 2, pp. 400–422, 1988. 41

    work page 1988

  30. [30]

    Athena: A New Code for Astrophysical MHD,

    J. M. Stone, T. A. Gardiner, P. Teuben, J. F. Hawley, and J. B. Simon, “Athena: A New Code for Astrophysical MHD,”The Astrophysical Journal Supplement Series, vol. 178, pp. 137–177, Sept. 2008

    work page 2008

  31. [31]

    Ideal and resistive edge stability calcu- lations with m3d-c1,

    N. M. Ferraro, S. C. Jardin, and P. B. Snyder, “Ideal and resistive edge stability calcu- lations with m3d-c1,”Physics of Plasmas, vol. 17, p. 102508, 10 2010

    work page 2010

  32. [32]

    An unsplit godunov method for ideal mhd via con- strained transport in three dimensions,

    T. A. Gardiner and J. M. Stone, “An unsplit godunov method for ideal mhd via con- strained transport in three dimensions,”Journal of Computational Physics, vol. 227, no. 8, pp. 4123–4141, 2008

    work page 2008

  33. [33]

    Comparative analysis of hall effect implementations in hall magnetohydrodynamics,

    K. Iwasaki and K. Tomida, “Comparative analysis of hall effect implementations in hall magnetohydrodynamics,”The Astrophysical Journal, vol. 984, p. 50, apr 2025

    work page 2025

  34. [34]

    P. D. Mullen,Magnetized models for the formation of the Moon. PhD thesis, University of Illinois at Urbana-Champaign, Urbana, IL, August 2021. Ph.D. Dissertation

    work page 2021

  35. [35]

    J. D. Huba,Hall Magnetohydrodynamics - A Tutorial, pp. 166–192. Berlin, Heidelberg: Springer Berlin Heidelberg, 2003

    work page 2003

  36. [36]

    Implicit highly-coupled single-ion hall-mhd formu- lation for hybrid particle-in-cell codes,

    C. Thoma, D. Welch, and D. Rose, “Implicit highly-coupled single-ion hall-mhd formu- lation for hybrid particle-in-cell codes,”Computer Physics Communications, vol. 261, p. 107823, 2021

    work page 2021

  37. [37]

    Divergence- and Curl-Preserving Prolongation and Restriction Formulas,

    G. Tóth and P. L. Roe, “Divergence- and Curl-Preserving Prolongation and Restriction Formulas,”Journal of Computational Physics, vol. 180, pp. 736–750, Aug. 2002

    work page 2002

  38. [38]

    University of Rochester, 2024

    FLASH Center for Computational Science, University of Rochester,FLASH User’s Guide (Development Version). University of Rochester, 2024. Development release, updated daily. 42

    work page 2024