Development of parallel direct simulation Monte Carlo method using a cut-cell Cartesian grid on a single graphics processor

被引:14
作者
Lo, M. -C. [1 ]
Su, C. -C. [1 ]
Wu, J-S [1 ]
Kuo, F. -A. [1 ,2 ]
机构
[1] Natl Chiao Tung Univ, Dept Mech Engn, Hsinchu 30010, Taiwan
[2] Natl Ctr High Performance Comp, Hsinchu 30076, Taiwan
关键词
DSMC; GPU; Cut-cell; Variable time-step; Transient adaptive sub-cell; WAVE REFLECTION TRANSITION; BOLTZMANN-EQUATION; SCHEME; IMPLEMENTATION; FLOWS;
D O I
10.1016/j.compfluid.2014.06.003
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
This study developed a parallel two-dimensional direct simulation Monte Carlo (DSMC) method using a cut-cell Cartesian grid for treating geometrically complex objects using a single graphics processing unit (GPU). Transient adaptive sub-cell (TAS) and variable time-step (VTS) approaches were implemented to reduce computation time without a loss in accuracy. The proposed method was validated using two benchmarks: 2D hypersonic flow of nitrogen over a ramp and 2D hypersonic flow of argon around a cylinder using various free-stream Knudsen numbers. We also detailed the influence of TAS and VTS on computational accuracy and efficiency. Our results demonstrate the efficacy of using TAS in combination with VTS in reducing computation times by more than 10x. Compared to the throughput of a single core Intel CPU, the proposed approach using a single GPU enables a 13-35x increase in speed, which varies according to the size of the problem and type of GPU used in the simulation. Finally, the transition from regular reflection to Mach reflection for supersonic flow through a channel was simulated to demonstrate the efficacy of the proposed approach in reproducing flow fields in challenging problems. (C) 2014 Elsevier Ltd. All rights reserved.
引用
收藏
页码:114 / 125
页数:12
相关论文
共 27 条
[1]  
Bird G., 1994, MOL GAS DYNAMICS DIR
[2]  
Bird G.A., 2006, P 25 INT S RAR GAS D
[3]  
Bird GA, 2005, AIP CONF PROC, V762, P541, DOI 10.1063/1.1941592
[4]   VECTORIZATION OF A MONTE-CARLO SIMULATION SCHEME FOR NONEQUILIBRIUM GAS-DYNAMICS [J].
BOYD, ID .
JOURNAL OF COMPUTATIONAL PHYSICS, 1991, 96 (02) :411-427
[5]   Scalar and parallel optimized implementation of the direct simulation Monte Carlo method [J].
Dietrich, S ;
Boyd, ID .
JOURNAL OF COMPUTATIONAL PHYSICS, 1996, 126 (02) :328-342
[6]  
Eric H., 1994, GRAPHICS GEMS, VIV, P24
[7]  
Gallis M., 2008, 40 THERM C, P3913
[8]   Graphics processing unit based direct simulation Monte Carlo [J].
Gladkov, Denis ;
Tapia, Jose-Juan ;
Alberts, Samuel ;
D'Souza, Roshan M. .
SIMULATION-TRANSACTIONS OF THE SOCIETY FOR MODELING AND SIMULATION INTERNATIONAL, 2012, 88 (06) :680-693
[9]  
Ivanov M, 1997, P PARALLEL CFD, V96, P485
[10]   Numerical analysis of shock wave reflection transition in steady flows [J].
Ivanov, MS ;
Markelov, GN ;
Kudryavtsev, AN ;
Gimelshein, SE .
AIAA JOURNAL, 1998, 36 (11) :2079-2086