摘要: 本文描述一个蒙特卡罗方法计算热中子通量和热利用因子的程序。利用全反射条件将问题转化为在直角三角形栅元内的求解问题。计算通量采用了通常的径迹长度方法。计算中考虑了热中子截面与能量相关、散射各向异性的情况。程序对栅元内的几何区域组成有很大的适应性。每区介质可由32种核素中的任何种核素组成。计算的结果是比较满意的。