Couche absorbante parfaitement adaptée

Schéma FDTD pour un problème de diffusion de la lumière. Les bordures rayées correspondent à des couches parfaitement adaptées, qui sont utilisées pour simuler des frontières ouvertes en absorbant les ondes sortantes.

Une couche absorbante parfaitement adaptée (en anglais Perfectly matched layer, PML) est une couche absorbante artificielle pour les équations d'ondes, couramment utilisée pour tronquer les domaines de calcul dans les méthodes numériques de simulation de problèmes à frontières ouvertes, particulièrement dans les méthodes FDTD et FEM.

La propriété essentielle d'une PML qui la distingue d'un matériau absorbant ordinaire est le fait qu'elle est conçue de telle sorte que les ondes incidentes l'atteignant depuis un matériau non PML ne se réfléchissent pas à l'interface. Cette propriété permet aux PML d'absorber fortement toutes les ondes sortant d'un domaine de calcul sans les renvoyer dans ce domaine.

Histoire

Les couches absorbantes parfaitement adaptées ont été initialement décrites par Jean-Pierre Bérenger en 1994 pour une utilisation dans les équations de Maxwell. Depuis, elles ont été reformulées plusieurs fois pour ces mêmes équations et pour d'autres équations d'ondes. La formulation initiale de Bérenger est appelée split-field PML, car elle divise le champ électromagnétique en deux champs non physiques dans la zone PML.

Une autre formulation qui est devenue plus connue à cause de sa simplicité et de son efficacité est l'uniaxial PML ou UPML (Gedney, 1996), dans laquelle la PML est décrite comme un matériau anisotrope absorbant artificiel. Bien que les formulations de Bérenger et les UPML aient été initialement dérivées par une construction manuelle des conditions d'incidence dans lesquelles les ondes planes ne se réfléchissent pas à l'interface de la PML avec un matériau homogène, les deux formulations se sont révélées plus tard équivalentes à une approche plus élégante et plus générale : la stretched-coordinate PML (Chew and Weedon, 1994; Teixeira and Chew, 1998). En particulier il a été montré que les PML correspondent à une transformation dans laquelle une ou plusieurs coordonnées sont attachées à des nombres complexes ; plus techniquement, il s'agit en réalité d'une prolongation analytique de l'équation d'onde dans le domaine complexe, où les ondes propagatives (oscillantes) sont remplacées par des ondes dont l'amplitude décroît exponentiellement. Ce point de vue permet aux PML d'être adaptées aux matériaux inhomogènes comme les guides d'ondes, ainsi qu'à d'autre systèmes de coordonnées et d'équations d'ondes.

Forme mathématique du problème

Le formalisme initial développé par Bérenger consiste en un passage dans le domaine complexe des coordonnées dans le domaine absorbant, sous la forme :

avec ω représentant la fréquence angulaire et σ est une fonction de x uniquement. Si σ est positive, les ondes se propageant sont atténuées le long de la direction x car, par le changement de coordonnées, on a :

Ainsi, les termes se propageant selon x sous la forme pour une constante de propagation k décroissent.

Notes et références

Voir aussi

Bibliographie

  • (en) J. Berenger, « A perfectly matched layer for the absorption of electromagnetic waves », Journal of Computational Physics, vol. 114, no 2,‎ , p. 185–200 (DOI 10.1006/jcph.1994.1159)
  • (en) S.D. Gedney, « An anisotropic perfectly matched layer absorbing media for the truncation of FDTD latices », Antennas and Propagation, IEEE Transactions on, vol. 44, no 12,‎ , p. 1630–1639 (DOI 10.1109/8.546249)
  • (en) W. C. Chew and W. H. Weedon, « A 3d perfectly matched medium from modified Maxwell's equations with stretched coordinates », Microwave Optical Tech. Letters, vol. 7, no 13,‎ , p. 599–604 (DOI 10.1002/mop.4650071304)
  • (en) F. L. Teixeira W. C. Chew, « General closed-form PML constitutive tensors to match arbitrary bianisotropic and dispersive linear media », IEEE Microwave and Guided Wave Letters, vol. 8, no 6,‎ , p. 223–225 (DOI 10.1109/75.678571)
  • (en) Allen Taflove et Susan C. Hagness, Computational Electrodynamics : The Finite-Difference Time-Domain Method, 3rd ed., Boston, Artech House Publishers, , 3e éd., 1006 p. (ISBN 978-1-58053-832-9, LCCN 2005045273)
  • S. G. Johnson, Notes on Perfectly Matched Layers, online MIT course notes (Aug. 2007).

Lien externe