预条件共轭梯度法在频域有限差分法分析二维柱体电磁散射问题中的应用.pdfVIP

第 11 卷第 1 期 　　　　　　　　　　　淮 阴 工 学 院 学 报 　　　　　　　　　　　Vol . 11 No . 1 2002 年 2 月 　　　　　　　　　Journal of Huaiyin Institute of Technology 　　　　　　　　　　　Feb . 2002 ( ) 文章编号 :1009 - 796 1 2002 0 1 - 0008 - 03 预条件共轭梯度法在频域有限差分法分析 二维柱体电磁散射问题中的应用 刘淑静 ,朱汉清 (淮阴师范学院物理系 ,江苏 　淮安 　22300 1) 摘 　要 :本文引入预条件共轭梯度法 ,提出了结合频域有限差分法结合分析二维柱体的电磁散射问题 。数值计算过程中利用 Mur 二阶吸收边界条件和场平均吸收条件截断网格 。作为算例 ,文中分析了一无限长理想导体柱对平面电磁波的散射 , 由于使用了预条件 共轭梯度法求解差分矩阵方程 ,从而减少了计算时间。数值结果表明了该方法的有效性 。 关键词 :预条件共轭梯度法 ;频域有限差分法 ;吸收边界条件 ; 电磁散射 中图分类号 : TN 817 ,O44 1 ;文献标识码 :A 1 　引　言 TMz 和 TEz 波分别求解 。下面我们将简要说明有限差 有限差分法是最早用来进行电磁场数值计算的方 分方程的建立和频域 Mur 二阶吸收边界条件的差分近 法 。在频域有限差分法中 ,无论是差分方程还是边界条 似构造 ,并引入 PC G 方法 。 件的差分近似都较时域法简单 、方便 ,而且在频域中不 必对时间作量化处理 ,可以减少迭代过程 ,提高计算精 2 . 1 内结点的频域有限差分方程 度[ 1] 。同时利用频域有限差分法分析二维柱体电磁散 由导体柱表面向外作法向射线并建立多层共形网格离 射 ,对不规则 、非均匀的散射 问题有着广泛的适用性 。 散求解区域 。用 u ( x ,y) 表示纵向场分量 ,显然每个结 在对差分方程的求解中 , 由于其系数矩阵非正定 ,经典 点上的总场 u 可分解为入射场 ui 与散射场 us 的和 ,即 i s ( ) 的迭代方法不能保证迭代过程中的收敛性 ,为此常采用 　　u = u + u 　　　　　　　　　 1 直接法求解 。但是 ,随着差分离散结点数的增加 ,矩阵 且 us 满足 Helmholtz 方程 ,设时间因子为 ,则有 2 S 2 s 方程阶数的提高 ,该方法求解方程的时间将会显著增 u u 2 ( ) 　　 + + k u = 0 2 x2 y2 0 大 ,从而极大地限制了频域有限差分法在分析许多实际 式中。在内结点上由环路积分方法可以得到关于 ( ) 问题中的应用 。共轭梯度法 简称为 C G 方法 作为求 散射场的有限差分方程