激波与火焰面相互作用数值模拟的GPU加速

计算物理 ›› 2016, Vol. 33 ›› Issue (1): 23-29.

激波与火焰面相互作用数值模拟的GPU加速

蒋华, 董刚, 陈霄

南京理工大学瞬态物理重点实验室, 江苏南京 210094

收稿日期:2014-12-04 修回日期:2015-04-03 出版日期:2016-01-25 发布日期:2016-01-25
作者简介:蒋华(1988-),男,博士生,主要从事计算流体力学研究,E-mail:18951895841@126.com;董刚(1970-),男,博士,研究员,主要从事计算流体力学研究,E-mail:dgvehicle@yahoo.com
基金资助:
国家自然科学基金（11372140）和江苏省研究生科研创新计划（KYLX_0328）资助项目

A Numerical Study of Interactions Between Shock Waves and Flame with GPU Acceleration

JIANG Hua, DONG Gang, CHEN Xiao

Key Laboratory of Transient Physics, Nanjing University of Science and Technology, Jiangsu Nanjing 210094, China

Received:2014-12-04 Revised:2015-04-03 Online:2016-01-25 Published:2016-01-25

摘要/Abstract

摘要： 为考察计算机图形处理器（GPU）在计算流体力学中的计算能力，采用基于CPU/GPU异构并行模式的方法对激波与火焰界面相互作用的典型可压缩反应流进行数值模拟，优化并行方案，考察不同网格精度对计算结果和计算加速性能的影响.结果表明，和传统的基于信息传递的MPI 8线程并行计算相比，GPU并行模拟结果与MPI并行模拟结果相同；两种计算方法的计算时间均随网格数量的增加呈线性增长趋势，但GPU的计算时间比MPI明显降低.当网格数量较小时（1.6×10⁴），GPU计算得到的单个时间步长平均时间的加速比为8.6；随着网格数量的增加，GPU的加速比有所下降，但对较大规模的网格数量（4.2×10⁶），GPU的加速比仍可达到5.9.基于GPU的异构并行加速算法为可压缩反应流的高分辨率大规模计算提供了较好的解决途径.

关键词: 激波, 火焰界面, GPU, 异构系统, 加速比

Abstract: To study performance of graphic processing unit (GPU) for computational fluid dynamics, simulation on interactions between shock waves and a flame interface, a typical compressible reactive flow, was carried out on CPU/GPU heterogeneous system. Several optimal strategies are taken to raise GPU code performance. Computational results and acceleration performance of GPU with different grid number were examined. It was found that computational results by parallel GPU are the same as those by traditional CPU based on MPI parallel of 8 threads. Computational times of two parallel methods linearly grow with increase of computational grid numbers. Compuational cost by GPU is less than cost by MPI. As grid number is small(1.6×10⁴), speedup ratio of 8. 6 was achieved on GPU. As grid number grows, speedup ratio decreases. However, a ratio of 5.9 still can be held on GPU when grid number is more (4.2×10⁶). Arithmetic on CPU/GPU heterogeneous system provides a good solution for large-scale computations of compressible reactive flows.

Key words: shock wave, flame interface, GPU, heterogeneous system, speedup ratio

中图分类号:

O354.5

蒋华, 董刚, 陈霄. 激波与火焰面相互作用数值模拟的GPU加速[J]. 计算物理, 2016, 33(1): 23-29.

JIANG Hua, DONG Gang, CHEN Xiao. A Numerical Study of Interactions Between Shock Waves and Flame with GPU Acceleration[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2016, 33(1): 23-29.

参考文献

[1] ELSEN E, LEGRESLY P, DARVE E. Large calculation of the flow over a hypersonic vehicle using a GPU[J]. Journal of Computational Physics,2008,227(24): 10148-10161.
[2] CPHEN J M, MOLEMAKER M J. A fast double precision CFD code using CUDA[J]. Journal of the Physical Society of Japan,2009,1(2): 237-341.
[3] BRANDVIK T, PULLAN G. An accelerated 3D Navier-Stokes solver for flows in turbomachines[J]. Journal of Turbomachinery,2011,133(2): 021-025.
[4] DONG T X, LI X L, LI S, et al. Acceleration of computational fluid dynamics codes on GPU[J]. Computer Systems & Applications,2011,20(1): 104-109
[5] LIU Feng, LI Hua, TIAN Zhengyu,et al. Heterogeneous parallel compressible flow solver based on MPI+CUDA[J]. Journal of National University of Defense Technology,2014,36(1): 6-10
[6] HUANG Shuo, XIAO Jinyou, HU Yucai, et al. GPU-accelerated boundary element method for Burton-Miller equation in acoustics[J]. Chinese J Comput Phys,2011,28(4): 481-487.
[7] KHOKHLOV A M, ORAN E S, CHTCHELKANOVA A Y,et al. Interaction of a shock with a sinusoidally perturbed flame[J]. Combustion and Flame,1999,117(1-2): 99-116.
[8] DONG G, JIN J M, YE J F, et al. A study of flame instability with wave propagation algorithm[J]. Chinese Journal of Computational Physics,2005,22(4): 371-376.
[9] DONG G, FAN B C,YE J F. Numerical investigation of ethylene flame bubble instability induced by shock waves[J]. Shock Waves,2008,17(6): 409-419.
[10] ZHU Yuejin, DONG Gang, FAN Baochun. Three-dimensional computation of the interactions between shock waves and flame in a confined space[J]. Journal of Propulsion Technology,2012,33(3): 405-411.
[11] YANG J, KUBOTA T, ZUKOSKI E E. Applications of shock-induced mixing to supersonic combustion[J]. AIAA Journal,1993,31(5): 854-862.
[12] ORAN E S, GAMEZO V N. Origins of the deflagration-to-detonation transition in gas-phase combustion[J]. Combustion and Flame,2007,148(1): 4-47.
[13] CHAKRAVARTHY S R, OSHER S. A new class of high accuracy TVD schemes for hyperbolic conservation laws[J]. AIAA Paper,1985,85-0363.
[14] JIANG Hua, DONG Gang, CHEN Xiao. Numerical study on the effects of small-amplitude initial perturbations on RM instability[J]. Chinese Journal of Theoretical and Applied Mechanics, 2014,46(4): 544-552.
[15] BROUILLETTE M. The Richtmyer-Meshkov instability[J]. Annual Review of Fluid Mechanics,2002,34(34): 445-468.
[16] IMAN Gohari S M, ESFAHANIAN V, MOQTADERI H. Coalesced computations of the incompressible Navier-Stokes equations over an airfoil using graphics processing units[J]. Computer and Fluids,2013,80(7): 102-115.