Parallel Optimization of Direct Simulation Monte Carlo Method Using Compute Unified Device Architecture Fortran

Abstract

Abstract:

Direct simulation Monte Carlo (DSMC) method is currently the most widely used numerical method for rarefied flow calculation. Based on the Fortran programming platform of GPU (Graphics Processing Unit) and CUDA (Compute Unified Device Architecture), the DSMC method was optimized. Taking the calculation of high supersonic aerodynamic heating as an example, the effect of serial and parallel computing speed and simulation molecules on parallel efficiency was studied. The comparison shows that the parallel results are consistent with each other and parallel program achieved an acceleration of 4 to 10 times. And accelerating performance is proportional to the size of calculation. The application of GPU parallel technology will greatly promote the development of DSMC method.

Key words: direct simulation Monte Carlo (DSMC), Fortran, graphic processing unit (GPU), compute unified device architecture (CUDA), aerodynamic heating

CLC Number:

YAN Li,DAI Xinyi,CHEN Jialuo,WANG Pingyang,OUYANG Hua. Parallel Optimization of Direct Simulation Monte Carlo Method Using Compute Unified Device Architecture Fortran[J]. Journal of Shanghai Jiaotong University, 2013, 47(08): 1198-1204.

References

［1］Bird G A. Molecular gas dynamics and the direct simulation of gas flows ［M］. Oxford: Clarendon Press, 1994.

［2］Wilmoth R G. Direct simulation monte carlo analysis of rarefied flows on parallel processors ［J］. Thermophysics, 1991, 5(3): 292300.

［3］Patnaik G, Obenschain K S. Using GPU on HPC applications to satisfy lowpower computational requirement［C］// 48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition. Orlando, Florida:AIAA,2010: AIAA2010524.

［4］Su C C. Parallel direct simulation Monte Carlo computation using CUDA on GPUs［C］// 27th International Symposium on Rarefied Gas Dynamics. California: American Institute of Physics, 2011: 343347.

［5］Phillips E H. Acceleration of 2D compressible flow solvers with graphics processing unit clusters ［J］. Journal of Aerospace Computing, 2011, 8: 237249.

［6］Lohner R, Corrigan A. Semiautomatic porting of a general Fortran CFD code to GPUs: The difficult modules ［J］. International Journal for Numerical Methods in Fluids, 2012, 69: 314331.

［7］黄昌盛,张文欢,侯志敏,等.基于CUDA的格子Boltzmann方法:算法设计与程序优化［J］. 科学通报, 2011, 56(2829): 24342444.

HUANG Changsheng, ZHANG Wenhuan, HOU Zhimin, et al. CUDA based lattice Boltzmann method: Algorithm design and program optimization ［J］. Chinese Sci Bull, 2011, 56(2829): 24342444.

［8］曾胜田,刘羽,马梦琦. 基于CUDA 的Prewitt算子并行实现［J］. 微计算机应用, 2011,32(11):7175.

ZENG Shengtian, LIU Yu, MA Mengqi. Parallel achieve prewitt operator based on CUDA ［J］. Microcomputer Applications, 2011, 32(11): 7175.

［9］Owens J D, Luebke D. A survey of general purpose computation on graphics hardware ［J］. Computer Graphics Forum, 2007, 26(1): 80113.

［10］王平阳,程慧尔,杨卫华. DSMC方法模拟双组元姿态发动机喷管流动［J］. 推进技术, 2003, 24(2): 111115.

WANG Pingyang, CHENG Huier, YANG Weihua. DSMC simulation of nozzle flow of bipropellant attitude control engine ［J］. Journal of Propulsion Technology, 2003, 24(2): 111115.

［11］Boyd I D, Penko P F, Meissner D L, et al. Experimental and numerical investigations of lowdensity nozzle and plume flows of nitrogen ［J］. AIAA Journal, 1992, 30(10): 24532461.

［12］Kirk D B, Hwu W W. Programming massively parallel processors—A handson approach ［M］. Burlington: Morgan, 2010.

［13］GladKov D. Graphics processing unit based direct simulation Monte Carlo ［J］. Simulation, 2012,88(6): 680693.

[1]	WANG Xiaoliang, YAO Xiaosong, GAO Shuang, LIU Guohua. Aerodynamic Drag Characteristics of Ultra-Low Orbit Satellites [J]. Journal of Shanghai Jiao Tong University, 2022, 56(8): 1089-1100.
[2]	WU Wanghao, DUAN Xu, ZhANG Xin, CHEN Dan, XU Zhendong. Research on the Integrated Calculation Method of Aerodynamics and Heat Transfer for Hypersonic Blunt Body [J]. Air & Space Defense, 2022, 5(3): 87-92.
[3]	CHEN Xiao-qinga (陈晓庆), WANG Yu-jieb (王宇杰), SUN Jian-qia*(孙建奇). Graphic Processing Unit Based Phase Retrieval and CT Reconstruction for Differential X-Ray Phase Contrast Imaging [J]. Journal of shanghai Jiaotong University (Science), 2014, 19(5): 550-554.
[4]	LIU Qiang,XIE Wei,QIU Liaoyuan,XIE Xueshen. Graphic Processing Unit Computing of Lattice Boltzmann Method on a Desktop Computer [J]. Journal of Shanghai Jiaotong University, 2014, 48(09): 1329-1333.
[5]	XU Xian，YANG Jie. High Degree of Accuracy Medical Image Fusion Algorithm Based on GPU Hardware Acceleration [J]. Journal of Shanghai Jiaotong University, 2010, 44(02): 248-0251.