[论文解读] Notes on Perfectly Matched Layers (PMLs)

Steven G. Johnson|arXiv (Cornell University)|Aug 4, 2021

Electromagnetic Simulation and Numerical Methods参考文献 20被引用 163

一句话总结

本 note 提供一个简洁、连贯的关于通过复坐标拉伸实现的完美匹配层（PML）吸收边界的视图，包含局限性与实际实现。

ABSTRACT

This note is intended as a brief introduction to the theory and practice of perfectly matched layer (PML) absorbing boundaries for wave equations, originally developed for MIT courses 18.369 and 18.336. It focuses on the complex stretched-coordinate viewpoint, and also discusses the limitations of PML.

研究动机与目标

• 动机：在不引入显著反射的前提下，对波动方程的计算域进行截断。
• 给出一个基于复坐标拉伸的统一、通用的PML框架。
• 展示PML如何在频域和时域中推导与实现。
• 讨论PML在离散化、角度 dependence、以及非均匀介质下的失效等局限性。

提出的方法

• 从解析延拓到复坐标、再回到实坐标的坐标变换推导PML。
• 将PML变换表征为在标量波方程中将 ∂/∂x 替换为 (1 + i σx/ω)^{-1} ∂/∂x。
• 解释当以变换后的本构关系表示时，PML如何表现为各向异性的吸收材料。
• 给出现时域实现，使用辅助微分方程 (ADEs) 实现 1/ω 衰减。
• 提供一维和二维标量波的示例来说明PML的构造与行为。

实验结果

研究问题

• RQ1 如何使用PML实现对任意入射角的波的无反射吸收？
• RQ2 复坐标拉伸与麦克斯韦方程组中变换后的材料特性之间有何关系？
• RQ3 在离散化、非均匀或衰减场景中，PML的局限性与失效模式有哪些？
• RQ4 如何在时域仿真中实现PML以维持频率独立的衰减？
• RQ5 在何种条件下PML会失效，以及在这些情况下存在哪些替代方案或变通方法？

主要发现

• 在理想连续情形下，PML将振荡解在感兴趣区域之外转化为指数衰减解且不产生反射。
• PML变换引入一种有效的各向异性、复值材料，通过坐标拉伸在各方向吸收波。
• 在时域中，ADEs可实现理想PML固有的1/ω衰减，使实际仿真成为可能。
• 离散化引入反射和与入射角相关的吸收；减慢开启速度和充足厚度可减少这些影响。
• 对于x方向不均匀介质、解析延拓假设失效的情况；或某些左手/右手材料或后向波结构，PML会失效；在这种情况下可能需要隔热/非PML吸收器等替代。
• PML的有效性随着更高的分辨率和适当选择的σx分布而提高，但过大σx会导致数值反射。