Abstract:A novel finite element scheme is studied for solving the time-dependent Maxwell's equations on unstructured grids efficiently. Similar to the traditional Yee scheme, the method has one degree of freedom for most edges and a sparse inverse mass matrix. This allows for an efficient realization by explicit time-stepping without solving linear systems. The method is constructed by algebraic reduction of another underlying finite element scheme which involves two degrees of freedom for every edge. Mass-lumping and additional modifications are used in the construction of this method to allow for the mentioned algebraic reduction in the presence of source terms and lossy media later on. A full error analysis of the underlying method is developed which by construction also carries over to the reduced scheme and allows to prove convergence rates for the latter. The efficiency and accuracy of both methods are illustrated by numerical tests. The proposed schemes and their analysis can be extended to structured grids and in special cases the reduced method turns out to be algebraically equivalent to the Yee scheme. The analysis of this paper highlights possible difficulties in extensions of the Yee scheme to non-orthogonal or unstructured grids, discontinuous material parameters, and non-smooth source terms, and also offers potential remedies.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop an efficient finite - element solution scheme for time - dependent Maxwell's equations applicable to unstructured grids. Specifically, the author proposes a new method similar to the traditional Yee scheme, which can effectively solve Maxwell's equations on unstructured grids and maintain the advantage of explicit time - stepping, that is, without the need to solve linear systems at each time step.
### Problem Background
Maxwell's equations describe the propagation behavior of electromagnetic waves in media. For the problem of electromagnetic wave propagation in linear, non - dispersive, and lossy media, the basic mathematical model is:
\[
\varepsilon \frac{\partial^2 E}{\partial t^2} + \sigma \frac{\partial E}{\partial t} + \nabla \times (\nu \nabla \times E) = f, \quad \text{in } \Omega, \; t > 0,
\]
\[
\mathbf{n} \times (\nu \nabla \times E) = g, \quad \text{on } \partial \Omega, \; t > 0.
\]
where:
- \(\varepsilon\) and \(\sigma\) represent the dielectric constant and conductivity of the dielectric respectively;
- \(\nu\) represents the magnetic permeability;
- \(E\) is the electric field strength;
- \(f\) and \(g\) describe the excitation of volume or surface currents;
- \(\Omega\) is the computational domain.
### Existing Methods and Their Limitations
The traditional Yee scheme (such as the FDTD method) performs well on orthogonal structured grids, but has difficulties in dealing with non - orthogonal or unstructured grids, discontinuous material parameters, and non - smooth source terms. In addition, these methods usually require implicit time - stepping, resulting in the need to solve linear systems at each time step, thus reducing efficiency.
### Main Contributions of the Paper
To solve the above problems, the author proposes the following improvements:
1. **New Finite - Element Method**: A new finite - element method is constructed through algebraic reduction. This method has only one degree of freedom on most edges and has a sparse inverse mass matrix.
2. **Mass - Lumping Technique**: Using the mass - lumping technique and appropriate numerical integration, the new method can be explicitly time - stepped, avoiding the complexity of solving linear systems at each time step.
3. **Error Analysis**: A strict error analysis is carried out on the new method, which proves its convergence and shows its equivalence to the traditional Yee scheme in specific cases.
4. **Scope of Application**: This method is not only applicable to unstructured grids, but can also be extended to structured grids and cases with discontinuous material parameters.
### Conclusion
This paper proposes an efficient finite - element solution scheme for time - dependent Maxwell's equations applicable to unstructured grids, which solves the limitations of existing methods in dealing with complex geometries and material properties. The effectiveness and accuracy of this method are verified through strict error analysis and numerical experiments.