定态中子输运方程本征值计算方法研究  

Numerical Computation of Eigenvalue for Steady Neutron Transport Equation

作　　者：洪振英[1] 袁光伟[1] 傅学东[1] 

机构地区：[1]北京应用物理与计算数学研究所,北京100094

出　　处：《核动力工程》2012年第3期6-11,共6页Nuclear Power Engineering

基　　金：中国工程物理研究院基金(2011B0202039);国家自然科学基金(11126135);国防基础科研计划(B1520110011);国家重点基础研究专项经费(2011CB309705)

摘　　要：散射源项各向异性展开阶数较大或者离散纵标方法的角度离散方向较多时,中子输运方程本征值的计算迭代容易失败。为了克服该问题,本文通过数学上的推导,构造了中子输运方程α本征值迭代求解的一种新方法,数值算例表明该方法提高了收敛速度,不收敛的问题也得到明显改善。对keff本征值的计算进行改进,改进后的方法不依赖迭代初值,数值算例表明改进方法的计算结果较好。The iteration progress will not succeed in the computation of α eigenvalue for the neutron transport equations, when considering more expanding number for anisotropic scattering source or more discrete direction for discrete ordinate method. In this paper, the numerical computation ofα eigenvalue is studied for spherical steady neutron transport equation and a new method is constructed. The new method canimprove the convergence rate and nonconvergent case. The modification of computation for keff eigenvalue will not dependent on the iterative initial value and the numerical results are stable.

关 键 词：α本征值 keff本征值 离散纵标 各向异性散射源 

分 类 号：O241.82[理学—计算数学]