A robust numerical algorithm for computing Maxwell's transmission eigenvalue problems

Tsung-Min Hwang, Wei Qiang Huang, Wen Wei Lin

Research output: Contribution to journalArticle

6 Citations (Scopus)

Abstract

We study a robust and efficient eigensolver for computing a few smallest positive eigenvalues of the three-dimensional Maxwell's transmission eigenvalue problem. The discretized governing equations by the Nédélec edge element result in a large-scale quadratic eigenvalue problem (QEP) for which the spectrum contains many zero eigenvalues and the coefficient matrices consist of patterns in the matrix form XY-1Z, both of which prevent existing eigenvalue solvers from being efficient. To remedy these difficulties, we rewrite the QEP as a particular nonlinear eigenvalue problem and develop a secant-type iteration, together with an indefinite locally optimal block preconditioned conjugate gradient (LOBPCG) method, to sequentially compute the desired positive eigenvalues. Furthermore, we propose a novel method to solve the linear systems in each iteration of LOBPCG. Intensive numerical experiments show that our proposed method is robust, although the desired real eigenvalues are surrounded by complex eigenvalues.

Original languageEnglish
Pages (from-to)A2403-A2423
JournalSIAM Journal on Scientific Computing
Volume37
Issue number5
DOIs
Publication statusPublished - 2015 Jan 1

Fingerprint

Transmission Problem
Robust Algorithm
Numerical Algorithms
Eigenvalue Problem
Eigenvalue
Conjugate gradient method
Computing
Quadratic Eigenvalue Problem
Linear systems
Iteration
Edge Elements
Preconditioned Conjugate Gradient
Preconditioned Conjugate Gradient Method
Nonlinear Eigenvalue Problem
Experiments
Chord or secant line
Governing equation
Linear Systems
Numerical Experiment
Three-dimensional

Keywords

  • LOBPCG
  • Maxwell's equations
  • Quadratic eigenvalue problems
  • Secant-type iteration
  • Transmission eigenvalues

ASJC Scopus subject areas

  • Computational Mathematics
  • Applied Mathematics

Cite this

A robust numerical algorithm for computing Maxwell's transmission eigenvalue problems. / Hwang, Tsung-Min; Huang, Wei Qiang; Lin, Wen Wei.

In: SIAM Journal on Scientific Computing, Vol. 37, No. 5, 01.01.2015, p. A2403-A2423.

Research output: Contribution to journalArticle

@article{9106533ede4846ca9ca2edf0f68153bb,
title = "A robust numerical algorithm for computing Maxwell's transmission eigenvalue problems",
abstract = "We study a robust and efficient eigensolver for computing a few smallest positive eigenvalues of the three-dimensional Maxwell's transmission eigenvalue problem. The discretized governing equations by the N{\'e}d{\'e}lec edge element result in a large-scale quadratic eigenvalue problem (QEP) for which the spectrum contains many zero eigenvalues and the coefficient matrices consist of patterns in the matrix form XY-1Z, both of which prevent existing eigenvalue solvers from being efficient. To remedy these difficulties, we rewrite the QEP as a particular nonlinear eigenvalue problem and develop a secant-type iteration, together with an indefinite locally optimal block preconditioned conjugate gradient (LOBPCG) method, to sequentially compute the desired positive eigenvalues. Furthermore, we propose a novel method to solve the linear systems in each iteration of LOBPCG. Intensive numerical experiments show that our proposed method is robust, although the desired real eigenvalues are surrounded by complex eigenvalues.",
keywords = "LOBPCG, Maxwell's equations, Quadratic eigenvalue problems, Secant-type iteration, Transmission eigenvalues",
author = "Tsung-Min Hwang and Huang, {Wei Qiang} and Lin, {Wen Wei}",
year = "2015",
month = "1",
day = "1",
doi = "10.1137/15M1018927",
language = "English",
volume = "37",
pages = "A2403--A2423",
journal = "SIAM Journal of Scientific Computing",
issn = "0036-1445",
publisher = "Society for Industrial and Applied Mathematics Publications",
number = "5",

}

TY - JOUR

T1 - A robust numerical algorithm for computing Maxwell's transmission eigenvalue problems

AU - Hwang, Tsung-Min

AU - Huang, Wei Qiang

AU - Lin, Wen Wei

PY - 2015/1/1

Y1 - 2015/1/1

N2 - We study a robust and efficient eigensolver for computing a few smallest positive eigenvalues of the three-dimensional Maxwell's transmission eigenvalue problem. The discretized governing equations by the Nédélec edge element result in a large-scale quadratic eigenvalue problem (QEP) for which the spectrum contains many zero eigenvalues and the coefficient matrices consist of patterns in the matrix form XY-1Z, both of which prevent existing eigenvalue solvers from being efficient. To remedy these difficulties, we rewrite the QEP as a particular nonlinear eigenvalue problem and develop a secant-type iteration, together with an indefinite locally optimal block preconditioned conjugate gradient (LOBPCG) method, to sequentially compute the desired positive eigenvalues. Furthermore, we propose a novel method to solve the linear systems in each iteration of LOBPCG. Intensive numerical experiments show that our proposed method is robust, although the desired real eigenvalues are surrounded by complex eigenvalues.

AB - We study a robust and efficient eigensolver for computing a few smallest positive eigenvalues of the three-dimensional Maxwell's transmission eigenvalue problem. The discretized governing equations by the Nédélec edge element result in a large-scale quadratic eigenvalue problem (QEP) for which the spectrum contains many zero eigenvalues and the coefficient matrices consist of patterns in the matrix form XY-1Z, both of which prevent existing eigenvalue solvers from being efficient. To remedy these difficulties, we rewrite the QEP as a particular nonlinear eigenvalue problem and develop a secant-type iteration, together with an indefinite locally optimal block preconditioned conjugate gradient (LOBPCG) method, to sequentially compute the desired positive eigenvalues. Furthermore, we propose a novel method to solve the linear systems in each iteration of LOBPCG. Intensive numerical experiments show that our proposed method is robust, although the desired real eigenvalues are surrounded by complex eigenvalues.

KW - LOBPCG

KW - Maxwell's equations

KW - Quadratic eigenvalue problems

KW - Secant-type iteration

KW - Transmission eigenvalues

UR - http://www.scopus.com/inward/record.url?scp=84945919789&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84945919789&partnerID=8YFLogxK

U2 - 10.1137/15M1018927

DO - 10.1137/15M1018927

M3 - Article

AN - SCOPUS:84945919789

VL - 37

SP - A2403-A2423

JO - SIAM Journal of Scientific Computing

JF - SIAM Journal of Scientific Computing

SN - 0036-1445

IS - 5

ER -