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.

Unconditional Convergence of the Euler Semi-Implicit Scheme for the Three-Dimensional Incompressible Mhd Equations

Stability and Error Analysis for the First-Order Euler Implicit/Explicit Scheme for the 3D MHD Equations

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

The almost unconditional convergence of the Euler implicit/explicit scheme for the three dimensional nonstationary Navier-Stokes equations

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

Decoupled Schemes for Unsteady MHD Equations. I. Time Discretization

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

An Efficient Hermite-Galerkin Spectral Scheme for Three-Dimensional Incompressible Hall-magnetohydrodynamic System on Infinite Domain

A Hybrid Scheme for Three-Dimensional Incompressible Two-Phase Flows

Convergence analysis of three finite element iterative methods for the 2D/3D stationary incompressible magnetohydrodynamics

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

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

Stability and Convergence of Semi-Implicit Time-Stepping Algorithm for Stationary Incompressible Magnetohydrodynamics

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

Decoupled schemes for unsteady MHD equations II: Finite element spatial discretization and numerical implementation.

Uniform Stability and Convergence with Respect to (ν , Μ , S, 1-Σ ) of the Three Iterative Finite Element Solutions for the 3D Steady MHD Equations

Convergence analysis of a fully discrete finite element method for thermally coupled incompressible MHD problems with temperature-dependent coefficients

Streamline Diffusion Finite Element Method for Stationary Incompressible Magnetohydrodynamics

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