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.

A Modular Grad-Div Stabilization for the 2D/3D Nonstationary Incompressible Magnetohydrodynamic Equations

Modular grad-div stabilization for the incompressible non-isothermal fluid flows

Analysis of the linearly extrapolated BDF2 fully discrete Modular Grad-div stabilization method for the micropolar Navier-Stokes equations

A staggered-grid multilevel incomplete LU for steady incompressible flows

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

Numerical Analysis and Computation of a Subgrid-Sparse-Grad-Div Stabilization Method for Incompressible Flow Problems

A Linear Stabilized Incompressible Magnetohydrodynamic Problem with Magnetic Pressure

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

A Parallel Finite Element Discretization Algorithm Based on Grad-Div Stabilization for the Navier–Stokes Equations

A second order partitioned method with grad-div stabilization for the non-stationary dual-porosity-Stokes model

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

Numerical analysis of a bdf2 modular grad-div Stabilization method for the Navier-Stokes equations

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

A stabilized Gauge-Uzawa discontinuous Galerkin method for the magneto-hydrodynamic equations

A Novel Stabilized Galerkin Meshless Method for Steady Incompressible Navier–Stokes Equations

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

Stability and optimal decay for the 3D magnetohydrodynamic equations with only horizontal dissipation

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

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

Gauge-Uzawa-based, highly efficient decoupled schemes for the diffuse interface model of two-phase magnetohydrodynamic

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