Umeå University's logo

In this paper, we discuss problems arising when computing resonances with a finite element method. In the pre-asymptotic regime, we detect for the one dimensional case, spurious solutions in finite element computations of resonances when the computational domain is truncated with a perfectly matched layer (PML) as well as with a Dirichlet-to-Neumann map (DtN). The new test is based on the Lippmann–Schwinger equation and we use computations of the pseudospectrum to show that this is a suitable choice. Numerical simulations indicate that the presented test can distinguish between spurious eigenvalues and true eigenvalues also in difficult cases.

We present an efficient procedure for computing resonances and resonant modes of Helmholtz problems posed in exterior domains. The problem is formulated as a nonlinear eigenvalue problem (NEP), where the nonlinearity arises from the use of a Dirichlet-to-Neumann map, which accounts for modeling unbounded domains. We consider a variational formulation and show that the spectrum consists of isolated eigenvalues of finite multiplicity that only can accumulate at infinity. The proposed method is based on a high order finite element discretization combined with a specialization of the Tensor Infinite Arnoldi method (TIAR). Using Toeplitz matrices, we show how to specialize this method to our specific structure. In particular we introduce a pole cancellation technique in order to increase the radius of convergence for computation of eigenvalues that lie close to the poles of the matrix-valued function. The solution scheme can be applied to multiple resonators with a varying refractive index that is not necessarily piecewise constant. We present two test cases to show stability, performance and numerical accuracy of the method. In particular the use of a high order finite element discretization together with TIAR results in an efficient and reliable method to compute resonances.

In this paper we consider scattering resonance computations in optics when the resonators consist of frequency dependent and lossy materials, such as metals at optical frequencies. The proposed computational approach combines a novel hp-FEM strategy, based on dispersion analysis for complex frequencies, with a fast implementation of the nonlinear eigenvalue solver NLEIGS. Numerical computations illustrate that the pre-asymptotic phase is significantly reduced compared to standard uniform h and p strategies. Moreover, the efficiency grows with the refractive index contrast, which makes the new strategy highly attractive for metal-dielectric structures. The hp-refinement strategy together with the efficient parallel code result in highly accurate approximations and short runtimes on multi processor platforms.

In this paper, we consider a sorting scheme for potentially spurious scattering resonant pairs in one- and two-dimensional electromagnetic problems and in three-dimensional acoustic problems. The novel sorting scheme is based on a Lippmann-Schwinger type of volume integral equation and can, therefore, be applied to structures with graded materials as well as to configurations including piece-wise constant material properties. For TM/TE polarized electromagnetic waves and for acoustic waves, we compute first approximations of scattering resonances with finite elements. Then, we apply the novel sorting scheme to the computed eigenpairs and use it to mark potentially spurious solutions in electromagnetic and acoustic scattering resonances computations at a low computational cost. Several test cases with Drude-Lorentz dielectric resonators as well as with graded material properties are considered.