基于辛算法的二维电磁场散射问题的研究

摘要

辛算法是保持Hamilton系统辛结构的一种新的数值方法,由于 Maxwell方程是一无穷维Hamilton系统,因此可将辛算法用于电磁场模拟中.本文提出一种基于辛分块Runge-Kutta(PRK)方法的显式辛算法,并将它成功应用于二维电磁散射问题的计算中.通过对金属方柱散射场的数值模拟,比较了FDTD法和低阶辛算法(一阶和二阶),结果表明低阶辛算法不仅与FDTD法精度相当,而且可以减少存储空间和计算时间,尤其是一阶辛算法节省了大约的CPU时间,提高了计算速度,体现了该算法的优越性.