【摘 要】
:
龙格库塔间断有限元(Runge-Kutta discontinuous Galerkin,RKDG)算法是一种高精度的流场算法,它具有能够捕捉间断,处理网格悬点等优点.然而由于DG算法的计算量大,即使是结构
【机 构】
:
北京航空航天大学数学与系统科学学院,北京,100191
【出 处】
:
第十六届全国流体力学数值方法研讨会
论文部分内容阅读
龙格库塔间断有限元(Runge-Kutta discontinuous Galerkin,RKDG)算法是一种高精度的流场算法,它具有能够捕捉间断,处理网格悬点等优点.然而由于DG算法的计算量大,即使是结构良好,编写高效的串行程序也要花费很长的时间才能完成流场求解.CUDA是近年来非常流行和活跃的高性能并行计算平台.CUDA具有编程简单,加速效果明显等优势,从而在科学计算领域迅速普及开来.本文的CUDA解法器是二维非结构网格上基于RKDG算法的GPU实现.文章使用NACA0012作为测试算例,串行程序使用Intel C++编译器编译(开启O2开关),运行在Intel Xeon X5690(主频3.47GHz)主机上;CUDA程序使用GCC4.4.3和NVCC5.0编译器编译(开启O2开关),运行在单块Nvidia Tesla M2070 GPU上.算例的网格规模为9340个三角形单元,马赫数0.4,攻角5度,从0开始计算流场到15秒.程序结果表明串行程序的墙上时间为24823秒(6.9小时),CUDA程序的墙上时间为995秒(16分钟),程序的加速比达到25倍.随着问题规模的增大,该加速比将进一步提升.
其他文献
目前CFD中应用广泛的数值方法都是基于网格离散流场,对于日益复杂的实际工程问题,高质量网格的生成已经成为数值研究面临的巨大挑战之一.对于超音速反应流,粘性的影响通常可以
炸药两点起爆情况下柱形爆轰驱动飞层问题同时涉及炸药爆轰驱动、多飞层相互作用、金属熔化、断裂等内容,是典型的高维多介质弹塑性大变形问题.张崇玉等人在文献中通过对铅材
本文发展了一套基于JASMIN框架的光滑粒子动力学并行程序.该程序被应用于液滴撞击液膜过程的数值模拟.不可压NS方程采用弱可压SPH方法进行求解.基于界面重构的算法被应用于表
本文基于叉树数据结构的自适应笛卡尔网格方法能够简化网格生成过程,并能实现流场信息剧烈变化区域的动态捕捉和提高流场解的分辨率。本文基于流场特征的网格自适应技术研究主
可压缩多介质流动和多相流问题在天体物理、惯性约束聚变(ICF)等很多科学与工程领域有着广泛的应用背景,由于存在强间断、大变形以及多相混合等复杂流动特征,如何对其有效模
在计算风工程中上风向的湍流问题是一个始终无法回避的问题.在模拟过程中很少考虑其对流场的影响.以方形柱体为研究对象,本文模拟了单个涡街的形成过程基于非稳态雷诺应力方
拉氏方法(Lagrangian)和ALE(arbitrary Lagrangian and Eulerian)方法具有自动处理界面运动、获得高清晰界面的优点,其在可压缩多介质领域如天体物理、惯性约束聚变(ICF)的数
在受限管道中,液滴遇到障碍物在石油工业、化学工业、微流体等领域有着广泛的应用是经常遇到的问题.本文采用基于Shan-Chen多相模型的(多松弛)格子Boltzmann方法对液滴穿过障
本文提出了一类新的基于相场理论的不可压格子Boltzmann模型.与已有的扩散界面模型相比,新模型可以模拟低粘性,有较大密度差流体的流动.新提出的模型利用了两类分布函数:一类
本文主要发展了一种新的数值求解球面上准地转方程的Fourier有限体积元方法.在球坐标下,经度方向用Fourier基近似,而纬度方向用有限体积元离散.在提出的数值方法中,为了得到