×

zbMATH — the first resource for mathematics

A spectral element method with transparent boundary condition for periodic layered media scattering. (English) Zbl 1373.78442
Summary: We present a high-order spectral element method for solving layered media scattering problems featuring an operator that can be used to transparently enforce the far-field boundary condition. The incorporation of this Dirichlet-to-Neumann (DtN) map into the spectral element framework is a novel aspect of this work, and the resulting method can accommodate plane-wave radiation of arbitrary angle of incidence. In order to achieve this, the governing Helmholtz equations subject to quasi-periodic boundary conditions are rewritten in terms of periodic unknowns. We construct a spectral element operator to approximate the DtN map, thus ensuring nonreflecting outgoing waves on the artificial boundaries introduced to truncate the computational domain. We present an explicit formula that accurately computes the Fourier coefficients of the solution in the spectral element discretization space projected onto the boundary which is required by the DtN map. Our solutions are represented by the tensor product basis of one-dimensional Legendre-Lagrange interpolation polynomials based on the Gauss-Lobatto-Legendre grids. We study the scattered field in singly and doubly layered media with smooth and nonsmooth interfaces. We consider rectangular, triangular, and sawtooth interfaces that are accurately represented by the body-fitted quadrilateral elements. We use GMRES iteration to solve the resulting linear system, and we validate our results by demonstrating spectral convergence in comparison with exact solutions and the results of an alternative computational method.

