Abstract:It has been extensively studied in the literature that solving Maxwell equations is very sensitive to the mesh structure, space conformity and solution regularity. Roughly speaking, for almost all the methods in the literature, optimal convergence for low-regularity solutions heavily relies on conforming spaces and highly-regular simplicial meshes. This can be a significant limitation for many popular methods based on polytopal meshes in the case of inhomogeneous media, as the discontinuity of electromagnetic parameters can lead to quite low regularity of solutions near media interfaces, and potentially worsened by geometric singularities, making many popular methods based on broken spaces, non-conforming or polytopal meshes particularly challenging to apply. In this article, we present a virtual element method for solving an indefinite time-harmonic Maxwell equation in 2D inhomogeneous media with quite arbitrary polytopal meshes, and the media interface is allowed to have geometric singularity to cause low regularity. There are two key novelties: (i) the proposed method is theoretically guaranteed to achieve robust optimal convergence for solutions with merely $\mathbf{H}^{\theta}$ regularity, $\theta\in(1/2,1]$; (ii) the polytopal element shape can be highly anisotropic and shrinking, and an explicit formula is established to describe the relationship between the shape regularity and solution regularity. Extensive numerical experiments will be given to demonstrate the effectiveness of the proposed method.
What problem does this paper attempt to address?
This paper attempts to solve the low - regularity problem encountered when solving the time - harmonic Maxwell equations in inhomogeneous media. Specifically, traditional methods often require meshes that conform to interfaces and have regular shapes to ensure optimal convergence when dealing with such problems. However, in practical applications, the discontinuity of electromagnetic parameters can lead to low regularity of the solution near the media interface, which makes many methods based on unstructured or polygonal meshes difficult to apply.
### Core problems of the paper
1. **Low - regularity problem**: In inhomogeneous media, due to the discontinuity of electromagnetic parameters (such as magnetic permeability, permittivity, and conductivity), the regularity of the solution near the media interface is low, usually only \(H^{\theta}\) regularity, where \(\theta\in(0, 1]\). This low regularity poses challenges to the convergence and accuracy of numerical methods.
2. **Difficulty in mesh generation**: To ensure the optimal convergence of traditional finite - element methods, it is usually necessary to generate meshes that conform to interfaces and have regular shapes. However, for complex media distributions and geometric shapes, generating such meshes is very difficult.
3. **Limitations of existing methods**: Many existing methods (such as the discontinuous Galerkin method, the hybrid DG method, etc.) rely on broken spaces when dealing with polygonal meshes, and these methods are difficult to maintain optimal convergence in the case of low regularity.
### Solutions
To solve the above problems, this paper proposes a new method based on the Virtual Element Method (VEM). The main contributions of this method are as follows:
1. **Theoretically guaranteed optimal convergence**: This method can theoretically guarantee optimal convergence when the solution has only \(H^{\theta}\) regularity (\(\theta\in(1/2, 1]\)). This is achieved by introducing a new Hodge mapping and optimized stabilization techniques.
2. **Support for highly anisotropic meshes**: This method can handle highly anisotropic polygonal meshes and establish an explicit relational formula between mesh - shape regularity and solution regularity. This makes the method more flexible when dealing with complex geometric shapes and media - interface problems.
3. **Verification by extensive numerical experiments**: Through a large number of numerical experiments, the effectiveness of this method in dealing with complex media interfaces (including complex topological structures, geometric singularities, and thin layers) has been verified.
### Summary
This paper aims to prove, through theoretical analysis and numerical experiments, that the Virtual Element Method can effectively solve the time - harmonic Maxwell equations in inhomogeneous media on polygonal meshes and highly anisotropic meshes, and maintain optimal convergence when the solution has only low regularity. This result provides a new and effective tool for dealing with complex electromagnetic phenomena.