Abstract:In this paper, we present exact divergence-free spectral method for solving the incompressible and resistive magneto-hydrodynamic (MHD) equations in two and three dimensions, as well as the efficient solution algorithm and unconditionally energy-stable fully-discretized numerical schemes. We introduce new ideas of constructing two families of exact divergence-free vectorial spectral basis functions on domains diffeomorphic to squares or cubes. These bases are obtained with the help of orthogonality and derivative relation of generalised Jacobi polynomials, several de Rham complexes, as well as the property of contravariant Piola transformation. They are well-suited for discretizing the velocity and magnetic fields, respectively, thereby ensuring point-wise preservation of the incompressibility condition and the magnetic Gauss's law. With the aid of these bases, we propose a family of exact divergence-free implicit-explicit $k$-step backward differentiation formula (DF-BDF-$k$) fully-discretized schemes for the MHD system. These schemes naturally decouple the pressure field from the velocity field. Consequently, the stability of the space-time fully-discretized numerical schemes based on these bases are significantly enhanced. These schemes exhibit unconditional stability for $k=1,2$, and demonstrate exceptional stability and accuracy for $k=3,4$, verified with extensive numerical results for long time simulations using large time step sizes. Furthermore, we present efficient solution algorithms for these two decoupled equations for the velocity and magnetic fields, respectively, by exploiting the sparsity and structure of the resultant linear algebraic systems. Ample numerical examples in two and three dimensions are provided to demonstrate the distinctive accuracy, efficiency and stability of our proposed method.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop an accurate divergence - free spectral method for solving the incompressible and resistive magnetohydrodynamic (MHD) equations in two - and three - dimensional cases. Specifically, the researchers proposed new construction methods to create two types of accurate divergence - free vector spectral basis functions, which are suitable for discretizing the velocity field and magnetic field on regions that are topologically homeomorphic to a square or a cube. By introducing the orthogonality and derivative relations of generalized Jacobi polynomials, several de Rham complexes, and the properties of the contravariant Piola transformation, these basis functions can ensure the point - by - point preservation of the incompressible condition and the magnetic Gauss law. ### Main Problems and Challenges 1. **Divergence - free Condition**: That is, $\nabla\cdot\mathbf{u} = 0$ and $\nabla\cdot\mathbf{B} = 0$. 2. **Properties of Saddle - point Systems**: The constraint $\nabla\cdot\mathbf{u} = 0$ introduced by the Lagrange multiplier $p$ causes the system to have saddle - point properties. 3. **Nonlinear Coupling**: The nonlinear coupling in the MHD system increases the complexity of numerical simulations. ### Solutions To address the above challenges, the author proposed a systematic method to construct accurate divergence - free H1 - and H(div) - conforming spectral basis functions using the de Rham complex. These basis functions not only ensure the point - by - point preservation of the divergence - free condition but also significantly enhance the stability of the fully space - time discrete numerical scheme. In addition, the author proposed a class of accurate divergence - free implicit - explicit k - step backward difference formula (DF - BDF - k) fully discrete schemes, which have unconditional stability for $k = 1,2$ and exhibit excellent stability and accuracy in long - time simulations for $k = 3,4$. ### Numerical Results A large number of numerical experiments show that the proposed algorithms exhibit excellent precision, efficiency, and stability in both two - and three - dimensional cases. In particular, these algorithms can perform long - time simulations with large time steps and only need to solve two symmetric positive - definite equations instead of the coupled saddle - point system, thereby improving computational efficiency. ### Summary The main contribution of this paper lies in providing an accurate divergence - free spectral method, which solves key problems in the numerical simulation of MHD equations, including the preservation of the divergence - free condition, the nonlinear coupling of the system, and the stability of long - time simulations.

An exact divergence-free spectral method for incompressible and resistive magneto-hydrodynamic equations in two and three dimensions

A well-balanced and exactly divergence-free staggered semi-implicit hybrid finite volume/finite element scheme for the incompressible MHD equations

A divergence-free hybrid finite volume / finite element scheme for the incompressible MHD equations based on compatible finite element spaces with a posteriori limiting

A Divergence-Free and $H(div)$-Conforming Embedded-Hybridized DG Method for the Incompressible Resistive MHD equations

A new discretely divergence-free positivity-preserving high-order finite volume method for ideal MHD equations

A high-fidelity finite volume scheme for ideal magnetohydrodynamics equations using boundary variation diminishing algorithm

A Locally Divergence-Free Oscillation-Eliminating Discontinuous Galerkin Method for Ideal Magnetohydrodynamic Equations

An Unstaggered Constrained Transport Method for the 3D Ideal Magnetohydrodynamic Equations

Stream Function-Velocity-magnetic Induction Compact Difference Method for the 2D Steady Incompressible Full Magnetohydrodynamic Equations

On Convergence of Numerical Solutions for the Compressible MHD System with Exactly Divergence-Free Magnetic Field

A novel divergence-free Finite Element Method for the MHD Kinematics equations using Vector-potential

An Adaptive Moving Mesh Method for Two-Dimensional Ideal Magnetohydrodynamics

Fully decoupled, linearized and stabilized finite volume method for the time-dependent incompressible MHD equations

A structure-preserving semi-implicit IMEX finite volume scheme for ideal magnetohydrodynamics at all Mach and Alfvén numbers

Two novel fully decoupled schemes for the two-phase MHD system with exactly divergence-free magnetic field

An arbitrary high-order Spectral Difference method for the induction equation

High-order Accurate Entropy Stable Adaptive Moving Mesh Finite Difference Schemes for Special Relativistic (Magneto)hydrodynamics.

Analysis of a semi-implicit and structure-preserving finite element method for the incompressible MHD equations with magnetic-current formulation

Second Order Unconditionally Convergent Fully Discrete Scheme for Incompressible Vector Potential MHD System

Unconditionally energy stable numerical schemes for the three-dimensional magneto-micropolar equations

Decoupled structure-preserving discretization of incompressible MHD equations with general boundary conditions