Adaptive finite element methods for computing band gaps in photonic crystals

S Giani, Ivan G Graham

Research output: Contribution to journalArticle

16 Citations (Scopus)

Abstract

In this paper we propose and analyse adaptive finite element methods for computing the band structure of 2D periodic photonic crystals. The problem can be reduced to the computation of the discrete spectra of each member of a family of periodic Hermitian eigenvalue problems on a unit cell, parametrised by a two-dimensional parameter - the quasimomentum. These eigenvalue problems involve non-coercive elliptic operators with generally discontinuous coefficients and are solved by adaptive finite elements. We propose an error estimator of residual type and show it is reliable and efficient for each eigenvalue problem in the family. In particular we prove that if the error estimator converges to zero then the distance of the computed eigenfunction from the true eigenspace also converges to zero and the computed eigenvalue converges to a true eigenvalue with double the rate. We also prove that if the distance of a computed sequence of approximate eigenfunctions from the true eigenspace approaches zero, then so must the error estimator. The results hold for eigenvalues of any multiplicity. We illustrate the benefits of the resulting adaptive method in practice, both for fully periodic structures and also for the computation of eigenvalues in the band gap of structures with defect, using the supercell method.
LanguageEnglish
Pages31-64
Number of pages34
JournalNumerische Mathematik
Volume121
Issue number1
DOIs
StatusPublished - 2012

Fingerprint

Adaptive Finite Element Method
Band Gap
Photonic crystals
Photonic Crystal
Error Estimator
Energy gap
Eigenvalue Problem
Eigenvalue
Finite element method
Eigenvalues and eigenfunctions
Computing
Eigenspace
Converge
Eigenfunctions
Zero
Periodic structures
Band structure
Adaptive Finite Elements
Discontinuous Coefficients
Band Structure

Cite this

Adaptive finite element methods for computing band gaps in photonic crystals. / Giani, S; Graham, Ivan G.

In: Numerische Mathematik, Vol. 121, No. 1, 2012, p. 31-64.

Research output: Contribution to journalArticle

@article{a6e55b2aab394cd0abfafc29c406caab,
title = "Adaptive finite element methods for computing band gaps in photonic crystals",
abstract = "In this paper we propose and analyse adaptive finite element methods for computing the band structure of 2D periodic photonic crystals. The problem can be reduced to the computation of the discrete spectra of each member of a family of periodic Hermitian eigenvalue problems on a unit cell, parametrised by a two-dimensional parameter - the quasimomentum. These eigenvalue problems involve non-coercive elliptic operators with generally discontinuous coefficients and are solved by adaptive finite elements. We propose an error estimator of residual type and show it is reliable and efficient for each eigenvalue problem in the family. In particular we prove that if the error estimator converges to zero then the distance of the computed eigenfunction from the true eigenspace also converges to zero and the computed eigenvalue converges to a true eigenvalue with double the rate. We also prove that if the distance of a computed sequence of approximate eigenfunctions from the true eigenspace approaches zero, then so must the error estimator. The results hold for eigenvalues of any multiplicity. We illustrate the benefits of the resulting adaptive method in practice, both for fully periodic structures and also for the computation of eigenvalues in the band gap of structures with defect, using the supercell method.",
author = "S Giani and Graham, {Ivan G}",
year = "2012",
doi = "10.1007/s00211-011-0425-9",
language = "English",
volume = "121",
pages = "31--64",
journal = "Numerische Mathematik",
issn = "0029-599X",
publisher = "Springer New York",
number = "1",

}

TY - JOUR

T1 - Adaptive finite element methods for computing band gaps in photonic crystals

AU - Giani, S

AU - Graham, Ivan G

PY - 2012

Y1 - 2012

N2 - In this paper we propose and analyse adaptive finite element methods for computing the band structure of 2D periodic photonic crystals. The problem can be reduced to the computation of the discrete spectra of each member of a family of periodic Hermitian eigenvalue problems on a unit cell, parametrised by a two-dimensional parameter - the quasimomentum. These eigenvalue problems involve non-coercive elliptic operators with generally discontinuous coefficients and are solved by adaptive finite elements. We propose an error estimator of residual type and show it is reliable and efficient for each eigenvalue problem in the family. In particular we prove that if the error estimator converges to zero then the distance of the computed eigenfunction from the true eigenspace also converges to zero and the computed eigenvalue converges to a true eigenvalue with double the rate. We also prove that if the distance of a computed sequence of approximate eigenfunctions from the true eigenspace approaches zero, then so must the error estimator. The results hold for eigenvalues of any multiplicity. We illustrate the benefits of the resulting adaptive method in practice, both for fully periodic structures and also for the computation of eigenvalues in the band gap of structures with defect, using the supercell method.

AB - In this paper we propose and analyse adaptive finite element methods for computing the band structure of 2D periodic photonic crystals. The problem can be reduced to the computation of the discrete spectra of each member of a family of periodic Hermitian eigenvalue problems on a unit cell, parametrised by a two-dimensional parameter - the quasimomentum. These eigenvalue problems involve non-coercive elliptic operators with generally discontinuous coefficients and are solved by adaptive finite elements. We propose an error estimator of residual type and show it is reliable and efficient for each eigenvalue problem in the family. In particular we prove that if the error estimator converges to zero then the distance of the computed eigenfunction from the true eigenspace also converges to zero and the computed eigenvalue converges to a true eigenvalue with double the rate. We also prove that if the distance of a computed sequence of approximate eigenfunctions from the true eigenspace approaches zero, then so must the error estimator. The results hold for eigenvalues of any multiplicity. We illustrate the benefits of the resulting adaptive method in practice, both for fully periodic structures and also for the computation of eigenvalues in the band gap of structures with defect, using the supercell method.

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

UR - http://dx.doi.org/10.1007/s00211-011-0425-9

U2 - 10.1007/s00211-011-0425-9

DO - 10.1007/s00211-011-0425-9

M3 - Article

VL - 121

SP - 31

EP - 64

JO - Numerische Mathematik

T2 - Numerische Mathematik

JF - Numerische Mathematik

SN - 0029-599X

IS - 1

ER -