中子输运方程源项反演的瑞利商迭代算法.pdfVIP

第29卷第3期 计 算 物 理 V01．29．No．3 2012年5月 CHINESE OFCOMPUTATIONALPHYSICS JOURNAL May，2012 文章编号：1001-246X(2012)03-0369—06 中子输运方程源项反演的瑞利商迭代算法 梅立泉， 方 叶， 原长琦 (西安交通大学数学与统计学院，西安710049) 摘 要：中子输运方程源项反演问题中，有效倍增因子的求解经过球谐函数展开和差分离散后，转化为求解大型矩 阵的特征值问题．针对矩阵的特点，利用Gerschgorin圆盘定理，给出反幂法迭代初值的选取方法，并进行相应的扰 动分析．针对有效倍增因子是矩阵最大特征值的特点，采用瑞利商迭代算法进行求解．理论和算例表明，该方法不 需选取特定初值，达到三阶收敛速度． 关键词：中子输运方程；Gerschgorin圆盘定理；扰动分析；瑞利商迭代