Abstract:In this paper, we consider mixed finite element semi-/full discretizations of the Rosensweig ferrofluid flow model. We first establish some regularity results for the model under several basic assumptions. Then we show that the energy stability of the weak solutions is preserved exactly for both the semi-discrete and fully discrete finite element solutions. Moreover, we prove the existence and uniqueness of the discrete solutions. We also derive optimal error estimates for the discrete schemes. Finally, we provide numerical experiments to verify the theoretical results.
What problem does this paper attempt to address?
This paper aims to solve the numerical simulation problem of the Rosensweig magnetohydrodynamic flow model. Specifically, it is to develop and analyze the energy - stable mixed finite - element method for this model. The following are the key problems that the paper attempts to solve:
1. **Establishing Regularity Results**:
- Firstly, the paper establishes regularity results for the Rosensweig magnetohydrodynamic flow model under some basic assumptions. These assumptions include the external magnetic field, initial conditions, and the properties of the solution domain.
2. **Maintaining Energy Stability**:
- The paper proves that the energy stability of the weak solution is exactly maintained in the semi - discrete and fully - discrete finite - element solutions. This is achieved by introducing appropriate energy estimates and regularity results.
3. **Existence and Uniqueness**:
- The paper proves the existence of the semi - discrete and fully - discrete finite - element solutions and, under certain conditions, proves their uniqueness. This ensures the effectiveness and reliability of the proposed numerical method.
4. **Error Estimation**:
- The paper derives the optimal error estimates of the discrete scheme. These estimates provide a theoretical guarantee for the accuracy of the numerical solution and verify the convergence of the method.
5. **Numerical Experiments**:
- In order to verify the theoretical results, the paper also conducts numerical experiments. These experiments show the performance and effectiveness of the proposed method in practical applications.
### Formula Summary
- **Energy Definition**:
\[
\mathcal{E}(t) := \rho \| \mathbf{u}(\cdot, t) \|^2 + \rho \kappa \| \boldsymbol{\omega}(\cdot, t) \|^2 + \| \mathbf{m}(\cdot, t) \|^2 + \mu_0 \| \mathbf{H}(\cdot, t) \|^2
\]
- **Dissipated Energy**:
\[
\mathcal{F}(t) := \eta \| \nabla \mathbf{u}(\cdot, t) \|^2 + \eta' \| \nabla \boldsymbol{\omega}(\cdot, t) \|^2 + (\eta' + \lambda') \| \text{div} \boldsymbol{\omega}(\cdot, t) \|^2 + \frac{1}{\tau} \| \mathbf{m}(\cdot, t) \|^2 + \chi_0 + (1 + \mu_0) \chi_0 \frac{1}{\tau} \| \mathbf{H}(\cdot, t) \|^2 + \sigma \| \text{div} \mathbf{m}(\cdot, t) \|^2 + \sigma \| \mathbf{k}(\cdot, t) \|^2 + \mu_0 \sigma \| \text{div} \mathbf{H}(\cdot, t) \|^2 + \zeta \| \text{curl} \mathbf{u}(\cdot, t) - 2 \boldsymbol{\omega}(\cdot, t) \|^2
\]
- **Energy Stability Inequality**:
\[
\mathcal{E}(t) + C_1 \int_0^t \mathcal{F}(s) \, ds \leq \mathcal{E}(0) + C_2 \int_0^t \left( \| \mathbf{H}_e(\cdot, s) \|^2_{H(\text{div})} + \| \partial_t \mathbf{H}_e(\cdot, s) \|^2 \right) ds
\]
### Conclusion
The paper solves the numerical simulation problem of the Rosensweig magnetohydrodynamic flow model through the above methods, ensuring energy stability, existence, and uniqueness.