波动数值模拟中边界稳定性.pdf

摘 要 摘 要 波动数值模拟是力学、地球物理和多个工程学科共同关注的领域，研 究者在基础和应用研究方面皆取得了丰硕成果。然而由边界引入局部失稳 问题，即使在线性范围亦尚未彻底解决。这里“边界”泛指人工边界、物理 边界或不同介质的分界面。本文从波动数值模拟稳定性分析方法着手展开 对局部失稳形成机理和消除方法的研究： 1. 从指导波动数值模拟稳定性分析的 Lax 等价定理出发，依据数值 解收敛方式的不同将定理中的稳定条件区分为强条件和弱条件，进而将稳 定性分析方法归纳为强稳定分析和弱稳定分析。对前者着重阐明了其对分 析边界引入局部失稳现象的重要性；对后者则阐明其与有限元方法相结合 应用于构建波动数值模拟稳定实现方案的价值。 2. 针对完美匹配层应用于正交各向异性介质波动数值模拟中存在的 数值失稳问题，基于一阶连续形式透射边界与有限元方法给出了稳定的数 值模拟实现方案，并采用弱稳定分析初步论证了该方案的稳定性。 3. 应用强稳定分析研究了透射边界引发局部失稳现象的机理及其消 除方法。首先完善了基于一维有限模型得到的透射边界高频振荡失稳机 理。然后从半无限模型出发阐明了透射边界与邻近内节点运动方程所激发 的高频局部失稳现象的机理，建议了一种