Perfectly matched layer

Perfectly matched layer

A perfectly matched layer (PML) is an artificial absorbing layer for wave equations, commonly used to truncate computational regions in numerical methods to simulate problems with open boundaries, especially in the FDTD and FEM methods. The key property of a PML that distinguishes it from an ordinary absorbing material is that it is designed so that waves incident upon the PML from a non-PML medium do not reflect at the interface—this property allows the PML to strongly absorb outgoing waves from the interior of a computational region without reflecting them back into the interior.

PML was originally formulated by Berenger in 1994 for use with Maxwell's equations, and since that time there have been several related reformulations of PML for both Maxwell's equations and for other wave equations. Berenger's original formulation is called a split-field PML, because it splits the electromagnetic fields into two unphysical fields in the PML region. A later formulation that has become more popular because of its simplicity and efficiency is called uniaxial PML or UPML (Gedney, 1996), in which the PML is described as an artificial anisotropic absorbing material. Although both Berenger's formulation and UPML were initially derived by manually constructing the conditions under which incident plane waves do not reflect from the PML interface from a homogeneous medium, "both" formulations were later shown to be equivalent to a much more elegant and general approach: stretched-coordinate PML (Chew and Weedon, 1994; Teixeira and Weedon, 1998). In particular, PMLs were shown to correspond to a coordinate transformation in which one (or more) coordinates are mapped to complex numbers; more technically, this is actually an analytic continuation of the wave equation into complex coordinates, replacing propagating (oscillating) waves by exponentially decaying waves. This viewpoint allows PMLs to be derived for inhomogeneous media such as waveguides, as well as for other coordinate systems and wave equations.

Technical description

