Numerical approximation based on a decoupled dimensionality reduction scheme for Maxwell eigenvalue problem

Abstract

We present a high‐accuracy numerical method based on a decoupled dimensionality reduction scheme for Maxwell eigenvalue problem in spherical domains. Using the orthogonality of vector spherical harmonics and the variable separation approach, we decompose the original problem into two classes of decoupled one‐dimensional TE mode and TM mode. For the TE mode, we establish a variational formulation and its discrete scheme and give the error estimations of the approximate eigenvalues and eigenfunctions. For the TM mode, it is different from TE mode which naturally meets the divergence‐free condition and will not generate some spurious eigenvalues. We design a numerical algorithm based on a parameterized method to filter out the spurious eigenvalues. Finally, some numerical results are presented to confirm the theoretical results and validate the algorithms.