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

Lechang Qin,Huiyuan Li,Zhiguo Yang
2023-12-19
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.
Numerical Analysis
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.