PML完美匹配层之中文翻译.docVIP

我们考虑在三维无界空间中传播的时变波场u，并假定该波场的源和初始扰动被限定在矩形区域Ω中，Ω表示为（。。。）。在Ω之外，我们进一步假定波场的传播速度c是一个常量；因此，所有在三维无界空间里传播的波都是由Ω中发出的。在Ω内部，波场u满足 我们希望截断无界外区，从而把计算区域限定在有限计算域Ω里。这种情况下，我们希望确保所有从Ω发出所传播的波没有恶意的反射（spurious reflection）。因此，我们用三个维度上厚度分别为L1 L2 L3的完美匹配层来包裹Ω，以此吸收从Ω里发出的波。在吸收层内部，u满足修正后的波动方程，该波动方程的解在吸收层内随着传播距离按指数规律进行衰减。 根据论文1,2，我们用u^来表示做拉普拉斯变换后的u，定义为 在Ω之外u^满足亥姆霍兹方程 接下来，我们引进坐标变换，把自然坐标x（i）变换到拉普拉斯变换域下的坐标x^（i） 在吸收层内部，阻尼曲线是正的，而在Ω为零。如果我们要求u做拉普拉斯变换后的u^满足坐标变换后的亥姆霍兹方程 那么，众所周知，此时在Ω内的波场u仍然保持不变，而在吸收层内的u将随距离呈指数规律快速衰减；因此吸收层是完美匹配的。事实上PML修正后的亥姆霍兹方程时标准的，但难点在于如何把（2.7）变回到时间域而不引入高阶导数或太多辅助量 从（2.6）坐标变换中可以发现，关于x（i）的偏导数和关于x^（i）的偏导数有关系，即如下所示 我们现