脉冲堆燃料元件运输临界安全计算剖析.docVIP

脉冲堆燃料元件运输临界安全计算剖析 　　摘要：脉冲堆燃料元件运输时需要进行临界安全计算分析。介绍了临界安全计算用蒙特卡罗方法和MCNP程序。利用当前MCNP程序进行脉冲堆燃料元件临界安全计算时出现可执行程序运行不了的问题，利用Fortran97和VC++6.0对出现的问题进行了编译处理，然后利用脉冲堆冷态临界实验时的实测临界棒位数据对编译后的MCNP程序进行了校核计算验证，最后利用经校核的MCNP程序完成了脉冲堆燃料元件运输过程中的临界安全分析，临界安全计算结果表明脉冲堆燃料元件运输过程中满足临界安全要求。 　　关键词：脉冲堆 MCNP 临界安全 　　中图分类号：TL943 文献标识码：A 文章编号：1007-9416（2016）03-0000-00 　　核临界安全是核科技工业中特有的、极其重要的安全问题，根据临界安全法律法规要求，必须对脉冲堆燃料元件运输过程中的临界安全问题进行计算分析，确保燃料元件运输过程中的临界安全[1]。 　　1 脉冲堆燃料运输临界安全问题 　　脉冲堆1987年开始建造，1990年首次冷态临界，采用铀氢锆燃料-慢化剂粗棒状元件。当前装载有89根燃料棒，燃料棒高度760mm[2]。 　　当前需要对脉冲堆乏燃料卸料、运输进行前期准备工作。脉冲堆燃料元件运输临界安全计算分析是需要提前进行的重要工作之一。拟利用MCNP程序对脉冲堆乏燃料运输过程中的临界安全问题进行计算分析，保证脉冲堆乏燃料运输过程中的临界安全。 　　2 蒙特卡罗方法简介 　　蒙特卡罗方法是计算数学的一个重要分支。它是通过随机模拟和统计试验方法来求数学、物理等方面问题近似解的数值方法，因而也称为随机抽样技巧或统计试验方法。 　　其基本思想是：为了求解数学、物理、工程技术以及生产管理等方面的问题，首先建立一个概率模型或随机过程，使它的参数等于问题的解，然后通过对模型或过程的观察、抽样来计算所求参数的统计特征；最后给出所求解的近似值，而解的精度可用估计值的标准方差来表示。 　　在中子输运问题的求解中，需要使用蒙特卡罗方法建立在特定几何结构中的单一中子真实运动历史，但是单一的中子运动带有很大的随机性，无法作为实验参考，因此需要使用大量的中子重复计算，对每中子的输运历史都进行跟踪，进而得到一个随机实验值[3]。接着使用统计平均的方法计算出随机变量某个数值特征的估计值，使用该估计值作为输运问题的解。由于蒙特卡罗方法的统计量较多，计算问题的解可以包括中子和光子的通量密度、堆芯的有效增殖系数keff、剂量率、粒子沉积能和反应率等。蒙特卡罗方法计算程序主要有MCNP程序、MORSE程序、KENO等程序，其中MCNP程序是应用最广泛的蒙特卡罗程序。 　　3 MCNP程序 　　MCNP是三维几何输运理论蒙特卡罗方法程序，全名为 Monte Carlo Neutron and Photon Transport Code，即蒙特卡罗中子-光子输运程序，它可用于计算光子、中子、电子或者耦合光子/中子/电子输运问题，也可以计算临界系统的本?椭滴侍猓?可用于脉冲堆燃料元件运输临界安全计算分析。 　　MCNP具有以下特点与功能：（1）灵活的几何处理能力；（2）采用较多的降低方差的技巧，如几何分裂与俄罗斯轮盘赌、能量分裂、重要性抽样、相关抽样以及权重窗等；（3）使用精细的点截面数据。可以采用能点线性插值法获得连续能量，也可以把所有截面压缩成240群，还可以在不需要太多能点的情况下，用较小能点的连续能量计算。因为不用进行并群处理，因而具有特殊优势，特别对于燃料或堆芯布置是非均匀的反应堆系统（如先进研究堆、加速器驱动的洁净能源系统或空间用核反应堆等）；（4）计算功能齐全。它适用于核科学和工程方面的各种计算，如粒子输运、屏蔽计算、武器试验的测试分析以及检验核截面数据等；（5）程序输入、输出以及绘图功能齐全。[4][5] 　　4 临界安全计算 　　4.1 程序计算时遇到问题的解决 　　利用现有版本的MCNP进行脉冲堆燃料元件运输临界安全计算，可执行程序运行时遇到“bad trouble in subroutine chgmem of imcn value of mdas is too small”错误提示[6]，导致程序不能运行，经分析主要原因是程序动态分配的内存空间不够。因此需要对“MCNP.exe”进行重新编译。 　　利用Fortran97和VC++6.0对MCNP程序进行编译，得到新的可执行程序“MCNP.exe”。新的“MCNP.exe”需要得到可靠实验数据的验证证明正确性之后才能继续开展计算分析。 　　利用脉冲堆首次冷态临界实验时的临界棒位数据对编译后的MCNP程序进行验证。首次冷态临界时的临界棒位为安全棒到顶390mm，调节棒在底，补偿