在多层无界区域中Helmholtz方程的数值解法及其运用.pdfVIP

摘要 本文主要提出了一种在无界区域上求解Helmholtz方程的有效数值解法。 首先，我们研究无界区域上Helmholtz方程的区域有界化问题：以往在无界 区域上求解Helmholtz方程的方法是设立假的边界条件将无界的求解区域有界 化，如设这些边界条件为第一、第二或第三类边界条件，这种做法只是将原方程 求解粗糙地近似化，所以无论是对此近似问题进行精确求解还是离散方程后进行 数值计算，解的偏差会变得很大，导致解的逼近效果变差。产生偏差的主要原因 在于这样的边界条件实际上并不符合无界区域中波传播的特性(产生了较强的反 Collino针对上述添加的PML，在 Layer，简称PML)的概念，1997年Francis 数学上相当于对坐标做了一个复的伸展变换量=x+i[xrr(r)dr，把方程转换为一 圳 个新的复方程，我们称之为改进的复Helmholtz方程。 其次，在有界区域上求解Helmholtz方程有许多直接的数值方法，如有限元 和有限差分法等。但若用有限差分或者有限元方法来处理这样一个水平区域很大 的Helmholtz方程时，产生的线性系统的阶数将非常的大，导致