西安电子科技大学
硕士学位论文
轴对称柱坐标FDTD网格剖分技术及其在天线辐射中的应用
姓名:朱湘琴
申请学位级别:硕士
专业:无线电物理
指导教师:葛德彪
20021201
摘要1
摘要
时域有限差分(FDTD)方法是计算天线辐射的一种有效方法。为了提高二维轴对称柱坐标系中FDTD方法的有效性和准确性,本文在网格剖分方面研究了亚网格技术和曲面共形技术。给出了粗细网格尺寸比为2、3的两种亚网格技术的理论分析,完成了相应程序并给出相应算例。为了更糖确的模拟二维轴对称柱坐标系中的曲面导体,本文研究了二维轴对称柱坐标系中的曲面共形技术并计算了球底面和平底面的圆锥天线的辐射特性。实际应用中等离子体和吸波材料等都具有色散介质特性。本文还应用移位算子方法研究二维柱坐标系中等离子体覆盖的单极子天线的辐射特性。
关键词:时域有限差分方法轴对称亚网格技术曲面共形技术等离子体
摘要
Abstract
TheFinite.-DifferenceTime-?Dommn(FDTD)techniquehasbeeneffectivelyappliedtotheanalysisofradiatingproblems.InordertoimprovethevalidityandtheaccuracyofFDTDmethodincylindricalcoordinatewithrotationalsymmetry,twokindsofsubgriddingmethodswithmeshrefinementbyafactoroftwoandthreearediscussed.Thesubgriddingtheoriesaregivenandthepro酽amsaresetup.Tosimulatethecurveperfectlyconductingobjectincylindricalcoordinatewithrotationalsymmetrymoreprecisely,wepresentasimpleyetaccurateconformalalgorithmfortheFDTDanalysisofantennaradiatingproblems.Accuracyofthismethodisillustratedbycomparingtheresultsderivedfromthismethodwiththeexperimentaldata.Finally,theradiantcharacteristicofmonopolecoveredwithplasma,whichisadispersivemedium,isdisposedbyashiftoperatormethod.
KeyWords:Finite-DifferenceTime—Domain(FDTD)methodAxialsymmetry
SubgriddingmethodConformaltechniquePlasma
创新性声明
本文声明所呈交的论文是我个人在导师指导下进行的研究工作及所取得的研究成果。尽我所知,除了文中特别加以标注和致谢中所罗列的内容以外,论文中不包含其他人已经发表或撰写过的研究或成果;也不包含为获得西安电子科技大学或其它教育机构的学位或证书使用过的材料。与我一同工作的同志对本研究所做的任何贡献均已在论文中做了明确的说明并表示了谢意。
申请学位论文与资料若有不实之处,本人承担一切相关责任。
本人签名:善豳,荔日期:2竺!:!!:耳
关于论文使用授权的说明
本人完全了解西安电子科技大学有关保留和使用学位论文的规定,即:研究生在校攻读学位期间论文工作的知识产权单位属于西安电子科技大学。本人保证毕业离校后,发表论文或使用论文工作成果时署名单位仍然为西安电子科技大学。学校有权保留送交论文的复印件,允许查阅和借阅论文;学校可以公布论文的全部内容或部分内容,可以允许采用影印、缩印或其它复制手段保存论文。(保密的论文在解密后遵循此规定)
本学位论文属于保密,在——年解密后适用本授权书
本人签名:塞捆差
翩虢觐鲤日期:21121丛12
日期:羔!!z:fz!丑
第一章前言!
第一章前言
1.1时域有限差分方法发展和现状
1966年K.S.Yee首次提出了一种电磁场数值计算的新方法——时域有限差分方法fFiniteDifferenceTimeDomain,简称FDTD)。在其著名论文“NumericalSolutionofInitialBoundaryValueProblemsInvolvingMaxwell’SEquationinIsotropicMedia”…中提出对E、H场分量在空间和时间上采取交替抽样的离散方式,应用这种离散方式将带时间变量的麦克斯韦旋度方程转化为一组差分方程,并在时间轴上逐步推进地求解空间电磁场。Yee提出的这种抽样方式后来被称为Yee氏网格。
FDTD方法经过三十多年的发展已成为一种成熟的数值方法【2¨”。近十几年来,随着计算机性能的提高和价格的下降,使FDTD方法更具吸引力,应用范围也越来越广,如计算天线辐射、分析周期性结构、研究微波器件和导行波结构的问题等。三十多年来国内、国际刊物上发表的论文数目可以充分说明这一点口l。前十年中,每年仅有几篇和FDTD有关的论文发表,此后十年中,每年都有十几篇相关的研究论文见诸各种刊物。近十几年来,每年发表的论文几乎按指数增长。
随着FDTD被人们广泛接受和应用,FDTD方法本身也在不断完善和发展,主要表现为:
——FDTD区。对于散射问题,通常在FDTD计算区域中引入总场边界(即连接边界),如图1.1所示。FDTD计算区域划分为总场区和散射场区。这样做的好处是:一、应用惠更斯(Huygens)原理,可以在连接边界处设置入射波,使入射波的加入变得简单易行;二、可以在截断边界(即吸收边界)处设置吸收边界条件,利用有限计算区域就能够模拟开域的电磁过程;三、根据等效原理,应用数据存储边界(即输出边界)处的近区场便可以实现远场的外推计算。对于辐射问题,激励源直接加到辐射天线上,整个FDTD计算区域为辐射场,如图1.2所示,不再区分总场区和散射场区。
——吸收边界条件(截断边界条件)吸收效果提高。为了在有限的计算区域模拟自由空间中的电磁问题,必须在计算区域的边界上设置吸收边界条件,将计算区域截断。吸收边界从开始简单的插值边界,到后来广泛采用的Mur吸收边界,以至近几年发展起来的完全匹配层(PML)吸收边界,其吸收效果越来越好。
凿1t散射计算时FDTD区域的划分
图12辐射计算时FDTD区域的划分
——近一远场变换。FDTD的模拟只能限于有限空间,为了获得计算域以外熬教鬈场,必须詹助等效琢理应翔诗算嚣域疆肉豹遥场散射数据实现诗舞区蠛以外远场的外推。
——潮格削分。在网格剖分方面新技术髂殿用提高了FDTD方法的肖效性和准确性。这些新技术包括:直角坐标系中的亚嘲格技术[8卜[121;直角坐标系中的姻面共形技术…3^[“j,广义正交坐标系中的差分格式和非正交变形网格技术等。这些薪技术匏应震大大辐宽了FDTD方法翁应鼹慈罄。
时域有限差分(FDTD)方法是计算天线辐射的一种有效方法。决定FDTD方法精发斡一个因素是掰取溺格豹尺寸。通常情况下,离散两格尺度矗鞭决于入射波长,即满足dsⅣ10。但在很多情况下。需疆对于感兴趣的部分进行细分,如染采弼均匀两格,列需要对整个计算区域都采用小尺寸网格,这样嚣受到计算机内存及运行时间的限制。鉴于这荦申情况,在直角坐标系孛,我嬲应用妻楚坐撂系中的亚嘲格技术{8j.f1“。即对于感兴趣的部分再进行细分,应用亚网格,而其它区域仍采腰粗网按:丽在二维辅霹称柱坐黎系中,与直焦坐标系中相似,巍可采用楗坐标系中的溉网格技术。这样既能保诞计算的精度,又节约内存及运行时间。键统FI)TD建模在应用台阶近儆模拟簸瑟时会葶i入误差,在煮焦坐标系中,出现
第一章裁富!
tt!ltⅫ目_目_目目_s——一Iutumt
术邮哪’来减少误差,其中有:非正交坐标系FDTD法,曲线坐标系法,回路积分法等。在二维轴对称柱坐标系中,也可以采用=维轴对称柱嫩标系中的曲丽共形技术来减少误差。色散分质是一类复杂介痰,等离子体、生物体、土壤等都具有色敖特性。在直角嫩标系中,色散介质时域递推方法已经商时域卷积方法m】,时域微分方程方法∽{僻},Z变换方法f2j}f。2}。在二维辘对称柱坐标系中,当西与盖的频瑶本构关系可以表示为有理分式函数时,将该本构关系首先过渡刘时域,再转换到时域的离散形式,进而导出从西到五的递推计算公式,即可实现色散介质在二维轴对称犍坐撼系中的时域遴接。褥该方法建来毳}}究等离予体覆盖熬单极子天线豹辗慰特性。
1.2论文的安排及主要贡献
本论文共分先五章。
第二章研究粗细网格尺寸比为2的二维轴对称柱坐标系中的FDTD娅网格技术。绘爨了二维辘对称壤况下耀缨网格尺寸比必2的溅网格技术蛉理论分析,完成了相应程序并进行了天线辐射的计算。本章所介绍的皿网格技术,是通过仿照直角攫檬系中耍隧捂技术靛捶擅方法愀”1来楚理耀缎露格速器。该方法密予特殊的阐捺排布,可以应用于亚网格边界穿越介质的情况。
第三章讨论粗细网格尺寸眈为3的二维辅对称聿主艇标系中的FDTD腰网格技术。给如了二维轴对称情况下粗细网格尺寸比为3的弧网格技术的理论分毫厅,完成了相应程序并进行了天线辐射的计算。本章所介绍的溉网格技术,是通过与直角坐振系中瓣插值方法m瓷织豹方法采处瑗耀细网格边界兹。该方法不需要特殊的耀格排布,但可以应用于粗细网格边界穿越金属介质的情况。
第疆牵FDTD曲蟊共形按术。本章绘出了=维轴对称坐标系中的萄面共形技术的理论分析,完成相应程序并计算了平底圆饿天线和球底两圆锥天线的反射电压。该筵形技术(MCFDTD方法)参照直角坐标系中采用局部紧缩交形网格的共形技术肿州7l,也采用局部紧缩变形阙格。农这些变形网接中,从积分形式麦克颠奄方程出发得到修改后的FDTD递推公式,丽在非变形网格中仍使用一般的递摧公式。越乡},还讨论了变形瘸格癍当满跫静筵形条穆。该技术操作琵较篱荜,计算稳定。该算法可用于模拟二维轴对称柱坐标系中的任意曲面形状的导体,第五章醑究了等离予俸覆盖酌筚极予天线的辐射祷往。本牵先介绍了二维轴对称柱坐标系中常用的一类色散介质的FDTD计算方法。该方法具有可行性。应耀此方法分析了等离子体覆盖的单极子天线的辐射特性。
。,一耋三坠翟些盐些蝥盥譬尘型型鏊尘—I—一T!__目目_目E|自__|s___E目__gE___I●●__Ej____155——
场重合。为了确保由粗网格到细网格的平稳过渡,应用了对时间外插值以及对时间和空间的内插值方法,这两种插值都具有二阶近似。
图2.2亚网格排布及粗细网格边界
图中灰色区域表示粗细网格的过渡区。对柱坐标中二维TM(TE),图
中+表示磁场(电场),既是粗网格的磁场(电场)又是亚网格的磁场
(电场),?表示粗网格的电场(磁场),0表示贬网格的电场(磁场)。2.1.3亚网格算法
亚网格算法的重点应该放在以下两个方面:
?保证整个亚网格区域的计算稳定性;
?将来自于亚网格边界的反射率减小到最低。在网格布局上应用一个过渡区域,细网格相对粗网格移动1/4个粗网格的大小。
研究轴对称柱坐标系中二维TM情形。记粗网格大小为占,:△,:位,时间步为At。:亚网格边长乃是粗网格的l/一,,时间步为△,r=At。肛r。用以”G,),B””2(f'/+1/2),丘”V2(f+l/2,,)表示粗网格的场值,用ⅣⅣ”0,g),勖m+WD,g+]/2),E矿m+W0+V2,9)表示亚网格的场值。各场分量之间关系如下:
.一。。I耋三茎』墅坠垦堡垒兰些望耋』型鐾璧毛22二窑窑垒耋奎。,一:(b)在时间上用Ⅳ耐一+%和Ⅳ∥”+%的插值,然后用空间插值来替换亚网格区域内相应的粗网格点的值Ⅳ。”+%,即
%“,4(i+1/2,J+1/2)=jllH∥“只(f+1,2,,+1/2)+%“%4(i+1/2,川/2)](2-4)fc)为了获得粗细网格边界(图2.2中,=‘,。一1/2行)上格点处因缩进亚网格而未计算的磁场值,先用ⅣⅣ“%和H,“+%的外插来获得N+3/4的场值,然后用空间内插来获得丢失的Ⅳd”+%。即对,=J。一1/2行上粗、细网格磁场重合的点,用时间外插值求得
%吲(f+I/zJo—I,2)=强弼“+I/2,.:o一1/2)一%喇(f+l/2,Jo一1/2)(2—5)
然后对不重合的点用空间内插值求得
qrN十必(f'矗-1/2)=弓IItt一”+%u-1/2,Jo-1/2)+H矿”+%(i+U2,矗-1/2)1(2—6)这样,=.,。一1/2行上小网格的磁场值均已求得。
5)t=N+4/4时,用FDTD来获得Erti*l'E:“1和亚网格的E矿“’,E矿“1
6)回到第~步,循环。
在此方法中插值扮演了一个非常重要的角色,在各种插值中选用了最简单的一种方式,即线性插值。
2.2粗细网格区域边界的非物理反射分析
在粗细网格边界,由于处理方法的不连续性,不可避免地要带来非物理性反射,这种反射是影响亚网格技术稳定性和精度的主要因素之一。一种好的亚网格技术应当使粗细网格边界的非物理反射尽量小。上节所描述的亚网格方法的这种反射非常小。
2.2.1计算区域的设置
为测试粗细网格边界的反射率,如图2f3所示,设计算区域为Q(含150x300个网格),四个角点为(0,一150)、(0,】50)、(150,150)、(150,一150)。作两次计算:第~次取亚网格区蜴,(含40×160个网格)的四个顶点坐标为(30,一80)、(30,80)、(70,80)、(70,一80):第二次取亚网格区Q2,(含90x160个网格)的四个顶点坐标为00,一80)、00,80)、020,80)、020,一80)。激励源设在原点。至少在某个时间步
8轴对称柱坐标FDTD网格割分技术及其在天线辐射中的应悄
内,当9I,右边界的反射波到达接收点P(68,o)时,绕,边界的反射波还未达到接收点。于是用两次计算得到的P点的场值相减,就可得到襁细网格边界的反射波。
图2.3计算襁细网格边界反射时的区域设爱
2.2.2粗细网格边界反射波的提取
对辐对称柱坐标系中豹TM波,高新脉冲没饕覆互:和,o)点,蕤形式为
E:陋“‘”%㈤
:e)=P/。2挖一71其中入射波上限频率为1.5x10”H:,上式中气=40At,At:3.33×10“sec,fo。0.8tm,如图2.4所示。图2.5中两条曲线分别怒亚网格区为奶,和Q2,时接收点P处的波形E,l(p)窝t2(P)。圈2.6曲线为
D(P)=疋。(P卜古:2(P)f2.81代表P点所接收蜀的边界的反射波。
对图2.6曲线备段进纷分析:
a)由O到A点(Timestep。144)是高斯脉冲由辐射源点到Q,右边界,及在该边赛产生黪反射波到达接收点P所零要的惑孵闯步;
b)由A点到B点(Timestep=224)是鸟,右边界的反射波:
c)垂c点(Timestep2344)到D点(Timestep:464)是g,右边界反射波及Q右边界反射波的迭加。
第二章粗细网格尺寸比为2的轴对称FDTD亚网格技术!
蔷—je
山“∑e
山”
10
08
O2
00
Timestep
a)入射高斯脉冲时域波形
f(GHz)
b)入射高斯脉冲频域特性
图2.4入射波的波形6
4
0
0
一o.。一时
塑轴对称柱坐标FDTD网格割分技术及其在天线辐射中的应用
o
o
—
Q-
o
.{5
—20
timestep
图2.5接收点P的场值
Timestep
圈2.6租细网格边界的反射波
婴上各时间段的划分与实际分析相一致,如A点(Timestep=144),经计算
可知该点是Q-,的右边界的反射波到达接收点的时间步:对于C点,可以用妥写
原理计算出Q2,右边界的反射波到达接收点的最短时间步是344。
在(2-8)式基础上,现在可通过时间加窗法提取0点到F点(Timestep=300)
的时域反射波,仅是9,右边界的反射波。
0
5
05
0
们¨∞
¨仙,。1.一d)“
第二章粗细网格尺寸比为2的轴对称FDTD亚网格技术旦
区场E。的时域分布。作为比较,图中也给出了采用均匀细网格占=lmm计算所得的结果。由图可见,二者符合很好。
图2.10球底面圆锥天线极化模型
O,O'lO
0005
时
0.000
—0,005
.O010
图2.1I接收角theta=60。时球底面圆锥天线的远区辐射分布图
2.4本章小结
本章介绍了一种既高效,反射系数又非常低的亚网格技术。数值计算测试结果表明,粗细网格界面反射系数可以达到(-90dB)以下。该方法可以用于亚网
格边界穿越介质界面的情形。采用本文的方法既可以获得精确的场分布,又节约
计算机内存和计算时间。
14
轴对称柱坐标FDTD网格剖分技术及其在天线辐射中的应用
第三章粗细网格尺寸比为3的轴对称FDTD亚网格技术
3.1亚网格技术原理
3.1.1亚网格布局
令粗细网格尺寸比值∥=3。网格布局如图3.1所示。由于粗细网格的尺寸比
是奇数,所以在细网格区的粗网格点均有对应的细网格重合点【12】。
-
-
。
I
Hl?
何
J
r
Ⅳ
J
10
日
J
一
^.
D
。
~F.“
丌
D
F‘
、芷’
。C。
E‘L’
‘‘g。
D
.
《
H岛《
;
I
’
I
!
I!
I
’
I
●
。1l
r
Ⅳ
i≥v
●
一
?一毛
H
j
啊?
r
一
?H占一.
一
V
一
一
一
}
’
●
’
I?I!
I
’
或2
丌
’
^
^
。E3tes占一l
u
一
V
U
一?
D
'
I
!
‘
j
I
!
I
,
e
一
-
●
●
—
V
●
’
’
胃
沁e
1日I与2
J
I
I
!Ⅳ.
?一罩2日
j
。
V
V
V
V
。
8
o
?
I
。
I
!
●
'
I
j
P
P
d
d
^.
p
D
’(‘’:凸’’‘C1
u
;
H
1
}H
』
’
圩
j:
日
』
P
图3.1亚网格排布及租细网格边界
对二维柱坐标系中的TM(TE),E表示粗网格的电(磁)场,或者表示既是粗网格的电
(磁)场又是亚网格的电(磁)场;Ⅳ表示粗网格的磁(电)场,或者表示既是粗网格的磁
(电)场又是亚网格的磁(电)场;e表示亚网格的电(磁)场;一~表示粗细网格边界
在该种布局中,粗细网格在空间无须有一定的偏移,对于柱坐标系中的TM极化,粗网格的磁场和相应的细网格磁场自然重合。同样对于TE极化,粗细网格的电场重合。为了确保由粗网格到细网格的平稳过渡,应用了加权、对对阔外插值及对时间和空间的内插值方法,其中两种内插都具有=阶近似。
第三章粗细网格尺寸比为3的轴对称FDTD亚网格技术竺3.1.2亚网格算法
亚网格算法的重点应该放在以下两个方面:
?保证整个亚网格区域的计算稳定性,次夕卜层的电场进行加权处理。
●将来自于粗细网格边界的反射率(主要是低频部分的)减小到最低。对细网
格区外层的磁场进行加权处理。计算证明这是一个很好的方法。
研究柱坐标系中二维TM情形。记粗网格大小为疋=Ar=△z,时间步为△f“亚网格边长D是粗网格的I/.,(”,=3),时间步为△‘,=△,。加r。用Ⅳp”O+1/2,,+1/2)、Ern+l/20+1/2,,)、Ezn+92(f,S+1/2)来表示粗网格的场值,用^p”(p+1/2,q+1/2)、erm+l/20+1/2,q)、ezra+l/2(p,q+l/2)表示亚网格的场值。各场分量之间关系如下
J孽pn+l/20+1/2,,+l/2)=百rn+l/20+1/2,J+1/2)
+FAI陋:”O+1,,+1/2)一E;”O,』+l/2)J
+FA2?p,”O+1/2,,)一E,”(f+1/2,J+1)J
Er”+1冬+i/2,,)=Er“+1(f十i/2,,)f3-1)
+FA2?睁∥O+1/2,,一1/2)一拂…I/O+l/2,,+l/2)l
Ezn+|(f,,+1/2)=Ez”1e,,+l/2)
+FAlh拂+i120+1/2,,+l/2)一‘一{口”20一1/2,,+l/2)I
其中
耻繁‘掰雌鳓姗抛删l-警p:,FA2:坐”…
FAl,FA2与离散网格尺寸有关。式(3-1)、(3-2)不仅适用于粗网格,同样适用于亚网格。当用于亚网格时,相应的匕,E,,Ez对应于^,,q,巳,瞑对应于拶,,At。对应于Atr。由上式可知,H,的计算需要邻近四个E,和Ez的格点值,而E,,E.的计算则需要邻近的二个爿0格点值。对于粗网格和亚网格我们可以用通常的FDTD递推公式进行计算。而在粗细网格边界上的电场无法计算。
对于柱坐标系中TM波,已知E”,∥-1和Ⅳ”%,亚网格方法的具体计算步骤如下:
(1)计算整个计算区域(包含粗网格区和细网格区)的Ⅳ”必和E“。
(2)用普通FDTD迭代来计算”+l/6时刻的细网格区内的磁场^。一+%。
(3)计算”+2/6时刻细网格区内的电场,即
(a)使用普通的FDTD迭代来计算亚网格区内部的。,一+%,e一+%:z
嘞唰哥粗纽网格边界士与粗丹椿电场榀垂会自争点冉争晦%‘坶酾钥掘崮涝。
I__—_——?_——_16轴对称柱坐标FDTD网格剖分技术及其在天线辐射中的应用
ii——_————●—●—●■_——__-—_———_IIII——
(b)对于粗细网格边界上与粗网格电场相重合的点的电场,计算所用公式为
E…=E“+Av+Bv2/2(3-3)式中
4=k“一E.-I/2,B:E““‘E“~一2E”(3—4)其中v=2/6:
(c)对于粗细网格边界上与粗网格的电场不相重合的点的郎“~,e:”%,使用线性内插来求出,如图3.1所示,具体公式为
%2;鲁+÷0(3-5)
…7
qz=号蜀+{也
式中的用到的E。、E。、E:均由(b)经时间内插得到。粗细网格边界上其它的与粗网格电场不重合的其它点的电场值可类推;
(d)对于细网格次外层上的电场值作加权处理,在图3.1中具体为
%=o.95e埘+0.05+兰≤鱼(3.6)
‘
其中次外层的岛,,及内层的P32由(a)步计算得,外层的岛由(a)或(b)步得到。细网格次外层的其它电场值可类似乎导到。
(4)计算”+3/6时刻细网格区内部的磁场,按照以下步骤:
(a)用普通FDTD迭代来计算细网格区的^。“邶;
(b)对细网格区磁场的最外层进行加权处理,如图3.1所示对粗细网格的磁场日I、h。计算公式分别为
Ⅳl2o?8蜀c+O?2啊7
n一∞
、~,
hl=0.2Hlf+O.8hl,
其中Hk表示由第(1)步所算出的n+l/2时刻的该点的粗网格场值,岛,表示由第(4)步所算出的该点的亚网格场值,啊及h。分别表示权重后的该点的粗、细网格场值。
(5)使用由(4)得到的加权后的粗网格磁场值日,“V2及历史上的粗网格电场值E一来重新计算粗细网格边界上的粗网格的电场值点”-。
(6)计算一+4/6时刻细网格区的电场。与第(3)步所述步骤相似。注意这里取(3b)步中v=4/6。
(7)用普通FDTD迭代来计算细网格区的hn+S]6。
(8)计算”+1时刻的细网格区的电场ern+l,%”1。重复第(3)步的步骤。注意这里取(3b)步中v=I。
回到第(1)步进行下一个循环。