Abstract:In this paper, we propose a conservative nonconforming virtual element method for the full stationary incompressible magnetohydrodynamics model. We leverage the virtual element satisfactory divergence-free property to ensure mass conservation for the velocity field. The condition of the well-posedness of the proposed method, as well as the stability are derived. We establish optimal error estimates in the discrete energy norm for both the velocity and magnetic field. Furthermore, by employing a new technique, we obtain the optimal error estimates in $L^2$-norm without any additional conditions. Finally, numerical experiments are presented to validate the theoretical analysis. In the implementation process, we adopt the effective Oseen iteration to handle the nonlinear system.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to propose a conservative non - conforming virtual element method (VEM) for the two - dimensional fully stationary incompressible magnetohydrodynamics (MHD) equations. Specifically, the paper focuses on how to maintain physical properties, especially mass conservation, during the numerical solution process and ensure that the divergence of the discrete velocity field is zero. In addition, the paper also aims to derive optimal error estimates and verify the effectiveness and stability of the proposed numerical method.
### Main problems and objectives
1. **Maintaining physical characteristics**:
- Ensure mass conservation of the velocity field by using the virtual element method with satisfactory divergence - free characteristics.
- Ensure that the discrete velocity field is piecewise divergence - free, that is, \( \text{div} \, u = 0 \).
2. **Stability and well - posedness of the numerical method**:
- Derive the well - posedness conditions and stability of the proposed method.
- Prove the uniqueness and stability of the discrete scheme.
3. **Error estimates**:
- Establish optimal error estimates in the discrete energy norm, applicable to the velocity field and the magnetic field.
- Obtain optimal error estimates in the \( L^2 \)-norm by introducing new techniques without additional conditions.
4. **Verification by numerical experiments**:
- Verify the correctness of the theoretical analysis through numerical experiments.
- Use the effective Oseen iteration method to handle the nonlinear system and ensure the accuracy of the numerical results.
### Mathematical model and method
The object of study in this paper is the two - dimensional fully stationary incompressible magnetohydrodynamics equations, and its mathematical model is as follows:
\[
\begin{cases}
-\frac{1}{R_\nu} \Delta u+(\nabla u)u + \nabla p - Sc \, \text{curl} \, b\times b = f, & \text{in } \Omega, \\
\frac{1}{R_m} Sc \, \text{curl} \, (\text{curl} \, b)-Sc \, \text{curl} \, (u\times b) = g, & \text{in } \Omega, \\
\text{div} \, u = 0, & \text{in } \Omega, \\
\text{div} \, b = 0, & \text{in } \Omega, \\
u = 0, & \text{on } \partial \Omega, \\
b\cdot n = 0, \quad n\times \text{curl} \, b = 0, & \text{on } \partial \Omega,
\end{cases}
\]
where \( u \) is the velocity field, \( b \) is the magnetic field, \( p \) is the pressure field; \( R_\nu \) is the hydrodynamic Reynolds number, \( R_m \) is the magnetic Reynolds number, \( Sc \) is the coupling coefficient; \( f \) and \( g \) are source - term functions.
### Method innovation points
Compared with the existing work, this research has the following two main differences:
1. **Model complexity**: The model discussed in this paper is more complex than the resistive - type MHD.
2. **Arbitrary - order virtual elements**: Existing literature [3, 7] mainly focuses on low - order virtual elements, while this paper develops an arbitrary - order (\( k\geq1 \)) scheme, which is suitable for shape - regular polygonal meshes.
Through these improvements, the paper not only improves the accuracy and efficiency of the numerical method but also ensures strict compliance with physical characteristics.