收稿日期:20180522; 修回日期:20180628 基金项目: 北京市科技创新服务能力建设—基本科研业务费 (科研类)资助项 目
(71E1810969);北京市科技创新服务能力建设—提升计划资助项目(PXM2017_014224_000009)
作者简介:刘硕(1992),男,河北磁县人,硕士研究生,主要研究方向为 GPU高性能计算(765375297@qq.com);朱希安(1962),男,教授,博
士,主要研究方向为智慧感知与信息处理、虚拟现实等;王占刚(1977),男,副教授,博士(后),主要研究方向为三维建模与可视化研究;厉夫兵
(1982),男,讲师,博士(后),主要研究方向为红外辐射特性建模、光线追踪算法等.
三维灰体辐射传递系数蒙特卡罗的 GPU计算
刘 硕,朱希安,王占刚,厉夫兵
(北京信息科技大学 信息与通信工程学院,北京 100101)
摘 要:为了计算漫反射灰体表面之间的辐射传递系数,建立了三维灰立方体模型,在所建立的三维立体空间
内表面选取发射表面,利用蒙特卡罗方法生成发射光线,并对光线进行光线追踪,计算发射表面与其他表面之间
的辐射传递系数。针对基于蒙特卡罗法计算三维灰体表面之间辐射传递系数时耗时过长的问题,进行了基于显
卡 GPU实现蒙特卡罗并行计算的研究。利用统一编程架构 CUDA对蒙特卡罗法进行了编程实现,采用 CUDA
随机数发生算法并行处理了求解辐射传递系数中的光线追踪过程,实现了计算辐射传递系数的并行方案。实验
结果表明,相对于
CPU串行计算方法,基于 CUDA架构的并行蒙特卡罗法求解辐射传递系数可以取得高达 80倍
的计算加速比。
关键词:蒙特卡罗方法;辐射传递系数;并行计算;GPU;三维灰体
中图分类号:TP301.4 文献标志码:A 文章编号:10013695(2019)11035335704
doi:10.19734/j.issn.10013695.2018.05.0358
MonteCarloforcalculatingradiativeheattransfercoefficientof
3DgraybodybasedonGPU
LiuShuo,ZhuXi’an,WangZhangang,LiFubing
(SchoolofInformation&CommunicationEngineering,BeijingInformationScience&TechnologyUniversity,Beijing100101,China)
Abstract:Inordertocalculatetheradiativeheattransfercoefficientbetweenthediffusesurfaceofthe3Dgraybody,thispa
persetupacellcubemodeltorepresentthegraysurface
,andgeneratedraysemittedfromonesurface(i.e.emissionsurface)
byusingMonteCarlomethod.Itwastracingthemultiplereflectionsforeachrayuntiltheraywastotallyabsorbed.Itrecorded
theabsorbedenergyateachrayfacetintersectiontocalculatetheheattransfercoefficientbetweentheemissionsurfaceandany
othersurface.TosolvetheoverconsumptionissueofradiativeheattransfercoefficientbasedonMonteCarlomethodbetween
thediffusesurfaceofthe3Dgraybody,thispaperusedaMonteCarlomethodbasedonGPUtorealizetheparallelcomputa
tion,andusedtheunifiedprogrammingarchitectureCUDAtogettherealizationofprogrammingtheMonteCarlomethod.Ac
cordingtotheCUDArandomnumbergenerationalgorithm,itparallellycomputedthetracingprocessoftheradiativetransfer
coefficient,andrealizedaparallelschemeforcalculatingtheradiativetransfercoefficient.Experimentalresultsshowthat,com
paredwiththeCPUserialcalculationmethod,theparallelMonteCarlomethodbasedonCUDAarchitecturecanobtainthecal
culatedaccelerationratioof80timesashighastheradiativeheattransfercoefficient.
Keywords:MonteCarlomethod;radiativeheattransfercoefficient;parallelcomputing;GPU;3Dgraybody
0 引言
众所周知,不同温度的物体之间的热量传递主要通过导
热、对流和热辐射三种方式来实现,尤其在真空中,唯一可以依
赖的换热方式就是辐射换热
[1]
。正因为如此,近三、四 十年
来,随着火箭技术、人造卫星、宇宙飞船和核工业的发展,辐射
换热的理论与实验研究受到了人们的极大重视。根据红外物
理理论,物体之间辐射换热的研究中最重要的是物体之间的辐
射传递系数的计算。蒙特卡罗方法是一种通过随机变量的统
计实验来求解数学物理或工程技术问题的数值方法,已有上百
年的历史,但早期由于模拟实验工具的限制,真正用于解决实
际问题的还比较少
[2]
。近二十年来,计算机的高速发展为红
外辐射特性计算的三维建模提供了有利条件;与此同时,蒙特
卡罗法也随之被逐步用于现实的应用中辐射换热问题的计算。
蒙特卡罗法求解表面辐射换热在近年来得到了广泛的应
用。文献[
3,4]介绍了在辐射换热领域,蒙特卡罗法求解辐射
换热是基于概率统计理论的随机抽样方法,利用重复的统计实
验对所研究的物理或数学问题进行直接或间接的描述,真实反
映了事物特点和物理实验过程。文献[
5]分析了在表面热辐射
问题中蒙特卡罗法求解辐射换热结果的精确度问题,在统计实
验特别大时,蒙特卡罗法的结果误差更小,给出了蒙特卡罗法是
求解辐射换热数值的金标准。文献[
6]根据辐射换热中的公式
法,对求解过程分析总结后,提出了漫反射辐射界面特性的单元
表面之间辐射传递系数的新型计算方法,并且与蒙特卡罗法计
算结果进行分析对比,间接证明了蒙特卡罗法的精确度和准确
性,为灰体表面之间的辐射传递系数的研究奠定了结果标准。
随着计算机技术的发展,高性能计算速度从每秒百万次浮
点运算提高到了每秒千万亿次浮点运算,使得三维建模仿真、
软硬件技术、算法实现和计算机应用等方面都得到了巨大的发
展。充分利用 GPU硬件资源对三维模型建立和运算进行并行
通用计 算是 一个 热门 研究 方向,所 以近 年来 很多 学者 利用
GPU实现了各种算法和三维模型仿真。文献[7,8]提出了蒙
第 36卷第 11期
2019年 11月
计 算 机 应 用 研 究
ApplicationResearchofComputers
Vol36No11
Nov.2019