Language
Country/Area
No result found
完美匹配层的神秘吸收技术：如何让波浪悄无声息地消失？
在波动方程的数值模拟中，完美匹配层（PML）作为一种人工吸收层，广泛应用于截断计算区域，特别是在有限差分时域（FDTD）和有限元素法（FE）中。 PML的主要特性是设计上能够使来自非PML介质的波浪在交界面上不反射，进而强效吸收来自计算区域内部的波浪。
PML的作用在于避免波浪的反射，这是其从一般吸收材料中的一大区别。
最早的PML由Berenger于1994年针对麦克斯韦方程组提出，自此以来，对于各种波动方程（如弹性动力学、线性欧拉方程组和Helmholtz方程）也出现了多种相关的重新表述。 Berenger提出的PML是所谓的「分裂场PML」，因为它将电磁场分隔为两个非物理场。而后来更受欢迎的是单轴PML（UPML），它被描述为一种人工各向异性的吸收材料，因为其设计简单且有效。
虽然Berenger的原始公式和UPML最初都是通过手动建构避免平面波在PML交界面上反射的条件来衍生，后来的研究表明这些公式其实可以用更优雅且一般化的方法来推导：伸缩坐标PML。这种观点使得PML能够被推导到不均匀介质中，如波导，并且适用于其他的坐标系和波动方程。
对于设计来吸收沿x方向传播波浪的PML，其在波动方程中的关键变换利用了复数坐标。具体来说，当波动方程中出现x导数时，会被一个带有衰减性质的复数项取代，这使得由于σ（与波的频率有关的吸收系数）为正时，波能量随着传播而迅速衰减系数随着坐标的变换而改变。
完美匹配层的核心机制就在于将波动方程延拓至复数坐标，从而有效地替代了传播波为指数衰减波。
然而，完美匹配层并不完美，它在某些情况下可能导致反射或什至指数增长。 PML在确切的连续波动方程中能够实现无反射，但一旦将波动方程离散化处理时，会出现微小的数值反射。为了解决这一问题，通常会在短距离内将PML吸收系数σ从零逐渐启用。
另外，存在着「反向波」解的情况，如在左手负指数超材料中，这些材料中的群速度和相速度的方向相反，会导致标准PML公式不稳定。幸运的是，对于左手材料，可以通过简单地颠倒σ的符号解决此问题，但这仍面临频率依赖性的挑战。
此外，PML需要介质在与边界垂直的方向上不变，这对于在周期性介质（如光子晶体或声子晶体）或波导进入边界的情况下使PML失效，因为无法保持反射性。
总体而言，完美匹配层作为波动模拟中的一项重要技术，无疑为我们的计算带来了极大的便利。然而，随着技术的进步和对波动行为的深入理解，我们是否能够找到新的方式来克服PML的限制，进一步提高模拟的准确性呢？
