Abstract:This paper develops an hybridizable discontinuous Galerkin (HDG) finite element method of arbitrary order for the steady thermally coupled incompressible Magnetohydrodynamics (MHD) flow. The HDG scheme uses piecewise polynomials of degrees $k(k\geq 1),k,k-1,k-1$, and $k$ respectively for the approximations of the velocity, the magnetic field, the pressure, the magnetic pseudo-pressure, and the temperature in the interior of elements, and uses piecewise polynomials of degree $k$ for their numerical traces on the interfaces of elements. The method is shown to yield globally divergence-free approximations of the velocity and magnetic fields. Existence and uniqueness results for the discrete scheme are given and optimal a priori error estimates are derived. Numerical experiments are provided to verify the obtained theoretical results.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop an arbitrary - order hybrid discontinuous Galerkin (HDG) finite element method applicable to steady - state thermally - coupled incompressible magnetohydrodynamics (MHD) flows. Specifically, the paper focuses on how to ensure that the proposed numerical method can provide globally divergence - free approximations of the velocity and magnetic fields, and guarantee the existence, uniqueness, and optimal a priori error estimates of the discrete scheme.
### Problem Background
The magnetohydrodynamics (MHD) equations describe the basic physical laws of the interaction between conductive fluids and magnetic fields and are widely used in the engineering field. This paper considers the steady - state thermally - coupled incompressible MHD model, which is a system coupled by the incompressible Navier - Stokes equations, Maxwell equations, and heat conduction equations. These equations need to satisfy two important divergence constraint conditions during the solution process:
1. **Mass Conservation**: The divergence of the velocity field is zero ($\nabla \cdot u = 0$).
2. **Magnetic Flux Conservation**: The divergence of the magnetic field is zero ($\nabla \cdot B = 0$).
If the numerical method cannot well preserve these physical properties, it may lead to numerical instability.
### Main Contributions of the Paper
The paper proposes a numerical scheme based on the HDG method, which has the following characteristics:
- **Globally Divergence - Free**: The approximations of the velocity and magnetic fields are globally divergence - free, which makes the method pressure - robust.
- **No Need for Large Parameters**: There is no need to introduce a "sufficiently large" stabilization parameter.
- **Existence and Uniqueness**: The existence and uniqueness results of the discrete scheme are established.
- **Optimal Error Estimates**: The optimal a priori error estimates are derived.
### Mathematical Model
Consider the following steady - state thermally - coupled incompressible MHD model:
\[
\begin{aligned}
-\frac{1}{H_a^2} \Delta u + \frac{1}{N} (u \cdot \nabla) u + \nabla p - \frac{1}{R_m} (\nabla \times B) \times B &= f_1 - \frac{Gr}{N R_e^2} T g, & \text{in } \Omega, \\
\nabla \cdot u &= 0, & \text{in } \Omega, \\
\frac{1}{R_m} \nabla \times (\nabla \times B) - \nabla \times (u \times B) + \nabla r &= f_2, & \text{in } \Omega, \\
\nabla \cdot B &= 0, & \text{in } \Omega, \\
-\frac{1}{Pr Re} \Delta T + (u \cdot \nabla) T &= f_3, & \text{in } \Omega,
\end{aligned}
\]
where:
- $H_a$ is the Hartmann number,
- $N$ is the interaction parameter,
- $Re$ is the Reynolds number,
- $Pr$ is the Prandtl number,
- $R_m$ is the magnetic Reynolds number,
- $Gr$ is the Grashof number,
- $g$ is the gravitational acceleration vector,
- $f_1$ and $f_2$ are external force terms,
- $f_3$ is the heat source term.
### Boundary Conditions
All variables satisfy homogeneous boundary conditions on the boundary:
\[
u|_{\partial \Omega} = 0, \quad (B \times n)|_{\partial \Omega} = 0, \quad T|_{\partial \Omega} = 0, \quad r|_{\partial \Omega} = 0.
\]