多孔介质中的波场传播与边界吸收条件

来源 :中国科学院研究生院 中国科学院大学 | 被引量 : 0次 | 上传用户:qiuqiuls
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
数值求解偏微分方程初值问题时,由于计算能力的约束,往往需要引入人工边界条件来限制计算区域的大小。至今,对声波和弹性波方程的吸收边界条件已经有了较为深入的研究,但对多孔介质中波动方程的吸收边界条件的研究却相对较少,特别是非PML型的边界吸收方程还未建立。本文重点研究了时域下三维多孔介质中波传播的无反射边界条件,并建立了相应的精确边界吸收方程。   为了更加全面的比较各种吸收边界条件,首先推导了多孔介质中波动方程的PML条件,并用交错网格有限差分法进行了数值计算。PML是现今较为广泛使用的吸收边界条件,但这是一种非精确的吸收边界条件。然后,在球坐标下,推导一个精确的吸收边界条件。这时的人工边界为半径R的球面B,并假设在B的外部介质是均匀,各向同性且无源的。在球坐标下,结合球面调和函数及向量球面调和函数的性质和数学物理背景知识,将波动方程分解为相互独立传播的压缩波和剪切波的形式,再引入微分算子和积分算子,得到压缩波和剪切波各分量的数量关系,变形处理最终得到一个时间局部,空间非局部的精确的吸收边界条件。由于这个结果中仪含一阶导数,因此可以和已有的求解B内部区域的方法,例如有限差分法,有限体积法很好的结合起来,避免了通常在边界上数值求解高阶偏微分方程的困难。   最后,对建立的新的多孔介质的边界吸收方程用有限差分法进行了数值计算。因为无反射边界条件是以球坐标的形式表示,所以在数值计算中采用球坐标下的差分格式。在球坐标下,波动方程和精确吸收边界条件的部分项的系数为分数形式并且分母上有因子Υ或sin。从而在Υ=0,π即极点时分母为零,差分格式无意义。因此极点处的差分格式的处理比较困难,论文中用控制体积法较好的解决了这个问题。  
其他文献
请下载后查看,本文暂不支持在线获取查看简介。 Please download to view, this article does not support online access to view profile.
期刊
在Hilbert空间中,逼近次微分理论已有很广泛的应用,获得了许多深入的成果,但在一般Banach空间中,这些理论研究最近也受到人们的关注.   在本文中,我主要是在p一致凸且q一致光滑
自1965年,美国控制论专家L.A.Zadeh教授在“Fuzzy sets”[1]一文中,提出模糊集合概念以来,许多学者便投身到模糊集合理论及应用的各个分支的研究之中。无论是理论探索还是实践应
本文研究并行超松弛(SOR)算法和GPU并行隐马尔可夫模型算法。   分布式并行SOR的目标是收敛速度快和通信时间短,本文用重复计算代替通信,将经典的PSOR算法的通信由五次减为
控制理论在以自动化、计算机、机器人等为代表的新技术革命中发挥着重要的作用,它的应用领域已经远远超出了工程技术的范畴。控制理论中的分布参数系统有着非常广泛的应用背景
随着大班额环境的发展,在大班额条件下如何提高学生学习效果的研究也急需发展.在大班额条件下,调动学生学习的积极性、有效性,就要提倡小组合作学习,这种教学方式的出现对提
序方法是处理非线性问题的一个基本的工具,其中锥上算子的不动点问题也是一个基本的问题。自从郭大钧教授和V.Lakskmikantham教授提出混合单调算子以来,它受到众多学者的关注,研
近年来,随着保险业的快速发展,团体寿险越来越受欢迎。团体寿险是一种多生命状态寿险,即多个被保人共有一张保单。目前团体寿险精算都是假定团体内单个生命间是独立的,但现实中,受
1、大力推进网络建设和电子商务 2000年在全国同行尚无成功先例的情况下,李振生同志勇于观念创新,大胆推行网上公开招标寻求合作伙伴和引进人才自主开发相结合的办法,建立了
本文的研究结果主要有两个方面的内容:讨论了全空间RN上一类非线性项在无穷远处为渐近线性的Laplace方程变号解的能量与该问题相应的最低能量之间的关系,通过移动平而的方法,给