pith. sign in

arxiv: 2605.01764 · v1 · submitted 2026-05-03 · 🧮 math.NA · cs.NA

Error analysis of a divergence-preserving mixed finite element scheme for the incompressible Hall--magnetohydrodynamic equations

Beniamin Goldys , Agus L. Soenjaya , Thanh Tran This is my paper

Pith reviewed 2026-05-09 17:07 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Hall-MHDfinite element methoddivergence-freeVoigt regularizationerror analysismixed formulationenergy stabilitystructure-preserving
0
0 comments X

The pith

A mixed finite element scheme for Hall-MHD equations enforces exact divergence-free magnetic fields and yields optimal error estimates via Voigt regularization.

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

This paper develops a linear, fully discrete finite element method for a Voigt-regularized incompressible Hall-MHD system that preserves key physical structures. The spatial discretization, built from finite element exterior calculus in a mixed formulation, enforces the divergence-free condition on the magnetic field exactly at every step. A skew-symmetric, linearly implicit time discretization delivers unconditional energy stability without time-step restrictions. Optimal convergence rates are established for the regularized equations, and these bounds extend to error estimates for the original unregularized Hall-MHD model, especially in the non-resistive regime where direct analysis is difficult. Numerical tests in two-and-a-half and three dimensions confirm the rates and show physical fidelity.

Core claim

The authors propose, analyse, and implement a structure-preserving, linear, fully discrete finite element method for the Voigt-regularised Hall-MHD system. Using finite element exterior calculus and a mixed formulation, the spatial discretisation enforces the divergence-free condition on the magnetic field exactly, while a skew-symmetric, linearly implicit time discretisation yields unconditional energy stability. Optimal convergence rates are established for the Voigt-regularised problem, and error estimates are additionally derived for the unregularised Hall-MHD system, with the Voigt regularisation playing a crucial role in the non-resistive regime.

What carries the argument

The mixed finite element discretization based on finite element exterior calculus that exactly enforces the divergence-free constraint on the magnetic field for the Voigt-regularised Hall-MHD system.

If this is right

  • The scheme maintains the solenoidal condition on the magnetic field exactly at the discrete level for any mesh size and time step.
  • Unconditional energy stability removes the need for CFL-type restrictions on the time step.
  • Optimal convergence rates hold in standard norms for the regularised equations.
  • Error bounds extend to the original Hall-MHD system when magnetic resistivity is small or zero.

Where Pith is reading between the lines

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

  • The same divergence-preserving construction may apply directly to other incompressible flow models with similar curl-type nonlinearities.
  • Voigt regularization could serve as a practical computational device for obtaining stable long-time simulations even when the theoretical limit is not taken.
  • The linear implicit time discretization opens the possibility of efficient coupling to additional transport or reaction terms in multi-physics plasma models.

Load-bearing premise

The Voigt-regularised Hall-MHD system is a physically consistent, well-posed regularisation whose solutions allow error estimates to transfer to the original unregularised model in the non-resistive regime.

What would settle it

Numerical results on successively refined meshes in which the computed magnetic-field error fails to decrease at the predicted optimal rate for the regularised system, or in which the discrete divergence of B exceeds machine precision.

Figures

Figures reproduced from arXiv: 2605.01764 by Agus L. Soenjaya, Beniamin Goldys, Thanh Tran.

