摘要

蒙特卡罗方法或离散纵标法模拟多群中子输运方程,散射角分布通常按Legendre级数展开,取L阶截断近似。当L较小时,近似角分布在[-1,1]定义域的多处出负,面积不归一,进而导致非物理解。文章针对这一问题,利用矩和Legendre系数等价关系,给出离散广义Gauss求积组,并用该求积组作节点构造等概率阶梯函数,逼近实际角分布,从而消除非物理解。计算取得了与实验和连续截面蒙特卡罗MCNP程序一致的结果。