Abstract:We present a divergence-free and $H(div)$-conforming hybridized discontinuous Galerkin (HDG) method and a computationally efficient variant called embedded-HDG (E-HDG) for solving stationary incompressible viso-resistive magnetohydrodynamic (MHD) equations. The proposed E-HDG approach uses continuous facet unknowns for the vector-valued solutions (velocity and magnetic fields) while it uses discontinuous facet unknowns for the scalar variable (pressure and magnetic pressure). This choice of function spaces makes E-HDG computationally far more advantageous, due to the much smaller number of degrees of freedom, compared to the HDG counterpart. The benefit is even more significant for three-dimensional/high-order/fine mesh scenarios. On simplicial meshes, the proposed methods with a specific choice of approximation spaces are well-posed for linear(ized) MHD equations. For nonlinear MHD problems, we present a simple approach exploiting the proposed linear discretizations by using a Picard iteration. The beauty of this approach is that the divergence-free and $H(div)$-conforming properties of the velocity and magnetic fields are automatically carried over for nonlinear MHD equations. We study the accuracy and convergence of our E-HDG method for both linear and nonlinear MHD cases through various numerical experiments, including two- and three-dimensional problems with smooth and singular solutions. The numerical examples show that the proposed methods are pressure robust, and the divergence of the resulting velocity and magnetic fields is machine zero for both smooth and singular problems.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop a new numerical method for solving the incompressible viscous resistive magnetohydrodynamic (MHD) equations. Specifically, the author proposes an Embedded - Hybridizable Discontinuous Galerkin (E - HDG) method that is divergence - free and conforms to the H(div) space. ### Main problems 1. **Multi - physics coupling**: Magnetohydrodynamics involves the coupling of fluid mechanics and electromagnetism, with different time and space scales. 2. **Non - linearity**: The MHD equations are non - linear, which increases the difficulty of solving them. 3. **Incompressibility**: It is necessary to precisely satisfy the mass conservation condition (i.e., the divergence of the velocity field is zero), which is crucial for pressure robustness. 4. **Solenoidal constraint of the magnetic field**: The divergence of the magnetic field must be zero to ensure the correct magnetic field topology and avoid introducing non - physical forces. 5. **Double saddle - point structure**: The discretization of the velocity - pressure system results in a large condition number and is difficult to solve. ### Solutions The author proposes two methods: - **Hybridizable Discontinuous Galerkin (HDG) method**: Reduce the computational cost by introducing surface variables and removing local degrees of freedom. - **Embedded - Hybridizable Discontinuous Galerkin (E - HDG) method**: Further reduce the computational cost, use continuous surface variables to approximate the velocity and magnetic field, and use discontinuous surface variables to approximate scalar variables (such as pressure and magnetic pressure). The key features of these two methods include: - **Divergence - free**: Ensure that the divergence of the velocity and magnetic field is zero. - **H(div) - conforming**: Ensure that the velocity and magnetic field conform to the H(div) space in the weak sense. - **High efficiency**: In particular, in three - dimensional, high - order, and fine - grid cases, the E - HDG method has significant computational advantages. ### Mathematical formulas The MHD equations can be expressed as: \[ -\frac{1}{Re}\Delta\mathbf{u}+\nabla p+(\mathbf{u}\cdot\nabla)\mathbf{u}+\kappa\mathbf{b}\times(\nabla\times\mathbf{b}) = \mathbf{g} \] \[ \nabla\cdot\mathbf{u} = 0 \] \[ \frac{\kappa}{Rm}\nabla\times(\nabla\times\mathbf{b})+\nabla r-\kappa\nabla\times(\mathbf{u}\times\mathbf{d})=\mathbf{f} \] \[ \nabla\cdot\mathbf{b} = 0 \] where: - $\mathbf{u}$ is the fluid velocity, - $\mathbf{b}$ is the magnetic field, - $p$ is the fluid pressure, - $r$ is the Lagrange multiplier related to the magnetic field divergence constraint, - $Re$ is the Reynolds number, - $Rm$ is the magnetic Reynolds number, - $\kappa$ is the coupling parameter, - $\mathbf{g}$ and $\mathbf{f}$ are external force terms. Through these methods, the author aims to provide an accurate and efficient numerical tool for solving complex magnetohydrodynamic problems.

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

High-order Discontinuous Galerkin Method for Vacuum Electron Device Modeling

An entropy stable nodal discontinuous Galerkin method for the resistive MHD equations: Continuous analysis and GLM divergence cleaning

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 Locally Divergence-Free Oscillation-Eliminating Discontinuous Galerkin Method for Ideal Magnetohydrodynamic Equations

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

Dual-robust iterative analysis of divergence-conforming IPDG FEM for thermally coupled inductionless MHD system

A hybridizable discontinuous Galerkin method for magnetic advection-diffusion problems

An exactly mass conserving space-time embedded-hybridized discontinuous Galerkin method for the Navier-Stokes equations on moving domains

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

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

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

A Divergence-Conforming Hybridized Discontinuous Galerkin Method for the Incompressible Reynolds Averaged Navier-Stokes Equations

Entropy-Stable Gauss Collocation Methods for Ideal Magneto-Hydrodynamics

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

eXtended Hybridizable Discontinous Galerkin (X-HDG) Method for Linear Convection-Diffusion Equations on Unfitted Domains

Hybridizable discontinuous Galerkin methods for solving the two-fluid plasma model

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

A robust solver for the finite element approximation of stationary incompressible MHD equations in 3D

A Multi-Dimensional High-Order Discontinuous Galerkin Method Based on Gas Kinetic Theory for Viscous Flow Computations

A DIVERGENCE-FREE ADAPTIVE MOVING MESH METHOD FOR IDEAL MAGNETOHYDRODYNAMICS