基于计算统一设备架物Fortran的直接模拟蒙特卡洛方法并行优化

上海交通大学学报（自然版） ›› 2013, Vol. 47 ›› Issue (08): 1198-1204.

基于计算统一设备架物Fortran的直接模拟蒙特卡洛方法并行优化

严立，戴欣怡，陈佳洛，王平阳，欧阳华

（上海交通大学机械与动力工程学院，上海 200240）

收稿日期:2012-10-30 出版日期:2013-08-29 发布日期:2013-08-29
基金资助:
上海市自然科学基金（12ZR1414700），上海交通大学研究生创新能力培养专项基金资助项目

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

YAN Li,DAI Xinyi,CHEN Jialuo,WANG Pingyang,OUYANG Hua

（School of Mechanical Engineering, Shanghai Jiaotong University， Shanghai 200240, China）

Received:2012-10-30 Online:2013-08-29 Published:2013-08-29

摘要/Abstract

摘要：

利用基于图形处理器（GPU）的计算统一设备架构(CUDA) Fortran编程平台，对直接模拟蒙特卡洛(DSMC)方法进行并行优化，并以高超声速气动热计算为例，考察了串行与并行计算速度以及不同仿真分子数对并行效率的影响.结果表明，在保证计算精度不变的情况下，程序取得了4～10倍的加速比，并且加速性能高低与计算规模大小成正比.

关键词: 直接模拟蒙特卡洛； Fortran；图形处理器；计算统一设备架构, 气动热

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

中图分类号:

严立，戴欣怡，陈佳洛，王平阳，欧阳华. 基于计算统一设备架物Fortran的直接模拟蒙特卡洛方法并行优化[J]. 上海交通大学学报（自然版）, 2013, 47(08): 1198-1204.

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.

参考文献

［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.