Mixed Finite Element Method for 2d Vector Maxwell's Eigenvalue Problem in Anisotropic Media

Abstract

Abstract-It is well known that the conventional edge element method in solving vector Maxwell's eigenvalue problem will lead to the presence of nonphysical zero eigenvalues. This paper uses the mixed finite element method to suppress the presence of these nonphysical zero eigenvalues for 2D vector Maxwell's eigenvalue problem in anisotropic media. We introduce a Lagrangian multiplier to deal with the constraint of divergence-free condition. Our method is based on employing the first-order edge element basis functions to expand the electric field and linear nodal element basis functions to expand the Lagrangian multiplier. Our numerical experiments show that this method can successfully remove all nonphysical zero and nonzero eigenvalues. We verify that when the cavity has a connected perfect electric boundary, then there is no physical zero eigenvalue. Otherwise, the number of physical zero eigenvalues is one less than the number of disconnected perfect electric boundaries.

Solving 2D Linear Isotropic Elastodynamics by Means of Scalar Potentials: A New Challenge for Finite Elements

“…As it is usual for treating equality constraints (see [19,20,21]), we are going to introduce a Lagrange multiplier in the space ker T . This leads us to write the following mixed problem…”

“…The last two equations of (146) show that it is equivalent to say that the solution ϕ h belongs to a subspace of V h of co-dimension dim M h + 3. This subspace is not a subspace of V N , in other words we realize a non-conform approximation of the space V N and that is why the multiplier η h is not 0 (see remark 9.1 in [20] and [21] for similar situations).…”

A Sufficient Condition of Having Independent TE and TM Modes in a Waveguide Filled with Homogenous Anisotropic Lossless Medium

“…. The transverse components of electromagnetic field are directly obtained from the equations ( 10)- (11). This is so called TEM mode.…”

“…t e t ) = 0. But we can not omit the divergence-free condition ∇ t • (µ −1 t e t ) = 0 in numerical computation, otherwise PDEs (25) will introduce spurious zero modes [11]. Because ∇ t • (µ −1 t e t ) = 0 and (18), we can get 1 µ+ab ∇ t • ( t e t ) = 0, then ∇ t • ( t e t ) = 0, which is just Gauss' law for electric field.…”

A Two-Level Preconditioned Helmholtz-Jacobi-Davidson Method for the Maxwell Eigenvalue Problem

“…In this paper, we develop an efficient numerical algorithm for solving the Maxwell eigenvalue problem which plays an important role in computational electromagnetism (see, e.g., [3,4,12,15,16,18]). The governing equations are…”