论文部分内容阅读
科学与工程中的许多问题,如时间调和声波遇到不可穿透障碍物的散射、海洋深处潜艇声呐探测系统及探测空间飞行器的电磁波的绕射与辐射等,都可归结为Helmholtz方程的外边值问题。因此,采用边界型数值方法求解此类问题要比域型方法,如有限差分法及有限元法等,更加的有效。目前,求解Helmholtz问题的边界型数值方法主要是基本解法和直接边界元法。基本解法通过虚拟边界避免奇异积分的计算,然而虚拟边界的合理选择是一个棘手的问题,通常凭借研究者的经验或误差实验来完成,因此缺乏客观性。传统的边界积分方程在求解Helmholtz外问题有时会出问题,亦即如果波数k是Dirichlet内问题或者是Neumann内问题的特征值时,传统的边界积分方程没有唯一解。反映在数值实施过程中是,当波数k等于或接近内问题的特征频率时,边界积分方程的影响系数矩阵具有非常大的条件数。值得注意的是,边界积分方程解的不唯一性纯粹是由数学公式本身的缺陷导致的,与原物理问题无关,即没有任何物理意思。因此,寻求求解外Helmholtz问题的避免内问题的虚拟特征频率的新的边界积分方程是该领域重要的研究课题。 避免虚拟频率最有效的方法是直接边界元法中的Burton-Miller公式和间接边界元法中的Le is公式。前者是将传统的直接变量边界积分方程与相应的梯度积分方程进行线性组合,后者是将传统的间接变量边界积分方程与相应的梯度积分方程进行线性组合。理论上,二者都能够在全频率上提供唯一解。然而, Burto n-M iller公式和Le is公式都包含超奇异积分,给数值计算带来了很大的困难。尽管已有一些处理超奇异积分的方法,但是它们的数值实施过程很复杂和麻烦,不便于大规模计算。 针对Le is公式中出现的超奇异积分,本文提出了新的规则化边界元的理论和方法。通过构造新的边值问题,将Le is公式中的超奇异积分转化为强奇异积分,进而对强奇异积分进行规则化,使得规则化Le is公式中排除了强奇异积分和超奇异积分的计算。二、三维Helmholtz问题的声辐射和散射的数值算例表明,本文方法具有程序设计容易,计算效率高,结果精度好的特点。因此,本文方法不仅能够求解全波数三维Helmholtz问题,更重要的是为超奇异积分的规则化计算提供了新的途径,是对边界元规则化理论的重要贡献。