c语言输运方程1001无标题,中子输运方程源项反演的瑞利商迭代算法.pdf

第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：)．

用球谐函数将中子角通量妒