资源描述:
《圆形槽波导的fdtd分析》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、圆形槽波导的FDTD分析
2、第1...摘要:圆形槽波导[1]是一种新型的宽频带、大功率,低损耗和大尺寸波导结构,在毫米波、亚毫米波波段有重要应用前景。利用时域有限差分法(FDTD)分析了圆形槽波导的传输特性,计算了圆形槽波导工作在主模500)this.style.ouseg(this)">情况下的色散曲线和截止频率,并与已知的理论值相比较,两者一致性好。关键词:圆形槽波导FDTD;Yee氏元胞;PML完全匹配层;截止波长时域有限差分法(FDTD)是一种近年来发展起来的计算各种复杂电磁问题的计算方法,由Ye
3、eK.S.[2]于1966年首次提出。它对电磁场E、H分量在空间和时间上采用交替抽样的离散方式,每个E(或H)场分量周围有4个H(或E)场分量环绕,应用这种方式将含时间变量的麦克斯韦旋度方程转化为一组差分方程,并在时间轴上逐步推进的求解空间电磁场。FDTD方法是求解麦克斯韦微分方程的直接时域方法。在计算中将空间某一网格点的电场(或磁场)与周围磁场(或电场)直接相关联,且介质参数已赋值给计算空间的每一个元胞,因此结合计算机技术能很方便地处理像具有复杂几何外形的物体。Choi和Hoefer[3](1986年
4、)首次用FDTD分析了波导腔体的谐振问题,计算了其谐振频率。相比于有限元法和表面积分法[4]分析任意截面波导TE和TM模式,FDTD不产生伪模式,且计算相对简单;使用FDTD的另一个优势是,通过一次时域计算,借助Fourier变换可以得到整个通带内的频率响应。圆形槽波导在1989被首次提出[1]。相比于其它波导结构,它有许多优势,如:低损耗、低色散、宽频带、大尺寸,单模工作,高功率容量等。本文利用FDTD分析了圆形槽波导的传输特性,计算了主模500)this.style.ouseg(this)">模的截止
5、频率和色散曲线,并将结果与已知的理论值做了对比,发现二者一致性好。1基本算法和公式FDTD方程是用有限差分式替代时域Maxouseg(this)">假设波导中心为空气,不考虑介质损耗,σ=0,σm=0。j=0,jm=0。则在直角坐标系中,根据FDTD的原理,采用Yee氏元胞,对麦克斯韦方程用中心差分离散,可以得到电磁场6个场分量的FDTD迭代公式,其中Ez和Hz的场分量表达式如下:500)this.style.ouseg(this)">500)this.style.ouseg(this)">其中:μ、ε
6、分别代表网格点的磁导率和介电常数,Δx、Δy、Δz分别代表x、y和z坐标方向的网格空间步长,i、j、k为整数,分别代表x、y和z坐标方向的网格标号或空间步长个数。时间步长为Δt,n表示时间步长个数。其它4个场分量Ex、Ey、Hx和Hy可以由类似方法推出。这样,电场和磁场在时间顺序上交替抽样,抽样时间间隔彼此相差半个时间步长,使得Maxouseg(this)">因为圆形槽波导是一种半开放的结构,而由于计算机容量的限制,FDTD计算只能是在有限区域进行,所以在计算区域的截断边界处必须给出吸收边界条件以限定计
7、算空间。本文采用Berenger完全匹配层PML。PML是一种特殊的介质层,该层的波阻抗与相邻介质波阻抗完全匹配,因而入射波将无反射的穿过分界面而进入PML层。并且,由于PML为有耗介质,进入PML层的透射波将迅速衰减,所以有限几层的PML介质能对入射波起到很好的吸收效果[5-7]。当在波导输入端某处加载激励脉冲以后,通过FDTD迭代,得到的时域上电场和磁场分量将自动满足Max=(i,j,k),其场分量Ψ(m,t)可以表示为所求空间模式的叠加:500)this.style.ouseg(this)">φs
8、为所研究模式k的空间分布函数,fs为其截止频率。采用FDTD法,场分量可以离散化为:Ψn(m,500)this.style.ouseg(this)">,n为时间步数。对Ψn(m,t)进行离散Fourier变换就可以得到频率的离散函数:500)this.style.ouseg(this)">这样由离散Fourier变换和连续Fourier变换的关系可得:500)this.style.ouseg(this)">由上式可以看出,频谱图上那些谐振峰所对应的频率,就是波导内电磁场在各个模式上的截止频率[3,6]
9、。2计算结果图1是圆形槽波导及其横截面示意图。具体尺寸2a=6mm,2c=3.9mm,c/a=0.65mm。计算用了50000步,激励源为置于波导中的一高斯脉冲点源:500)this.style.ouseg(this)">500)this.style.ouseg(this)">采用点源可以激励起多数模式。计算结果如表1,可以发现使用FDTD计算得到圆形槽波导截止波长与文献[1]给出的理论值相符,误差<1%。500)t