反应堆α本征值蒙特卡罗计算方法研究

在蒙特卡罗(MCNP)模拟系统α本征值计算中常采用(k,α)回归法。该方法在超临界和浅度次临界系统计算中,可以得到较好的结果,但由于算法本身的局限,在深度次临界问题模拟时并不能给出正确的结果。为了解决蒙特卡罗方法计算深度次临界系统α本征值的问题,本文提出了基于蒙特卡罗瞬态程序(TMCC)计算α本征值的衰变曲线拟合法。该方法解决了计算α本征值本征分布的难点,可以有效地计算动态系统的α本征值。文中利用超临界和浅次临界问题对MCNP中(k,α)回归法和衰变曲线拟合法的计算结果进行了对比验证,两者结果符合较好,证明衰变曲线拟合法可以有效地处理超临界和浅次临界问题。在此基础上,利用实验测得的深次临界情况下的α本征值结果,对衰变曲线拟合法进行了验证。结果显示,衰变曲线拟合法可以有效地计算反应堆α本征值,并可以在深次临界情况下得到准确的计算结果。