收藏本站
《山东大学》 2011年
收藏 | 手机打开
二维码
手机客户端打开本文

基于多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模型的辐射光功率,并利用顶部光子晶体提高了其辐射功率。
【学位授予单位】:山东大学
【学位级别】:博士
【学位授予年份】:2011
【分类号】:TN011;TP391.41

【引证文献】
中国期刊全文数据库 前2条
1 李春荣;袁奇伦;王新政;吕怀武;;非均匀微带线电磁脉冲耦合特性分析[J];电讯技术;2012年10期
2 赖积保;孟圆;余涛;王玉璟;林英豪;吕天然;;一种基于Dual-GPU的三次卷积插值并行算法研究[J];计算机科学;2013年08期
中国博士学位论文全文数据库 前2条
1 陈锐;CSAMT三维交错采样有限差分数值模拟并行算法研究[D];中国地质大学(北京);2012年
2 陈新莲;微纳米结构增强蓝光LED发光效率的研究[D];山东大学;2013年
【参考文献】
中国期刊全文数据库 前1条
1 廖振鹏,黄孔亮,杨柏坡,袁一凡;A TRANSMITTING BOUNDARY FOR TRANSIENT WAVE ANALYSES[J];Science in China,Ser.A;1984年10期
【共引文献】
中国期刊全文数据库 前10条
1 顾媛媛;吴先良;;时域有限差分方法分析微带天线散射参数和方向图[J];合肥师范学院学报;2008年03期
2 丁让箭,吴先良,张玉梅,赵谨;四阶指数差分及其在FDTD中的应用[J];安徽大学学报(自然科学版);2003年02期
3 尹玉,吴先良;基于时域多分辨分析方法的散射场计算[J];安徽大学学报(自然科学版);2005年03期
4 王建永;陆雪平;杨建设;陈秉岩;;一种新的角点和边棱处理方法[J];安徽大学学报(自然科学版);2006年02期
5 赵佳;吴先良;程媛;;MLC-FDTD在波导带通滤波器设计中的应用[J];安徽大学学报(自然科学版);2008年02期
6 ;A 1D time-domain method for in-plane wave motions in a layered half-space[J];Acta Mechanica Sinica;2007年06期
7 蔡二平 ,蒋全兴;用时域有限差分法模拟千兆赫兹横电磁波室的馈电接头[J];安全与电磁兼容;2004年06期
8 陈鑫;胡军;;电磁脉冲对埋地电缆耦合的FDTD仿真分析[J];安全与电磁兼容;2006年06期
9 曾繁春;陈加清;梁妙元;王盼;;冲击泄流时接地体周围电势分布的FDTD分析[J];安全与电磁兼容;2010年01期
10 于大鹏;宋航;;PCB的FDTD-PSPICE混合分析方法研究[J];安全与电磁兼容;2010年05期
中国重要会议论文全文数据库 前10条
1 吴琼;;人工神经网络在基片集成波导滤波器设计中的应用[A];安徽省第五届“兴皖之光”青年学术年会论文集(工科卷)[C];2005年
2 ;Comparative study of the free-surface boundary condition in two-dimensional finite-difference elastic wave field simulation[A];中国科学院地质与地球物理研究所第11届(2011年度)学术年会论文集(上)[C];2012年
3 曹明;于小利;罗中涌;公勋;章德;;声表面波器件的有限时域差分模拟[A];中国声学学会第九届青年学术会议论文集[C];2011年
4 赵永久;梁昌洪;;一种基于场量Taylor展开的吸收边界条件[A];1997年全国微波会议论文集(下册)[C];1997年
5 王春;;一种新型宽带HF单极子天线[A];2005年海峡两岸三地无线科技学术会论文集[C];2005年
6 杨凌霞;葛德彪;魏兵;;金属表面搭接缝隙电磁耦合效应分析[A];2005年海峡两岸三地无线科技学术会论文集[C];2005年
7 胡晓娟;葛德彪;杨利霞;;MC-FDTD分析低频电磁脉冲对钢筋混凝土的穿透特性[A];2005年海峡两岸三地无线科技学术会论文集[C];2005年
8 郑宏兴;郁道银;;添加各向异性参数修正ADI-FDTD算法中的数值色散[A];2005年海峡两岸三地无线科技学术会论文集[C];2005年
9 王均宏;;有耗媒质中目标探测的数值实验[A];信息科学与微电子技术:中国科协第三届青年学术年会论文集[C];1998年
10 王梦麟;张学庆;王新鲁;;战场电磁环境分析模型体系结构研究[A];节能环保 和谐发展——2007中国科协年会论文集(一)[C];2007年
中国博士学位论文全文数据库 前10条
1 刘少刚;基于单一平面电容传感器的木材含水率检测系统研究[D];哈尔滨工程大学;2009年
2 徐刚;不规则波中浮体二阶水动力时域数值模拟[D];哈尔滨工程大学;2010年
3 易煦农;人工电磁材料的光学特性及应用研究[D];华中科技大学;2010年
4 王海涛;月球车热控关键问题研究及材料光谱性质的选择性调节[D];中国科学技术大学;2010年
5 李学识;基于表面波驱动的等离子体天线及其阵列技术研究[D];华南理工大学;2010年
6 赵立华;纳米尺度光波导及微结构若干性质研究[D];南开大学;2010年
7 廖臣;三维电磁粒子模拟并行算法及其应用研究[D];电子科技大学;2010年
8 荆燕飞;线性方程组迭代法与预条件技术及在电磁散射计算中的应用[D];电子科技大学;2010年
9 任志刚;预条件算法及在电磁场数值模拟中的应用[D];电子科技大学;2010年
10 赖生建;计算电磁学中的径向基无网格法[D];电子科技大学;2010年
中国硕士学位论文全文数据库 前10条
1 田禹;煤岩双相介质弹性波场数值模拟[D];山东科技大学;2010年
2 滕巍;左手介质基本电磁特性的时域分析[D];哈尔滨工程大学;2010年
3 刘兴鹏;基于微波同轴谐振腔物料密度测量技术的研究[D];哈尔滨工程大学;2010年
4 王化宇;基于时域有限差分法的微带结构电磁兼容研究[D];哈尔滨工程大学;2010年
5 白彦丽;电磁脉冲对开孔屏蔽腔耦合效应的仿真分析[D];沈阳理工大学;2010年
6 王富明;探地雷达的FDTD建模与仿真[D];南昌大学;2010年
7 卢怡;高速动车组车厢屏蔽效能研究[D];北京交通大学;2010年
8 蔡云秀;基于ANSYS的普洱茶微波干燥过程温度场的数值分析[D];昆明理工大学;2009年
9 夏灿;介质目标电磁散射特性的矩量法分析[D];安徽大学;2010年
10 蒋征;通信系统中有源和无源器件的高功率电磁脉冲效应研究[D];浙江大学;2011年
【同被引文献】
中国期刊全文数据库 前10条
1 童孝忠;柳建新;郭荣文;;复杂二维/三维大地电磁的有限单元法正演模拟策略[J];CT理论与应用研究;2009年01期
2 陈建华;牛中奇;;基于FDTD方法的微带线耦合噪声研究[J];电波科学学报;2009年02期
3 厉天威;阮江军;黄道春;刘忠杰;;大规模电磁场数值计算中并行迭代方法的比较[J];电工技术学报;2007年08期
4 赵国泽;陈小斌;汤吉;;中国地球电磁法新进展和发展趋势[J];地球物理学进展;2007年04期
5 王志刚;何展翔;魏文博;;井中垂直双极源体积分方程法三维模拟研究[J];地球物理学进展;2007年06期
6 吴小平,徐果明,李时灿;利用不完全Cholesky共轭梯度法求解点源三维地电场[J];地球物理学报;1998年06期
7 沈金松;用交错网格有限差分法计算三维频率域电磁响应[J];地球物理学报;2003年02期
8 谭捍东,余钦范,John Booker,魏文博;大地电磁法三维交错采样有限差分数值模拟[J];地球物理学报;2003年05期
9 谭捍东,余钦范,John Booker,魏文博;大地电磁法三维快速松弛反演[J];地球物理学报;2003年06期
10 陈明生,闫述;CSAMT勘探中场区、记录规则、阴影及场源复印效应的解析研究[J];地球物理学报;2005年04期
中国博士学位论文全文数据库 前10条
1 鲁晶津;地球电磁三维数值模拟的多重网格方法及其应用研究[D];中国科学技术大学;2010年
2 仲妍;大型稀疏线性方程组并行求解及预处理技术研究[D];国防科学技术大学;2010年
3 李小康;基于MPI的频率域航空电磁法有限元二维正演并行计算研究[D];中国地质大学(北京);2011年
4 邓居智;可控源音频大地电磁法三维交错采样有限差分数值模拟研究[D];中国地质大学(北京);2011年
5 程汤培;地下水流动数值模拟的高效并行计算研究[D];中国地质大学(北京);2011年
6 李月卉;电磁场数值求解中迭代方法与预条件技术研究[D];电子科技大学;2011年
7 曹国强;基于并行计算技术的离心式叶片泵流场数值模拟及三维设计软件开发[D];辽宁工程技术大学;2005年
8 平学伟;电磁场中的快速有限元分析[D];南京理工大学;2007年
9 林昌洪;大地电磁张量阻抗三维共轭梯度反演研究[D];中国地质大学(北京);2009年
10 张卉;非均匀传输线场路特性研究[D];北京交通大学;2008年
中国硕士学位论文全文数据库 前10条
1 魏卓群;GPU并行免疫算法在冷轧生产调度中的应用[D];大连理工大学;2010年
2 李熙铭;基于GPU的高性能有限元方法研究[D];吉林大学;2011年
3 秦晋;Krylov子空间方法的GPU加速算法研究[D];国防科学技术大学;2010年
4 陈波;基于CPU-GPU异构平台的性能优化及多核并行编程模型的研究[D];中国科学技术大学;2011年
5 黄均才;网络环境下的并行生物计算研究与实现[D];电子科技大学;2004年
6 洪钊岗;MPI-2的分析与实现研究[D];北京科技大学;2005年
7 郝立杰;分布式并行计算中MPI标准的初步实现[D];北京化工大学;2005年
8 谢昊;并行计算在CFD问题中的应用[D];四川大学;2005年
9 张永杰;大型稀疏方程组预处理迭代快速求解技术研究[D];西北工业大学;2006年
10 司华陆;大地电磁三维正反演算法对比研究[D];中国地质大学(北京);2006年
【相似文献】
中国期刊全文数据库 前10条
1 郑士明;关于多项式求根的一个并行算法的收敛性[J];数学研究与评论;1987年04期
2 康立山,陈毓屏;并行算法简介(续)[J];数值计算与计算机应用;1988年04期
3 蒋昌俊;阎春钢;;由数组元素存贮地址求其下标的快速并行算法[J];山东科技大学学报(自然科学版);1989年02期
4 张汝清,胡宁;弹塑性分析的子结构并行算法[J];力学与实践;1991年06期
5 武震东;一类并行算法的提出[J];苏州大学学报(工科版);1993年02期
6 成礼智;任意长—维DFT的MIMD并行算法[J];应用数学学报;1996年04期
7 黄国满;等值线三维显示及其并行算法[J];测绘科学;2001年02期
8 张宝琳,陆金甫,陶应学,杜正平;非线性热传导方程的分组隐式解法及数值结果[J];计算物理;2002年01期
9 杨明辉;矩阵方程AX+XB=C的一种数值算法[J];南京林业大学学报(自然科学版);2002年06期
10 朱旭,吴慧卓;一阶常微分方程系统的并行算法及收敛性[J];西安交通大学学报;2002年04期
中国重要会议论文全文数据库 前10条
1 王元元;曾建潮;谭瑛;;基于并行计算模型的并行微粒群算法的性能分析[A];2007年全国开放式分布与并行计算机学术会议论文集(上册)[C];2007年
2 王轶;席裕庚;;并行的极点配置算法及其新构想[A];1994年中国控制会议论文集[C];1994年
3 杜云飞;王攀峰;富弘毅;周海芳;杨学军;;矩阵LU分解的容错并行算法设计与实现[A];2008年全国开放式分布与并行计算机学术会议论文集(下册)[C];2008年
4 曹淑瑛;郭圣权;樊水康;;并行极点配置算法探讨[A];1995年中国控制会议论文集(下)[C];1995年
5 王韶娟;曾国荪;;分形维数的一个并行算法[A];2005年全国开放式分布与并行计算学术会议论文集[C];2005年
6 陈宁涛;王能超;施保昌;;生物多序列比对的并行算法[A];2005年全国开放式分布与并行计算学术会议论文集[C];2005年
7 孙明珠;王鹏;邱枫;杨玉良;;自洽场理论的实空间求解方法的并行算法以及在嵌段共聚物相分离中的应用[A];2005年全国高分子学术论文报告会论文摘要集[C];2005年
8 李朝鹏;李肯立;;基于分层聚类的并行数据预处理算法[A];2007年全国开放式分布与并行计算机学术会议论文集(下册)[C];2007年
9 姜弘道;余天堂;;有限元并行分析的进展——第九届全国结构工程学术会议特邀报告[A];第九届全国结构工程学术会议论文集第Ⅰ卷[C];2000年
10 张红虎;吴硕贤;;声学辐射度方法的并行实现[A];第十届全国噪声与振动控制工程学术会议论文集[C];2005年
中国重要报纸全文数据库 前10条
1 ;并行算法研究进展[N];中国计算机报;2004年
2 李讯 特约记者 刘程 本报记者 于莘明;银河有颗闪亮的星[N];科技日报;2005年
3 ;名师上讲台:每堂课都是精心编排的演出[N];科技日报;2006年
4 清华大学计算机系 薛巍;电网仿真考验高性能计算[N];计算机世界;2006年
5 陈文光 郑纬民;高性能计算的三大研究领域[N];计算机世界;2006年
6 曙光信息产业有限公司研发中心 温鑫;并行计算软件开发概述[N];中国计算机报;2007年
7 英特尔微处理器技术实验室主任工程师 Anwar Ghuloum;并行编程为什么这么难?[N];中国计算机报;2007年
8 本报记者 赵波 本报特约记者 刘程 本报通讯员 章立丁;一位巾帼学者的无悔征程[N];解放军报;2005年
9 新华社记者 奚启新 本报记者 廖文根;三次选择 无怨无悔[N];人民日报;2005年
10 孙;"生物信息处理专用计算机与算法研究" 通过验收[N];计算机世界;2005年
中国博士学位论文全文数据库 前10条
1 戚晶晶;热物性反问题高效并行算法研究[D];武汉理工大学;2013年
2 付红;量子动力学中并行算法的发展与应用[D];中国科学院研究生院(武汉物理与数学研究所);2013年
3 吴超;信息检索中top-k问题的并行算法及优化研究[D];中国科学技术大学;2011年
4 韩丛英;若干优化问题的并行算法研究[D];上海交通大学;2008年
5 周杰;合成孔径雷达数据处理应用的细粒度并行算法与结构[D];国防科学技术大学;2010年
6 白洪涛;基于GPU的高性能并行算法研究[D];吉林大学;2010年
7 佘春东;数据挖掘算法分析及其并行模式研究[D];电子科技大学;2004年
8 谭国真;时变、随机网络最优路径算法及其应用研究[D];大连理工大学;2002年
9 王建;并行最优化算法与软件设计及数值软件移植[D];中国科学院研究生院(软件研究所);2005年
10 郭荣祥;基于VEGA网格的Dixon结式分布式计算[D];中国科学院研究生院(成都计算机应用研究所);2006年
中国硕士学位论文全文数据库 前10条
1 辛师勃;并行信号处理系统研究[D];南京理工大学;2004年
2 吴青;基于地理信息系统的配电网优化规划的研究[D];广西大学;2004年
3 刘淼;基于平衡分割的并行序列比对[D];西安电子科技大学;2007年
4 潘海琳;高拱坝体形优化决策系统研究[D];河海大学;2006年
5 焦建英;求解微分方程的区域分解两重网格算法[D];西安理工大学;2007年
6 游聪伟;基于GPU平台的KLU并行算法的研究:对角线块的LU分解[D];内蒙古大学;2011年
7 于战华;跨声速欧拉方程并行算法研究及应用[D];南京航空航天大学;2004年
8 平晓慧;最短路径问题的并行算法研究[D];大连理工大学;2006年
9 孙超;快速分形图像压缩并行算法研究与实现[D];曲阜师范大学;2011年
10 彭强;基于并行Boost图库的单源最短路径并行算法的研究[D];华南理工大学;2010年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62791813
  • 010-62985026