高阶辛算法的稳定性与数值色散性分析(1)

＼８－１男３。３
（４）
０
男＝
ａ
０ｚ
ａ Ｏｚ
０
ａ ａｙ
ａ ａｘ
（５）
其中占和ｐ是媒质的介电常数和磁导率． 式（３）从ｔ＝０到ｔ＝△。的时间演化为
ａ
ａ
０
０ｙ
ａｘ
（：）（△。）＝ｅｘｐ（△。（Ｕ＋矿））（：）（０），
（６）
其中ｅｘｐ（Ａ。（Ｕ＋ｙ））为麦克斯韦方程的时间演化矩阵． １．１时间演化矩阵的近似
（９）
ａ孝Ｉｆ：船， 鲁’
△ｆ
。‘
其中手＝茹，Ｆ，Ｚ，皿ｉ，ｊ，ｋ，Ｗ，是空间差分系数，其定义如下：
ｑ／２
Ｅ＝２∑Ｉ形，Ｉ ｒ。１
（１０）
表１给出了常见的矽，和Ｅ之值．
Ｔａｂｌｅ １
表１ ｑ阶交错差分系数 Ｓｔａｇｇｅｒｅｄ ｄｉｆｆｅｒｅｎｃｅ ｃｏｅｆｆｉｃｉｅｎｔｓ ｏｆ ｏｒｄｅｒ口
本文主要研究四阶交错差分，对于四阶以上差分算法分析方法类似．经过空间ｑ阶交错差分、时间Ｐ阶 辛积分离散后，我们建立了Ｍａｘｗｅｌｌ方程的离散辛框架，即ＳＦＤＴＤ（ｐ，ｑ）算法．
１２０
１６０
母严
（ｂ）相对相速度误差随球面角妒的变化
图１高阶ＳＦＤＴＤ（４，４），ＦＤＴＤ（２，４），ＦＤＴＤ（２，２）方法的色散性比较 Ｆｉｇ．１ Ｄｉｓｐｅｒｓｉｏｎｓ ｉｎ ＳＦＤＴＤ（４，４），ＦＤＴＤ（２，４）ａｎｄ ＦＤＴＤ（２，２）
其次，我们将ＳＦＤＴＤ（４，４）方法同高阶Ｒ—Ｋ（４，４）方法和高阶Ｊ－Ｆａｎｇ（４，４）方法相比较．其中， Ｒ—Ｋ（４，４）方法使用４阶龙格库塔时间步进和４阶空间交错差分．从图２（ａ），（ｂ）可以看出，ＳＦＤＴＤ（４，４）算 法的色散性优于Ｒ－Ｋ（４，４）方法但比Ｊ－Ｆａｎｇ（４，４）方法略差．虽然如此，一方面由于Ｊ－Ｆａｎｇ（４，４）方法有幅度 耗散，另一方面该方法使用了空间高阶偏导数，极难处理复杂边界，因而在电磁场工程中很少应用．
（１．安徽大学计算智能与信号处理教育部重点实验室，安徽合肥２３００３９； ２．香港大学电机电子工程系，香港；３．合肥师范学院物理电子工程系，安徽合肥２３０６０１）
摘要：利用Ｍａｘｗｅｌｌ方程的哈密尔顿函数，导出对应的欧拉一哈密尔顿方程．利用辛积分技术与高阶交错差分
技术，建立求解三维时域Ｍａｘｗｅｌｌ方程的高阶辛算法；结合电磁场中的物理概念，借助矩阵分析和张量分析理
次相乘，可得
ｔｒ（Ｓ）＝２＋∑ｇｌ（Ｆ。２△Ａ ２。仉２，Ｘ‘，
（２１）
反２，。＾。ｌ｝２。Ｚ秒。Ｃｉｔ嘭ｔｃｉ２ｄｊ２＂ｃｉｚ叱＋。ｌ。Ｊ。；‘ Ｚ。…。＾‘ 。ｄｑｉＩ ＾‘ｌ｝２‘丘‘…‘ｉｒＵｈ。 ｍ
其中ｖｏ＝圭为介质中的波速． ｑ雌
ｌ‘ ｉＩ ‘ｊｌ‘ １２＜１２ ‘…‘ｉ‘ ｌ Ｊｌ‘ｍ
第２７卷第１期 ２０１０年１月
计
算
物
理
ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＣＯＭＰＵＴＡＴｌ０ＮＡＬ ＰＨＹＳＩＣＳ
文章编号：１００１—２４６Ｘ（２０１０）Ｏｌ－００８２－０７
Ｖ０１．２７，Ｎｏ．１ Ｊａｎ．。２０１０
高阶辛算法的稳定性与数值色散性分析
黄志祥１， 沙 威２， 吴先良１’３， 陈明生３， 况晓静１
响，并且比较了不同高阶算法的数值色散特性．从结果
图３ 区域中心的Ｅ，场迭代２００个周期后的波形
中可以看出： Ｉ）ＳＦＤＴＤ（４，４）与ＦＤＴＤ相比稳定性与色散性都
Ｆｉｇ．３
Ｗａｖｅｆｏｒｍ ｏｆ Ｅ＃ａｔ ｃｅｎｔｒｅ ｏｆ ｔｈｅ ｃｏｍｐｕｔａｔｉｏｎａｌ ｄｏｍａｉｎ ｗｉｔｈ ２００ ｐｅｒｉｏｄｓ ｉｔｅｒａｔｉｏｎ
教育厅重点项目（ＫＪ２００８Ａ１００＆ＫＪ２００８Ａ０３６）资助项目 作者简介：黄志祥（１９７９一），男，安徽广德，教授，博士，主要从事计算电磁学、散射与逆散射、雷达成像的研究．
万方数据
第１期
黄志祥等：高阶辛算法的稳定性与数值色散性分析
８３
Ｕ＝
一ｐ’１男３ｘ３、
ｔ０１，。，）６ｘ６ ９
ｙ：ｆ Ｉ ｏ‰
行进一步的离散．在此，我们基于交错差分近似算符Ｒ中一阶偏导数，引入记号厂州”（ｉ，．『，ｋ）＝
八讼，，ｊＡ，，ｋＡ：；（ｒｔ＋下。）△；）为近似离散点（１ａｘ，ｊＡｙ，ｋＡｚ）在第儿时间步的第Ｚ级真实解，这里，每个时间步需
要ｍ级时间步进，且第ｚ级步进对应的时间增量为下。△。．
Ｉ：兰Ｅ坐型生土盟≠竺暨业＋ｏ（△ｍ 了ａＦｎ＋ｔｌｍ 利用ｑ阶交错差分离散空间一阶偏导数
ＦＤＴＤ（４，４）方法的最大稳定度．这里，所有的空间差分均
采用显式交错差分．其中ＳＦＤＴＤ（４，４）方法使用的辛算子
可参考文［４］．
２．２色散性分析
由于ＩＡ¨Ｉ－Ｉ，该算法是无耗散的．根据反余弦求出
特征值的相角，则色散方程可表述为
ｔｏＡｔ＝ａｒｃｃｏｓ［ｔｒ（Ｓ）／２］，
（３０）
裘２各种算法的最大稳定度比较 Ｔａｂｌｅ ２ ＭａｘｉｍｉｚａｔｌＯｎ ｓｔａｂｉｌｉｔｉｅｓ
论，获得高阶时域方法及高阶辛算法的稳定性和数值色散性的统一处理新方法．用数值结果证实方法的正确性，
与ＦＤＴＤ算法和其它时域高阶方法相比，高阶辛算法具有较大的计算优势，为电磁计算提供了新的途径．
关键词：哈密尔顿函数；辛积分技术；稳定性和数值色散性；高阶辛算法
中图分类号：ＴＮ８
文献标识码：Ａ
Ｏ 引言
时域有限差分…（ＦＤＴＤ）法以其简单直观的特点已成为电磁学数值计算的一种常用方法，它直接求解 依赖时间的Ｍａｘｗｅｌｌ旋度方程，利用二阶精度的中心差分近似旋度方程中的时间及空间微分算符，极易处理 非均匀媒质的情形．然而，它的计算精度相对较低，计算的时间步长与空间离散网格的大小必须满足 Ｃｏｕｒａｎｔ．Ｆｒｉｅｄｒｉｅｈｓ．Ｌｅｖｙ（ＣＦＬ）稳定性条件，并且随着计算时间步的增加，误差将会累积．为克服ＦＤＴＤ方法 的这些缺点，国内外学者进行了一系列的改进，主要分为时间离散方式和空间差分格式．时间离散方式主要 有交替方向隐式时间步进法（ＡＤＩ）［２１、四阶龙格库塔法（Ｒ—Ｋ）∞１及辛积分技术Ｈ圳；空间差分格式的改进主 要为高阶交错差分技术＂‘８】、四阶隐式紧差分技术”］，多分辨扩展方法ｐ’及伪谱近似方法¨叫等．
ｅｘｐ（Ａ。（Ｕ＋ｙ））＝兀（，＋ｄｌＡ。Ｖ）（，＋ｃａＡ。Ｕ）＋Ｄ（∥１），
（８）
可以看出，该方法具有显式逐级递推的特点．因而，同Ｒ－Ｋ和其它一些隐式算法相比，大大节约了内存，提高
了计算的效率．
１．２空间近似
由于算符ｕ，ｙ中含有三维旋度算符男，为了得到Ｍａｘｗｅｌｌ方程的数值解，必须在空间方向上对方程进
ｓｓ＝ｎ［一２了叩１，Ｈｄ，△。。ｌｌ一／：了一石叩１ ，ｄ，△。 ： ｐ：７：，７，‘ｆ△△‘‘］』．
ｃ‘·１７，’
ｓ：ｐ ｓｔ：１，
（１８）
万方数据
第１期
其解Ａｍ：型
黄志祥等：高阶辛算法的稳定性与数值色散性分析 稳定性条件要求ＩＡ，．：Ｉ＿１，所以只需ｔｔｒ（Ｓ）Ｉ≤２即可．将（１７）式中矩阵依
Ａｌｇｏｒｉｔｈｍ
ｍ似 阳似
＂
ｍ
辱；哪脚珊但＠””水ｍ朋朋”∽
０．５７７ ０．４９５ ０．５７７
Ｏ．７∞
０．７４３
定义相对相速度误差
万方数据
这里％２若为数值相速度·
计
算
物
理
肪蚴·ｇ Ｉ警』，
第２７卷 （３１）
３ 数值结果
首先，我们将ＳＦＤＴＤ（４，４）方法同高阶ＦＤＴＤ（２，４）方法和传统的ＦＤＴＤ（２，２）方法相比较．从图１（ａ）， （ｂ）可以看出，与这些算法相比，ＳＦＤＴＤ（４，４）算法有明显的优势．
ＳＦＤＴＤ（４，４）方法．从图中可以看出，高阶辛算法较传
甲
统的ＦＤＴＤ法有明显的优势．
昌
之
‘一
４ 结论
本文导出了基于辛积分技术的高阶时域计算方法
的统一辛框架，利用矩阵分析，首次系统探讨了常见高
阶时域方法的稳定性与数值色散性方程，详细地讨论 了ＣＦＬ、时间步长以及传播方向对数值色散误差的影
Ｒｅａｌｔｉｖｅ ｔｉｍｅ ｓｔｅｐ（ｎ）
接下来，我们用１０ ＧＨｚ的余弦源进行试验，分析区域中心（２０，２０，２０）处的时域场值．余弦波经过２００个
万方数据
第１期
黄志祥等：高阶辛算法的稳定性与数值色散性分析
８７
周期的迭代后，我们记录其一个周期的波形如图３所
示．计算中稳定度常数ＣＦＬ＝０．５，空间分辨率每波长
６个采样点，算法采用ＦＤＴＤ（２，２）方法和
鼍
晕 ５
窘
量
喜
喜
＞
莹
蓍 甚
耋
譬
要 蛊
ＰＰｗ （ａ）相对相速度误差随空间分辨率的变化
ｍｆ （ｂ）相对相速度误差随球面角矿的变化
图２高阶ＳＦＤＴＤ（４，４），Ｒ—Ｋ（４，４）方法，高阶Ｊ－Ｆａｎｇ（４，４）方法的色散性比较 Ｆｉｇ．２ Ｄｉｓｐｅｒｓｉｏｎｓ ｉｎ ＳＦＤＴＤ（４，４），Ｒ—Ｋ（４，４）ａｎｄ Ｊ－Ｆａｎｇ（４，４）
可以通过空间稳定度因子和Ｉｔｒ（Ｓ）ｌ≤２的约束求出．这种时空分离的形式，对分析任意ＳＦＤＴＤ（Ｐ，ｇ）算法
的稳定度是很方便的，尤其适用于时间方向或空间方向其中一个离散方式固定而另外一种改变的情况．
表２列出了传统的ＦＤＴＤ（２，２）方法，高阶ＦＤＴＤ（２，４）
方法，Ｊ－Ｆａｎｇ（４，４）方法¨４１，龙格库塔（４，４）方法和高阶辛
有无可比拟的优势，ＳＦＤＴＤ（４，４）方法在每波长取６个
点时，误差可达到一６０ ｄＢ的效果；ＦＤＴＤ（２，２）在每波长取１２个点时，误差仅达到一４０ ｄＢ的效果；
ＦＤＴＤ（２，４）在每波长取１２个点时，误差仅达到一５０ ｄＢ的效果；
２）ＳＦＤＴＤ（４，４）与Ｒ．Ｋ（４，４）及Ｊ－Ｆａｎｇ（４，４）相比具有更好的稳定性，ＳＦＤＴＤ（４，４）较Ｒ－Ｋ（４，４）的数值
利用ｍ级Ｐ阶显式辛积分技术Ｈ’５１近似时间演化矩阵ｅｘｐ（ｄ，（Ｕ＋ｙ）），即
ｅｘｐ（Ａ。（Ｕ＋Ｖ））＝丌ｅｘｐ（ｄｌＡ。Ｖ）ｅｘｐ（ｃｔＡ。Ｕ）＋ｏ（△∥），
（７）
这里Ｃ，和ｄ，为辛算子或者辛传播子．注意到无论对于连续的微分矩阵或者离散后的差分矩阵，总可以证明 Ｕ７＝∥＝Ｏ（∥≥２）．使用泰勒展开技术，（７）式可以写成
¨ｚ％…‰，
（２２）
去（：）＝［÷。叼，口，＋二，＋叩，口，，×吾（，７，Ｐ＃＋：Ｐｙ＋叩ｚＰｊ）×］（：）， （２３）
一如型丛型型等—业生堕巡，
（２４）
＂如巫丛型型笔—业生丝幽，
（２５）
二［豆ｊ匾
［亟ｊ匾 ｏ
旦严＼一 ＯｔＩＥＪ＿
ｐ
（２６）
。
占
这里，露是由球面角定义的张量矩阵…１．虽然形式上式（２６）是６ Ｘ６的矩阵，但它只有两个独立的本征值．
本文利用辛积分技术与高阶交错差分技术，建立了求解时域Ｍａｘｗｅｌｌ方程的高阶辛时域有限差分算法， 结合电磁场中的物理概念，借助矩阵分析和张量分析理论，将文［１１］中二维高阶辛算法的稳定性分析推广 至三维情形，得到了常见高阶时域方法及高阶辛时域差分算法的数值稳定性和色散性的统一处理新方法，为 电磁计算提供了新的理论与计算方法．
１ 理论
无源、均匀、无耗媒质中的Ｍａｘｗｅｌｌ方程的哈密尔顿函数可定义如下：
宣，（日，日）：了１－【—１—日．ｖ×Ｈ＋ｊ二Ｅ．ｖ×层），
厶８
斗
Байду номын сангаас
其对应的欧拉一哈密尔顿方程为
Ｏ～ＯＨｔ ＝一一ａＯＥ．ｇ’
Ｏ一ＯＥｔ ＝——ＯａＨ毫．’
（２）
通过变分法，式（２）可以改写为
击（日）刈川㈣
（３）
收藕日期：２００８—０９—０４；修回日期：２００９一０２—１８ 基金项目：国家自然科学基金重点项目（Ｎｏ．６０９３１００２），国家自然科学基金（Ｎｏ．６０６７１０５１），高校博士点基金（Ｎｏ．２００６０５５７００４）及安徽省
ｏ５
霍
呻Ｏ 。
＋ＳＦＤＴＤ（４，４）
＋ＦＤＴＤ（２，∞
。５
·－·一ＦＤＴＤ（２，２）
耄．
重
巧 Ｏ ．ｒ ｒ ｒ ｒ，ｒ ｒ ｒ ｒ ｒ
ｒ
ｒ
ｒ ｒ—ｒ
ｒｒｒｒ
ｒ
弓５
麦
粤毛ｇＱ茸８一竽骘Ｂｔｌ厶爹薯一星 咱Ｏ
ｅ＝６｛妒。ＰＰｗ一８．ｃＦｌ商．４
ＰＰＷ （ａ）相对相速度误差随空间分辨率的变化
吒５
０
４０
８０
万方数据
２数值稳定性、色散性分析
计算
物
理
第２７卷
警：一芝茎ＥＥ弛型世丛上生堕丝产坠业＃三业业止型Ｆ
其中叩，：兰Ｅ竺丛二丛生二旦型生垒丢二＿竺丛型三二型．
击（≥）：ｆ一三旦一言啬］（≥），
ｃ－５，
击（≥）＝［一圣仇一言仇］（≥），
ｃ－６，
记ｕ＝［三一言仇］，ｙ＝［一；仇三］，时间上采用高阶辛积分技术逼近Ｍａｘｗｅｕ的时间演化矩阵，则稳定
因此，式（１６）和（２６）本质上是同构的．采用类似于一维的分析方法，可求出
ｔｒ（Ｓ）＝２＋∑ｇ。｛移。２△２。（’７：＋叼；＋田：）｝１．
（２７）
一般而言，时域算法的最大稳定度ＣＦＬ。。可以写成时空分离的形式
ＣＦＬ…＝竿，
（２８）
＾Ｓ
这里Ａ。为空间稳定度因子，可以表示为
Ａｓ＝ｄｖ厮Ｗ，，
（２９）
ｄｉｍ＝Ｉ，２，３为空间维数，形，的定义已经在式（１０）中给出，为所有空间差分系数之和．Ａ，为时间稳定度因子，