波动数值模拟中透射边界的稳定实现-中国地震信息网.PDF

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