时域有限差分法(姚伟)讲解.doc

伊犁师范学院硕士研究生 ————期末考核 科目：电磁波有限时域差分方法 姓名：姚伟 学号：1076411203009 学院：电子与信息工程学院 专业：无线电物理 时域有限差分法 选题背景 在多种可用的数值方法中，时域有限差分法(FDTD)是一种新近发展起来的可选方法。 2 原理分析 2.1 FDTD的Yee元胞 E,H场分量取样节点在空间和时间上采取交替排布，利用电生磁，磁生电的原理 图1 Yee模型 如图1所示，Yee单元有以下特点[2]： 1）E与H分量在空间交叉放置，相互垂直；每一坐标平面上的E分量四周由H分量环绕，H分量的四周由E分量环绕；场分量均与坐标轴方向一致。 2）每一个Yee元胞有8个节点，12条棱边，6个面。棱边上电场分量近似相等，用棱边的中心节点表示，平面上的磁场分量近似相等，用面的中心节点表示。 3）每一场分量自身相距一个空间步长，E和H相距半个空间步长 4）每一场分量自身相距一个时间步长，E和H相距半个时间步长，电场取n时刻的值，磁场取n+0.5时刻的值；即：电场n时刻的值由n-1时刻的值得到，磁场n+0.5时刻的值由n-0.5时刻的值得到；电场n时刻的旋度对应n+0.5时刻的磁场值，磁场n+0.5时刻的旋度对应(n+0.5)+0.5时刻的电场值，逐步外推。 5）3个空间方向上的时间步长相等， 以保证均匀介质中场量的空间变量与时间变量完全对称。应用这种离散方式，将含时间变量的Maxwell方程转化为一组差分方程，并在时间轴上逐步推进地求解空间电磁场。由电磁问题的初值和边界条件，就可以逐步推进地求解以后各时刻空间电磁场分布。 2.2 Maxwell方程FDTD的差分格式 麦克斯韦第一、二方程 （1） 式中，时电流密度，反映电损耗，是磁流密度，单位，反映磁损耗。主要与上式对应。各向同性介质中的本构关系： （2） 其中是磁阻率，计算磁损耗的。 以为变量，在直角坐标中，展开麦克斯韦第一、二方程，分别为 （3） （4） 令代表在直角坐标中的任何一个分量，离散符号取为 （5） 关于时间和空间的一阶偏导数取中心差分近似为 （6） 可以看出，每一节点上沿某一方向场分量的一阶偏微分可以用在该方向上相邻两点的一阶中心差商来描述，将式(1)用一阶中心差商方程取代，整理后便得到一阶差分方程，它具有二阶精度[3]。 Yee元胞如图1所示，规定为 1）剖分节点与场分量所在棱边中点不同，场分量的位置，即节点是Yee元胞节点的相对位置，不需要单独编码； 2）当空间存在媒质分界面时，场量自动满足场的连续性条件， 电磁分量的取样方式不仅符合法拉第电磁感应定律和安培环路定律的自然结构，也符合麦克斯韦方程的差分计算。其次，时间步长可以取为电磁波传播一个空间步长所需时间的一半，因此与在时间顺序上交替抽样，时间间隔相差半个时间步长。 2.3 一维问题 均匀平面波（TEM波）是一维问题，设电磁波沿z轴方向传播，则，场量和介质参数均与x，y无关，即，麦克斯韦方程为 （7） 和 （8） 旋转坐标轴后可以只保留一组公式[4]，设保留（7） Yee元胞如图2所示 图2 一维Yee元胞 差分格式为 （9） （10） 如果介质无损耗，则 2.4 二维问题 三维通常是散射问题，二维是TE、TM波问题，一维是TEM波问题。 在二维场中，所有物理量与Z坐标无关，既。于是在TE和TM波的表达式分别为 TE波（） （11） TM波（） （12） 图3分别给出了TM波和TE波的Yee元胞图 图3 TM波的Yee元胞 图4 TE波的Yee元胞 对于TE波，只要令，在上，不随z变化，m中去掉k即可得到： 式中：