透射边界与无穷远辐射条件.pdf

维普资讯 中 橱 研 学 (E 辑) 第 31卷 第 3期 SCIENCEINCHINA fSeriesE) 2001年 6月 透射边界与无穷远辐射条件 廖振鹏 (哈尔滨工程大学船舶与海洋工程 系．哈尔滨 150001 摘要 研究了无穷远辐射条件和用于近场波 动数值模拟 的透射边界之 间的关系． 结论为：透射边界近似地等价于一大类无限舟质 中的无穷远辐射条件，这 一近似 的误 差和用有限元或有限差分模拟波动的数值误差具有相 同的数量级．这一结果为用于 近场波动数值模拟的透射边界提供 了健垒的理论基础，并对将透射边界应用于虼数 值模拟所积累的主要经验提供 了垒面的解释 ． 关键词 波动数值模拟 无穷远辐射条件 局部人工边界条件 多次透射公式 在无 限介质 内波动现象的解析研究中，出于物理考虑需要对无穷远处 的波动行为施加一 定 限制 ，这种限制称为无穷远辐射条件_l(简称辐射条件)，它们对保证物理上适定的波动问题 的解析解的惟一性至关重要．与波动解析研究中的辐射条件相对应 ，在利用有 限元或有限差 分进行相应 的数值研究时 ，需要在有限计算区的虚拟人工边界上加上所谓人工边界条件 (arti． ficialbo,md~-ycondition，简记为ABC) 模拟外部无 限域对计算区内波动的影响并保证数值解 的惟一性 ．一类 ABC满足计算 区之外无穷域 内的所有数学物理条件 ，因而它们是和辐射条件 精确地等价的．在许多领域中已对这类 ABC作了大量研究，并在很多情况下已成功地用于解 决有关近场波动 问题_2“J．遗憾 的是 ，就时步积分解法而言，此类 ABC在每一时步 中不仅将 所有边界点的运动耦联起来 (空间耦联)，而且将从激励开始到所考虑的时刻之 间的所有时刻 的运动耦联起来 (时间耦联)．这一全局特征使此类 ABC在某些情形下不实用 ．例如 ，当计算 区内出现非线性而外区仍保持线性和计算模型的自由度总数极大 的情形L5J．这类 ABC的缺 点导致寻求另一类具有局部特征的ABC．它们在时间和空间上都是局部的，即它们仅涉及 与 所考虑边界点邻近 的点在一个短时间窗 内的运动_6“J．局部 ABC的吸引力不仅在于它们在 保证数值计算效率方面至关重要 ，而且在适当选择和正确使用的条件下它们也具有足够的精 度 ．例如 ，对于二维和三雏基准情形的一系列数值试验结果令人信服地表 明，利用透射边界_l一 所获得的数值解的精度与满足辐射条件的数值解的精度难分高下_8 …．但是 ，透射边界及其 他常见的局部 ABC，从物理观点看来均建立在平面波假定之上 ，即其导出和精度分析在物理上 均基于假定射 向人工边界的波动为平面波_3l4J．这一假定显然和辐射条件所要求的波 函数在 无穷远处为零 不相 容．如何在理论上 理解透射边 界在 近场波 动数值 模拟 中的 出色表现 20004~．17收稿 ，2000-07-24收修改稿 *地震科学联台基叠 资助项 目 * 目前通信地址 ：中国地震局工程力学研究所，哈尔滨 150080 维普资讯 第 3期 廖振鹏 ：透射边界与无穷远辐射条件 255 就成了一个有趣的疑 问．此外 ，按照平面波假定 ，局部 ABC的误差应与人工边界至波源 区的 距离无关，但是 ，人工边界却必须经验地设置在源区的一定距离之外 ，尽管这一距离不大 ．为 了回答上述疑 问，需放弃平面波假定 ，进一步研究局部 ABC的理论基础 ．本文试 图通过阐明 透射边界与辐射条件之间的关系探讨这一理论基础 ． 1 射 向人工边界的波动 本节将从无穷远辐射条件 的经典理论出发导出射 向人工边界的波动的一维表达 ． 1．1 无穷远辐射条件的近场表达 考虑在均匀、各 向同性和线弹性介质 中从分布源区 射出的波动 ，假定人工边界 B 是凸 的，且包古区域 (图 1)． 首先考虑