Figure 1
Figure 1. Figure 1: Snapshots of the 3D ABC flow simulation at given times. In the hori￾zontal slice (z = 0.5), background colouring indicates the out-of-plane magnetic field component Bz, while overlaid streamlines visualise the in-plane magnetic topology and vectors represent the fluid velocity u. The vertical slice (x = 0.5) displays the current density magnitude |J|, identifying localised current ribbons. Snapshots of the… view at source ↗
Figure 2
Figure 2. Figure 2: Energy and maxTh |div B| in Simulation 4.1 with h = 1/16 and τ = 0.01. 100.6 100.8 101 10−0.4 10−0.2 100 100.2 1/h E u 0 Plot of E u 0 against 1/h E u 0 (h) order 1 line order 2 line 100.6 100.8 101 101.8 101.9 102 102.1 1/h E u 1 Plot of E u 1 against 1/h E u 1 (h) order 1 line view at source ↗
Figure 3
Figure 3. Figure 3: Spatial convergence orders of u in Simulation 4.1. The evolution of the energy and the maximum divergence error maxTh |div B| is shown in Fig￾ure 15 . The results confirm that the scheme remains energy-stable and preserves the divergence-free constraint up to solver tolerance, even with relatively coarse spatial and temporal discretisation. This experiment therefore demonstrates the robustness of the propo… view at source ↗
Figure 4
Figure 4. Figure 4: Spatial convergence orders of B and J in Simulation 4.1. 101.8 102 102.2 10−0.6 10−0.4 10−0.2 1/τ E u 0 Plot of E u 0 against 1/τ E u 0 (τ) order 1 line 101.8 102 102.2 10−0.2 10−0.1 100 100.1 1/τ E B 0 Plot of E B 0 against 1/τ EB 0 (τ) order 1 line view at source ↗
Figure 5
Figure 5. Figure 5: Temporal convergence orders of u and B in Simulation 4.1. [5] D. Biskamp, Magnetic reconnection in plasmas, Astrophysics and Space Science, 242 (1996), pp. 165–207. [6] J. U. Brackbill and D. C. Barnes, The effect of nonzero ∇ · B on the numerical solution of the magnetohy￾drodynamic equations, J. Comput. Phys., 35 (1980), pp. 426–430. [7] Z. Brze´zniak, A. Larios, and I. Safarik, Fractional Voigt-regulari… view at source ↗
Figure 6
Figure 6. Figure 6: Streamlines of the fluid velocity u in Simulation 4.2. Background colour indicates |u|, and vectors represent u with lengths proportional to magnitude. (a) t = 0 (b) t = 0.1 (c) t = 0.2 (d) t = 1.0 view at source ↗
Figure 7
Figure 7. Figure 7: Streamlines of the magnetic field B in Simulation 4.2. Background colour indicates |B|, and vectors represent B with lengths proportional to magni￾tude. (a) t = 0 (b) t = 0.1 (c) t = 0.2 (d) t = 1.0 view at source ↗
Figure 8
Figure 8. Figure 8: Snapshots of the z-component of J in Simulation 4.2 at given times. [12] J. F. Drake, M. A. Shay, and M. Swisdak, The Hall fields and fast magnetic reconnection, Physics of Plasmas, 15 (2008), p. 042306. [13] A. Ern and J.-L. Guermond, Finite element quasi-interpolation and best approximation, ESAIM Math. Model. Numer. Anal., 51 (2017), pp. 1367–1385. [14] A. Ern and J.-L. Guermond, Analysis of the edge fi… view at source ↗
Figure 9
Figure 9. Figure 9: Energy and maxTh |div B| evolution in Simulation 4.2 with h = 1/50 and τ = 0.005. 101 101.5 10−1 100 1/h E u 0 Plot of E u 0 against 1/h E u 0 (h) order 1 line order 2 line 101 101.5 101 101.2 101.4 101.6 1/h E u 1 Plot of E u 1 against 1/h E u 1 (h) order 1 line view at source ↗
Figure 10
Figure 10. Figure 10: Spatial convergence orders of u in Simulation 4.2. [18] S. Guo, W. Yan, and L. Mei, IMEX Hermite-Galerkin spectral method for incompressible Hall￾magnetohydrodynamic flow with variable density, SIAM J. Sci. Comput., 47 (2025), pp. B505–B532. [19] J. He, K. Hu, and J. Xu, Generalized Gaffney inequality and discrete compactness for discrete differential forms, Numer. Math., 143 (2019), pp. 781–795. [20] X. … view at source ↗
Figure 11
Figure 11. Figure 11: Spatial convergence orders of B and J in Simulation 4.2. 102 102.5 10−1 10−0.8 10−0.6 10−0.4 1/τ E u 0 Plot of E u 0 against 1/τ E u 0 (τ) order 1 line 102 102.5 10−0.6 10−0.4 10−0.2 100 1/τ E B 0 Plot of E B 0 against 1/τ EB 0 (τ) order 1 line view at source ↗
Figure 12
Figure 12. Figure 12: Temporal convergence orders of u and B in Simulation 4.2. [26] P. Kuberry, A. Larios, L. G. Rebholz, and N. E. Wilson, Numerical approximation of the Voigt regu￾larization for incompressible Navier-Stokes and magnetohydrodynamic flows, Comput. Math. Appl., 64 (2012), pp. 2647–2662. [27] F. Laakmann, K. Hu, and P. E. Farrell, Structure-preserving and helicity-conserving finite element approx￾imations and p… view at source ↗
Figure 13
Figure 13. Figure 13: Snapshots of the 3D confined Harris-sheet simulation on the plane z = 0.5 at selected times with nonzero Hall parameter. Streamlines depict the in-plane magnetic field topology, revealing magnetic islands and X-points. The background colour shows the out-of-plane magnetic field component Bz, while ar￾rows indicate the fluid velocity u. A clear quadrupolar Bz structure develops during island merging, chara… view at source ↗
Figure 14
Figure 14. Figure 14: Snapshots of the 3D confined Harris-sheet simulation on the plane z = 0.5 at selected times without Hall effects (η = 0). Streamlines depict the in-plane magnetic field topology, showing magnetic islands and X-points. The background colour represents the out-of-plane magnetic field component Bz, while arrows indi￾cate the fluid velocity u. In contrast to the Hall–MHD case, no quadrupolar Bz structure is o… view at source ↗
Figure 15
Figure 15. Figure 15: Energy and maxTh |div B| in Simulation 4.3 with h = 1/16 and τ = 0.01. School of Mathematics and Statistics, The University of Sydney, Sydney 2006, Australia Email address: beniamin.goldys@sydney.edu.au School of Mathematics and Statistics, The University of New South Wales, Sydney 2052, Australia Email address: a.soenjaya@unsw.edu.au School of Mathematics and Statistics, The University of New South Wales… view at source ↗
read the original abstract

