波动数值模拟中的外推型人工边界条件

当前波动数值模拟中的人工边界条件(artificial boundary condition,ABC)数量繁多,但用于串联它们的理论及公式体系还需进一步完善,以便在复杂波动问题模拟中更准确地选取ABC并评估其性能.本工作发展一种外推型人工边界条件理论,将一系列利用临近边界一组节点在前若干时刻的运动来外推人工边界节点运动的经典ABC纳入一个体系.这些ABC包括廖氏透射边界(multi-transmittig formula,MTF)、旁轴近似边界、Higdon边界以及Givoli-Neta、Hagstrom-Warburton、AWWE辅助变量边界等.针对现有边界公式存在的不足,分别提出一种引入多个人工波速进行优化的MTF公式(离散公式)和一种定义在统一局部坐标之上并采用多个人工波速作为控制参数的统一Higdon边界公式(连续公式或微分方程形式ABC),作为外推型ABC的两个基本公式.这二者是最简单实用的外推型ABC,其他同类ABC大多可以由它们转化得到,或者通过某种等价的中间形式与之相关联.数值试验证实了理论的正确性,并初步展示了多人工波速ABC比传统单一人工波速ABC所具有的优势.研究结果不仅具有重要理论价值,还为更好地解决具有差异较大的多种物理波速的复杂波动,如大纵横波速比的软土介质中波动或海洋声学、气象学中的频散波动等的ABC问题提供了实用方法.