电磁波场数值计算中时间域有限差分法与时间域伪谱法的-现代地质.PDFVIP

第２６卷　第６期 现　代　地　质 Ｖｏｌ２６　Ｎｏ６ ２０１２年１２月 ＧＥＯＳＣＩＥＮＣＥ Ｄｅｃ２０１２  电磁波场数值计算中时间域有限差分法与 时间域伪谱法的对比研究 林树海，王伟利 （中国地质大学 地球物理与信息技术学院，北京　１０００８３） 摘要：为提高电磁波场数值计算精度，对时间域有限差分法与时间域伪谱法进行了对比研究。时间域有限差分法是一 种目前流行的电磁场时域数值计算方法，已被广泛应用于求解与时间有关的偏导数方程。对于大规模数值计算，时间 域有限差分法需要较多的内存空间。时间域伪谱法基于时间域有限差分法，该方法使用快速傅里叶度换来计算麦克斯 韦方程中的空间导数。由于傅里叶变换的准确性，时间域伪谱法使计算精度提高，数值计算时所需的格子数大大减少， 这极大地节省了计算机内存空间，适合于大规模正反演问题的数值计算。 关键词：数值计算；时间域有限差分法；时间域伪谱法；快速傅里叶变换；精度 中图分类号：Ｐ６３１３　　　　　文献标志码：Ａ　　　　　文章编号：１０００－８５２７（２０１２）０６－１１９３－０６ ＳｔｕｄｙｏｎｔｈｅＤｉｆｆｅｒｅｎｃｅｓＢｅｔｗｅｅｎＦｉｎｉｔｅＤｉｆｆｅｒｅｎｃｅＴｉｍｅＤｏｍａｉｎＭｅｔｈｏｄ ａｎｄＰｓｅｕｄｏｓｐｅｃｔｒａｌＴｉｍｅＤｏｍａｉｎＡｌｇｏｒｉｔｈｍｉｎＮｕｍｅｒｉｃａｌ ＣａｌｃｕｌａｔｉｏｎｏｆＥｌｅｃｔｒｏｍａｇｎｅｔｉｃＷａｖｅＦｉｅｌｄｓ ＬＩＮＳｈｕｈａｉ，ＷＡＮＧＷｅｉｌｉ （ＳｃｈｏｏｌｏｆＧｅｏｐｈｙｓｉｃｓａｎｄＩｎｆｏｒｍａｔｉｏｎＴｅｃｈｎｏｌｏｇｙ，ＣｈｉｎａＵｎｉｖｅｒｓｉｔｙｏｆＧｅｏｓｃｉｅｎｃｅｓ，Ｂｅｉｊｉｎｇ　１０００８３，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：Ｉｎｏｒｄｅｒｔｏｉｍｐｒｏｖｅｔｈｅｒｅｓｏｌｕｔｉｏｎｏｆｎｕｍｅｒｉｃａｌｃａｌｃｕｌａｔｉｏｎｏｆｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｗａｖｅｆｉｅｌｄｓ，ｔｈｅｄｉｆｆｅｒ ｅｎｃｅｓｂｅｔｗｅｅｎｆｉｎｉｔｅｄｉｆｆｅｒｅｎｃｅｔｉｍｅｄｏｍａｉｎ（ＦＤＴＤ）ｍｅｔｈｏｄａｎｄｐｓｅｄｏｓｐａｃｔｒａｌｔｉｍｅｄｏｍａｉｎ（ＰＳＴＤ）ａｌｇｏｒｉｔｈｍ ｗｅｒｅｓｔｕｄｉｅｄ．ＴｈｅＦＤＴＤｍｅｔｈｏｄｉｓａｃｏｍｍｏｎｎｕｍｅｒｉｃａｌｃａｌｃｕｌａｔｉｏｎｍｅｔｈｏｄｏｆｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｗａｖｅｆｉｅｌｄｓ，ａｎｄ ｉｔｉｓｗｉｄｅｌｙａｐｐｌｉｅｄｔｏｔｈｅｓｏｌｕｔｉｏｎｏｆｔｈｅｐａｒｔｉａｌｄｅｒｉｖａｔｉｖｅｅｑｕａｔｉｏｎｓｒｅｌａｔｅｄｔｏｔｉｍｅ．Ｍｏｒｅｍｅｍｏｒｙｉｓｎｅｅｄｅｄ ｕｓｉｎｇＦＤＴＤｍｅｔｈｏｄｉｎｌａｒｇｅｓｃａｌｅｎｕｍｅｒｉｃａｌｃａｌｃｕｌａｔｉｏｎ．ＴｈｅＰＳＴＤａｌｇｏｒｉｔｈｍｉｓｂａｓｅｄｏｎＦＤＴＤｍｅｔｈｏｄａｎｄ ｅｍｐｌｏｙｓｆａｓｔＦｏｕｒｉｅｒｔｒａｎｓｆｏｒｍ（ＦＦＴ）ｔｏｃａｌｃｕｌａｔｅｔｈｅｓｐａｔｉａｌｄｅｒｉｖａｔｉｖｅｓｉｎＭａｘｗｅｌｌｓｅｑｕａｔｉｏｎｓ．Ｄｕｅｔｏａｃｃｕｒａ ｃｙｂｙＦＦＴ，ｔｈｅｃａｌｃｕｌａｔｉｏｎｐｒｅｃｉｓｉｏｎｉｓｉｍｐｒｏｖｅｄａｎｄｇｒｉｄｎｕｍｂｅｒｓｎｅｅｄｅｄｉｓｄｅｃｒｅａｓｅｄｇｒｅａｔｌｙｂｙｕｓｉｎｇＰＳＴＤ ａｌｇｏｒｉｔｈｍ．ＴｈｅＰＳＴＤａｌｇｏｒｉｔｈｍｈａｓｔｈｅａｄｖａｎｔａｇｅｏｆｓａｖｉｎｇｃｏｍｐｕｔａｔｉｏｎａｌｍｅｍｏｒｙｓｐａｃｅａｎｄｉｓｓｕｉｔａｂｌｅｆｏｒ ｌａｒｇｅｓｃａｌｅｎｕｍｅｒｉｃａｌｃａｌｃｕｌａｔｉｏｎｏｆｆｏｒｗａｒｄｍｏｄｅｌｉｎｇａｎｄｉｎｖｅｒｓｉｏｎ． Ｋｅｙｗｏｒｄｓ：ｎｕｍｅｒｉｃａｌｃａｌｃｕｌａｔｉｏｎ；ｆｉｎｉｔｅｄｉｆｆｅｒｅｎｃｅｔｉｍｅｄｏｍａｉｎｍｅｔｈｏｄ；ｐｓｅｕｄｏｓｐｅｃｔｒａｌｔｉｍｅｄｏｍａｉｎａｌｇｏ ｒｉｔｈｍ；ｆａｓｔＦｏｕｒｉｅｒｔｒａ