Abstract:In this paper, we introduce a novel unfitted finite element method to solve the quad-curl interface problem. We adapt Nitsche's method for curlcurl-conforming elements and double the degrees of freedom on interface elements. To ensure stability, we incorporate ghost penalty terms and a discrete divergence-free term. We establish the well-posedness of our method and demonstrate an optimal error bound in the discrete energy norm. We also analyze the stiffness matrix's condition number. Our numerical tests back up our theory on convergence rates and condition numbers.
What problem does this paper attempt to address?
### What problem does this paper attempt to solve?
This paper aims to solve the **quad - curl interface problem**. Specifically, the author proposes a new **unfitted finite element method** to solve this problem, and ensures the stability of the numerical method by introducing the Nitsche method and doubling the degrees of freedom on the interface elements. In addition, in order to further enhance the stability and ensure the discrete divergence - free condition, the author also introduces **ghost penalty terms**.
#### Background of the quad - curl interface problem
The quad - curl problem is of great significance in multi - physical field simulations, especially in simulating magnetized plasma in magnetohydrodynamics (MHD). It also plays a key role in approximating the Maxwell transmission eigenvalue problem, especially in the inverse electromagnetic scattering theory. Although many numerical methods have been developed to solve the quad - curl problem, there are still few studies on the quad - curl problem with interfaces.
#### Main challenges of the research
1. **Stability problem**: To ensure the stability of the numerical method, the author introduces two ghost penalty terms, and establishes interpolation estimates and weak Galerkin orthogonal relations in the discrete energy norm.
2. **Error estimate**: By combining the estimates of the ghost penalty terms, the author derives the error estimate between the numerical solution and the exact solution.
3. **Condition number estimate**: The author also studies the condition number estimate of the stiffness matrix, which is crucial for evaluating the performance of the numerical method.
#### Method overview
The method proposed by the author is based on the Cut Finite Element Method (CutFEM), which is a non - fitting method that allows the interface to pass through the grid. The specific steps include:
- **Introducing the Nitsche method**: Used to handle the interface conditions and ensure that the transmission conditions are weakly satisfied on the interface.
- **Increasing degrees of freedom**: Doubling the degrees of freedom on the interface elements to improve the accuracy of the method.
- **Introducing ghost penalty terms**: To improve the stability of the numerical method and ensure the discrete divergence - free condition.
- **Discretizing the weak form**: Defining the discrete weak form and the lowest - order curl - curl conforming finite element space.
- **Existence and uniqueness**: Proving the existence and uniqueness of the discrete scheme.
- **Error analysis**: Deriving a priori error estimates in the discrete norm.
- **Condition number estimate**: Providing a condition number estimate of the stiffness matrix.
### Summary
This paper solves the quad - curl interface problem by introducing a new unfitted finite element method, and theoretically proves the existence, uniqueness and convergence of this method. Through a series of numerical tests, the effectiveness and accuracy of this method are verified.