多媒质区域中电磁波传播的TDCB模拟.pdfVIP

1 一^f 1{t： 多媒质区域中电磁波传播的TDCB模拟 褚庆昕 胡福刚 商安电子科技大学710071 eduCn email：qxchu@xidian 摘要：本文将时域特征波法j{j丁模拟多媒质区域的电磁波传播。’o时域有限筹分法 相比，时域特征波法在处理媒质交界面和截断边界上具有优越性。 关键词：时域特征波法，多媒质区域，电磁波传播。 I．引言 多媒质区域的电磁波传播问题有着广泛的实际背景，例如：介质波导，光纤，微 带电路，复杂地层等中的电磁波传播。由r边界复杂、区域媒质不均匀，模拟多媒质 区域中任意波形的电磁波传播具有难度。解析方法不仅复杂．两置只能模拟一些简单 此在媒质交界面以厦计算域截断边界上是不成立的。通常，在处理媒质交界面时采用 平均介电常数和磁导率近似，这在某种程度上影响r计算精度。对，截断边界则需要 寻找新的方程及差分格式。时域特征波法(TDCB){3．4]平q用系数矩阵的特征分解将三 数步时问算法依次沿各空闯坐际方向采崩迎风差分格式求解。由于迎风差分格式是按 照一维方程的特征波方向构造的，因丽在媒质交界两上可以得到严格的差分公式。在 截断边界上只需令相应的反射波为零就可以得到吸收边界条件，无需寻我新方程。本 文将时域特征波法用于模拟多媒质区域中的电磁波传插。作为例子，模拟了电磁波垂 直入射到媒质交界面的问题。与精确嬲和FD’FD模拟结果比较表明，时域特征波法在 媒质交界面上不会引入新的误著。 II．时域特征波算法 在无源、各向同性媒质中，Maxwell方程 V。厅：s望 优一 (1) vx￡：一Ⅱ一cgH e} 螽 区角 坐 际 系 kP 以 丧 ～小 ， 鼍警 十 4 十 4 A = 盟＆ 塑砂 盟瑟 一49一 式中，衫=tE，B，t，日，，Ⅳ，，也-7 4=[一警≥，三协‘蠢，]，，=t，2，。 000} 00—1 B。=0oI l，B2 o00{，B= 0．】0 1 00 湖 米崩分数步舁绣，l2)趵群口3以由备^瞿杯鞫刀圆嗣待J仳设群狱得【4j 口”‘==LxL，上：三：t，三，疗” (3) 式中，算子上。，￡，，￡：分别袤示沿x，Y，z方向的～维方程解： 上式中的备矢量方程都是由F列形式的标量方程组成的 |盟一三鹎：o 。新 {曼 (5)… l堡一三强：o 【df ∥讹 式中，f，．，，k分别代表坐标x，y，z中的一个。于是，求解(1)的最基本问题是如何求 解(5)。利用矩阵的特征分解，(5)可以分解为两个独立的特征波方程： 型一 V--0：二蟹宴 幅¨(i。) 百一p