论文部分内容阅读
声波散射问题一直受到人们广泛关注,也是数学和物理学中的重要研究领域。相关的研究不仅有理论意义,还有更重要的实用价值,在医学、地球物理、信号图像处理等,尤其是大波数问题。研究大波数散射问题面临的困难主要有两个:解的剧烈震荡性和计算区域的无界性。对于该类问题的模拟逼近,已经存在很多数值方法。但是对于大波数问题产生的“污染效果”,在高维空间中还是没有消除。本文中,利用(Wang et al,2015)提出的新型有限差分方法,我们对一些特殊区域上具有大波数的Helmholtz方程进行了研究,并希望这些工作能促进声波散射问题数值算法研究的发展。 主要工作有: 1.对于具有不可穿透的圆柱形障碍物的三维浅水波中的散射问题,首先,应用 DtN算子将无界区域转化到有界区域。其次,在柱坐标系下对问题进行变量分离将三维问题转化为一系列一维问题,然后利用泰勒展式及方程固有的特性,考虑泰勒展式中无限多项的影响,构造无污染的新型有限差分格式,对问题进行模拟逼近,并给出了一系列的数值算例来验证该格式的可行性。 2.基于上述问题的结果,研究具有可穿透的障碍物的声波传输问题。这类问题的数值模拟面临的困难主要有两个:一是波在传输界面上的突变性,另一个是区域的无界性。对于此类无界传输问题,应用 DtN算子,我们将无界区域转化为有界区域,得到一个包含障碍物的具有人工边界的散射体问题,它包含两个区域,即障碍区域与外散射区域。我们对不同区域采用不同剖分,这样既可以降低计算量,同时又能保证数值解在传输界面上的精度。 3.对于大波数问题,“污染效果”只能在一维方程中彻底消除,对一般高维维问题却不能完全消除。但是随着有限差分格式收敛阶的提高,在一定程度上可以提高模拟效果。以环形区域上一般Helmholtz方程的大波数问题为研究对象,运用方程解的光滑性与波数的关系,构造新型高阶有限差分方法,提高大波数问题的模拟结果。