Abstract:We propose a new mixed finite element method for the three-dimensional steady magnetohydrodynamic (MHD) kinematics equations for which the velocity of the fluid is given. Although prescribing the velocity field leads to a simpler model than full MHD equations, its conservative and efficient numerical methods are still active research topic. The distinctive feature of our discrete scheme is that the divergence-free conditions for current density and magnetic induction are both satisfied. To reach this goal, we use magnetic vector potential to represent magnetic induction and resort to H(div)-conforming element to discretize the current density. We develop an preconditioned iterative solver based on a block preconditioner for the algebraic systems arising from the discretization. Several numerical experiments are implemented to verify the divergence-free properties, the convergence rate of the finite element scheme and the robustness of the preconditioner.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a new hybrid finite element method for solving the three - dimensional steady - state magnetohydrodynamic (MHD) equations of motion. Specifically, the paper aims to ensure that the discrete approximations of the current density $\mathbf{J}$ and the magnetic field induction $\mathbf{B}$ strictly satisfy the divergence - free conditions (i.e., $\nabla\cdot\mathbf{J} = 0$ and $\nabla\cdot\mathbf{B} = 0$). To achieve this goal, the author adopts the following strategies:
1. **Use of magnetic vector potential**: By introducing the magnetic vector potential $\mathbf{A}$, such that $\mathbf{B}=\nabla\times\mathbf{A}$, thus naturally ensuring that $\nabla\cdot\mathbf{B} = 0$.
2. **H(div)-conforming elements**: Use H(div)-conforming elements to discretize the current density $\mathbf{J}$ to ensure that $\nabla\cdot\mathbf{J} = 0$.
3. **Pre - conditioned iterative solver**: Develop an iterative solver based on block pre - conditioning for solving the algebraic system generated by discretization to improve the solution efficiency.
The focus of the paper is on proposing a finite element method that can accurately maintain the divergence - free property and designing an effective pre - conditioned iterative solver. This method has been verified in numerical experiments, demonstrating its good convergence and robustness.
### Key issues
- **Divergence - free conditions**: Ensure that the discretized current density and magnetic field induction strictly satisfy the divergence - free conditions.
- **Efficient solution**: Develop an efficient pre - conditioned iterative solver to deal with complex saddle - point systems.
### Application background
This research is of great significance for the application of magnetohydrodynamics, especially in fields such as generators, solar dynamos, salt water and geodynamos, where the Lorentz force $\mathbf{J}\times\mathbf{B}$ is relatively small and the electromagnetic effect on the fluid can be ignored. Therefore, the change of the magnetic field can be studied through the given velocity field $\mathbf{w}$.
### Mathematical model
The main equations involved in the paper are:
\[
\begin{aligned}
&\frac{\partial\mathbf{B}}{\partial t}+\nabla\times\mathbf{E} = 0,\\
&\nabla\times\mathbf{H}=\mathbf{J},\\
&\sigma(\mathbf{E}+\mathbf{w}\times\mathbf{B})=\mathbf{J},\\
&\nabla\cdot\mathbf{J} = 0,\quad\nabla\cdot\mathbf{B} = 0.
\end{aligned}
\]
By introducing the magnetic vector potential $\mathbf{A}$ and the scalar electric potential $\phi$, and using the Coulomb gauge $\nabla\cdot\mathbf{A} = 0$, the above equations are transformed into a new PDE system, thus ensuring the divergence - free property.
### Conclusion
The method proposed in the paper performs well in numerical experiments, demonstrating its advantages in maintaining the divergence - free property and solution efficiency. Future research can be further extended to transient cases, more robust pre - processing methods and coupling with the Navier - Stokes equations.