三维非平面断层破裂数值模拟的GPU实现

三维非平面断层破裂数值模拟的GPU实现

论文摘要

在地震动力学研究中,实际发生地震的断层往往不仅具有复杂的表面地形环境,断层本身的形状分布也非常复杂,没有有效的分析解法,一直以数值分析模拟地方式进行研究。三维非平面断层破裂模拟采用曲线有限差分方法,能够模拟复杂地形表面,倾斜、具有不规则几何形状的断层破裂过程,拓宽了有限差分方法所能适用的模型。但在继承了有限差分方法精度高、直观的优点的同时,也同时有有限差分方法计算密集,对计算资源要求较高的缺点。随着摩尔定律的失效,CPU(Central Process Unit)性能提升缓慢,而图像处理单元GPU(Graphics Process Unit)在单机游戏图像处理的市场需求下迅速发展,渲染像素的图像处理过程需要大规模并行计算,为满足这一需要,GPU在硬件设计即为大规模并行运算架构。近年英伟达公司发展的的CUDA编程模型,为将通用计算程序移植到GPU设备上提供了方便。使用GPU并行处理加速计算任务、缩短运算时间,迅速在学界和业界得到大规模应用,并成为新的超算基准。本文从介绍GPU异构编程和CUDA并行编程模型谈起,介绍了三维曲线有限差分方法在描述地震波波场传播和滑动弱化准则下断层自发破裂过程的理论基础。并结合两者利用CUDA并行编程模型进行程序设计,完整实现了三维非平面断层破裂模拟的曲线有限差分程序,即曲线坐标下地震波波场传播过程,地震断层自发破裂作为波场内边界条件与波场的耦合,牵引力镜像法实现地表自由表面的模拟,和CPML吸收边界条件。最后对比了本研究编写的GPU程序和三维非平面断层破裂数值模拟的CPU端多核心MPI程序在全空间均匀介质平面断层和带自由表面的倾向断层的模拟结果,在精度范围内两者具有很好的一致性。并应本用GPU程序带自由表面的粗糙断层的模拟上验证其对复杂断层模拟的普适性。在保证精确性和普适性的基础上,对比了GPU相对CPU在使用单精度浮点数模拟的加速效果,相比于1/4/12/24处理器最高分别获得413.6/181.7/76.6/30.7倍的加速效果。截至目前,本研究为使用GPU加速三维非平面断层破裂数值模拟的首次实现。

