Giani, S. and Graham, I.G. (2012) 'Adaptive finite element methods for computing band gaps in photonic crystals.', Numerische Mathematik., 121 (1). pp. 31-64.
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.
|Keywords:||65M50, 65M60, 65F15.|
|Full text:||(AM) Accepted Manuscript|
Download PDF (907Kb)
|Publisher Web site:||http://dx.doi.org/10.1007/s00211-011-0425-9|
|Publisher statement:||The final publication is available at Springer via http://dx.doi.org/10.1007/s00211-011-0425-9.|
|Record Created:||20 Apr 2015 09:35|
|Last Modified:||20 Apr 2015 10:39|
|Social bookmarking:||Export: EndNote, Zotero | BibTex|
|Look up in GoogleScholar | Find in a UK Library|