Specifically, for a PML designed to absorb waves propagating in the "x" direction, the following transformation is included in the wave equation. Wherever an "x" derivative partial/partial x appears in the wave equation, it is replaced by::frac{partial}{partial x} o frac{1}{1 + frac{isigma(x)}{omega frac{partial}{partial x}where ω is the angular frequency and σ is some function of "x". Wherever σ is positive, propagating waves are attenuated because::e^{i(kx - omega t)} o e^{i(kx - omega t) - frac{k}{omega} int^x sigma(x') dx'} ,where we have taken a planewave propagating in the +"x" direction (for k > 0) and applied the transformation (analytic continuation) to complex coordinates: x o x + frac{i}{omega} int^x sigma(x') dx', or equivalently dx o dx (1 + isigma/omega). The same coordinate transformation causes waves to attenuate whenever their "x" dependence is in the form e^{ikx} for some propagation constant "k": this includes planewaves propagating at some angle with the "x" axis and also transverse modes of a waveguide.

The above coordinate transformation can be left as-is in the transformed wave equations, or can be combined with the material description (e.g. the permittivity and permeability in Maxwell's equations) to form a UPML description. Note also that the coefficient σ/ω depends upon frequency—this is so the attenuation rate is proportional to "k"/ω, which is independent of frequency in a homogeneous material (not including material dispersion, e.g. for vacuum) because of the dispersion relation between ω and "k". However, this frequency-dependence means that a time domain implementation of PML, e.g. in the FDTD method, is more complicated than for a frequency-independent absorber, and involves the auxiliary differential equation (ADE) approach (equivalently, "i"/ω appears as an integral or convolution in time domain).

Perfectly matched layers, in their original form, only attenuate propagating waves; purely evanescent waves (exponentially decaying fields) oscillate in the PML but do not decay more quickly. However, the attenuation of evanescent waves can also be accelerated by including a real coordinate stretching in the PML: this corresponds to making σ in the above expression a complex number, where the imaginary part yields a real coordinate stretching that causes evanescent waves to decay more quickly.

One caveat with perfectly matched layers is that they are only reflectionless for the "exact" wave equation. Once the wave equation is discretized for simulation on a computer, some small numerical reflections appear. For this reason, the PML absorption coefficient σ is typically turned on gradually from zero (e.g. quadratically) over a short distance on the scale of the wavelength of the wave.

References

*cite journal | author= J. Berenger | title= A perfectly matched layer for the absorption of electromagnetic waves | journal= Journal of Computational Physics | year= 1994 | volume= 114 | pages= 185–200 | doi= 10.1006/jcph.1994.1159
*cite journal | author= S.D. Gedney | title= An anisotropic perfectly matched layer absorbing media for the truncation of FDTD latices| journal= Antennas and Propagation, IEEE Transactions on | year= 1996 | volume= 44 | pages= 1630–1639 | doi= 10.1109/8.546249
*cite journal | author= W. C. Chew and W. H. Weedon | title= A 3d perfectly matched medium from modified Maxwell's equations with stretched coordinates| journal= Microwave Optical Tech. Letters | year= 1994 | volume= 7 | pages= 590–604
*cite journal | author= F. L. Teixeira W. C. Chew | title= General closed-form PML constitutive tensors to match arbitrary bianisotropic and dispersive linear media| journal= IEEE Microwave and Guided Wave Letters | year= 1998 | volume= 8 | pages= 223–225 | doi= 10.1109/75.678571
* A. Taflove and S. C. Hageness, "Computational Electrodynamics: The Finite-Difference Time-Domain Method" (Norwood, MA: Artech, 2000).
* S. G. Johnson, [http://math.mit.edu/~stevenj/18.369/pml.pdf Notes on Perfectly Matched Layers] , online MIT course notes (Aug. 2007).


Wikimedia Foundation. 2010.

Игры ⚽ Поможем сделать НИР

Look at other dictionaries:

  • Zone absorbante parfaitement adaptée — Une zone absorbante parfaitement adaptée (en anglais Perfectly matched layer, PML) est une zone absorbante artificielle pour les équations d ondes, couramment utilisée pour tronquer les domaines de calcul dans les méthodes numériques de… …   Wikipédia en Français

  • Метод конечных разностей во временной области — (англ. Finite Difference Time Domain, FDTD)  один из наиболее популярных методов численной электродинамики, основанный на дискретизации уравнений Максвелла, записанных в дифференциальной форме. Содержание 1 Описание 2 Алгоритм Йи …   Википедия

  • Finite-difference time-domain method — Finite difference time domain (FDTD) is a popular computational electrodynamics modeling technique. It is considered easy to understand and easy to implement in software. Since it is a time domain method, solutions can cover a wide frequency… …   Wikipedia

  • Charles E. Smith (engineer) — Dr. Charles E. Smith was professor and Chair in the Department of Electrical Engineering at the University of Mississippi. His focus research areas included: application of electromagnetic theory to microwave circuits, application of numerical… …   Wikipedia

  • List of numerical analysis topics — This is a list of numerical analysis topics, by Wikipedia page. Contents 1 General 2 Error 3 Elementary and special functions 4 Numerical linear algebra …   Wikipedia

  • PML — may refer to: * Pakistan Muslim League, faction of the Pakistani political parties with similar names * Partido Movimiento Libertario, libertarian political party in Costa Rica * PartnerML, Vodafone s proprietary replacement for Wireless Markup… …   Wikipedia

  • List of mathematics articles (P) — NOTOC P P = NP problem P adic analysis P adic number P adic order P compact group P group P² irreducible P Laplacian P matrix P rep P value P vector P y method Pacific Journal of Mathematics Package merge algorithm Packed storage matrix Packing… …   Wikipedia

  • UPML — may refer to: * Ukrainian Lyceum of Physics and Mathematics of T. Schevchenko Kiev University * Unaxial Perfectly Matched Layer, numerical truncation methodology …   Wikipedia

  • Electromagnetic field solver — Electromagnetic field solvers (or sometimes just field solvers) are specialized programs that solve (a subset of) Maxwell s equations directly. They form a part of the field of electronic design automation, or EDA, and are commonly used in the… …   Wikipedia

  • PML — Die Abkürzung PML steht für: IATA Code von Port Moller, Arkansas, USA Pakistan Muslim Liga, zwei pakistanische Volksparteien (PML N und PML Q) Palm Markup Language Perfectly Matched Layer, ein Verfahren in der numerischen Elektrodynamik zur… …   Deutsch Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”