MSC:
78M22 Spectral, collocation and related methods applied to problems in optics and electromagnetic theory
65N35 Spectral, collocation and related methods for boundary value problems involving PDEs
78A45 Diffraction, scattering
78M30 Variational methods applied to problems in optics and electromagnetic theory
65F10 Iterative numerical methods for linear systems
78A48 Composite media; random media in optics and electromagnetic theory
35Q60 PDEs in connection with optics and electromagnetic theory
35J05 Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation
PDF BibTeX XML Cite
Full Text: DOI
References:
[1] Banerjee, P.K., Butterfield, R.: Boundary element methods in engineering science. McGraw-Hill, London (1981) · Zbl 0499.73070
[2] Bonnet, M.: Boundary integral equation methods for solids and fluids. Wiley, New Jersey (1999)
[3] Greengard, L; Rokhlin, V, A fast algorithm for particle simulations, J. Comput. Phys., 73, 325-348, (1987) · Zbl 0629.65005
[4] Bruno, OP; Reitich, F, Numerical solution of diffraction problems: a method of variation of boundaries, J. Opt. Soc. Am. A, 10, 1168-1175, (1993)
[5] Bruno, OP; Reitich, F, Numerical solution of diffraction problems: a method of variation of boundaries. II. finitely conducting gratings, Padé approximants, and singularities, J. Opt. Soc. Am. A, 10, 2307-2316, (1993)
[6] Bruno, OP; Reitich, F, Numerical solution of diffraction problems: a method of variation of boundaries. III. doubly periodic gratings, J. Opt. Soc. Am. A, 10, 2551-2562, (1993)
[7] Milder, D.M.: An improved formalism for rough-surface scattering of acoustic and electromagnetic waves. In: Proceedings of SPIE—the international society for optical engineering (San Diego, 1991), 1558, pp. 213-221. Bellingham, WA (1991)
[8] Milder, DM, An improved formalism for wave scattering from rough surfaces, J. Acoust. Soc. Am., 89, 529-541, (1991)
[9] Milder, DM, Role of the admittance operator in rough-surface scattering, J. Acoust. Soc. Am., 100, 759-768, (1996)
[10] Milder, DM, An improved formalism for electromagnetic scattering from a perfectly conducting rough surface, Radio Sci., 31, 1369-1376, (1996)
[11] Milder, D.M., Sharp, H.: Efficient computation of rough surface scattering. In: Mathematical and numerical aspects of wave propagation phenomena (Strasbourg, 1991), pp. 314-322, SIAM, Philadelphia, PA (1991) · Zbl 0853.73066
[12] Milder, DM; Sharp, HT, An improved formalism for rough surface scattering. II: numerical trials in three dimensions, J. Acoust. Soc. Am., 91, 2620-2626, (1992)
[13] Nicholls, DP; Shen, J, A rigorous numerical analysis of the transformed field expansion method, SIAM J. Numer. Anal., 47, 2708-2734, (2009) · Zbl 1410.65465
[14] Nicholls, DP; Shen, J, A stable, high-order method for two-dimensional bounded-obstacle scattering, SIAM J. Sci. Comput., 28, 1398-1419, (2006) · Zbl 1130.78007
[15] Fang, Q; Nicholls, DP; Shen, J, A stable, high-order method for three-dimensional bounded-obstacle scattering, J. Comput. Phys., 224, 1145-1169, (2007) · Zbl 1123.76049
[16] Petit, R.: Electromagnetic theory of gratings. Springer, Berlin (1980)
[17] Nicholls, DP; Reitich, F, A new approach toanalyticity of Dirichlet-Neumann operators, Proc. R. Soc. Edinb. Sect. A. Math., 131, 1411-1433, (2001) · Zbl 1016.35030
[18] Nicholls, DP; Reitich, F, Stability of high-order perturbative methods for the computation of Dirichlet-Neumann operators, J. Comput. Phys., 170, 276-298, (2001) · Zbl 0983.65115
[19] Nicholls, DP; Reitich, F, Analytic continuation of Dirichlet-Neumann operators, Numer. Math., 94, 107-146, (2003) · Zbl 1030.65109
[20] He, Y; Nicholls, DP; Shen, J, An efficient and stable spectral method for electromagnetic scattering from a layered periodic structure, J. Comput. Phys., 231, 3007-3022, (2012) · Zbl 1243.78025
[21] Nicholls, DP, A method of field expansions for vector electromagnetic scattering by layered periodic crossed gratings, J. Opt. Soc. Am. A, 32, 701-709, (2015)
[22] Deville, M.O., Fischer, P.F., Mund, E.H.: High-order methods for incompressible fluid flow. Cambridge University Press, Cambridge (2002) · Zbl 1007.76001
[23] Han, HD; Wu, XN, Approximation of infinite boundary condition and its application to finite element methods, J. Comput. Math., 3, 179-192, (1985) · Zbl 0579.65111
[24] Keller, J; Givoli, D, Exact nonreflecting boundary conditions, J. Comput. Phys., 82, 172-192, (1989) · Zbl 0671.65094
[25] Givoli, D, Nonreflecting boundary conditions, J. Comput. Phys., 94, 1-29, (1991) · Zbl 0731.65109
[26] Givoli, D; Keller, J, Special finite elements for use with high-order boundary conditions, Comput. Methods Appl. Mech. Eng., 119, 199-213, (1994) · Zbl 0853.73066
[27] Givoli, D.: Numerical methods for problems in infinite domains. Elsevier Scientific Publishing Co., Amsterdam (1992) · Zbl 0788.76001
[28] Grote, M; Keller, J, On nonreflecting boundary conditions, J. Comput. Phys., 122, 231-243, (1995) · Zbl 0841.65099
[29] Givoli, D, Recent advances in the dtn FE method, Arch. Comput. Methods Eng., 6, 71-116, (1999)
[30] Nicholls, DP; Nigam, N, Exact non-reflecting boundary conditions on general domains, J. Comput. Phys., 194, 278-303, (2004) · Zbl 1049.65132
[31] Nicholls, DP; Nigam, N, Error analysis of a coupled finite element/dtn map algorithm on general domains, Numer. Math., 105, 267-298, (2006) · Zbl 1162.65402
[32] Binford, TL; Nicholls, DP; Nigam, N; Warburton, T, Exact non-reflecting boundary conditions on general domains and hp-finite elements, J. Sci. Comput., 39, 265-292, (2009) · Zbl 1203.76135
[33] Strutt, JW; Rayleigh, L, On the manufacture and theory of diffraction gratings, Philos. Mag., 47, 193-205, (1874)
[34] Bao, G, Finite element approximation of time harmonic waves in periodic structures, SIAM J. Numer. Anal., 32, 1155-1169, (1995) · Zbl 0853.65134
[35] Barnett, A; Greengard, L, A new integral representation for quasi-periodic fields and its application to two-dimensional band structure calculations, J. Comput. Phys., 229, 6898-6914, (2010) · Zbl 1197.78025
[36] Saad, Y; Schultz, MH, A generalized minimal residual algorithm for solving nonsymmetric linear systems, SIAM J. Sci. Stat. Comput., 7, 856-869, (1986) · Zbl 0599.65018
[37] Nicholls, DP, Efficient enforcement of far-field boundary conditions in the transformed field expansions method, J. Comput. Phys., 230, 8290-8303, (2011) · Zbl 1408.76405
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.