第29卷第3期 计 算 物 理 V01.29.No.3
2012年5月 CHINESE OFCOMPUTATIONALPHYSICS
JOURNAL May,2012
文章编号:1001-246X(2012)03-0369—06
中子输运方程源项反演的瑞利商迭代算法
梅立泉, 方 叶, 原长琦
(西安交通大学数学与统计学院,西安710049)
摘 要:中子输运方程源项反演问题中,有效倍增因子的求解经过球谐函数展开和差分离散后,转化为求解大型矩
阵的特征值问题.针对矩阵的特点,利用Gerschgorin圆盘定理,给出反幂法迭代初值的选取方法,并进行相应的扰
动分析.针对有效倍增因子是矩阵最大特征值的特点,采用瑞利商迭代算法进行求解.理论和算例表明,该方法不
需选取特定初值,达到三阶收敛速度.
关键词:中子输运方程;Gerschgorin圆盘定理;扰动分析;瑞利商迭代
中图分类号:0175.2 文献标识码:A
0 引言
.一
中子输运方程源项反演问题关注求解有效倍增因子k.||}的重要性基于以下物理事实¨’23:对于给定系
统参数(即给定材料成分和性质),只有在一定的几何尺寸下,系统才能达到临界状态:反之,对于给定的几
何条件就要求一定的材料特性与之配合.对于任意给定的系统,一般情况下方程式并不一定有解.但是在物
理上可以借助引进一个特征参数值来调节使它达到临界.裂变项前的这个系数就是k.因此,源项反演问题
在临界和次临界试验研究中有重要意义”3,可用来对实验装置设计进行理论估计,近年来,并行算法在该领
域中也得到广泛应用,与试验测量值进行参照,从而校验试验方法和计算方法.
传统源项反演特征值问题求解使用的数值计算方法是乘幂法.乘幂法简单易行,但是收敛速度较慢.文
[4]使用反幂法计算.如果事先可以对倍增因子做出较好的估计,该方法可以显著提高计算速度.但文[4]中
取方法,并做了相应的扰动分析.
根据矩阵稀疏和有效倍增因子是最大特征值这些特点,使用了瑞利商迭代算法.该方法适用于求解最大
特征值和相应特征向量,一般情况下收敛速度为二阶,矩阵对称正定时达到三阶口圳.将该方法与乘幂法、反
幂法做比较,结果表明,瑞利商迭代方法在收敛速度上与反幂法相同,远高于乘幂法.而且不用选取特定的迭
代初值,是求解此类问题更出色的一种方法.
1数学模型
1.1 中子输运方程与离散化
裂变介质核反应系统满足中子输运方程
1-2.v妒(戈,门,E)+盯。(耳,E)妒(菇,D,E)一f∑矿,(戈,E7)f(石,力,E,一n,E)妒(茁,力,,E,)dO,dE,:
÷』石1出,玑E)“x,E『)9(戈,n,E’)dO恤’, (1)
F1p(戈,E’,E)为裂变反应指示函数.
09—30
收稿日期:2011—06—22;修回日期:2011
基金项目:国家自然科学基金及中央高校基本科研业务费专项资金(xjj20looll2)资助项目
作者简介:梅立泉(1969一),男,博士,教授,蕊事偏微分数值解研究,E—mial:lqmei@mail.xjtu.edu,c”
万方数据
370 计 算 物 理 第29卷
采用迭代方法求解有效倍增因子k.设多。x,三,E)为初始中子角通量,i=1,2,…,有
7dE’, (2)
i1J石1王,(菇,E’,E)盯f(戈,E’)‰(z,力,E’)dn
则露=lim(1l疹;ll:/ll多。ll:).
用球谐函数将中子角通量妒