负折射率媒质模型的时域FDTD分析研究.pdf

200 5’全国微波毫米波会议论文集 负折射率媒质模型的时域FDTD分析 林振梁昌洪 linzhen@people．com．ctl (西安电子科技大学天线与微波技术国家重点实验室 陕西西安710071) 摘要负折射率媒质是一种人造媒质，在某一微波波段内其介电常数和磁导率同时为负值。本文对TE模电磁 场满足的麦克斯韦方程组进行数值处理，解决了直接利用FDTD方法在Yee’s网格下计算负折射率媒质内部 场分布时的数值发散问题。通过数值算例，验证了方法的可行性，模拟了负折射率媒质内部及周围的电磁 场分布。 关键词负折射率媒质，TE模，FDTD方法 1．引言 近年来，左手性媒质的研究正在成为科学研究的热点之一。电磁学中，介电常数占和磁 导率F是描述均匀媒质中电磁场性质最基本的两个物理量。对于自然界所有已知的材料而言， 介电常数s和磁导率F都是非负的常数。由麦克斯韦方程组可知，在P和∥都为正值的物质 handed 中，电场、磁场和波矢三者构成右手关系，我们称这样的物质为右手性媒质(right handed medium，LHM)的物理思想，该理论认为平面电磁波在一个同时具有负介电常数和负导磁率 的媒质中传播时，将射向与Snell定律不同的方向，即发生了微波异常传播的现象。此外该 媒质的折射率n取负数，即当媒质同时具有￡，o，∥，o时，媒质的折射率为F／：～、／占，Ⅳ一 直到1999年，由于Pendry等人预言利用某种特定的人造复合材料结合金属线阵列可以制作出 在某一频率区间满足介电常数s和磁导率∥同时为负的物质，即左手性媒质或负媒质晗’31。 能用计算机模拟出左手媒质奇特的传播特性，将大大有利于人们对它的理解和研究工 作。时域有限差分法(FDTD)是直接由麦克斯韦方程组导出的数值计算方法。它没有对计算 区域中反射系数的正负以及波的传播方向做任何假设，因此用作计算左手媒质的传播特性非 常实用。 但是，直接简单的在传统FDTD计算式中设左手媒质区域的介电常数s和磁导率∥为负 数，将会导致计算的不稳定，模拟数值发散。本文对传统的FDTD计算公式进行相应的数值 处理，解决了数值计算中的不稳定性，模拟出左手媒质区域TE模电磁场的分布情况。 2．左手媒质中电磁场的FDTD计算式 由于在利用FDTD方法数值计算负媒质区域电磁场传播时，在材料边界处应该满足匹配 条件。对于TE模式的电磁波，其满足的麦克斯韦方程为： 占—i。，2{一√， 占冬：譬一以 占—o=一—一三一√．占堡：一堡一，． jjaHz—aExOEy “—o=——一——一 采用Yee’S的标准网格，在媒质分界面上，场分量计算中需要分别取边界两边媒质的介 电常数￡和磁导率Ⅳ的平均值作为该点的介电常数和磁导率。以ex分量为例， 1467 2005’全国微波毫米波会议论文集 训Ein+lz√2 E川n亿，一丧慨‰一H川in+亿l／2川：)(2) 其中，上标n为时问步，下标i，j为空INx，y方向坐标，占。。为媒质界面两边占l和占2的平均 值(占。。=(占l+占2)／2)。 这样，当分界面一边为右手媒质，一边为左手媒质，且满足以=一,a2，q=一岛时，将会 导致界面上的占。为o，使计算数值发散。对于E。和致有同样的结果。即使‰不为O(如 率∥为负数，仍然会出现数值发散现象。 为了避免这样的数值发散，采用2DTLMH山1方法，间接地来设置负媒质区域的负的介电