论文目录

  • 摘要
  • ABSTRACT
  • 第1章 绪论
  •   1.1 课题背景及研究的目的和意义
  •   1.2 数值模拟平台简介
  •     1.2.1 CPU集群
  •     1.2.2 GPU通用计算
  •   1.3 本文的主要研究内容
  • 第2章 三维非平面断层破裂数值基本原理
  •   2.1 引言
  •   2.2 断层破裂动力学基本方程
  •   2.3 曲线网格坐标系和曲线有限差分方法
  •   2.4 吸收边界条件
  •     2.4.1 指数吸收边界条件
  •     2.4.2 CPML吸收边界条件
  •   2.5 曲线坐标系下地震波波场断层自发破裂内边界条件
  •   2.6 断层分裂节点及临近断层格点的速度更新
  •     2.6.1 断层内边界条件和断层分裂节点的速度更新
  •     2.6.2 临近断层格点的速度更新
  •   2.7 断层分裂节点及临近断层格点的应力更新
  •     2.7.1 临近断层格点的应力更新
  •     2.7.2 断层分裂节点应力更新
  •   2.8 断层分裂节点的荣格库塔积分
  •   2.9 本章小结
  • 第3章 GPU异构计算和CUDA并行编程模型
  •   3.1 引言
  •   3.2 并行编程与GPU并行编程
  •     3.2.1 并行编程
  •     3.2.2 从顺序编程到并行编程
  •     3.2.3 CPU平台并行处理协议简介
  •     3.2.4 GPU与 CPU在并行运算上的对比
  •   3.3 CUDA编程模型和基本概念
  •     3.3.1 内核函数
  •     3.3.2 线程层次结构
  •     3.3.3 存储层次结构
  •     3.3.4 GPU异构编程
  •   3.4 GPU计算能力
  •   3.5 本章小结
  • 第4章 CUDA编程模型实现GPU程序
  •   4.1 引言
  •   4.2 初始化常量及数据结构
  •     4.2.1 系数及常量
  •     4.2.2 波场和断层参数数据结构
  •   4.3 地震波波场
  •     4.3.1 内核函数设计
  •     4.3.2 线程块维度及大小
  •   4.4 滑动摩擦准则及断层分裂节点
  •     4.4.1 内核函数设计
  •     4.4.2 线程块维度及大小
  •   4.5 荣格库塔时间积分
  •   4.6 本章小结
  • 第5章 数值验证
  •   5.1 引言
  •   5.2 全空间均匀介质平面断层
  •     5.2.1 模型简介
  •     5.2.2 波形对比
  •   5.3 带自由表面的倾向断层
  •   5.4 粗糙断层模型模拟结果
  •   5.5 GPU程序加速效果分析
  •   5.6 本章小结
  • 结论
  • 参考文献
  • 致谢
  • 文章来源

    类型: 硕士论文

    作者: 李孟阳

    导师: 张振国

    关键词: 数值模拟,有限差分,非平面断层,断层破裂动力学,计算地震学

    来源: 哈尔滨工业大学

    年度: 2019

    分类: 基础科学

    专业: 数学,地质学,地球物理学

    单位: 哈尔滨工业大学

    分类号: P315;O241.3

    DOI: 10.27061/d.cnki.ghgdu.2019.004867

    总页数: 72

    文件大小: 6146K

    下载量: 35

    相关论文文献

    • [1].非平面元素在视觉传达中的探索与应用[J]. 美术文献 2020(06)
    • [2].平面广告的“非平面”表达[J]. 中国市场 2012(45)
    • [3].平面设计中非平面要素的扩展[J]. 大舞台 2014(01)
    • [4].平面设计的非平面化构建策略探究[J]. 美术大观 2014(07)
    • [5].平面设计中非平面元素的运用探索[J]. 美术教育研究 2013(11)
    • [6].浅析VI设计中非平面元素的表现形式[J]. 美术大观 2010(12)
    • [7].非平面直写式光刻电控系统的研究[J]. 电子工业专用设备 2015(03)
    • [8].面向指路标志系统的非平面交通路网模型[J]. 测绘科学技术学报 2011(06)
    • [9].横向补给高架索系统非平面振动数值研究[J]. 兵工学报 2010(10)
    • [10].现代平面设计的非平面拓展研究[J]. 艺术科技 2016(02)
    • [11].含有一个非平面杂环胺配体的新型反铂抗癌药物的水解机理(英文)[J]. 物理化学学报 2009(11)
    • [12].非平面态屏蔽金属网建模的等效算法研究[J]. 电子测量技术 2008(12)
    • [13].流场非均匀性对非平面激波诱导的Richtmyer-Meshkov不稳定性影响的数值研究[J]. 爆炸与冲击 2019(04)
    • [14].基于非平面机动的弹炮对抗仿真[J]. 火力与指挥控制 2015(06)
    • [15].“非平面阴极”槽生产过程的控制策略[J]. 甘肃冶金 2012(06)
    • [16].一种面向应用的非平面路网数据模型及其实现[J]. 测绘科学 2008(06)
    • [17].材料非平面部位组成与状态的原位统计分布分析表征技术[J]. 冶金分析 2013(04)
    • [18].无线传感器网络中非平面监测区域的分簇路由算法[J]. 计算机与数字工程 2018(04)
    • [19].大视角下非平面场景的图像特征匹配[J]. 清华大学学报(自然科学版) 2010(04)
    • [20].活用面积法解非平面几何问题[J]. 中学生数学 2011(13)
    • [21].非理想情况下通电螺线管的磁矩[J]. 大学物理 2013(11)
    • [22].非平面U型管件制造工艺改进方法[J]. 锅炉技术 2015(02)
    • [23].非平面齿PDC钻头研制与应用[J]. 石化技术 2019(01)
    • [24].非平面9杆巴氏桁架的位移分析[J]. 北京邮电大学学报 2008(04)
    • [25].非平面弯管成形过程的回弹补偿研究[J]. 锻压技术 2009(01)
    • [26].复合材料非平面连接接头强度预测[J]. 宇航材料工艺 2015(04)
    • [27].基于灰度共生矩阵的非平面表面粗糙度的图像纹理研究[J]. 光学与光电技术 2008(06)
    • [28].一类非平面图的色性(英文)[J]. 数学杂志 2018(05)
    • [29].非平面芯片及碲化太阳能电池的研发与分析[J]. 节能 2009(04)
    • [30].典型场景·典型事件·典型细节——电视人物报道人物典型化的表现手法[J]. 新闻与写作 2008(03)

    标签:;  ;  ;  ;  ;  

    三维非平面断层破裂数值模拟的GPU实现
    下载Doc文档

    猜你喜欢