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.

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

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

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

Error estimates for the highly efficient and energy stable schemes for the 2D/3D two-phase MHD

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

Stability and error analysis of the SAV schemes for the inductionless MHD equations

The Euler Implicit/explicit Scheme for the 2D Time-Dependent Navier-Stokes Equations with Smooth or Non-Smooth Initial Data

Optimal error estimates of a decoupled finite element scheme for the unsteady inductionless MHD equations

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

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

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

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

Analysis of a time filtered finite element method for the unsteady inductionless MHD equations

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

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

Decoupled Schemes for Unsteady MHD Equations. I. Time Discretization

A Diffuse Interface Model and Semi-Implicit Energy Stable Finite Element Method for Two-Phase Magnetohydrodynamic Flows

An optimal error estimate of the BDF-Galerkin FEM for the incompressible MHD system

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

New error analysis of charge-conservative finite element methods for stationary inductionless MHD equations

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