二维非光滑散射体声波散射场数值解.docVIP

二维非光滑散射体声波散射场数值解 赵海峰，刘继军 东南大学数学系，南京, 210096 摘要： 在实际应用领域, 给定入射波和不可穿透散射体, 声波散射问题的研究有着重要意义. 如果散射体的边界不光滑, 采用标准的算法, 如有限元 (FEM) 或边界元 (BEM), 边界角点会影 响算法的精度. 考虑散射体边界带有角点的声学散射问题. 在二维情况, 此模型可以用 Helmholtz 方程外问题描述. 散射体角点的奇异性会影响算法的数值结果. 本文通过引入人工 边界条件避免角点附近散射场的计算, 并引入完全匹配层 (PML) 将散射问题截断为有界域问 题, 其中 PML 层的外边界赋予零边值 Dirichlet 条件. 最后用间断 Galerkin 有限元求解此边值 问题. 数值结果表明算法是有效的. 关键词：散射波, 非光滑散射体, 间断 Galerkin 有限元, 完全匹配层, 数值算法. 中图分类号： O241.8 On numerical solution to acoustic wave scattering by a 2-dimensional non-smooth obstacle Zhao Haifeng, Liu Jijun Department of Mathematics, Southeast University, Nanjing, 210096 Abstract: For given incident wave and an impenetrable obstacle, the wave scattering problems are of great importance in applied area. In case of the obstacle boundary being nonsmooth, the numerical solution of scattered wave outside of the obstacle by standard schemes such as FEM and BEM will be contaminated due to the singularity of scattered wave at the obstacle corner. Consider the acoustic wave scattering by an obstacle with a corner on its boundary. In 2-dimensional case, this problem is governed by an exterior problem for the Helmholtz equation, which is much complicated for stable numerical computation due to the existence of obstacle corner. We introduce an arti?cial boundary condition (ABC) to eliminate the singularity of scattered wave at the corner and introduce a perfectly matched layer (PML) to truncate the unbounded domain where the zero Dirichlet boundary condition is setup on the outer boundary of PML. Then the discontinuous Galerkin method is applied to solve this boundary value problem. Finally, some numerical examples are presented to show the validity of the proposed scheme. 基 金 项 目： This work is supported by Research Found for the Doctoral Program of Higher Education of China (No.20110092110018) 作者简介： Correspondence author：Liu Jijun（1965-），male，Professor，major research direction：inverse scatterings, medical imaging, industry mathematics. Email: jjliu@seu.edu.cn.