时域有限差分方法林志立研讨.pptVIP

FDTD的应用实例 模拟线光源进入平板玻璃光波导的过程 玻璃平板波导 线光源 光在波导中的传播 光开始进入光波导 光波长:　0.8微米； 平波厚度：2.4微米； 折射率: n=1.6 运行程序: FDTD_Slab.m 参考书目： [1] A. Taflove, and S. C. Hagness, Computational Electrodynamics: The Finite- Difference Time-Domain Method, 3 ed. Artech House: Boston, MA, 2005. [2] K. S. Kunz and R. J. Luebbers, The Finite Difference Time Domain Method for Electromagnetics, CRC Press: Boca Raton, FL, 1993. [3] D. M. Sullivan, Electromagnetic Simulation Using the FDTD Method, New York: IEEE Press, 2000. [4] A. Elsherbeni and V. Demir, The Finite Difference Time Domain Method for Electromagnetics: With MATLAB Simulations, SciTech Publishing, 2009. [5] 王长清，祝西里，电磁场计算中的时域有限差分方法，北京：北京大学出版社，1994. [6] 高本庆，时域有限差分法，北京：国防工业出版社，1995. [7] 葛德彪，闫玉波，电磁波时域有限差分方法，第二版，西安：西安电子科技大学出版社，2005. 英文书： 中文书： Beihang University Beihang University 时域有限差分方法 林志立 zllin2008@ 北航仪器光电学院光电工程系 The Finite-Difference Time-Domain Method for Electromagnetics 计算电磁学中的 电磁学中几种重要的数值计算方法 有限差分法 Finite Difference Method 静电场、静磁场的有限差分法；　 时域行波的电磁场的时域有限差分法； 有限元法 (Finite Element Method) 数值求解各类独立的偏微分方程；（电磁学、材料力学、工程热力学、声学等等） 矩量法 (Method of Moments) 适合于细线、平面形状结构的电磁场问题 电磁学基本方程 麦克斯韦方程组( Maxwell Equations) James Clerk Maxwell (1831–1879) (法拉第感应定律) (安培环路定律) (高斯定律－电场) (高斯定律－磁场) 物质本构关系 (Constitutive Relations) （欧姆定律） （电极化） （磁化） 一切电磁场和电磁波问题均可由以上方程，以及各类具体的边界条件所决定！ 麦克斯韦方程组中的运算符 散度（Divergence) 旋度（Curl) 连续函数的偏微分运算 设有一连续函数 , 现欲求　　　。 FDTD的基本思想 －连续偏微分的有限阶近似 －时域和空间域的离散化 二阶中心差分近似表达式： 当　　越小时，上式的近似程度越高。 实际上： FDTD空间域的离散化 (1)空间域的分割离散化 Ex分量的空间离散分布图 Yee元胞（Δx, Δy, Δz) Hx分量的空间离散分布图 节点 FDTD空间域的离散化 YEE 元胞 各电磁场分量在元胞中的位置 K.S. Yee, “Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media,”IEEE Trans. Antennas Propagat., vol. 14, 1966, pp. 302-307. 例如： Hz为Ex和Ey所环绕。 FDTD空间偏微分的近似 以Hz为例： 类似地，可实现各电磁场分量的 空间偏微分计算。 FDTD时间偏微分的近似 以Hz为例： t=(n+1/2)Δt t=nΔt t=(n-1/2)Δt 时间上的推移 (1)电场在时间上取整数倍的Δt; t=n *Δt; (2)磁场在时间上取(整数＋1/2)倍的Δt; t=(n +1/2)*Δt; 麦克斯韦方程的离散化近似 以Hz为例： 上式即为Hz的更新方程，由前一时刻的磁场和