论文部分内容阅读
本文对代数多重网格法(AMG)在求解由偏微分方程有限元离散得到的大型稀疏对称线性方程组上的应用进行了研究。特别针对信号完整性(SI)分析系统中,对Maxwell方程在时间域和频率域进行边型有限元(Edge FiniteElement)离散后得到的大型稀疏对称线性方程组上代数多重网格法的应用进行了研究,针对S1分析系统中特殊的三棱柱网格结构提出了经过修改的AMG算法(SI-AMG),并和Krylov子空间法结合来求解S1分析系统中碰到的超大规模线性方程组。SI-AMG算法中,由于三维有限元网格是由多层具有完全相同结构的二维有限元网格相互连接组成,因此只需对一层网格进行粗化,其中辅助矩阵的传递算子的构造与普通AMG法一样,系数矩阵的传递算子按照普通AMG中方法构造后再沿拓到整个三维有限元网格,粗网格上辅助矩阵和系数矩阵的计算仍采用Galerkin方法。本文也对经典的迭代算法进行了描述,其中有Gauss-Seidel法,实系数共轭梯度法(CG)、复系数共轭梯度法(Complex CG)、BiCG(BiConjugate Gradient)法、BiCGSTAB(BiConiugate Gradient Stabilized)法,并且推导了实系数预处理共轭梯度法(PCG)和预处理BiCGSTAB法。针对S1分析系统在时间域生成的实对称正定(SPD)线性方程组,我们将SI-AMG作为预处理共轭梯度法(PCG)中的预处理算子;针对频率域生成的复系数稀疏对称线性方程组,将SI-AMG作为预处理BiCGSTAB法中的预处理算子。当然,后者也可用于求解SPD线性方程组。最后,我们给出了S1分析系统碰到的一些实例,给出了网格粗化结果,并对不同的迭代法进行了比较,相对于其它迭代方法,SI-AMG预处理方法具有非常快的收敛速度。