第十一讲表面阻抗边界条件.doc

表面阻抗边界条件和近-远场变换 第9讲 褚庆昕 PAGE 9 第9讲 表面阻抗边界条件和近-远场变换 9.1 表面阻抗边界条件 早在40年代，表面阻抗边界条件（SIBC）就已被应用于电磁波相互作用问题。在FDTD模拟中，SIBC可以用于计算有耗介质或导电结构外部的场，而不需要模拟结构内部。本讲介绍模拟SIBC的FDTD方法。 9.1.1 Maloney-Smith法 表面阻抗本质上是频域概念，在自由空间和导电媒质的交界面上，切向电场与表面电流满足: (9-1) 式中，为复值表面阻抗，是表面的单位法向矢量。 变换到时域 (9-2) 考虑在自由空间中(媒质1)时谐平面波以关于有耗半空间(媒质2)边界外法向成角度的方向传播。有耗媒质具有本构关系，，，损耗角正切。假设两介质的相对折射率足够大，即 (9-3) 这时对于自由空间任何入射角， 在有耗介质半空间的折射角都接近0，则表面阻抗具有简单表达式： (9-4) 上式对任何入射平面波角度都成立。 根据(9-4)的Laplace反演，得到适用于(9-2)的等效表面阻抗冲击响应： (9-5) 式中,和分别是第一类0阶和一阶修正Bassel函数。是Heaviside单位阶跃函数，是Dirac单位冲击函数。 将(9-5)代入(9-2)，得 (9-6) 设在时间上是分段线性的，离散上式中的积分， (9-7) 式中， 显然，(9-7)的卷积的数值估计涉及一个非常大的计算机负担，因为交界上所有点在所有时间值都必须存在计算机中，而且在每一时间步都要进行处理。 为此，采用Prony方法，将展开为指数和形式。 (9-8) 式中，，和通过在最小二乘意义上误差最小来确定。 将(9-8)代入(9-7)，得 (9-9) 可以证明，满足递推关系 (9-10) 所以，(9-6)简化为 (9-11) 上式是非常有效的，在每一边界点上所需的计算机存储只是步求和。而且每一时间步求和只需进行两次乘法运算。另外，计算时域阻抗函数(9-5)及其指数级数近似(9-8)所需的时间与FDTD本身的负担相比可以忽略。 在FDTD网格中执行这种SIBC的困难是，需要知道在结构边界上与电场相同时-空点上的。这与Yee网格中电场磁场分量相互分离是冲突的。通常，在边界上所需的可由离边界远、时间提前的平行分量近似。 9.1.2 Beggs法 考虑自由空间和有耗介质半空间的平面交界面位于处,单色场有复分量沿交界面在点处，(9-1)为 (9-12) 式中，表面阻抗通过严格求解电磁场边值问题而得。 如果有耗介质是良导体，Leontovich（里奥托维奇）阻抗边界条件提供了(9-12)一阶近似 (9-13) 式中, (9-14) 由(9-13),(9-14)得 (9-15) 式中，，。 假设，在下面的处理中保持常数，即 (9-16) 由Laplace反演，得时域条件 (9-17) 上式需要在相同空-时点的电场和磁场。处理中，用近似，于是, (9-18) 考虑Faraday定理及围线积分，有 (9-19) 将(9-18) 代入(9-19)，得 于是，利用导数的差分近似及时间平均算法，最后得 (9-20) 9.2 近—远场变换 在计算天线辐射特性和雷达双站散射截面时，常要用到远区电磁场分布。而FDTD限于计算机容量和计算时间，只可能计算近区场分布。为此需要利用近一远场变换从近场求出远场。 9.2.1 基本公式 首先在FDTD计算域中（包括了所有结构）建立一虚拟立方体。设在该立方体表面上用FDTD计算得到的场为。 根据等效原理，当立方体区域中为零场时立方体外的场是表面上的等效电磁流 (9-21) 产生的。利用辅助位函数法，在频域位函数 (9-22) 式中，为立方体表面积。 ， 电磁场为 (9-23) 在远区，球坐标系下各场分量为