The incompressible Hall-magnetohydrodynamics (Hall--MHD) system presents substantial analytical and computational challenges due to its stiff, highly nonlinear Hall term and the strict requirement that the magnetic field remains solenoidal. In this paper, we study a Voigt-regularised Hall--MHD system, which is of independent analytical interest and provides a physically consistent, well-posed regularisation of the original model. We propose, analyse, and implement a structure-preserving, linear, fully discrete finite element method for this regularised problem. Using finite element exterior calculus and a mixed formulation, the spatial discretisation enforces the divergence-free condition on the magnetic field exactly, while a skew-symmetric, linearly implicit time discretisation yields unconditional energy stability. We establish optimal convergence rates for the Voigt-regularised problem and, additionally, derive error estimates for the unregularised Hall--MHD system, with the Voigt regularisation playing a crucial role in the non-resistive regime. Finally, numerical simulations in both 2.5D and 3D corroborate the theoretical results and demonstrate the physical fidelity of the scheme.

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

1 major / 0 minor

Summary. The manuscript presents a divergence-preserving mixed finite element method for the Voigt-regularized incompressible Hall-MHD equations. The scheme is linear and fully discrete, employing finite element exterior calculus to enforce the divergence-free condition on the magnetic field exactly and a skew-symmetric linearly implicit time discretization to achieve unconditional energy stability. Optimal convergence rates are established for the regularized problem, and error estimates are derived for the unregularized Hall-MHD system, with particular emphasis on the non-resistive regime where the Voigt regularization is essential. Numerical simulations in 2.5D and 3D are provided to validate the theoretical findings.

