Choice of the perfectly matched layer boundary condition for iterative solvers of the frequency-domain Maxwell's equations

We show that the performance of iterative solvers of the frequency-domain Maxwell's equations is greatly affected by the kind of the perfectly matched layer (PML) used. In particular, we demonstrate that using the stretchedcoordinate PML (SC-PML) results in significantly faster convergence speed as compared with using the uniaxial PML (UPML). Such a difference in convergence behavior is explained by an analysis of the condition number of the coefficient matrices. Additionally, we develop a diagonal preconditioning scheme that significantly improves solver performance when UPML is used.

[1]  Stephen D. Gedney,et al.  A parallel finite-element tearing and interconnecting algorithm for solution of the vector wave equation with PML absorbing medium , 2000 .

[2]  F.-J. Schmuckle,et al.  Optimizing the FDFD Method in Order to Minimize PML-Related Numerical Problems , 2007, 2007 IEEE/MTT-S International Microwave Symposium.

[3]  Shanhui Fan,et al.  Choice of the perfectly matched layer boundary condition for iterative solvers of the frequency-domain Maxwell's equations , 2012, Other Conferences.

[4]  E. Palik.  Handbook of Optical Constants of Solids , 1997 .

[5]  G. Veronis,et al.  Modes of Subwavelength Plasmonic Slot Waveguides , 2007, Journal of Lightwave Technology.

[6]  Shanhui Fan,et al.  Theoretical investigation of compact couplers between dielectric slab waveguides and two-dimensional metal-dielectric-metal plasmonic waveguides. , 2007, Optics express.

[7]  S..  A Perfectly Matched Anisotropic Absorber for Use as an Absorbing Boundary Condition-Antennas and Propagation, IEEE Transactions on , 2004 .

[8]  R. Freund,et al.  QMR: a quasi-minimal residual method for non-Hermitian linear systems , 1991 .

[9]  S. Gedney.  An anisotropic perfectly matched layer-absorbing medium for the truncation of FDTD lattices , 1996 .

[10]  Jin-Fa Lee,et al.  A perfectly matched anisotropic absorber for use as an absorbing boundary condition , 1995 .

[11]  Steven G. Johnson,et al.  Perturbation theory for anisotropic dielectric interfaces, and application to subpixel smoothing of discretized numerical methods. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[12]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[13]  R. W. Christy,et al.  Optical Constants of the Noble Metals , 1972 .

[14]  J. L. Volakis,et al.  A robust iterative scheme for FEM applications terminated by the perfectly matched layer (PML) absorbers , 1998, Proceedings of the Fifteenth National Radio Science Conference. NRSC '98 (Cat. No.98EX109).

[15]  Weng Cho Chew,et al.  A 3D perfectly matched medium from modified maxwell's equations with stretched coordinates , 1994 .

[16]  Carey M. Rappaport,et al.  Perfectly matched absorbing boundary conditions based on anisotropic lossy mapping of space , 1995 .

[17]  Raj Mittra,et al.  A new look at the perfectly matched layer (PML) concept for the reflectionless absorption of electromagnetic waves , 1995 .

[18]  Fernando L. Teixeira,et al.  General closed-form PML constitutive tensors to match arbitrary bianisotropic and dispersive linear media , 1998 .

[19]  James G. Berryman,et al.  FDFD: a 3D finite-difference frequency-domain code for electromagnetic induction tomography , 2001 .

[20]  Lloyd N. Trefethen,et al.  How Fast are Nonsymmetric Matrix Iterations? , 1992, SIAM J. Matrix Anal. Appl..

[21]  B. Stupfel,et al.  A study of the condition number of various finite element matrices involved in the numerical solution of Maxwell's equations , 2004, IEEE Transactions on Antennas and Propagation.

[22]  Jian-Ming Jin,et al.  Combining PML and ABC for the finite-element analysis of scattering problems , 1996 .

[23]  Jin-Fa Lee,et al.  A comparison of anisotropic PML to Berenger's PML and its application to the finite-element method for EM scattering , 1997 .

[24]  J. T. Smith.  Conservative modeling of 3-D electromagnetic fields, Part I: Properties and error analysis , 1996 .

[25]  John L. Volakis,et al.  Preconditioned generalized minimal residual iterative scheme for perfectly matched layer terminated applications , 1999 .

[26]  Zongfu Yu,et al.  Phase front design with metallic pillar arrays , 2010, CLEO/QELS: 2010 Laser Science to Photonic Applications.

[27]  D. A. H. Jacobs,et al.  A Generalization of the Conjugate-Gradient Method to Solve Complex Systems , 1986 .