修正辛算法在波动方程计算中的应用研究

修正辛算法在波动方程计算中的应用研究

论文摘要

发展高效的波动方程数值求解方法具有重要的现实意义,它影响到生产与生活的方方面面。用计算机求解偏微分方程同时涉及空间和时间两方面的数值离散。目前,针对空间部分的离散得到了广泛的研究。相比空间离散,对时间离散的研究还比较少。时间的离散也会引入数值频散,研究与空间离散相匹配的时间离散格式具有重要的价值。本文关注波动方程的时间离散方法,在传统辛算法中加入额外的空间离散算子,用来弥补数值计算过程中由于时间离散带来的数值频散,该系列算法被命名为“修正辛算法”。为了验证修正辛时间离散算法与各种常用的空间离散方法相结合的适应性能力,本文将修正辛时间离散与准粒子法空间离散相结合用于求解准粒子体系下的波动方程--修正辛准粒子法波动方程计算;将修正辛时间离散与三角形有限元法空间离散相结合用于求解表面不规则模型--修正辛有限元法波动方程计算。论文从理论分析和数值计算两个方面对修正辛算法进行了研究,主要涉及以下四部分内容:1.研究并总结了现有波动方程数值求解方法的优缺点以及适应场合,通过综合对比各种时间离散方法,发现用二阶辛算法对偏微分方程时间部分进行离散时,具有时间迭代步数少,计算精度高的优点。然而传统二阶辛算法在计算非均匀介质物理现象时,会出现大时间步长计算能力差的问题,其具有进一步改进的空间。本文拟在二阶辛算法基础之上对算法作进一步的发展,并在此算法基础之上推导出修正辛算法。2.研究了准粒子法空间离散体系。粒子法适合解决涉及单元结构、电磁及其耦合等问题的计算。基于Hamilton力学,从分子动力学角度建立了针对地震波传播的空间离散的准粒子体系,给出了粒子间的相互作用系数。3.针对有限元空间离散内存消耗巨大的缺点,研究了在三角形离散情形下的刚度矩阵存储策略,以便实现在计算过程中降低内存的需求。4.将本文构造的修正辛算法与不同的空间离散格式进行时空结合,通过数学方法从理论上研究了它们的数值频散性能以及数值稳定性情况。为了在有限的模型中模拟计算无限大的空间,研究了完美匹配层(PML)吸收边界条件。通过理论分析和具体数值算例进行测试,证实了本文提出的各类修正辛格式在求解波动方程上具有数值频散弱、稳定性强与常见空间离散时空结合性能好等优点。

