- 1、本文档共5页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
一种基于波动方程的新高阶FDTD方法.pdf
第 19卷 第l2期 强 激 光 与 粒 子 束 Vo1.19,NO.12
2007年 l2月 HIGH POWER LASER AND PARTICLE BEAMS Dec.。2007
文章编号: i001-4322(2007)12—2069—05
一 种基于波动方程的新高阶FDTD方法
王 烁 , 邵振海 , 文光俊 , 李世贵。, 谢 康
(1.电子科技大学通信与信息工程学院射频集成电路实验室,成都610054; 2.Micro Design Technology Pte Ltd,Singapore 117674
3.重庆科技学院数理学院,重庆400042; 4.电子科技大学 宽带光纤传输与通信网技术教育部重点实验室,成都 610054)
摘 要: 提出了一种基于二阶波动方程的(2M,4)高阶时域有限差分(FDTD)方法,通过使用辛积分传播
子(SIP)在时域上获得 4阶精度,使用离散奇异卷积(DSC)方法在空域上达到 2M阶精度。与已有的(2M,4)
阶FDTD方法相比,虽然两者都采用SIP和DSC方法,但是此二者的不同点在于:第一,新方法基于二阶波动
方程;第二,在离散计算空间时使用单一网格而不是传统的Yee网格;第三,单独计算某一场分量从而节约内存
并减少计算量。数值计算结果表明,与传统高阶算法相比,基于波动方程的高阶FDTD方法耗费的机时只有
它的5O ,内存消耗下降 1O ,而两者的计算结果之间相对误差小于5‰。
关键词: 时域有限差分; 波动方程; 离散奇异卷积; 拉格朗日多项式; 辛积分传播子
中图分类号: TN0l1; TP3l1.11 文献标识码: A
传统时域有限差分(FDTD)法[1 在解决实际电磁场问题时,需要将6个场分量同时进行交叉迭代运算。
即使在2维TM/TE问题中,基于Yee网格的传统FDTD法至少需要3个分量。而从二阶波动方程出发的
2D-FDTDc。 只需要计算3个场分量中的一个,在导波结构场分析等应用中具有较好的灵活性。另一方面,尽
管FDTD方法有着广泛的应用,但其精度却无法与以格林函数为基础的积分方法相比。近年来,学者们提出
了许多方案来提高FDTD的精度,例如具有(4,4)阶精度的基于龙格一库塔方法的FDTD[3]、基于辛积分传播子
(SIP)方法的FDTD[4 ]、基于拉格朗日多项式插值与SIP方法[E 以及基于离散奇异卷积(DSC)与龙格一库塔方
法的高阶FDTD【7 等。然而,从已发表的文献来看,以上这些方案都是Yee网格的FDTD法。
本文提出一种基于波动方程的具有(2M,4)阶精度的FDTD算法。其特点如下:首先,波动方程是一个场
分量的方程,故该方案在计算时只需要考虑单个场分量,减少了计算场量,从而节约了计算量和存储量。又因
为在离散计算区域时只使用单一网格而不需要引入形如Yee空问的交错网格,因此空间结构得以简化。其
次,该(2M,4)方案在空域和时域上都采用高阶算法,其精度远远高于标准的FDTD。在空域上,采用基于拉格
朗日 核的DSC方法,时域上采用4阶精度的SIP方法。选择适当的计算带宽M值能够大大提高计算精度。
本文介绍基于波动方程的高阶FDTD方法的基本理论,以导波系统为例对所提出的方法进行了计算精度
验证并讨论了该方法的优势。
1 基本原理
1.1 空域2M阶DSC方法
文献[7]提出了一种DSC时域方法,其在计算电磁场问题时,引用DSC方法离散计算空问,构造高阶插值
方程从而提高计算精度。
假设 T表示系统的特征函数,厂(z)为一个测试函数,则.厂(z)的奇异卷积表示为
广+∞
F(£)=T(£)*厂(£)一I 了、(£一z)厂(z)dx (1)
J 一 ∞
其离散的奇异卷积(DSC)形式为
F。(£)一 T。(£一zt)厂(zt)
您可能关注的文档
最近下载
- 2025年浙江安防职业技术学院单招职业技能测试题库(名校卷).docx
- 2024年广州中考英语听说考试——“信息询问” 情景训练课件.pptx
- 四年级下册单元作文《热爱生命》.ppt
- 骆驼奶营养价值课件.pptx
- 统编版高中语文教材中杜甫四首诗《登高》《登岳阳楼》《蜀相》《客至》联读.docx
- 附件3-2:12千伏一二次融合环网柜(箱)及配电自动化终端(DTU)标准化设计方案(2021版).pdf VIP
- 第6课 戊戌变法 课件(共34张PPT).pptx
- 2023年上海杉达学院工程管理专业《管理学》科目期末试卷A(有答案).docx VIP
- 分包考核管理办法和考核表.pdf VIP
- 送教上门康复教案.docx
文档评论(0)