Abstract:We study a system of Maxwell's equations that describes the time evolution of electromagnetic fields with an additional electric scalar variable to make the system amenable to a mixed finite element spatial discretization. We demonstrate stability and energy conservation for the variational formulation of this Maxwell's system. We then discuss two implicit, energy conserving schemes for its temporal discretization: the classical Crank-Nicholson scheme and an implicit leapfrog scheme. We next show discrete stability and discrete energy conservation for the semi-discretization using these two time integration methods. We complete our discussion by showing that the error for the full discretization of the Maxwell's system with each of the two implicit time discretization schemes and with spatial discretization through a conforming sequence of de Rham finite element spaces converges quadratically in the step size of the time discretization and as an appropriate polynomial power of the mesh parameter in accordance with the choice of approximating polynomial spaces. Our results for the Crank-Nicholson method are generally well known but have not been demonstrated for this Maxwell's system. Our implicit leapfrog scheme is a new method to the best of our knowledge and we provide a complete error analysis for it. Finally, we show computational results to validate our theoretical claims using linear and quadratic Whitney forms for the finite element discretization for some model problems in two and three spatial dimensions.
What problem does this paper attempt to address?
### What problem does this paper attempt to solve?
This paper aims to solve the problem of numerical discretization of Maxwell's equations in the time - evolution process of electromagnetic fields, especially how to achieve high - order accurate time and space discretization while maintaining energy conservation and stability. Specifically, the author has studied a system of Maxwell's equations with an additional electric scalar variable \(p\) to make it suitable for mixed finite - element space discretization.
#### Main problems include:
1. **Introduction of additional variables**:
- Introduce a physically fictitious electric scalar variable \(p\), which is related to the time - varying divergence of the electric field, so that Maxwell's equations can be discretized by mixed finite elements.
- This additional variable \(p\) helps to directly approximate \(\nabla\cdot(\varepsilon E)\) and makes the system easier to handle.
2. **Stability and energy conservation**:
- The stability of the variational form and its energy conservation properties are studied.
- Two implicit time - discretization methods are proposed: the classical Crank - Nicholson scheme and the implicit leapfrog scheme, and their discrete stability and discrete energy conservation are proved.
3. **Error analysis**:
- For full discretization (i.e., both time and space are discretized), the de Rham finite - element space sequence is used for space discretization, and an error analysis is carried out for each implicit time - discretization scheme.
- It is proved that the error converges as the square of the time step and converges as an appropriate polynomial power of the mesh parameter.
4. **Verification of theoretical results**:
- Through the finite - element discretization of linear and quadratic Whitney forms, the computational results are shown in two - dimensional and three - dimensional model problems to verify the validity of the theoretical analysis.
#### Summary of mathematical formulas:
- Maxwell's equations:
\[
\begin{aligned}
&\frac{\partial p}{\partial t}+\nabla\cdot(\varepsilon E)=f_p, \\
&\nabla p+\varepsilon\frac{\partial E}{\partial t}-\nabla\times H = f_E, \\
&\mu\frac{\partial H}{\partial t}+\nabla\times E=f_H.
\end{aligned}
\]
- Boundary conditions:
\[
p = 0,\quad E\times n=0,\quad H\cdot n = 0\quad\text{on}\quad\partial\Omega\times(0, T].
\]
- Initial conditions:
\[
p(x, 0)=p_0(x),\quad E(x, 0)=E_0(x),\quad H(x, 0)=H_0(x)\quad\text{for}\quad x\in\Omega.
\]
#### Conclusion:
The main contribution of this paper is to propose a new implicit leapfrog scheme and carry out a detailed error analysis, filling some gaps in this field. In addition, the paper also verifies the effectiveness and feasibility of these methods in actual calculations.