论文目录

  • 摘要
  • abstract
  • 第一章 绪论
  •   1.1 研究意义及目的
  •   1.2 国内外关于求波动方程数值解的方法现状
  •   1.3 研究内容
  •   1.4 本文的组织结构
  • 第二章 波动方程数值模拟基础
  •   2.1 震源设置
  •     2.1.1 爆炸源
  •     2.1.2 集中力源
  •     2.1.3 剪切力源
  •   2.2 吸收边界设置
  •   2.3 本章小结
  • 第三章 修正辛算法的构造
  •   3.1 哈密尔顿系统与辛几何
  •   3.2 辛Runge-Kutta方法及其相关方法
  •     3.2.1 辛RK方法
  •     3.2.2 辛PRK方法
  •     3.2.3 辛RKN方法
  •   3.3 波动方程的多级辛格式
  •   3.4 修正辛算法的构造
  •     3.4.1 第一类修正辛
  •     3.4.2 第二类修正辛
  •     3.4.3 第三类修正辛
  •     3.4.4 第四类修正辛
  •     3.4.5 第五类修正辛
  •   3.5 本章小结
  • 第四章 修正辛有限差分法波动方程计算
  •   4.1 弹性波方程的Hamilton构造
  •   4.2 稳定性分析
  •   4.3 数值频散分析
  •   4.4 数值计算
  •     4.4.1 二维均匀介质
  •     4.4.2 非均匀介质的长时程计算
  •   4.5 本章小结
  • 第五章 修正辛准粒子法波动方程计算
  •   5.1 弹性波方程的准粒子法空间离散
  •   5.2 时间离散-修正辛格式
  •     5.2.1 稳定性分析
  •     5.2.2 频散分析
  •     5.2.3 数值实验
  •   5.3 声波方程的Hamilton构造
  •     5.3.1 声波方程的准粒子法空间离散
  •     5.3.2 频散分析
  •     5.3.3 数值实验
  •   5.4 本章小结
  • 第六章 修正辛有限元法波动方程计算
  •   6.1 有限元基础
  •   6.2 弹性波方程的有限元化
  •   6.3 第二类修正辛FEM
  •     6.3.1 数值稳定性分析
  •     6.3.2 数值频散分析
  •     6.3.3 数值实验
  •   6.4 第五类修正辛FEM
  •     6.4.1 数值模拟与分析
  •   6.5 本章小结
  • 第七章 总结与展望
  •   7.1 研究成果总结
  •   7.2 创新点
  •   7.3 后续工作及展望
  • 致谢
  • 参考文献
  • 附录A 英文缩略词
  • 附录B 博士研究生期间发表的学术论文
  • 附录C 博士研究生期间参与的科研项目
  • 文章来源

    类型: 博士论文

    作者: 苏波

    导师: 刘知贵

    关键词: 波动方程,修正辛算法,有限差分,准粒子离散,有限元

    来源: 中国工程物理研究院

    年度: 2019

    分类: 基础科学

    专业: 数学,数学

    单位: 中国工程物理研究院

    分类号: O241.82

    DOI: 10.27498/d.cnki.gzgwy.2019.000025

    总页数: 128

    文件大小: 9526K

    下载量: 45

    相关论文文献

    • [1].桥梁在移动荷载作用下动力学响应的广义多辛算法[J]. 振动与冲击 2008(04)
    • [2].辛算法相位误差特性[J]. 上海第二工业大学学报 2015(04)
    • [3].气动声学中声传播问题的辛算法研究[J]. 工程热物理学报 2018(06)
    • [4].带约束动力学辛算法的符号计算[J]. 系统科学与数学 2010(09)
    • [5].基于保辛算法的绳系卫星系统闭环反馈控制问题求解[J]. 计算机辅助工程 2014(02)
    • [6].大阻尼杆振动的广义多辛算法[J]. 动力学与控制学报 2013(01)
    • [7].基于哈密顿动力系统新变分原理的保辛算法之一:变分原理和算法构造[J]. 计算力学学报 2013(04)
    • [8].高阶辛算法的稳定性与数值色散性分析[J]. 计算物理 2010(01)
    • [9].显辛算法在电力系统暂态稳定性计算中的应用[J]. 陕西电力 2011(04)
    • [10].高阶辛算法在典型量子器件模拟中的应用[J]. 光子学报 2015(04)
    • [11].基于辛算法模拟探地雷达在复杂地电模型中的传播[J]. 地球物理学报 2013(02)
    • [12].哈密尔顿偏微分方程多辛算法(英文)[J]. 应用数学与计算数学学报 2013(02)
    • [13].Hamilton体系的辛算法[J]. 赤峰学院学报(自然科学版) 2013(12)
    • [14].含力梯度显辛算法拓广于摄动二体问题[J]. 江西科学 2009(04)
    • [15].Maxwell方程组的多辛算法[J]. 微波学报 2015(01)
    • [16].基于离散变分原理的Birkhoff辛算法[J]. 科技创新导报 2013(19)
    • [17].Hamilton系统下基于相位误差的精细辛算法[J]. 应用数学和力学 2019(06)
    • [18].线性阻尼杆振动问题基于变量变换的多辛算法[J]. 噪声与振动控制 2017(01)
    • [19].非完整系统Boltzmann-Hamel方程的Birkhoff化及其广义辛算法[J]. 物理学报 2012(23)
    • [20].带三次非线性项的四阶Schrdinger方程的分裂多辛算法(英文)[J]. 计算物理 2011(05)
    • [21].辛和非辛算法求解薛定谔方程误差分析[J]. 工程热物理学报 2014(02)
    • [22].用辛算法模拟单摆的振荡(英文)[J]. 中国科学院研究生院学报 2012(06)
    • [23].广义Benjamin-Bona-Mahoney方程的多辛算法[J]. 西北工业大学学报 2008(06)
    • [24].基于哈密顿动力系统新变分原理的保辛算法之三:数值算例[J]. 计算力学学报 2013(04)
    • [25].不平水底浅水波问题的位移法[J]. 水动力学研究与进展(A辑) 2016(05)
    • [26].短波近似的保辛算法[J]. 计算力学学报 2008(01)
    • [27].长波近似水波问题中哈密顿系统的辛几何算法(英文)[J]. 中国科学院大学学报 2013(05)
    • [28].动力学平衡方程的辛两步求解算法[J]. 计算力学学报 2008(06)
    • [29].Kdv方程的辛-谱算法[J]. 湖南工程学院学报(自然科学版) 2009(02)
    • [30].基于哈密顿动力系统新变分原理的保辛算法之二:算法保辛性质证明[J]. 计算力学学报 2013(04)

    标签:;  ;  ;  ;  ;  

    修正辛算法在波动方程计算中的应用研究
    下载Doc文档

    猜你喜欢