定态中子输运方程本征值计算方法研究 numerical computation of eigenvalue for steady neutron transport equation.pdfVIP

第33卷第3期 核动力工程 V01．33．No，3 2012正6月 NuclearPower JUn．20I2 Engineering 文章编号：0258-．0926(2012)03-0006-06 定态中子输运方程本征值计算方法研究 洪振英，袁光伟，傅学东 北京应用物理与计算数学研究所．北京。100094 摘要：散射源项各向异性展开阶数较大或者离散纵标方法的角度离散方向较多时，中子输运方程碎征 值的计算迭代容易失败。为了克服该问题，本文通过数学上的推导，构造了中子输运方程钵征值迭代求解 的一种新方法，数值算例表明该方法提高了收敛速度，不收敛的问题也得到明显改善。对t甜本征值的计算 进行改进，改进后的方法不依赖迭代初值，数值算例表明改进方法的计算结果较好。 关键词：a本征值；七《本征值；离散纵标；各向异性散射源 中图分类号：0241．82文献标志码：A 1 引 言 2中子输运方程D本征值求解过程 中子输运方程描述中子群体在介质内的输运 口本征值可以通过解非定态中子输运方程求 过程。除了特殊例子外，它不能解析求解而只能 得，也可以通过解定态中子输运方程来求解。求 数值求解(1吲。在很多应用领域中，数值求解此类 解盘本征值的定态中子输运方程为： 方程是数值模拟的关键。中子输运方程的数值求 (1) 解分为确定论方法和随机模拟方法。确定论方法 钒+专噍2‰善+Qf’s 根据角变量的不同处理方式分为离散纵标方法 对于一维球几何问题： (Sn)【4叫和球谐函数展开方法(Pn)【7JJ。在中 子输运方程中，精确计算各类本征值[1,9,101对工程轮笋掣t掣峨；噍 应用有很大影响。 实际运用中，当口本征值给出的中子通量密 (2) 度谱非常接近实际情况，即处理高超临界和深次 源项绞包含裂变源、散射源： 临界的问题时，口本征值较符合实际，而对于临 界问题，口刁t征值的计算结果往往有一定偏差。 纬=去妇|g+G名) (3) Sn方法中，在对定态问题求解口本征值时，当散 射源项中勒让德展开采用高阶形式(如n，P5) 外边界条件：珐以，∥)=0；∥s0。 或sn中角度离散方向较多(如S16，S知)的情况通常不直接求口本征值，而是采用尝试差值 下，有时迭代过程会失败。针对这些问题，本文 法，通过计算本征值丘啦来确定口【l】： 对一类迭代求解方法进行分析，初步给出口本征 值收敛速度缓慢甚至不收敛的一些原因，对原求 k妨+詈缸=‰+甓 解方法进行改进，提出新的方法。在一维球几何 (4) 多群定态问题中对一些算例进行对比计算与分 本文并不利用盘疗本征值对口本征值进行求 析；对于恕fF本征值本文也对其数值求解方法进行 解，而是直接利用式(1)，通过内外迭代得出口本 改进，使得计算结果更加稳定。 征值： 收稿日期：2011-03-24；修回日期：2011．-12．23 研究专硬经费(2011CB309705) 万方数据 洪振英等：定态中