谱元方法求解含吸收边界的声学问题

推导并实现了带吸收边界条件的二维波动方程的切比雪夫谱元解法．该解法引入一阶Clayton-Engquist—Majda吸收边界条件，在空间上利用谱元方法，在时间上利用中心差分的积分方法得到波动方程的离散形式，并给出具体算例进行了验证．结果表明：该解法在空间上具有谱精度，在时间上达到二阶精度；与第一类边界条件相比，吸收边界有效地削弱了边界上的数值反射，避免了解的失真；使用中心差分的时间积分方法同隐式积分方法相比，适合波的传播特性，避免了矩阵求逆运算，并且占用的计算机内存小．

Chebychev Spectral Elements Method for Wave Equation with Absorbing Boundary Conditions

A Chebychev spectral elements approximation of the acoustic wave equation with firstorder Clayton-Engquist-Majda absorbing boundary conditions was derived. Its discretization is based on spectral elements in space and central differences method in time. The numerical result, shows that this approximation has spectral accuracy in space and up to second-order in time. Compared with the same wave problems with conventional Dirichlet boundary conditions, the absorbing boundary conditions can reduce the numerical reflection on boundaries and avoid solution distortion. In addition, the central differences method is economic in memory requirement and computing time compared with implicit integral method.