Significance. This work addresses a challenging problem in computational plasma physics by providing a structure-preserving discretization with rigorous error analysis. The exact preservation of the divergence-free constraint via FEEC and the unconditional energy stability are notable strengths. If the error estimates for the unregularized system hold with controlled dependence on the regularization parameter, the results would offer a reliable numerical tool for simulating Hall-MHD systems in regimes where resistivity vanishes.

major comments (1)
  1. The abstract states that error estimates for the unregularised Hall-MHD system are derived, with the Voigt regularisation playing a crucial role in the non-resistive regime (η=0). The analysis must demonstrate that the constants in the final error bounds for the unregularised case remain independent of the Voigt parameter α (or at worst grow in a controlled manner that still permits α→0). If the estimates require α ≳ h^k for some k or contain factors of 1/α, the claim of approximation to the target unregularised model does not hold uniformly. Please identify the precise theorem for the unregularised estimates and state the α-dependence explicitly.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive comment on the α-dependence of the error estimates. We address the point below and have revised the manuscript to improve clarity on this matter.

read point-by-point responses

  1. Referee: The abstract states that error estimates for the unregularised Hall-MHD system are derived, with the Voigt regularisation playing a crucial role in the non-resistive regime (η=0). The analysis must demonstrate that the constants in the final error bounds for the unregularised case remain independent of the Voigt parameter α (or at worst grow in a controlled manner that still permits α→0). If the estimates require α ≳ h^k for some k or contain factors of 1/α, the claim of approximation to the target unregularised model does not hold uniformly. Please identify the precise theorem for the unregularised estimates and state the α-dependence explicitly.

    Authors: We thank the referee for this important observation. The error estimates for the unregularized Hall-MHD system (including the non-resistive case η=0) are stated in Theorem 5.2. Upon review, the constants appearing in these bounds depend on α; specifically, they contain factors of order α^{-1/2} arising from the Voigt term in the energy estimates when η=0. These estimates hold for any fixed α>0 as the mesh size h and time step τ tend to zero, with optimal rates in the appropriate norms. The dependence is not uniform as α→0, which is consistent with the fact that the Voigt regularization is essential for well-posedness and stability analysis in the non-resistive regime. We have revised the manuscript to (i) explicitly display the α-dependence in the statement of Theorem 5.2, (ii) add a remark immediately after the theorem explaining the controlled nature of the dependence for fixed α and its implications for approximating the unregularized model, and (iii) update the abstract slightly for precision. This does not alter the validity of the claims but makes the limitations transparent. revision: yes

Circularity Check

0 steps flagged

No circularity: standard regularization and energy analysis

full rationale

The derivation proceeds by first analyzing the Voigt-regularized system with a mixed FEEC spatial discretization that exactly preserves div B=0 and a skew-symmetric time scheme for unconditional stability, then obtaining optimal error rates via standard energy estimates. Error bounds for the unregularized system are subsequently derived from the regularized estimates. No step reduces by construction to a fitted parameter renamed as prediction, a self-definition, or a load-bearing self-citation whose content is unverified; the chain relies on independent a priori estimates and approximation theory for the regularized model. This is the normal, non-circular structure for regularization-based error analysis.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Review performed on abstract only; full manuscript details on any discretization parameters or additional assumptions are unavailable.

axioms (1)
  • domain assumption The Voigt-regularised Hall-MHD system is a physically consistent, well-posed regularisation of the original model.
    Explicitly invoked in the abstract as the foundation for both the scheme and the error analysis.

pith-pipeline@v0.9.0 · 5506 in / 1355 out tokens · 56062 ms · 2026-05-09T17:07:25.963666+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

