论文部分内容阅读
自1873年,麦克斯韦建立了描述电场强度和磁场强度随着时间和空间变化关系的电磁波方程以来,电磁波理论发展已经有一百多年的历史。目前,电磁波的研究和应用已经深入到各个领域。然而,由于实际环境的复杂性,通常只有一些经典问题有解析解,因此往往需要通过数值解得到具体环境下的电磁波特征。辛算法是由我国科学家冯康先生在上个世纪80年代提出并系统总结出来的一类拥有稳定性和长时间计算准确性的数值方法。近年来,Bridges和Marsden等人推广了辛算法,提出了多辛算法的概念。随着多辛理论的逐步建立,许多学者尝试了用多辛算法来计算比较著名的偏微分方程,比如:非线性波方程,schr(o|¨)dinger方程,Korteweg-de Vries方程,Zakharov-Kuznetsov方程,Kadomtsev-Petviashvili方程并取得了满意的效果。目前,多辛算法计算麦克斯韦方程还没有人尝试过。本文的出发点是发展麦克斯韦方程的多辛理论和多辛算法并检验了多辛算法求解2+1维电磁波方程的有效性。基于麦克斯韦方程的Bridges多辛结构,我们得到了三个求解麦克斯韦方程的多辛格式。在第二章,我们推得一个求解麦克斯韦方程的多辛Euler-box格式,并且证明这个格式在时间方向是自共轭的。接着利用这个格式得到了关于时间方向的组合格式,而且这个组合格式也是多辛的。第三章,我们在时间和空间方向同时用辛中点离散规则离散麦克斯韦方程的Bridges形式,消去参量从而得到了一个求解麦克斯韦方程的两层102点格式。通过推导,我们证明了这个多辛中点格式能够精确保持离散的局部能量。在这章中,我们对中点格式做了后误差分析,得到了修正方程。通过修正方程,得到了这个格式的截断误差。自共轭格式和组合格式进行的两个数值实验结果表明,这两个多辛格式用来计算2+1维的麦克斯韦方程是相当有效的。实验结果也表明了组合思想的是切实可行的。Preissman格式与Yee格式的数值实验比较发现,中点格式能够得到更好的数值结果。数值结果也反应了Preissman格式相当好的保持了离散的局部能量和解的波形。因此,我们可以说多辛算法是计算求解麦克斯韦方程的一个有力的工具。