Fast transform based preconditioners for 2D finite-difference frequency-domain - Waveguides and periodic structures

The fields scattered by dielectric objects placed inside parallel-plate waveguides and periodic structures in two dimensions may efficiently be computed via a finite-difference frequency-domain (FDFD) method. This involves large, sparse linear systems of equations that may be solved using preconditioned Krylov subspace methods. Our preconditioners involve fast discrete trigonometric transforms and are based on a physical approximation. Simulations show significant gain in terms of computation time and iteration count in comparison with results obtained with preconditioners based on incomplete LU (ILU) factorization. Moreover, with the new preconditioners, the required number of iterations is independent of the grid size.

[1]  Elisabeth Larsson,et al.  A Domain Decomposition Method for the Helmholtz Equation in a Multilayer Domain , 1999, SIAM J. Sci. Comput..

[2]  D. R. Fokkema,et al.  BiCGstab(ell) for Linear Equations involving Unsymmetric Matrices with Complex Spectrum , 1993 .

[3]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[4]  Raymond H. Chan,et al.  Conjugate Gradient Methods for Toeplitz Systems , 1996, SIAM Rev..

[5]  Dianne P. O'Leary,et al.  Efficient Iterative Solution of the Three-Dimensional Helmholtz Equation , 1998 .

[6]  Samuel P. Marin,et al.  Variational methods for underwater acoustic problems , 1978 .

[7]  Thomas Huckle,et al.  Fast transforms for tridiagonal linear equations , 1994 .

[8]  Cornelis Vuik,et al.  A Novel Multigrid Based Preconditioner For Heterogeneous Helmholtz Problems , 2005, SIAM J. Sci. Comput..

[9]  Dianne P. O'Leary,et al.  A Multigrid Method Enhanced by Krylov Subspace Iteration for Discrete Helmholtz Equations , 2001, SIAM J. Sci. Comput..

[10]  Elisabeth Larsson,et al.  Iterative Solution of the Helmholtz Equation by a Second-Order Method , 1999, SIAM J. Matrix Anal. Appl..

[11]  Dianne P. O'Leary,et al.  Eigenanalysis of some preconditioned Helmholtz problems , 1999, Numerische Mathematik.

[12]  Elisabeth Larsson,et al.  Parallel Solution of the Helmholtz Equation in a Multilayer Domain , 2003 .

[13]  Alexander Graham,et al.  Kronecker Products and Matrix Calculus: With Applications , 1981 .