arrow
Volume 2, Issue 4
Comparison of Finite Difference and Mixed Finite Element Methods for Perfectly Matched Layer Models

V. A. Bokil & M. W. Buksas

Commun. Comput. Phys., 2 (2007), pp. 806-826.

Published online: 2007-02

Export citation
  • Abstract

We consider the anisotropic uniaxial formulation of the perfectly matched layer (UPML) model for Maxwell's equations in the time domain. We present and analyze a mixed finite element method for the discretization of the UPML in the time domain to simulate wave propagation on unbounded domains in two dimensions. On rectangles the spatial discretization uses bilinear finite elements for the electric field and the lowest order Raviart-Thomas divergence conforming elements for the magnetic field. We use a centered finite difference method for the time discretization. We compare the finite element technique presented to the finite difference time domain method (FDTD) via a numerical reflection coefficient analysis. We derive the numerical reflection coefficient for the case of a semi-infinite PML layer to show consistency between the numerical and continuous models, and in the case of a finite PML to study the effects of terminating the absorbing layer. Finally, we demonstrate the effectiveness of the mixed finite element scheme for the UPML by a numerical example and provide comparisons with the split field PML discretized by the FDTD method. In conclusion, we observe that the mixed finite element scheme for the UPML model has absorbing properties that are comparable to the FDTD method.

  • Keywords

  • AMS Subject Headings

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address
  • BibTex
  • RIS
  • TXT
@Article{CiCP-2-806, author = {V. A. Bokil and M. W. Buksas}, title = {Comparison of Finite Difference and Mixed Finite Element Methods for Perfectly Matched Layer Models}, journal = {Communications in Computational Physics}, year = {2007}, volume = {2}, number = {4}, pages = {806--826}, abstract = {

We consider the anisotropic uniaxial formulation of the perfectly matched layer (UPML) model for Maxwell's equations in the time domain. We present and analyze a mixed finite element method for the discretization of the UPML in the time domain to simulate wave propagation on unbounded domains in two dimensions. On rectangles the spatial discretization uses bilinear finite elements for the electric field and the lowest order Raviart-Thomas divergence conforming elements for the magnetic field. We use a centered finite difference method for the time discretization. We compare the finite element technique presented to the finite difference time domain method (FDTD) via a numerical reflection coefficient analysis. We derive the numerical reflection coefficient for the case of a semi-infinite PML layer to show consistency between the numerical and continuous models, and in the case of a finite PML to study the effects of terminating the absorbing layer. Finally, we demonstrate the effectiveness of the mixed finite element scheme for the UPML by a numerical example and provide comparisons with the split field PML discretized by the FDTD method. In conclusion, we observe that the mixed finite element scheme for the UPML model has absorbing properties that are comparable to the FDTD method.

}, issn = {1991-7120}, doi = {https://doi.org/}, url = {http://global-sci.org/intro/article_detail/cicp/7929.html} }
TY - JOUR T1 - Comparison of Finite Difference and Mixed Finite Element Methods for Perfectly Matched Layer Models AU - V. A. Bokil & M. W. Buksas JO - Communications in Computational Physics VL - 4 SP - 806 EP - 826 PY - 2007 DA - 2007/02 SN - 2 DO - http://doi.org/ UR - https://global-sci.org/intro/article_detail/cicp/7929.html KW - AB -

We consider the anisotropic uniaxial formulation of the perfectly matched layer (UPML) model for Maxwell's equations in the time domain. We present and analyze a mixed finite element method for the discretization of the UPML in the time domain to simulate wave propagation on unbounded domains in two dimensions. On rectangles the spatial discretization uses bilinear finite elements for the electric field and the lowest order Raviart-Thomas divergence conforming elements for the magnetic field. We use a centered finite difference method for the time discretization. We compare the finite element technique presented to the finite difference time domain method (FDTD) via a numerical reflection coefficient analysis. We derive the numerical reflection coefficient for the case of a semi-infinite PML layer to show consistency between the numerical and continuous models, and in the case of a finite PML to study the effects of terminating the absorbing layer. Finally, we demonstrate the effectiveness of the mixed finite element scheme for the UPML by a numerical example and provide comparisons with the split field PML discretized by the FDTD method. In conclusion, we observe that the mixed finite element scheme for the UPML model has absorbing properties that are comparable to the FDTD method.

V. A. Bokil and M. W. Buksas. (2007). Comparison of Finite Difference and Mixed Finite Element Methods for Perfectly Matched Layer Models. Communications in Computational Physics. 2 (4). 806-826. doi:
Copy to clipboard
The citation has been copied to your clipboard