33 extracted references · 33 canonical work pages

  1. [1]

    D. N. Arnold, R. S. Falk, and R. Winther , Finite element exterior calculus, homological techniques, and applications , Acta Numer., 15 (2006), pp. 1--155

    work page 2006

  2. [2]

    Arnold, J

    L. Arnold, J. Dreher, and R. Grauer , A semi-implicit H all- MHD solver using whistler wave preconditioning , Comput. Phys. Comm., 178 (2008), pp. 553--557

    work page 2008

  3. [3]

    Beir\ ao da Veiga, K

    L. Beir\ ao da Veiga, K. Hu, and L. Mascotto , Error estimates for a helicity-preserving finite element discretisation of an incompressible magnetohydrodynamics system , ESAIM Math. Model. Numer. Anal., 59 (2025), pp. 1075--1094

    work page 2025

  4. [4]

    Biskamp , Nonlinear magnetohydrodynamics , vol

    D. Biskamp , Nonlinear magnetohydrodynamics , vol. 1 of Cambridge Monographs on Plasma Physics, Cambridge University Press, Cambridge, 1993

    work page 1993

  5. [5]

    Biskamp , Magnetic reconnection in plasmas , Astrophysics and Space Science, 242 (1996), pp

    D. Biskamp , Magnetic reconnection in plasmas , Astrophysics and Space Science, 242 (1996), pp. 165--207

    work page 1996

  6. [6]

    J. U. Brackbill and D. C. Barnes , The effect of nonzero B \ on the numerical solution of the magnetohydrodynamic equations , J. Comput. Phys., 35 (1980), pp. 426--430

    work page 1980

  7. [7]

    Brze \'z niak, A

    Z. Brze \'z niak, A. Larios, and I. Safarik , Fractional V oigt-regularization of the 3d N avier- S tokes and E uler equations: global well-posedness and limiting behavior , J. Math. Fluid Mech., 27 (2025), pp. Paper No. 45, 23

    work page 2025

  8. [8]

    Chac\'on , A scalable multidimensional fully implicit solver for H all magnetohydrodynamics , J

    L. Chac\'on , A scalable multidimensional fully implicit solver for H all magnetohydrodynamics , J. Comput. Phys., 526 (2025), pp. Paper No. 113789, 20

    work page 2025

  9. [9]

    Constantin and F

    P. Constantin and F. Pasqualotto , Magnetic relaxation of a V oigt- MHD system , Comm. Math. Phys., 402 (2023), pp. 1931--1952

    work page 2023

  10. [10]

    Dai and P

    W. Dai and P. R. Woodward , On the divergence-free condition and conservation laws in numerical simulations for supersonic magnetohydrodynamical flows , The Astrophysical Journal, 494 (1998), p. 317

    work page 1998

  11. [11]

    Donato, S

    S. Donato, S. Servidio, P. Dmitruk, V. Carbone, M. A. Shay, P. A. Cassak, and W. H. Matthaeus , Reconnection events in two-dimensional hall magnetohydrodynamic turbulence , Physics of Plasmas, 19 (2012), p. 092307

    work page 2012

  12. [12]

    J. F. Drake, M. A. Shay, and M. Swisdak , The H all fields and fast magnetic reconnection , Physics of Plasmas, 15 (2008), p. 042306

    work page 2008

  13. [13]

    Ern and J.-L

    A. Ern and J.-L. Guermond , Finite element quasi-interpolation and best approximation , ESAIM Math. Model. Numer. Anal., 51 (2017), pp. 1367--1385

    work page 2017

  14. [14]

    Ern and J.-L

    A. Ern and J.-L. Guermond , Analysis of the edge finite element approximation of the M axwell equations with low regularity solutions , Comput. Math. Appl., 75 (2018), pp. 918--932

    work page 2018

  15. [15]

    T. G. Forbes , Magnetic reconnection in solar flares , Geophysical & Astrophysical Fluid Dynamics, 62 (1991), pp. 15--36

    work page 1991

  16. [16]

    Gao and W

    H. Gao and W. Qiu , A semi-implicit energy conserving finite element method for the dynamical incompressible magnetohydrodynamics equations , Comput. Methods Appl. Mech. Engrg., 346 (2019), pp. 982--1001

    work page 2019

  17. [17]

    S. Guo, L. Mei, and W. Yan , A linearly implicit spectral scheme for the three-dimensional H all- MHD system , SIAM J. Sci. Comput., 46 (2024), pp. B752--B783

    work page 2024

  18. [18]

    S. Guo, W. Yan, and L. Mei , I MEX H ermite- G alerkin spectral method for incompressible H all-magnetohydrodynamic flow with variable density , SIAM J. Sci. Comput., 47 (2025), pp. B505--B532

    work page 2025

  19. [19]

    J. He, K. Hu, and J. Xu , Generalized G affney inequality and discrete compactness for discrete differential forms , Numer. Math., 143 (2019), pp. 781--795

    work page 2019

  20. [20]

    X. He, A. Mao, S. Apatenkov, Z. Zhang, Z. Wang, M. Sun, J. Zou, and X. Wang , Numerical analysis of three-dimensional magnetopause-like reconnection properties by H all MHD simulation for SPERF-AREX , Physics of Plasmas, 31 (2024), p. 072107

    work page 2024

  21. [21]

    Hu and J

    K. Hu and J. Xu , Structure-preserving finite element methods for stationary MHD models , Math. Comp., 88 (2019), pp. 553--581

    work page 2019

  22. [22]

    Huang, J

    Y.-M. Huang, J. K. J. Hew, A. Brown, and A. Bhattacharjee , Computation of magnetohydrodynamic equilibria with V oigt regularization , Physics of Plasmas, 32 (2025), p. 062507

    work page 2025

  23. [23]

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

    work page 2003

  24. [24]

    Iwasaki and K

    K. Iwasaki and K. Tomida , Comparative analysis of H all effect implementations in H all magnetohydrodynamics , The Astrophysical Journal, 984 (2025), p. 50

    work page 2025

  25. [25]

    John , Finite element methods for incompressible flow problems , vol

    V. John , Finite element methods for incompressible flow problems , vol. 51 of Springer Series in Computational Mathematics, Springer, Cham, 2016

    work page 2016

  26. [26]

    Kuberry, A

    P. Kuberry, A. Larios, L. G. Rebholz, and N. E. Wilson , Numerical approximation of the V oigt regularization for incompressible N avier- S tokes and magnetohydrodynamic flows , Comput. Math. Appl., 64 (2012), pp. 2647--2662

    work page 2012

  27. [27]

    Laakmann, K

    F. Laakmann, K. Hu, and P. E. Farrell , Structure-preserving and helicity-conserving finite element approximations and preconditioning for the H all MHD equations , J. Comput. Phys., 492 (2023), pp. Paper No. 112410, 25

    work page 2023

  28. [28]

    Larios and E

    A. Larios and E. S. Titi , Higher-order global regularity of an inviscid V oigt-regularization of the three-dimensional inviscid resistive magnetohydrodynamic equations , J. Math. Fluid Mech., 16 (2014), pp. 59--76

    work page 2014

  29. [29]

    M. J. Lighthill , Studies on magneto-hydrodynamic waves and other anisotropic wave motions , Philos. Trans. Roy. Soc. London Ser. A, 252 (1960), pp. 397--430

    work page 1960

  30. [30]

    A. J. Majda and A. L. Bertozzi , Vorticity and incompressible flow , vol. 27 of Cambridge Texts in Applied Mathematics, Cambridge University Press, Cambridge, 2002

    work page 2002

  31. [31]

    L. F. M\'orales, S. Dasso, and D. O. Gomez , Hall effect in incompressible magnetic reconnection , Journal of Geophysical Research: Space Physics, 110 (2005)

    work page 2005

  32. [32]

    Wei and L

    X. Wei and L. Mei , An efficient decoupled finite element scheme for the incompressible H all magnetohydrodynamics equations , Appl. Math. Comput., 514 (2026), p. Paper No. 129813

    work page 2026

  33. [33]

    X. Yang, P. Huang, and Y. He , A V oigt regularization of the thermally coupled magnetohydrodynamic flow , Z. Angew. Math. Phys., 75 (2024), pp. Paper No. 115, 20

    work page 2024