时域有限差分法用类库的研制.pdf

摘 要* 时域有限差分法(FDTD)是计算电磁场中最常用的一种数值计算方法 它把时域麦 克斯韦旋度方程转化为差分方程来求解电磁场各分量 FDTD 广泛应用于研究电磁波 与各种物体的相互作用 如微波 光镊 近场光学 细胞中光场分布的数值求解 然而 由于FDTD 的建模问题具有其特殊性 如激励源的设置 介质参数 吸收 边界条件各不相同 传统的计算程序基于结构化程序设计方法设计 往往只能解决某 一特定问题 可重用性差 对不同问题必须重新编制程序 设计通用计算软件较为困 难 本文通过对 FDTD 算法进行面向对象分析 设计了一套用于 FDTD 数值计算的 C++通用类库(KFDTD) 类库设计的主要思想是将上述各特殊问题抽象成单独的类 与迭代过程分离 本文对类库及其核心类做了详细介绍 分析了每个类的内部结构和 函数接口 在此成功建立类库的基础上 进一步编制了高度聚焦高斯光束的 FDTD 计算程 序 所得数值计算结果与精确解进行了比较 二者吻合良好 该类库可用于确定性电 磁场问题的FDTD 数值求解 由于在设计中充分考虑了后续扩充的要求 代码重用率 高 易于维护和扩充 关键词 时域有限差分法 面向对象 激励源 高斯波束 *本课题受到国家自然科学基金和国家高技术研究发展计划项目(2003AA231011) 的资 助 I Abstract Finite-difference time-domain (FDTD) method provides numerical solution to the electromagnetic field. It transfers the Maxwell differential equations with time variant to the difference equations to find solve the electromagnetic fields. FDTD is widely used in the analysis of electromagnetic problems, such as microwave, Optical tweezers, near field optics and so on. Unfortunately, due to the particularity of the problems, which is caused by variety of conditions, such as the coupling of incident wave, the placement of scatterer, the absorbing boundary condition and so on, program codes for FDTD are usually model-specified. It is hard to make the code to use repeatedly. Usually, we have no choice but modify the source code to solve different problem., Till now, there is few common-used scientific computing software for FDTD. Object oriented Programming method brings us the possibility of realizing general FDTD software. The absorbing boundary condition, source excitation and all other special transactions are abstracted into individual classes, and are se