Plane-Wave Analysis and Comparison of Split-Field, Biaxial, and Uniaxial PML Methods as ABCs for Pseudospectral Electromagnetic Wave Simulations in Curvilinear Coordinates

In this paper, we discuss and compare split-field, biaxial, and uniaxial perfectly matched layer (PML) methods for absorbing outgoing vector waves in cylindrical and spherical coordinates. We first extend Berenger's split-field formulation into spherical and cylindrical coordinates in such a way that it maintains all the desirable properties it exhibits in rectangular coordinates. Then we discuss the biaxial and the uniaxial medium PML methods in Cartesian coordinates and extend them to spherical and cylindrical coordinates. Properties of plane-wave solutions of the PML methods are analyzed. In particular, the decay and boundness properties of the solutions are considered in order to provide further insight into the different formulations presented herein. Moreover, we propose a set of symmetric hyperbolic equations for both the biaxial and the uniaxial PML methods in the time-domain, which is fine-tuned in numerical experiments and very suitable for time-domain problems. All three types of spherical and cylindrical PML methods are applied in simulations of plane wave scattering as well as radiating dipole problems. We use a multidomain pseudospectral (Chebyshev) numerical scheme, and the effectiveness of the PML methods is demonstrated through the accurate numerical results obtained. The order of outer-boundary reflection is as low as 0.1% of the exact solution.

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

[2]  J. Bérenger.  Three-Dimensional Perfectly Matched Layer for the Absorption of Electromagnetic Waves , 1996 .

[3]  Raj Mittra,et al.  The use of the FFT for the efficient solution of the problem of electromagnetic scattering by a body of revolution , 1990 .

[4]  Thomas L. Geers,et al.  Evaluation of the Perfectly Matched Layer for Computational Acoustics , 1998 .

[5]  Marcus J. Grote,et al.  Nonreflecting Boundary Conditions for Maxwell's Equations , 1998 .

[6]  David Gottlieb,et al.  A Mathematical Analysis of the PML Method , 1997 .

[7]  Chen Wu,et al.  Application of PML superabsorbing boundary condition to non-orthogonal FDTD method , 1994 .

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

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

[10]  Weng Cho Chew,et al.  PML-FDTD in cylindrical and spherical grids , 1997 .

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

[12]  Raj Mittra,et al.  Frequency dependence of the constitutive parameters of causal perfectly matched anisotropic absorbers , 1996 .

[13]  Raj Mittra,et al.  Investigation of nonplanar perfectly matched absorbers for finite-element mesh truncation , 1997 .

[14]  Jan S. Hesthaven,et al.  Spectral Simulations of Electromagnetic Wave Scattering , 1997 .

[15]  W. Chew,et al.  Systematic derivation of anisotropic PML absorbing media in cylindrical and spherical coordinates , 1997 .

[16]  Andreas C. Cangellaris,et al.  GT-PML: generalized theory of perfectly matched layers and its application to the reflectionless truncation of finite-difference time-domain grids , 1996, IMS 1996.