A Gradient Recovery Method for Electron Magnetohydrodynamics with Fractional Dissipation

Abstract: We propose and analyze a structure-preserving numerical method for the $2\tfrac{1}{2}$-dimensional (2.5D) electron magnetohydrodynamics system with fractional dissipation on the periodic torus. The method works directly with the magnetic field components and combines this component formulation with the gradient recovery operator of [T. Chu, H. Guo, and Z. Zhang, SIAM J. Numer. Anal., 63 (2025), pp. 23--53]. We establish discrete energy stability for a semi-implicit structure-preserving formulation and use an explicit-Hall integrating-factor implementation for efficient computation on periodic grids. The fractional dissipation is treated exactly in Fourier space, and the in-plane divergence constraint is enforced by a spectral Hodge projection. Numerical experiments demonstrate second-order spatial convergence and stable Hall-driven dynamics across several benchmark tests.

• The paper introduces a structure-preserving method for 2.5D EMHD that bypasses high-order elements through direct discretization and gradient recovery.
• The paper employs a semi-implicit integrating-factor scheme to linearize Hall terms and ensure unconditional energy stability with second-order spatial convergence.
• The paper validates its approach with benchmarks such as Kelvin–Helmholtz shear layers and magnetic island mergers, accurately capturing Hall effects and fractional dissipation.

Gradient Recovery for Electron Magnetohydrodynamics with Fractional Dissipation

EMHD Model and Physical Context

The paper introduces a structure-preserving numerical method for the analysis and simulation of the $2.5$D electron magnetohydrodynamics (EMHD) system with fractional magnetic dissipation. The EMHD model governs the magnetic field evolution in plasmas where electron dynamics dominate and is critically relevant for understanding electron-scale magnetic reconnection, whistler wave propagation, and current-sheet dynamics. The Hall term, $((\operatorname{curl} B)\times B)$, introduces a third-order spatial nonlinearity, substantially increasing both modeling and computational complexity compared to standard MHD.

The fractional Laplacian $(-\Delta)^{\alpha/2}$ with $1 < \alpha \leq 2$ is incorporated to model anomalous, scale-dependent dissipation, thereby capturing regimes where standard resistive models are not physically or mathematically adequate. This model aligns with local well-posedness results under fractional dissipation established by Chae--Wan--Wu [ChaeWanWu15].

Component Formulation and Gradient Recovery

A salient methodological choice is the direct discretization of the magnetic field components $(B_x, B_y, B_z)$, rather than relying on scalar flux-function representations. This enables direct enforcement of divergence-free conditions for the in-plane components and yields a formulation amenable to low-order finite element discretization, avoiding third-order Hall terms present in scalar forms.

The gradient recovery operator $G_h$, as developed in [CGZ25], is pivotal. It reconstructs superconvergent and smoother gradients from $C^0$ finite element solutions, thus enabling stable and accurate computation of second-order discrete Hall terms on linear elements. The patch-based recovery strategy achieves second-order spatial convergence and exhibits stability equivalent to discrete gradient norms, which is critical for handling the derivative-intensive Hall coupling.

Discrete Structure and Algorithmic Implementation

The spatial discretization employs $C^0$ linear finite element spaces on shape-regular triangulations complemented with gradient recovery. The approach reduces the derivative count in Hall terms, thus making the discretization feasible with lower-order elements. The discrete operators are separately constructed for in-plane curl (current), Hall nonlinearity, and the fractional Laplacian, the latter of which is treated exactly in Fourier space.

The time-stepping scheme is semi-implicit: one factor in each Hall term is treated explicitly to linearize the update, while fractional diffusion is handled unconditionally stably via an integrating-factor method. Spectral Hodge projection enforces the in-plane divergence constraint, maintaining physical consistency to round-off accuracy. The resultant computational loop leverages FFTs for all differential operators, reducing complexity and enhancing scalability.

Notably, the paper establishes rigorous discrete energy stability for the semi-implicit, structure-preserving formulation, where Hall cancellation is preserved at the discrete level. The $L^2$ stability theorem guarantees unconditional nonlinear stability with respect to magnetic energy, while a high-regularity $H_h^s$ estimate follows via discrete Littlewood–Paley theory, provided $((\operatorname{curl} B)\times B)$0.

Numerical Results and Accuracy

Numerical experiments validate the spatial convergence, energy stability, and physical fidelity of the scheme. Second-order convergence is observed in both $((\operatorname{curl} B)\times B)$1 and $((\operatorname{curl} B)\times B)$2 norms, irrespective of fractional exponent $((\operatorname{curl} B)\times B)$3. Benchmark tests include:

• Kelvin–Helmholtz Shear Layer: Hall-driven perturbations induce undulations in sharp-gradient current sheets, with persistent layer coherence and energy decay consistent with physical expectations.
• Magnetic Island Merging: Hall and fractional dissipation mechanisms drive reconnection and coalescence of flux islands, with fine-scale current sheet formation reflecting real-world plasma behavior.
• Helical Initial Data: The integrating-factor scheme produces energy decay in alignment with analytical predictions; Hall-induced nonlinear corrections were quantified at sub-3% levels for moderate amplitude initial data.

The divergence-free constraint is maintained to round-off accuracy throughout, confirming the effectiveness of the spectral Hodge projection.

Theoretical and Practical Implications

The method demonstrates that structure-preserving discretization—retaining discrete Hall cancellation and exact treatment of fractional diffusion—is essential for robust EMHD simulation, especially in regimes where higher-order nonlinear terms intensify numerical instability and physical ill-posedness. The component-based formulation and recovery operator together allow accurate and stable simulation on low-order finite element spaces, circumventing the need for complex high-order elements.

The theoretical framework, especially the discrete $((\operatorname{curl} B)\times B)$4 stability and Hall cancellation mechanism, provides a blueprint for future EMHD algorithms. The approach can be extended to more general boundary conditions, variable coefficients, and stochastic EMHD models [HuPengYang26]. The unconditional stability under fractional dissipation suggests utility in practical plasma simulation scenarios with anomalous transport.

Future Directions

Potential avenues include the adaptation to fully three-dimensional EMHD systems, integration with adaptive mesh refinement targeting current sheets, and extension to stochastic or turbulent plasma regimes where fractional dissipation is physically motivated. Further mathematical analysis could address rigorous discrete paraproduct estimates and well-posedness on non-uniform meshes, as well as coupling with ion dynamics.

Conclusion

This work rigorously develops and analyzes a gradient recovery-based, structure-preserving method for $((\operatorname{curl} B)\times B)$5D EMHD with fractional dissipation. By combining component-form discretization, spectral gradient recovery, integrating-factor time stepping, and spectral divergence projection, the scheme attains second-order spatial accuracy, unconditional energy stability, and faithful reproduction of Hall-driven and dissipative plasma dynamics. The theoretical guarantees and practical results position this approach as a robust foundation for advanced numerical EMHD studies (2606.19716).