时域有限差分方法在层状介质空间中的处理技术毕设外文翻译译文1Word格式文档下载.doc
《时域有限差分方法在层状介质空间中的处理技术毕设外文翻译译文1Word格式文档下载.doc》由会员分享,可在线阅读,更多相关《时域有限差分方法在层状介质空间中的处理技术毕设外文翻译译文1Word格式文档下载.doc(19页珍藏版)》请在冰豆网上搜索。
![时域有限差分方法在层状介质空间中的处理技术毕设外文翻译译文1Word格式文档下载.doc](https://file1.bdocx.com/fileroot1/2022-10/30/0a47ad38-8523-42c7-9d21-6702dc12f903/0a47ad38-8523-42c7-9d21-6702dc12f9031.gif)
在这一章中,我们提出,演示了在第二章(近场到远场的变换)和第三章(总场/散射场边界条件)在层状介质有限时域差分分析方法的技术应用实例。
在前面的章节中,我们考虑的例子,涉及的方法是各自章的主题。
例如,我们只考虑从放置在半空间和第2章中接地板结构的远场辐射。
同样,我们只研究了第三章中的TF/SF(总场/散射场)的边界性能(即,抑制入射平面波,散射波的透射度)在这一章中,我们使用第2-3章中发展的两种方法;
也就是说,利用TF/SF的边界和近场到远场变换这两个例子。
4.1层状介质的散射
在本节中,我们考虑两个问题,当入射场为平面波放置的对象,涉及放置在多层介质中理想导电物体的散射。
引入到FDTD的网格入射平面波的过程,有必要使用TF/SF多层介质边界,而从散射物体远区场的散射计算是由NFFFT完成。
使用这种方法获得成功,在图6的NFFFT的字段中必须只包括散射场。
这确实是的TF/SF的边界技术的基本前提:
只允许散射场退出边界。
在3.4节中提到散射波TF/SF的边界透明度。
在多层介质中散射问题的一般配置如图33所示。
在多层介质中TF/SF的边界包含的缺陷(或散射),记为A集体应用于现场改正的TF/SF的边界包含一个平面波,用TF/SF的边界左上角箭头表示。
散射结构创建了一个散射波,用TF/SF的边界右上角箭头表示。
这波完全透明退出边界,并达到近场到远场的变换,如图33用S表示。
由于多层介质中的散射结构,在S中切向电场Et和磁场Ht属于散射波。
以下内容是这里的重点:
在图33方框S区域中NFFFT,不包括入射平面波。
这是因为入射平面波几乎完全包含在TF/SF的边界里面,即NFFFT的方框S里面。
图33:
在多层介质中平面电磁波散射问题所采用的配置
作为我们第一个例子,我们重新考虑在第3.4节的情况。
这个例子如几何图31所示。
在这个例子中,一个PEC直角棱镜的尺寸是5cm×
0.4cm×
5cm,被放置在一个1厘米深的半空间介质(或地面)中,一个在半空间中的高斯包络平面波是入射的。
在图32(b)中PEC棱镜的散射波是可见的,但没有讨论电场的辐射。
我们现在量化几何辐射电场采用第2章中研究的NFFFT。
在图34和图35中,从PEC棱镜辐射电场(被1/r规一化)的散射,在xz平面和YZ平面上以各自不同的角度显示。
这两个数字,第一子图用ψ表示该部分,第二子图用φ表示辐射电场。
时间被τL=L/v归一化,其中L=5厘米是最大尺寸的PEC棱镜。
是在介质中的地面传播速度。
电场的辐射波形在时间轴的正半轴,点远离该轴(逆时针方向)。
图34(a),图35(a)和图35(b)中的波形绘制在同一尺度中;
而图34(b)被放大了10倍。
入射平面波的高斯波形为参考原点。
在图34(a)–图35(b)振幅ψ,电场辐射分量φ的相对差异可以通过简单的参数定性解释。
例如,在图34(b)中,我们观察到在xz平面φ子源的辐射电场具有最小的平均幅度。
这完全是因为在PEC棱镜的表面电流大多在与xz平面平行的平面上,这是一个无效的平面,由电流辐射的电场分量φ组成。
在xz平面的辐射电场φ分量是由平行于YZ面和xy平面中的表面电流产生的。
与平行于xz平面的两个面相比,这个面又较小的面积。
注意图35(a)中,ψ接近90°
时ψ的辐射电场的振幅将减小。
这是因为在平行于xz平面的表面上的电流ψ分量,Z轴是零位轴。
正如我们的第二个例子,我们考虑如图36所示的几何形状。
在这几何图中,一个直角棱镜PEC再次放置在两层介质中,类似前面图31所示的例子。
PEC的棱镜的尺寸和前面的一样:
长=5厘米,宽=0.4厘米,高=4厘米。
在这里,上半部空间比下半部空间有较高的相对介电常数:
r0=2.5>
r1=2。
这将允许非均匀平面波为在3.3.3章节的例子做创造。
如图24,图26和图30中,入射平面波的传播方向在xz平面和z轴角度θ=70°
上。
入射平面波的极化如图24(a)所示。
模拟时域有限差分和入射平面波电场波形的参数与图30中使用的相同。
正弦调制高斯波形的参考点为原点。
图34:
平面波从埋在地下介质中的PEC矩形棱镜的散射(a)ψ部分,在xz平面的散射电场(b)φ部分(由10个扩增)。
图35:
平面波从埋在地下介质中的PEC矩形棱镜的散射(a)ψ部分,在yz平面的散射电场(b)φ部分。
图36:
非均匀平面波散射问题的几何形状。
(a)截面在xz平面。
(b)截面在yz平面。
没有PEC的散射,入射平面波完全在分界面反射,一个在标准的分界面方向衰减的非均匀平面波,在较低的半空间形成。
该例子是先前在第3.3.3节图30(a)所示。
在这里,我们探讨了放置在下半部分的空间PEC散射体的效果,如图36。
在图37中,散射过程说明了FDTD几何模拟一个场快照级数。
在图37(a)中,非均匀平面波尚未与下面的PEC的棱镜界面接触。
在图37(b)中,非均匀平面波穿过PEC棱镜,散射波就开始出现。
在图37(c)中,非均匀平面波通过PEC棱镜,球面散射波在分界面的上方和下方清晰可见。
在图37(d)中,散射波的足迹远离PEC棱镜,以及大部分散射能量向下传播作为一种在较低半空间的非均匀平面波的延伸。
在图30(b)中,平面波和非均匀散射波之间的区域看起来像平面波和非均匀平面波传播之间的界面。
散射过程如图37所示。
总结以上,本质上是一个非均匀平面波到传播波的转换。
在3.3.3节中,这现象是指抑制内部的全反射环境(红外),其中一个非均匀平面波产生了一个在较低层传播的波。
图30(b)中说明这种现象。
正如红外的例子,创作从一个携带能量的非均匀平面波到远离没有能量承载的界面,带来了明显的矛盾。
在3.3.3章节中,已经给出了这个问题的答案:
两个不均匀波携带一个方向的净能量,在相反的方向衰减。
利用第2章讨论的NFFFT,我们可以量化辐射电场,如在前面的例子。
在图38中,辐射电场ψ分量显示在xz平面上的不同的角度。
在较低半空间的辐射电场(π<
ψ<
2π),使用NFFFT扩展,在2.3节得到解释。
时间轴被τL=L/v归一化,其中L=5厘米是PEC棱镜的最大尺寸,和是在较低半空间的传播速度。
电场的辐射波形对ψ点远离的时间轴(反时针方向)是确定的。
在上半空间中的短线的ψ=153.43°
和ψ=26.56°
是指在上半空间中定义的辐射电场角度范围的边界。
这个范围包含所有可能的方向,即一个平面波从下半空间折射传播到上半部分空间。
图37:
非均匀平面波从PEC直角棱镜转换到散射波的传播。
在xz平面上的电场强度在不同的瞬间显示(a)-(d)。
图38:
非均匀平面波从在两个层状介质的中的PEC直角棱的镜散射。
电场散射ψ分量在xz平面上显示不同的角度。
附录A
一种在多层介质中的时域辐射电场
在这个附录中,我们将获得2.2章节中在上半空间辐射电场和时域表面电流Jt,Mt之间的关系,如图图6。
研究结果将提交到一个更一般的环境下,即一个任意容量的电流分布J、M是存在的。
专业化的表面电流Jt,Mt是那么微不足道。
让我们首先考虑图6对应的频域。
其中局部源J(r'
),M(r'
)在无损多层介质中辐射。
设z=z0的任何平面上的局部源。
在随后的分析中,我们将使用这平面上的平面波频谱,推导出辐射电场。
电场的二维傅立叶变换ρ=(x,y)→kρ=(kx,ky)(或频谱分析)在z=z0的平面记为:
如果没有z分量,(239)通常被称为该电场的平面波频谱。
利用图6中的几何不变性,在XY平面,(239)可表示频谱叠加的积分:
其中J(kρ,z),M(kρ,z)是源电流J(r'
),M(r'
)的频谱分解。
在(240),G-EJ,EM(kρ,z|z'
)记为光谱格林函数,在附录B,(308)–(309)中给出。
在(240)中取代二维傅立叶积分Ĵ(kρ,z'
)和M(kρ,z'
),在其中我们得到体积V包围的所有来源。
推导的关键点是,上述平面z=z0的电场可以直接表示为平面波的频谱范围(239)。
为此,平面波频谱(239)被插入到傅立叶重积分当中[1,eq.(3.83a)]:
其中ρ=(x2+y2)1/2,及kz=k02−kρ2)1/2,Im(kz)<
0。
由此产生的场(242)在z>
z0区域满足麦克斯韦方程组和平面z=z0的边界条件,因此是解决该电场唯一的方法。
但是,注意这一构想需要一些详细的分析。
由于整个空间是无损的,多层介质可以支持本征模而不是绝对可积(例如表面波以ρ−1/2衰减,平面波以r−1衰减),因此,实数kx,ky传统的二维傅里叶变换是不存在的。
出于这个原因,在(242)中沿着实数kρ积分是无效的。
为了规避这一问题并且不放弃传统的二维傅立叶变换,通常认为,在最上面的半空间有一个微乎其微的损失,使多层介质在各个方向的所有本征模指数衰减。
这是为了确保绝对可积性和传统的二维傅里叶变换对存在。
在(242)中r→∞时,最上面半空间的辐射电场可以得到渐近积分评估。
将(241)插入到(242)和评估V积分,我们得到不同阶的贝塞尔函数。
剩下的半无限kρ积分变换成一个使用汉克尔函数的无穷积分[99]。
通过路径积分变换到经过鞍点kρ=k0sinθ的最陡降路径,最速下降法[13,99]可以应用到无穷积分kρ当中。
对于足够大的k0r,这个贡献来自于在各个方向上鞍点附近占主导地位的领域。
因此,我们忽视了来自kρ平面(引起面波)奇点的贡献,并只考虑鞍点的贡献。
在这种情况下,没有必要详细的进行以上繁琐的步骤;
我们不是简单套用固定相来直接讨论双向积分(242)[1,3.6.2章节]。
用给出的这种方法获得电场的辐射:
其中k0=ω/c,α=k0cosФsinθ,β=ksinφsinθ.
以(243)和(244)–(245)和(241)的逆傅立叶时间变换,我们得到
其中,tr=t+(x'
cosφsinθ+y'
sinφsinθ)/c是横向延迟时间,J(r'
t),M(r'
t)是时域源电流,且GEJ,EMθ,φ(θ,φ,t|z'
)各自的逆傅立叶时间变换
在附录B中(308)–(309)取代kρ=(kx,ky)=(α,β),在(248)–(249)中应用点乘,最后逆傅立叶时间变换,我们得到以下GEJ,EMθ,φ(θ,φ,t|z'
)的表达式:
在这里,ε'
r表示源点的相对介电常数,Z0=1/Y0=是自由空间的波阻抗。
Vpv(t|z'
),Vpi(t|z'
)in(250)–(253)的函数给出:
其中Vpv(z0,t|z'
),Vpi(z0,t|z'
)是源坐标(z,t)=(z'
0)的时域TL格林函数和观察坐标(z,t)=(z0,t)。
这些TL格林函数是附录B(285)–(288)中传输线方程的基本解,替代kρ=(kx,ky)=(α,β)后逆时域: