收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

基于多GPU的FDTD并行算法及其在电磁仿真中的应用

杜刘革  
【摘要】:理论、实验与计算相结合已成为科学研究的基本模式,在电磁科学与工程领域中,时域有限差分(FDTD)算法已成为进行电磁场分析的重要方法。FDTD算法是一种麦克斯韦(Maxwell)方程组的时域求解方法,直接将电磁场按照Yee网格的方式进行离散,在空间及时域上利用中心差分近似Maxwell旋度方程中的偏微分,就可以实现电磁场在时域的交替递推。其实现简洁,易于理解,对各种形状以及各种材料的介质有着广泛的适应性;因为FDTD方法直接求解Maxwell方程组,所以各种电磁现象均隐含其中,因此其适用于求解电磁场的辐射、传输及散射等各种问题。自从1966年FDTD由Yee提出以来,也在不断地发展并已广泛地应用于各频段的电磁场仿真领域。 作为一种差分方法,受到数值色散及数值稳定性的影响,为保证FDTD算法的精度,对网格划分有着较为严格的限制。一般其空间步长要小于波长的1/10,当物体结构更为复杂时,空间取样点更要足够多以尽可能真实地模拟物体,而时间步长要满足Courant稳定性条件,与空间步长相关。因此进行电大问题或者精细结构问题的计算时,FDTD方法往往是十分耗时的。 FDTD算法具有天然可并行优势,因此进行并行计算可有效地减少计算时间,加速仿真设计进度。FDTD并行计算主要集中在基于网络设备的并行算法上,如超级计算机以及个人计算机集群,但由于成本及网络速度影响,这种并行方式的性价比并不高;基于可编程器件的FDTD并行算法也得到部分研究者关注,不过由于可编程器件的复杂性以及器件发展问题也并未得到广泛应用。 近年来,图形处理器(GPU)受到游戏市场需求的带动以超过摩尔定律的速度发展,而且其浮点运算能力远高于同时期CPU的运算能力,所以GPU在通用科学计算领域中的应用也逐渐受到关注,如今随着通用图形处理器(GPGPU)技术的迅速发展,GPU已广泛应用于各种通用算法以及各领域的科学计算中,在电磁计算方面特别是FDTD算法上的应用得到了研究者的广泛关注。计算统一设备架构(CUDA)模型出现以后,使得通用图形处理器并行程序的开发更为快速高效,受到科学研究者的欢迎并迅速应用于各学科的计算领域。 本论文研究课题来源于国家重点基础研究发展计划项目:金属/介质纳米异质结构中的局域耦合效应及其在光电转换器件中的应用,本论文研究内容为其中的应用GPU技术进行发光二极管(LED)并行仿真计算系统研究部分,主要研究了基于GPU的FDTD并行算法,最终实现了多GPU平台上的FDTD混合并行运算,极大地提高了利用FDTD算法进行电磁仿真的运算速度,已应用于LED的仿真设计中,进行了LED发光增强研究。论文主要分为以下几个部分: 首先,本论文对研究相关的基础做了介绍,包括电磁计算以及并行计算基础,说明了本文的研究意义以及主要内容,然后对并行计算技术进行了研究,分析了各种并行方法的特点,并对GPU以及通用图形处理器技术的发展应用作深入探讨,研究了CUDA模型的软硬件基础以及编程模型,最终选择CUDA模型作为研究FDTD并行算法的基础。 其次,本文研究了基本FDTD算法原理以及相关知识,如数值色散、边界条件以及激励源等,然后讨论了并行FDTD计算的发展现状,引出本文所要研究的具体内容。 论文提出了一种在CUDA架构下二维及三维FDTD并行算法的实现方式,并实现了二维FDTD算法的各向异性完全匹配层(UPML)吸收边界条件,以及三维FDTD算法的UPML和卷积完全匹配层(CPML)吸收边界条件,实现的入射源包括二维线电流源,三维偶极子源以及平面波入射源,并且在平面波入射源的加入中也实现了一维Mur吸收边界条件的FDTD并行算法。本文提出利用二维线程组织控制电磁场的递推的方式处理二维问题,并提出了多种存储器访问优化方案,包括共享存储器的两种访问方式以及纹理存储器的使用等。在处理三维问题时,本文提出并实现了两种线程组织方案,并对两种方案进行了优化,对比了其计算速度,相对于传统CPU串行算法均达到了10倍以上的加速比。针对UPML和CPML的不同特点,本文采取了扩展PML以及分立计算的不同处理方式,并采取了相应的优化方式,在保证计算精度的前提下,均实现了较高的计算速度,与串行算法相比普遍达到20倍以上的速度提升,最高达到了58倍的加速比。 在单GPU并行计算的基础上,本文将并行算法扩展到多GPU平台。采用FDTD区域分解以及合理的边界交换方案,并利用GPU与CPU内存之间的同步数据传输方案实现了FDTD算法的多GPU并行,为降低数据传输的影响,本文针对多GPU的FDTD算法提出了异步数据传输方案,经验证本方案能够有效地提升多GPU的并行效率。首次实现了GPU内部并行计算,GPU之间并行计算以及数据传输与计算之间的任务并行的FDTD混合并行计算。本文对多GPU算法进行性能测试,包含10层CPML的FDTD算法,在8块GTX295组成的计算平台上达到了4000Mcells/s以上的运算速度。 本文利用GPU运算平台研究了三维FDTD算法中CPML各参数对其吸收效果的影响,进行了微带天线以及滤波器的仿真分析。本文提出了利用FDTD算法计算偶极子辐射功率的方法,在多GPU平台上进行了验证,并利用此方法计算了LED模型的辐射光功率,并利用顶部光子晶体提高了其辐射功率。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 关福宏;陈彬;陈海林;方大纲;俞文明;张力军;;脉冲平面波斜入射旋转对称体电磁散射的FDTD算法(英文)[J];计算物理;2007年03期
2 刘亮元;文伟;;〔A(BC)〕~m结构光子晶体光纤传输特性[J];激光杂志;2011年04期
3 黄重庆,刘靖,安伟科,罗文华,元秀华;波导电路中用数值算法建模时源的研究[J];湖南理工学院学报(自然科学版);2004年03期
4 柴立群,葛德彪,许乔,朱湘琴;光学元件表面微结构的局域场分布特性研究[J];光学技术;2005年01期
5 T.Wang ,晓尘;三维介质中电磁波传播的FDTD模拟[J];勘探地球物理进展;1997年01期
6 张清河;时域有限差分(FDTD)法中的吸收边界条件[J];三峡大学学报(自然科学版);2004年05期
7 简国树,柏菲,潘石;PSTM成像的FDTD数值模拟[J];光学技术;2005年01期
8 周永江;陈朝辉;程海峰;曹义;刘海韬;;用FDTD方法研究颗粒型复合材料微波等效介电常数[J];材料科学与工程学报;2006年06期
9 王飞;葛德彪;魏兵;;磁化铁氧体电磁散射的移位算子FDTD分析[J];物理学报;2009年09期
10 林福明;郑勤红;姚斌;王建敏;李升辉;胥万松;李琳;;用FDTD法分析填充介质金属同轴线的截止频率[J];云南师范大学学报(自然科学版);2011年01期
11 朱湘琴,葛德彪,闫玉波,田春明;轴对称柱坐标系中的FDTD亚网格技术及粗细网格边界反射特性研究[J];汉中师范学院学报;2002年06期
12 蓝朝晖;胡希伟;江中和;李建胜;;基于Haar小波包的时域多分辨分析法[J];电波科学学报;2008年01期
13 冯彦谦;王银;刘四新;;基于FDTD的探地雷达数值仿真与成像研究[J];铁道工程学报;2009年06期
14 董慧;闫玉波;李清亮;;FDTD模拟SLF\ELF水平电偶极子在非均匀地-电离层波导中的场[J];电波科学学报;2010年02期
15 周建;姜洪舟;张林;王琳;曾春年;陈戟;;FDTD方法对二维光子晶体的数值模拟[J];武汉理工大学学报;2008年03期
16 余燕忠;林梅;;基于FDTD方法研究一维孔径菲涅耳衍射[J];泉州师范学院学报;2008年06期
17 李金霞;常明;;理想介质中电磁导波特性的FDTD分析[J];天津理工大学学报;2009年02期
18 李静;曾昭发;黄玲;吴丰收;;三维探地雷达数值模拟中UPML边界研究[J];物探化探计算技术;2010年01期
19 郭斌;曾进启;张焕德;;基于FDTD法等离子体光子晶体的禁带特性研究[J];武汉理工大学学报(交通科学与工程版);2011年04期
20 史林兴,李相银;随机介质中产生随机激光振荡特性的半经典理论分析方法[J];江西科学;2005年04期
中国重要会议论文全文数据库 前10条
1 ;Computation of Scattering from Anisotropically Coated Bodies Using Conformal FDTD[A];2002海峡两岸三地无线科技研讨会论文集[C];2002年
2 ;Three-Dimensional FDTD Analysis of Chiral Discontinuities in the Waveguide[A];2002海峡两岸三地无线科技研讨会论文集[C];2002年
3 曹佳;薛正辉;任武;李伟明;徐晓文;;超宽带脉冲激励下电磁场分布的并行FDTD仿真[A];第二十届全国电磁兼容学术会议论文集[C];2010年
4 张子俊;高本庆;;用FDTD法研究敷有吸波薄层导体目标的RCS[A];1995年全国微波会议论文集(上册)[C];1995年
5 刘小利;孙忠良;周健义;;FDTD法研究复杂结构中的介质谐振电路[A];1997年全国微波会议论文集(下册)[C];1997年
6 肖建明;朱之伟;梁昌洪;;表面分形柱体的电磁散射研究—FDTD方法[A];1995年全国微波会议论文集(上册)[C];1995年
7 马积福;高本庆;;FDTD/表面积分法求金属散射体的远场[A];1995年全国微波会议论文集(上册)[C];1995年
8 兰婧;戴婧怡;牛帅;杨阳;;基于弱条件稳定的FDTD方法的复杂微带滤波器的电磁特性分析[A];2011年全国微波毫米波会议论文集(下册)[C];2011年
9 胡淑欣;窦文斌;;几种T结波导环行器的FDTD法分析[A];2002海峡两岸三地无线科技研讨会论文集[C];2002年
10 袁伟良;梁昌洪;;FDTD求解平面波照射下的导体散射问题[A];1997年全国微波会议论文集(下册)[C];1997年
中国博士学位论文全文数据库 前10条
1 杜刘革;基于多GPU的FDTD并行算法及其在电磁仿真中的应用[D];山东大学;2011年
2 陈智慧;扩展FDTD算法研究[D];西安电子科技大学;2007年
3 孔永丹;基于分裂步长的无条件稳定FDTD算法研究[D];华南理工大学;2011年
4 姚斌;FDTD方法在电磁场本征值问题中的应用研究[D];云南大学;2011年
5 王健;电磁脉冲与舰船典型结构之间电磁相互作用的混合FDTD方法研究[D];上海交通大学;2012年
6 雷继兆;PC和服务器集群下的并行FDTD算法及其应用研究[D];西安电子科技大学;2009年
7 李娟;粗糙面及其与目标复合电磁散射的FDTD方法研究[D];西安电子科技大学;2010年
8 丁海;突变介质区域高精度FDTD算法研究[D];西安电子科技大学;2008年
9 金涛斌;不同结构DGS微波滤波器的传输特性研究[D];天津大学;2011年
10 曹德安;手征介质的FDTD算法研究[D];华南理工大学;2011年
中国硕士学位论文全文数据库 前10条
1 李鹏飞;基于ADI-FDTD的微波加热模型研究[D];电子科技大学;2010年
2 黄明红;FDTD方法在半空间散射问题中的应用[D];西安电子科技大学;2011年
3 李新硕;并行FDTD算法和并行遗传算法及其在复杂天线设计中的应用研究[D];四川大学;2004年
4 王祎君;等离子体电磁散射CDLT-FDTD算法及M-UPML吸收边界研究[D];江苏大学;2010年
5 曾浩;FDTD及其并行算法在粗糙面和目标复合电磁散射中的应用[D];西安电子科技大学;2010年
6 吴锋涛;高效率微带天线阵研究[D];国防科学技术大学;2003年
7 沈琛;基于GPU加速的FDTD算法对电磁辐射与散射问题的研究[D];安徽大学;2010年
8 宋健;微波光子晶体的共形FDTD仿真研究[D];西安电子科技大学;2004年
9 李菲菲;三维FDTD并行算法的研究及应用[D];西安电子科技大学;2011年
10 荆敏娟;单光纤光镊的FDTD模拟分析及实验设计[D];燕山大学;2010年
中国重要报纸全文数据库 前10条
1 记者 严天秀;昆明新机场全球首安“GPU”[N];昆明日报;2010年
2 ;GPU市场Q2创八年最大增幅[N];计算机世界;2009年
3 本报记者 宋家雨;GPU还是CPU? 高性能计算的一个问题![N];网络世界;2010年
4 ;并行算法研究进展[N];中国计算机报;2004年
5 电脑报 朱文利;GPU助力,高性能计算机“异构化”[N];电脑报;2009年
6 本报记者 霍光;GPU计算逐渐大规模商用[N];中国计算机报;2009年
7 游讯;图形处理器GPU[N];人民邮电;2011年
8 记者 刘欢;全国首个GPU计算应用测试中心在京成立[N];北京日报;2011年
9 李大庆;首套分布式GPU超级计算系统启用[N];中国技术市场报;2010年
10 本报记者 霍光;CPU集成GPU挑战主板设计[N];中国计算机报;2010年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978