[论文解读] Efficient PML for the wave equation

Marcus J. Grote, Imbo Sim|arXiv (Cornell University)|Jan 2, 2010

Electromagnetic Simulation and Numerical Methods参考文献 19被引用 52

一句话总结

本文提出了一种新颖且高效的二维和三维二阶波动方程完美匹配层（PML）公式，二维仅需两个辅助变量，三维仅需四个辅助变量，避免了将波动方程重写为一阶系统的需求。该方法在使用标准有限差分或有限元方法进行数值模拟时，确保了强稳定性与长时间精度，且在非均质介质和复杂几何结构中表现出良好的鲁棒性。

ABSTRACT

In the last decade, the perfectly matched layer (PML) approach has proved a flexible and accurate method for the simulation of waves in unbounded media. Most PML formulations, however, usually require wave equations stated in their standard second-order form to be reformulated as first-order systems, thereby introducing many additional unknowns. To circumvent this cumbersome and somewhat expensive step, we instead propose a simple PML formulation directly for the wave equation in its second-order form. Inside the absorbing layer, our formulation requires only two auxiliary variables in two space dimensions and four auxiliary variables in three space dimensions; hence it is cheap to implement. Since our formulation requires no higher derivatives, it is also easily coupled with standard finite difference or finite element methods. Strong stability is proved while numerical examples in two and three space dimensions illustrate the accuracy and long time stability of our PML formulation.

研究动机与目标

• 开发一种针对二阶波动方程的稳定且高效的PML公式，避免一阶系统重构带来的计算开销。
• 最小化吸收层中所需的辅助变量数量——二维两个，三维四个，从而降低实现复杂度与内存成本。
• 通过拉普拉斯变换分析，在连续设置下确保PML公式的强稳定性和适定性。
• 通过避免高阶导数，实现与标准有限差分法和有限元法的直接耦合。
• 在二维和三维配置中，包括非均质介质和点源，通过数值模拟展示长时间稳定性和精度。

提出的方法

• 通过在二阶形式的波动方程上应用复坐标拉伸，将PML公式推导于拉普拉斯变换域中。
• 在变换域中引入精心选择的辅助变量，以保持二阶结构并最小化变量数量。
• 仅在垂直于吸收层的方向上应用该变换，确保出射波的指数衰减且无反射。
• 利用标准能量法与半群理论证明稳定性，建立具有无界PML的连续柯西问题的适定性。
• 采用空间上的二阶中心有限差分与时间上的蛙跳格式对PML修改后的波动方程进行空间离散化。
• 通过数值实验调优最优阻尼参数，将阻尼轮廓ζi设为计算域内为零、吸收层内为正的函数。

实验结果

研究问题

• RQ1 能否构建一种针对二阶波动方程的PML公式，避免将其系统转换为一阶形式？
• RQ2 在二维和三维中，为保持稳定性和准确性，PML层所需的最少辅助变量数量是多少？
• RQ3 所提出的PML公式在长时间模拟中是否仍保持稳定，即使在非均质介质或复杂波前存在的情况下？
• RQ4 与现有PML公式相比，该方法在反射误差和长时间稳定性方面表现如何？
• RQ5 该公式能否轻松扩展至有限元法及其他弱形式，而无需引入高阶导数？

主要发现

• 所提出的PML公式在二维仅需两个辅助变量，三维仅需四个辅助变量，与以往的二阶PML公式相比，未知数数量显著减少。
• 二维和三维的数值结果表明，L²误差在时间t=8内持续降低七个数量级，表现出长时间稳定性。
• 该方法在不同阻尼系数（ζ̄i）下均保持鲁棒性，即使在长时间模拟中也未观察到不稳定性。
• 在二维中，点源产生圆形波并向外传播，仅伴随极小的虚假反射，通过快照和误差衰减得到验证。
• 在波速变化的非均质介质中，PML能有效吸收斜入射波，未产生反射或不稳定性。
• 三维点源模拟显示，球形波向外传播，即使在t=1时也未出现可见反射或数值发散。