Fast integral equation solver for Maxwell’s equations in layered media with FMM for Bessel functions. (English) Zbl 1306.78010

Summary: The paper presents a new fast integral equation solver for Maxwell’s equations in 3-D layered media. First, the spectral domain dyadic Green’s function is derived, and the 0-th and the 1-st order Hankel transforms or Sommerfeld-type integrals are used to recover all components of the dyadic Green’s function in real space. The Hankel transforms are performed with the adaptive generalized Gaussian quadrature points and window functions to minimize the computational cost. Subsequently, a fast integral equation solver with \(O(N_z^2N_xN_y\log(N_xN_y))\) in layered media is developed by rewriting the layered media integral operator in terms of Hankel transforms and using the new fast multipole method for the \(n\)-th order Bessel function in 2-D. Computational cost and parallel efficiency of the new algorithm are presented.


78M16 Multipole methods applied to problems in optics and electromagnetic theory
65F10 Iterative numerical methods for linear systems
65R20 Numerical methods for integral equations
65N22 Numerical solution of discretized equations for boundary value problems involving PDEs
35J05 Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation
65T50 Numerical methods for discrete and fast Fourier transforms
65N80 Fundamental solutions, Green’s function methods, etc. for boundary value problems involving PDEs
78A48 Composite media; random media in optics and electromagnetic theory
65D30 Numerical integration


Full Text: DOI


[1] Abramowitz M, Stegun I A. Handbook of Mathematical Functions. New York: Dover, 1970 · Zbl 0171.38503
[2] Atwater, H A; Polman, A, Plasmonics for improved photovoltaic devices, Nature Materials, 9, 205-213, (2010)
[3] Barnes, J; Hut, P, A hierarchical \(O\)(\(N\) log\(N\)) force-calculation algorithm, Nature, 324, 446-449, (1986)
[4] Bulldozer for Servers: Testing AMD’s “Interlagos” Opteron 6200 Series. http://www.anandtech.com/show/5058/amdsopteron-interlagos-6200
[5] Cai W. Computational Methods for Electromagnetic Phenomena: Electrostatics in Solvation, Scattering, and Electron Transport. Cambridge: Cambridge University Press, 2013 · Zbl 06136469
[6] Cai, W; Wang, J, Adapative multiresolution collocation methods for initial boundary value problems of nonlinear pdes, SIAM J Numer Anal, 33, 937-970, (1996) · Zbl 0856.65115
[7] Cai, W; Yu, T, Fast calculation of dyadic green’s functions for electromagnetic scattering in a multilayered medium, J Comp Phys, 165, 1-21, (2000) · Zbl 1040.78015
[8] Cai, W; Zhang, W, An adapative spline wavelet ADI (SW-ADI) method for two-dimensional reaction-diffusion equations, J Comp Phys, 139, 92-126, (1998) · Zbl 0905.65103
[9] Chakraborty, U; Chatterjee, S; Chowdhury, S K; etal., A compact microstrip patch antenna for wireless communication, PIER C, 18, 211-220, (2011)
[10] Cho, M H; Cai, W, A parallel fast algorithm for computing the Helmholtz integral operator in 3-D layered media, J Comp Phys, 231, 5910-5925, (2012) · Zbl 1277.35143
[11] Erteza, A; Park, B K, Nonuniqueness of resolution of Hertz vector in presence of boundary, and the horizontal problem, IEEE Trans Ant Prop, 17, 376-378, (1969)
[12] Felsen L B, Marcuvitz N. Radiation and Scattering of Waves. Englewood Cliffs, NJ: Prentice-Hall, 1973 · Zbl 0050.41802
[13] Greengard L F. The Rapid Evaluation of Potential Fields in Particle Systems. Cambridge: MIT Press, 1988 · Zbl 1001.31500
[14] Harrington R F. Field Computation by Moment Method. New York: Macmillian, 1968
[15] Kupec, J; Witzigmann, B, Computational electromagnetics for nanowire solar cells, J Comput Electron, 11, 153-165, (2012)
[16] Lam, V D; Kim, J B; Tung, N T; etal., Dependence of the distance between cut-wire-pair layers on resonance frequencies, Optics Express, 16, 5934-5941, (2008)
[17] Lam, V D; Tung, N T; Cho, M H; etal., Influence of lattice parameters on the resonance frequencies of a cut-wire-pair medium, J Appl Phys, 105, 113102, (2009)
[18] Li, P; Johnston, H; Krasny, R, A Cartesian treecode for screened Coulomb interactions, J Comp Phys, 228, 3858-3868, (2009) · Zbl 1165.78304
[19] Ma, J; Rokhlin, V; Wandzura, S, Generalized Gaussian quadrature rules for systems of arbitrary functions, SIAM J Numer Anal, 33, 971-996, (1996) · Zbl 0858.65015
[20] Michalski, K; Mosig, J, Multilayered media green’s functions in integral equation formulations, IEEE Trans Antennas Propagation, 45, 508-519, (1997)
[21] Michalski, K A; Zhang, D, Electromagnetic scattering and radiation by surfaces of arbitrary shape in layered media, I: theory, IEEE Trans Antennas Propagation, 38, 335-344, (1990)
[22] Mosig, J R; Itoh, T (ed.), Integral equation technique, 133-213, (1989), New York
[23] OpenMP. http://openmp.org/wp/
[24] Pthreads. https://computing.llnl.gov/tutorials/pthreads/ · Zbl 0905.65103
[25] Rokhlin, V, Rapid solution of integral equations of scattering theory in two dimensions, J Comp Phys, 86, 414-439, (1990) · Zbl 0686.65079
[26] Saad Y. Iterative Methods for Sparse Linear Systems. Boston, MA: PWS Publishing Company, 1996 · Zbl 1031.65047
[27] Sommerfeld A. Partial Differential Equations in Physics. New York: Academic Press, 1949 · Zbl 0034.35702
[28] Stratton J A. Electromagnetic Theory. Hoboken, NJ: John Wiley & Sons, 2007
[29] Vorst, H A, BI-CGSTAB: A fast and smoothly converging variant of bi-CG for the solution of nonsymmetric linear systems, SIAM J Sci Stat Comput, 13, 631-644, (1992) · Zbl 0761.65023
[30] Wait J R. Electromagnetic Waves in Stratified Media. Elmsford, NY: Pergamon, 1970 · Zbl 0122.45005
[31] Wang, S; Koleilat, G I; Tang, J; etal., Tandem colloidal quantum dot solar cells employing a graded recombination layer, Nature Photonics, 5, 480-484, (2011)
[32] Xia, L; Wang, C F; Li, L W; etal., Resonant behaviors of microstrip antenna in multilayered media: an efficient full-wave analysis, Prog Electromagnetics Res, 31, 55-67, (2001)
[33] Yu, T; Cai, W, High-order window functions and fast algorithms for calculating dyadic electromagnetic green’s functions in multilayered media, Radio Sci, 36, 559-569, (2001)
[34] Yu, T; Cai, W, FIFA-fast interpolations and filtering algorithm for calculating dyadic green’s function in the electromagnetic scattering of multi-layered structures, Commun Comput Phys, 1, 229-260, (2006) · Zbl 1114.78361
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.