Remote Access Mathematics of Computation
Green Open Access

Mathematics of Computation

ISSN 1088-6842(online) ISSN 0025-5718(print)

 
 

 

Analysis of a finite PML approximation for the three dimensional time-harmonic Maxwell and acoustic scattering problems


Authors: James H. Bramble and Joseph E. Pasciak
Journal: Math. Comp. 76 (2007), 597-614
MSC (2000): Primary 78M10, 65F10, 65N30
DOI: https://doi.org/10.1090/S0025-5718-06-01930-2
Published electronically: December 7, 2006
MathSciNet review: 2291829
Full-text PDF Free Access

Abstract | References | Similar Articles | Additional Information

Abstract: We consider the approximation of the frequency domain three-dimensional Maxwell scattering problem using a truncated domain perfectly matched layer (PML). We also treat the time-harmonic PML approximation to the acoustic scattering problem. Following work of Lassas and Somersalo in 1998, a transitional layer based on spherical geometry is defined, which results in a constant coefficient problem outside the transition. A truncated (computational) domain is then defined, which covers the transition region. The truncated domain need only have a minimally smooth outer boundary (e.g., Lipschitz continuous). We consider the truncated PML problem which results when a perfectly conducting boundary condition is imposed on the outer boundary of the truncated domain. The existence and uniqueness of solutions to the truncated PML problem will be shown provided that the truncated domain is sufficiently large, e.g., contains a sphere of radius $ R_t$. We also show exponential (in the parameter $ R_t$) convergence of the truncated PML solution to the solution of the original scattering problem inside the transition layer.

Our results are important in that they are the first to show that the truncated PML problem can be posed on a domain with nonsmooth outer boundary. This allows the use of approximation based on polygonal meshes. In addition, even though the transition coefficients depend on spherical geometry, they can be made arbitrarily smooth and hence the resulting problems are amenable to numerical quadrature. Approximation schemes based on our analysis are the focus of future research.


References [Enhancements On Off] (What's this?)

  • 1. C. Amrouche, C. Bernardi, M. Dauge, and V. Girault.
    Vector potentials in three-dimensional non-smooth domains.
    Math. Methods Appl. Sci., 21(9):823-864, 1998. MR 1626990 (99e:35037)
  • 2. J.-P. Bérenger.
    A perfectly matched layer for the absorption of electromagnetic waves.
    J. Comput. Phys., 114(2):185-200, 1994. MR 1294924 (95e:78002)
  • 3. J.-P. Bérenger.
    Three-dimensional perfectly matched layer for the absorption of electromagnetic waves.
    J. Comput. Phys., 127(2):363-379, 1996. MR 1412240 (97h:78001)
  • 4. W. Chew and W. Weedon.
    A 3d perfectly matched medium for modified Maxwell's equations with streched coordinates.
    Microwave Opt. Techno. Lett., 13(7):599-604, 1994.
  • 5. F. Collino and P. Monk.
    The perfectly matched layer in curvilinear coordinates.
    SIAM J. Sci. Comp., 19(6):2061-2090, 1998. MR 1638033 (99e:78029)
  • 6. V. Girault and P. Raviart.
    Finite Element Approximation of the Navier-Stokes Equations.
    Lecture Notes in Math. 749, Springer-Verlag, New York, 1981. MR 0548867 (83b:65122)
  • 7. C. I. Goldstein.
    The finite element method with nonuniform mesh sizes applied to the exterior Helmholtz problem.
    Numer. Math., 38(1):61-82, 1981/82. MR 0634753 (82k:65062)
  • 8. M. Grote and J. Keller.
    Nonreflecting boundary conditions for Maxwell's equations.
    J. Comput. Phys., 139:327-342, 1998. MR 1614098 (99a:78026)
  • 9. N. Kantartzis, P. Petropoulis, and T. Tsiboukis.
    A comparison of the Grote-Keller exact ABC and the well posed PML for Maxwell's equations in spherical coordinates.
    IEEE Trans. on Magnetics, 35:1418-1422, 1999.
  • 10. M. Lassas and E. Somersalo.
    On the existence and convergence of the solution of PML equations.
    Computing, 60(3):229-241, 1998. MR 1621305 (99a:65133)
  • 11. M. Lassas and E. Somersalo.
    Analysis of the PML equations in general convex geometry.
    Proc. Roy. Soc. Edinburgh Sect. A, 131(5):1183-1207, 2001. MR 1862449 (2002k:35020)
  • 12. P. Monk.
    Finite Element Methods for Maxwell's Equations.
    Numerical Mathematics and Scientific Computation. Oxford University Press, Oxford, UK, 2003. MR 2059447 (2005d:65003)
  • 13. J. Peetre.
    Espaces d'interpolation et théorèm de Soboleff.
    Ann. Inst. Fourier, 16:279-317, 1966. MR 0221282 (36:4334)
  • 14. P. G. Petropoulos.
    Reflectionless sponge layers as absorbing boundary conditions for the numerical solution of Maxwell equations in rectangular, cylindrical, and spherical coordinates.
    SIAM J. Appl. Math., 60(3):1037-1058 (electronic), 2000. MR 1750090 (2001c:35230)
  • 15. L. Tartar.
    Topics in Nonlinear Analysis.
    Math. d'Orsay, Univ. Paris-Sud, 1978. MR 0532371 (81b:35001)

Similar Articles

Retrieve articles in Mathematics of Computation with MSC (2000): 78M10, 65F10, 65N30

Retrieve articles in all journals with MSC (2000): 78M10, 65F10, 65N30


Additional Information

James H. Bramble
Affiliation: Department of Mathematics, Texas A&M University, College Station, Texas 77843-3368
Email: bramble@math.tamu.edu

Joseph E. Pasciak
Affiliation: Department of Mathematics, Texas A&M University, College Station, Texas 77843-3368
Email: pasciak@math.tamu.edu

DOI: https://doi.org/10.1090/S0025-5718-06-01930-2
Keywords: Maxwell's equations, Helmholtz equation, time-harmonic acoustic and electromagnetic scattering, div-curl systems, perfectly matched layer, PML
Received by editor(s): August 9, 2005
Received by editor(s) in revised form: April 20, 2006
Published electronically: December 7, 2006
Additional Notes: This work was supported in part by the National Science Foundation through grant No. 0311902.
Article copyright: © Copyright 2006 American Mathematical Society
The copyright for this article reverts to public domain 28 years after publication.

American Mathematical Society