Krylov子空间方法在本征值问题中应用研究.pdfVIP

第十三届反应堆数值计算与粒子输运学术会议 Krylov 子空间方法在本征值问题中的应用 研究 邓志红 孙玉良 李富 王登营 郭炯 （清华大学 核能与新能源技术研究院，北京 10084） 摘要：本征值问题是反应堆物理计算中的基本问题，快速求解本征值对实际工程意义重大。本文研究 Krylov 子空间方法在求解本征值问题中的应用。在扩散方程的有限差分近似下，首先研究了共轭梯度方法 代替超松弛块迭代；其次应用隐式再启动Arnoldi 算法代替幂迭代方法求解本征值问题。基于IAEA-3D 基 准问题，将应用Krylov 子空间方法的程序计算结果同有限差分扩散计算程序CITATION 的计算结果进行了 对比。结果表明，程序计算结果是一致的，说明实现是正确的；限制迭代次数的预处理共轭梯度方法能够 提高幂迭代方法的计算速度；在计算精度相当的情况下，隐式再启动Arnoldi 算法计算速度为不带Chebyshev 加速的幂迭代方法的4~5 倍，并且该算法能快速求解多个高阶本征值和高阶谐波。 关键词：本征值问题；Krylov 子空间方法；隐式再启动Arnoldi ； 中图分类号：TL325 文献标识码： A 1 引 言 CG ）方法代替超松弛块迭代对求解临界本 在反应堆物理计算中，求解多群中子扩 征值问题的改善情况；其次采用隐式再启动 散方程的本征值问题非常重要、非常普遍。 Arnoldi 算法（Implicit Restarted Arnoldi 对于本征值问题，采用有限差分近似的情况 Method, IRAM ）代替幂迭代方法，考察其计 下，对形成的大型线性方程组通常采用超松 算效率。以 IAEA-3D 基准问题为对象，将 弛块迭代作为内迭代进行空间上的求解，在 计算结果同扩散计算程序 CITATION 的结 外迭代采用幂迭代方法来求解最大本征值 果进行比较，研究 Krylov 子空间方法在本 （有效增值系数keff ）以及相应的本征函数 征值问题中的应用前景。 （中子通量）。幂迭代方法的收敛速度由占 2 基本方程 优比——次大特征值与最大特征值的比值 多群中子扩散方程的一般形式如下： 决定[1] 。当问题的占优比接近于1时，幂迭 −∇iD ∇φ +∑ φ S r g g R ,g g g ( ) 代方法收敛缓慢，甚至可能不收敛[2] 。通常， g −1 S (r) =∑ φ r + (1) 大型反应堆的本征值占优比非常接近于1， g ∑ g →g g ( ) 传统幂迭代方法收敛速度较慢，有必要研究 g 1 更有效的算法。 χg G ν ∑ φ r g 1, 2,...G k ∑( f )g g ( ) 为此，国外学者进行了一系列的研究。 g 1 Verdu等发展了子空间迭代与对称 为计算多维多群中子扩散方程，将方程 Raleigh-Ritz投影结合的方法[3] 。E.J. Allen等 在空间上进行有限差分离散，离散后得