时域有限差分法论文

《时域有限差分法论文》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、时域有限差分法1选题背景在多种可用的数值方法中，时域有限差分法(FDTD)是一种新近发展起来的可选方法。1966年，K.S.Yee首次提出电磁场数值计算的新方法—时域有限差分法(FiniteDifference-TimeDomain，简称FDTD)。经历了二十年的发展FDTD法才逐渐走向成熟。上世纪80年代后期以来FDTD法进入了一个新的发展阶段，即由成熟转为被广泛接受和应用的阶段。FDTD法是解决复杂问题的有效方法之一，是一种直接基于时域电磁场微分方程的数值算法，它直接在时域将Maxwell旋度方程用

2、二阶精度的中心差分近似，从而将时域微分方程的求解转换为差分方程的迭代求解。是电磁场和电磁波运动规律和运动过程的计算机模拟。原则上可以求解任意形式的电磁场和电磁波的技术和工程问题，并且对计算机内存容量要求较低、计算速度较快、尤其适用于并行算法。现在FDTD法己被广泛应用于天线的分析与设计、目标电磁散射、电磁兼容、微波电路和光路时域分析、生物电磁剂量学、瞬态电磁场研究等多个领域[1]。经过了近四十年的发展，FDTD法在计算方法和应用上取得了大量成果。近几年来，讨论FDTD法的深入发展和实际应用的文章几乎按指

3、数增长，目前FDTD法的许多重要问题得到了很好的解决，已经发展成为一种成熟的数值计算方法。随着计算机数据处理性能的快速提高和计算机价格的下降，使得FDTD法的应用范围越来越广，而FDTD法本身在应用中又有新的发展.2原理分析2.1FDTD的Yee元胞E,H场分量取样节点在空间和时间上采取交替排布，利用电生磁，磁生电的原理图1Yee模型如图1所示，Yee单元有以下特点[2]：1）E与H分量在空间交叉放置，相互垂直；每一坐标平面上的E分量四周由H分量环绕，H分量的四周由E分量环绕；场分量均与坐标轴方向一致。

4、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、推。5）3个空间方向上的时间步长相等，以保证均匀介质中场量的空间变量与时间变量完全对称。应用这种离散方式，将含时间变量的Maxwell方程转化为一组差分方程，并在时间轴上逐步推进地求解空间电磁场。由电磁问题的初值和边界条件，就可以逐步推进地求解以后各时刻空间电磁场分布。2.2Maxwell方程FDTD的差分格式麦克斯韦第一、二方程（1）式中，时电流密度，反映电损耗，是磁流密度，单位，反映磁损耗。主要与上式对应。各向同性介质中的本构关系：（2）其中是磁阻率，计算磁损耗的。以为变量，在直角坐标中，展开麦克斯

6、韦第一、二方程，分别为（3）（4）令代表在直角坐标中的任何一个分量，离散符号取为（5）关于时间和空间的一阶偏导数取中心差分近似为（6）可以看出，每一节点上沿某一方向场分量的一阶偏微分可以用在该方向上相邻两点的一阶中心差商来描述，将式(1)用一阶中心差商方程取代，整理后便得到一阶差分方程，它具有二阶精度[3]。Yee元胞如图1所示，规定为1）剖分节点与场分量所在棱边中点不同，场分量的位置，即节点是Yee元胞节点的相对位置，不需要单独编码；2）当空间存在媒质分界面时，场量自动满足场的连续性条件，电磁分量的取

7、样方式不仅符合法拉第电磁感应定律和安培环路定律的自然结构，也符合麦克斯韦方程的差分计算。其次，时间步长可以取为电磁波传播一个空间步长所需时间的一半，因此与在时间顺序上交替抽样，时间间隔相差半个时间步长。2.3一维问题均匀平面波（TEM波）是一维问题，设电磁波沿z轴方向传播，则，场量和介质参数均与x，y无关，即，麦克斯韦方程为（7）和（8）旋转坐标轴后可以只保留一组公式[4]，设保留（7）Yee元胞如图2所示图2一维Yee元胞差分格式为（9）（10）如果介质无损耗，则2.4二维问题三维通常是散射问题，二维

8、是TE、TM波问题，一维是TEM波问题。在二维场中，所有物理量与Z坐标无关，既。于是在TE和TM波的表达式分别为TE波（）（11）TM波（）（12）图3分别给出了TM波和TE波的Yee元胞图图3TM波的Yee元胞图4TE波的Yee元胞对于TE波，只要令，在上，不随z变化，m中去掉k即可得到：式中：（13）式中：（14）式中，（15）对TM波，只要令，在上，不随z变化，m中去掉k，即可得到：式中，（16）式中，（17）式中：（18）为了编写统