Efficient implementation of multi-pole UPML using trapezoidal approximation for general media


Journal Article

© 2014 Elsevier B.V. Based on the uniaxial anisotropic perfectly matched layer (UPML) with multi-poles, unsplit-field implementation of the higher-order PML using the trapezoidal approximation (TA) method is proposed to terminate the finite-difference time-domain (FDTD) computational domains. From the point of view of the Courant-Friedrichs-Levy (CFL) condition, to the best of our knowledge, time step based on the TA only needs to meet CFL condition, whereas time step based on the matched Z-transform (MZT) method has to make it smaller for retaining stability and desirable accuracy. Moreover, these formulations are completely independent of the material properties of the FDTD domains and hence can be applied to truncate arbitrary media without any modification because of the D-B constitutive relations used in Maxwell's equations. Four numerical examples have been carried out in three dimensional (3D) FDTD computational domains to validate these formulations. It is shown that the proposed UPML formulations with two poles are effective in terms of attenuating both the low-frequency propagating waves and evanescent waves and reducing late-time reflections, and also can produce results as accurate as the published MZT-UPML but with fewer number of steps and less CPU time.

Full Text

Duke Authors

Cited Authors

  • Feng, N; Liu, QH

Published Date

  • December 1, 2014

Published In

Volume / Issue

  • 111 /

Start / End Page

  • 59 - 65

International Standard Serial Number (ISSN)

  • 0926-9851

Digital Object Identifier (DOI)

  • 10.1016/j.jappgeo.2014.09.020

Citation Source

  • Scopus