波动数值模拟中透射边界的稳定实现.pdfVIP

第 45 卷 第 4 期 地 球 物 理 学 报 Vol . 45 , No . 4 　　　　　　 　　　　　　 　2002 年 7 月 CHINESE JOURNAL OF GEOPHYSICS 　 July , 2 0 0 2 ( ) [文章编号] 　0001 - 5733 2002 04 - 0533 - 13 [ 中图分类号] 　P631 波动数值模拟中透射边界的稳定实现 1 ,2) 1) 3) 廖振鹏 　周正华 　张艳红 1) 中国地震局工程力学研究所 ,哈尔滨 　150080 2) 哈尔滨工程大学船舶与海洋工程系 ,哈尔滨 　150001 3) 中国水利水电科学研究院 ,北京 　100044 [ 摘 　要] 　从波动能量在计算区内累积增大的观点出发 , 通过简单的一维弹性波模型, 系统 地阐明了在近场波动数值模拟中透射边界两类数值失稳 ———“高频振荡”和“零频飘移”的机 理 :前者源于对波动数值模拟无意义的高频波动在人工边界上的放大和波在有限计算区内多 次反射产生的反复放大; 后者则源于透射边界允许零频和接近零频的分量不断进入计算区. 由此提出了稳定实现透射边界的完整方案包括两项简单措施 :第一, 在全部计算区内按文中 建议的方法注入小阻尼, 以消除高频振荡; 第二, 给出一种具有明确物理意义的消除零频飘移 的算子算法. 最后, 提供了三维波源问题和散射问题的详细数值试验结果. [ 关键词] 　波动, 数值模拟, 近场, 透射边界, 数值稳定性. 1 　引　　言 在地球物理 、土木工程和电磁学等领域内, 只要涉及对局部区域内某种波动现象的解 释, 并在此解释中需要考虑环境介质的影响, 就会涉及近场波动问题. 由于局部区域内的 介质可能具有非均匀组成和不规则形状, 这类问题一般需要通过数值模拟加以研究. 为 了实现波动的数值模拟, 需要引入一个包围所考虑局部区域的有限计算区. 在计算区内 ( ) 用有限元 或有限差分 技术将场方程和物理边界条件转化为离散的运动方程, 从而实现 计算区内波动的数值模拟; 而在计算区的边界上则需建立适当边界条件以模拟外部无限 域对区内波动过程的影响 。由于计算区的边界是人为引入的, 其上的边界条件可称为人 工边界条件. 人工边界条件是近场波动数值模拟的一个特有的重要问题,20 世纪 60 年代 末以来, 关于这一问题已进行广泛而深入的研究, 发表了一系列研究结果[ 1 —8] . 在已提出 的各种人工边界条件中, 我国学者近 20 余年提出和发展的透射理论和在此基础上建立的 多次透射公式简称 MTF 不失为一种可供广泛使用的方案[ 7 , 8] . MTF 在近场波动数值模拟中的稳定实现是实际应用这一人工边界条件的关键问题. 一系列模拟近场波动的数值试验结果表明:MTF 有时能稳定地实现, 并获得正确的数值结 [ 收稿 日期] 　2001 - 03 - 28 收到 ,2001 - 12 - 15 收到修定稿. ( ) [ 基金项 目] 　地震科学联合基金 95 - 07 - 442 . [作者简介] 　廖振鹏 ,男 ,1961 年毕业于清华大学土木系 , 中国工程院院士 ,研究员 , 哈尔滨工业大学船舶与海洋工程 学科教授 ,主要从事工程波动理论和地震工程研究. Email :iemliao @public . hr . hl . cn