SPH Method Based Partly on Riemann Solution in Hydrodynamics and It’s Programing on GPU.docVIP

分类号 O352 学号 ＵＤＣ　 密级 公 开 工学硕士学位论文 水动力学中不完全黎曼解的SPH方法及其在GPU上的实现 硕士生姓名 姜 伟 学科专业 力 学 研究方向 计算流体力学与应用 指导教师 汤文辉 教 授 徐志宏 副教授 国防科学技术大学研究生院 二〇一三年十一月   SPH Method Based Partly on Riemann Solution in Hydrodynamics and It’s Programing on GPU Candidate： Advisor：Prof. Tang Wen-hui A dissertation Submitted in partial fulfillment of the requirements for the degree of Master of Engineering in Mechanics Graduate School of National University of Defense Technology Changsha，Hunan，P.R.China November，2013  目 录 摘 要 i ABSTRACT iii 第一章 绪论 1 1.1 研究背景及意义 1 1.1.1 研究背景 1 1.1.2 研究意义 3 1.2 国内外研究现状 3 1.3 本文的主要内容 6 第二章 水动力学中不完全黎曼解的SPH方法 9 2.1 SPH基本原理 9 2.1.1 函数的核函数近似法 9 2.1.2 函数导数的核函数近似法 10 2.1.3 粒子近似法 11 2.1.4 常用的光滑函数 13 2.2 Navier-Stokes方程组的SPH表达式 14 2.2.1 Navier-Stokes方程组 14 2.2.2 密度的粒子近似法 15 2.2.3 动量守恒方程的粒子近似法 16 2.2.4 能量守恒方程的粒子近似法 17 2.3 SPH方法计算的实施 18 2.3.1 状态方程 18 2.3.2 边界条件 19 2.3.3 近邻粒子搜索 20 2.3.4 时间积分 22 2.4 水动力学中传统带人工黏性项的SPH方法 23 2.5 水动力学中完全黎曼解的SPH方法 24 2.5.1 Godunov格式基本思想 24 2.5.2 HLLC一阶黎曼求解器 26 2.5.3 HLLC二阶黎曼求解器 27 2.5.4 完全黎曼解的SPH方法 29 2.6 水动力学中不完全黎曼解的SPH方法 31 2.7 本章小结 32 第三章 水动力学中的SPH方法在GPU上的实现 33 3.1 从GPGPU到CUDA 33 3.1.1 早期的GPGPU开发 33 3.1.2 CUDA架构及其发展简介 33 3.2 CUDA编程模型 34 3.2.1 主机与设备 34 3.2.2 CUDA线程模型 35 3.2.3 CUDA存储器模型 36 3.3 CUDA软件体系 38 3.3.1 CUDA软件栈 38 3.3.2 CUDA C语言简介 39 3.4 水动力学中含不同修正的SPH方法的GPU程序流程 39 3.5 本章小结 40 第四章 水动力学中不同SPH方法数值试验及对比分析 41 4.1 不完全黎曼解SPH方法的适用性验证 41 4.1.1 模拟参数设置 41 4.1.2 二维溃坝问题 41 4.1.3 孤立波造波理论 43 4.1.4 二维孤立波爬坡问题 44 4.2 不完全黎曼解的SPH方法对压力和数值黏性的改善 45 4.2.1 三维静态问题 45 4.2.2 三维溃坝问题 53 4.2.3 三维孤立波爬坡问题 61 4.3 本章小结 69 第五章 结束语 71 5.1 全文总结 71 5.2 对本课题未来工作的展望 72 致 谢 73 参考文献 75 作者在学期间取得的学术成果 81 表 目 录 表3.1 GPU上各种存储器比较 37 表4.1 二维溃坝、孤立波爬坡问题数值模拟参数设置 41 表4.2 二维溃坝问题粒子数 42 表4.3 二维孤立波爬坡问题粒子数 45 表4.4 三维静态问题不同初始水深粒子数 46 表4.5 三维溃坝问题不同初始水深粒子数 54 表4.6 三维孤立波爬坡问题不同初始水深粒子数 61 图 目 录 图1.1 CPU与GPU算术逻辑单元（ALU）数量对比 1 图1