高效拉盖尔基时域有限差分算法研究
本文关键词:高效拉盖尔基时域有限差分算法研究 出处:《南京邮电大学》2014年博士论文 论文类型:学位论文
更多相关文章: 时域有限差分 加权拉盖尔多项式 无条件稳定算法 吸收边界条件
【摘要】:时域有限差分算法(FDTD)是目前应用最为广泛的电磁场数值计算方法之一,无条件稳定算法、计算效率和精度的提高以及吸收边界条件等课题一直是该算法研究的热点。拉盖尔基时域有限差分算法(WLP-FDTD)是一种以加权拉盖尔多项式作为时域基函数的FDTD法,此方法是无条件稳定的。该算法需要求解大型稀疏矩阵,对计算机内存需求大,计算效率低。近年来提出的高效WLP-FDTD算法,将大型稀疏矩阵转化成三对角矩阵,利用追赶法高效求解,减小了算法的内存需求,提高了计算效率。但是该算法由于引入了微扰项,产生了较大的计算误差。本文重点研究了新的高效WLP-FDTD算法,旨在提高WLP-FDTD算法的计算效率和计算精度。主要创新工作如下: 1.提出了一种新的二维高效WLP-FDTD算法。该算法在传统的WLP-FDTD算法中引入了新的微扰项,将大型稀疏矩阵方程转化成三对角方程求解。和现有的二维高效WLP-FDTD算法相比,减少了微扰项的项数,降低了微扰项中场微分的阶数,并且取消了两步算法中的非物理中间变量,,有效的减少了分裂误差,节省了计算机内存,提高了计算效率。在此基础上,又提出了基于该算法的迭代算法,进一步提高了算法的计算精度。 2.提出了一种新的三维高效WLP-FDTD算法。该算法是本文提出的新的二维高效WLP-FDTD算法向三维算法的推广,和现有的三维高效WLP-FDTD算法的区别在于引入的微扰项不同,可以有效的减少微扰项引入的分裂误差。算法中还引入了迭代算法,进一步提高了算法的计算精度。 3.从理论上证明了新的二维、三维高效WLP-FDTD迭代算法的收敛性;推导出了新的三维高效WLP-FDTD迭代算法的时域分段计算差分方程;并采用时域分段技术以谐振腔为例,验证了三维高效WLP-FDTD迭代算法的数值稳定性。 4.提出了适用于新的二维高效WLP-FDTD算法的PML、UPML吸收边界条件以及适用于新的三维高效WLP-FDTD算法的PML、CPML吸收边界条件,研究了匹配层参数对反射误差的影响规律。
[Abstract]:The finite difference time domain method (FDTD) is currently the most widely used electromagnetic field numerical calculation method of unconditionally stable algorithm, improve the computational efficiency and accuracy and absorbing boundary condition has been a hot topic in the research of algorithm. La Gail Ki finite-difference time-domain method (WLP-FDTD) is a kind of weighted Laguerre polynomials as FDTD method of time domain basis function, this method is unconditionally stable. This algorithm requires the solution of large sparse matrix, demand for computer memory and low computing efficiency. The efficient WLP-FDTD algorithm proposed in recent years, the transformation of a large sparse matrix into three diagonal matrix, by chasing method efficient solution, reduces the memory requirement however, the computational efficiency is improved. The algorithm due to the perturbation, produced a large calculation error. This paper focuses on the efficiency of the new WLP-FDTD algorithm, aiming at improving WLP-FDTD The calculation efficiency and accuracy of the algorithm. The main innovations are as follows:
1. this paper puts forward a new two-dimensional efficient WLP-FDTD algorithm. The algorithm in the traditional WLP-FDTD algorithm is introduced in the new perturbation, the large sparse matrix equation is transformed into solving three diagonal equations. Compared with two-dimensional, the existing WLP-FDTD algorithm, reduce the number of perturbation, the perturbation of low order reduction a midfielder differential, and the abolition of the non physical variables two step algorithm, effectively reduce the splitting errors, saving the computer memory and improve the computational efficiency. On this basis, and the iterative algorithm is proposed based on the algorithm, further improve the precision of the algorithm.
2. this paper puts forward a new efficient three-dimensional WLP-FDTD algorithm. The algorithm is a generalization of new efficient two-dimensional WLP-FDTD algorithm proposed in this paper to the three-dimensional algorithm, and the difference between the existing 3D WLP-FDTD algorithm is introduced, the different perturbation terms, can effectively reduce the perturbation into the split algorithm has introduced error. The iterative algorithm, to further improve the accuracy of the algorithm.
3. theoretically proved that the new two-dimensional, three-dimensional, convergence of the WLP-FDTD iterative algorithm is derived; subsection time domain WLP-FDTD iterative algorithm is a new efficient 3D calculation of the difference equation; and using time domain segmentation technology to the resonant cavity as an example, to verify the numerical stability of three-dimensional WLP-FDTD efficient iterative algorithm.
4., we proposed the PML, UPML absorbing boundary conditions suitable for the new two-dimensional efficient WLP-FDTD algorithm and the PML and CPML absorbing boundary conditions suitable for the new three-dimensional efficient WLP-FDTD algorithm. We also studied the influence of matching layer parameters on the reflection error.
【学位授予单位】:南京邮电大学
【学位级别】:博士
【学位授予年份】:2014
【分类号】:TM15
【相似文献】
相关期刊论文 前10条
1 邢丽;;地震声波数值模拟中的吸收边界条件[J];上海第二工业大学学报;2006年04期
2 赵永久,张新军,梁昌洪;一种简明的频域二阶吸收边界条件[J];西安电子科技大学学报;1998年05期
3 邢丽;;地震声波数值模拟中的透射边界条件[J];科学技术与工程;2009年13期
4 郭立伟;刘晓宇;李世海;;波动数值模拟中的吸收边界条件[J];工程力学;2011年09期
5 彭军;鞠智芹;;一种改进的吸收边界条件[J];弹箭与制导学报;2004年S1期
6 常建梅;冯怀平;冯文杰;;吸收边界条件下散射波场模拟及其误差分析研究[J];振动与冲击;2005年06期
7 袁伟良,梁昌洪;时域有限差分法的优化吸收边界条件[J];西安电子科技大学学报;1998年05期
8 李永明,俞集辉,黄键;时域有限差分法中的吸收边界条件与角点处理[J];重庆大学学报(自然科学版);2001年03期
9 李凤;高益庆;朱泉水;;基于时域有限差分法平面光波的一维吸收边界条件[J];南昌航空大学学报(自然科学版);2010年01期
10 陈志夫;文桂林;卿宏军;王明;王艳广;;用于气动声学计算的完全耦合层吸收边界条件[J];中国机械工程;2014年17期
相关会议论文 前10条
1 赵永久;梁昌洪;任获荣;;高阶频域吸收边界条件的研究[A];1999年全国微波毫米波会议论文集(上册)[C];1999年
2 袁伟良;史小卫;梁昌洪;;吸收边界条件的等波纹优化[A];1999年全国微波毫米波会议论文集(上册)[C];1999年
3 陈军;洪伟;;一种有效的数值吸收边界条件及其在散射问题中的应用[A];1995年全国微波会议论文集(下册)[C];1995年
4 喻志远;薛良金;林为干;;时域有限差分算法中的混合吸收边界条件[A];1995年全国微波会议论文集(下册)[C];1995年
5 贾洪祥;;时域有限差分法电磁场计算采用吸收边界条件的一种二维角点处理[A];第五届全国青年计算物理学术交流会论文摘要[C];2008年
6 魏素花;;声波方程吸收边界条件[A];中国工程物理研究院科技年报(1999)[C];1999年
7 胡驰;;三维修正的完全匹配层吸收边界条件[A];电波科学学报[C];2011年
8 王耘;杨慧珠;;各向异性介质吸收边界条件[A];1994年中国地球物理学会第十届学术年会论文集[C];1994年
9 马延路;周蕙兰;;SEM计算中联合吸收边界条件的初步尝试[A];中国地球物理学会第二十届年会论文集[C];2004年
10 薛东川;戴世坤;;2.5维电磁测深有限元模拟中的吸收边界条件[A];中国地球物理学会第二十四届年会论文集[C];2008年
相关博士学位论文 前3条
1 赵永久;吸收边界条件的研究及其应用[D];西安电子科技大学;1998年
2 陈志夫;基于色散关系分析的高精度气动声学计算方法研究[D];湖南大学;2013年
3 陈正;高效拉盖尔基时域有限差分算法研究[D];南京邮电大学;2014年
相关硕士学位论文 前10条
1 周庆华;吸收边界条件的研究及其在电磁散射中的应用[D];湖南大学;2009年
2 宋鹏;高精度吸收边界条件及声波方程无反射递推算法逆时偏移[D];中国海洋大学;2005年
3 张铭;地震波传播的弹性波方程数值模拟的吸收边界条件[D];同济大学;2004年
4 周丹凤;微波电路的时域谱元快速分析[D];南京理工大学;2010年
5 史杰民;时域有限差分法中积分匹配吸收边界条件的研究[D];华北电力大学;2013年
6 龚亚军;高阶吸收边界条件及有限元方法解电磁边值问题[D];电子科技大学;2011年
7 陶贵生;全人体电磁模型及其在电磁场生物效应中的初步应用研究[D];浙江大学;2004年
8 汪宁;时域有限元法及其二阶和自适应ABC的研究[D];西安电子科技大学;2014年
9 朱长江;射频电磁波有限公式计算技术[D];河海大学;2007年
10 闫玲玲;采用局域吸收边界条件求解激光与原子相互作用的含时Schr(?)dinger方程[D];吉林大学;2009年
本文编号:1401022
本文链接:https://www.wllwen.com/kejilunwen/dianlilw/1401022.html