声场一阶与二阶吸收边界条件下有限元实现 Absorbing Boundary Condtion for Acoustic Computation

在处理声场问题时，有一类重要的问题就是吸收边界问题。它主要解决声波传播时避免反射波影响原始波形和波速。 通常数值模拟吸收边界条件都会采用 Engquist and Majda 吸收边界并使用经典的Mur有限差分法。 使用有限差分法处理已有十几年时间了，差分步骤繁琐并且不适用于任意形状的边界。 笔者使用有限元法（FEM）结合拉格朗日乘子（Lagrange Multiplier）来耦合声波方程和Engquist and Majda 一阶和二阶吸收边界条件，并且在时域和频域比较一阶和二阶吸收效果。 这篇文章基于我的论文 Absorbing Boundary Condtion for Acoustic Computation， 有兴趣的读者可以与我联系。以下是论文第一章简介部分。

 

Introduction

 

Many computational problems in acoustics, as well as in other elds of application like electromagnetic and geophysics, involve waves in an unbounded medium. However, the solution of such problems using the Finite Element (FE) method and Finite Difference (FD) method or other domain-type methods usually requires the use of a fnite computational domain in which the entire calculation is to be done. Thus, one has to introduce an articial boundary that encloses the region of interest . We want such a boundary to have the property of transparency whi