有效增殖因子的动态模拟计算.pdfVIP

第 31卷 第 4期 核 动 力 工 程 VO1．31．NO．4 2 O 10 年 8月 NuclearPowerEngineering Aug．2 0 10 文章编号：0258-0926(2010)04-0001—04 有效增殖因子的动态模拟计算 谢奇林 ，刘汉刚，杨成德 ，贺仁辅，范晓强 (中国工程物理研究院核物理与化学研究所，四川绵阳，621900) 摘要：讨论了有效增殖因子(fr)的动态计算方法 ，建立了基于Geant4toolkit平台的 fr蒙特卡罗 (MC) 动态计算模型。分别采用动态方法和 K静态方法计算了GodivaI反应堆基准实验和偏离临界较远系统时的 fr。 结果表明，K 静态方法仅适用于靠近临界的系统，而动态方法无论在临界附近还是偏离临界较远均可 得到正确的结果。 关键词 ：蒙特卡罗；K静态方法 ；动态方法 ；有效增殖因子 中图分类号：TL329+．2 文献标识码：A 1 引 言 = ，ESf，EJ+S+ ，J 运用蒙特卡罗 (MC)方法进行临界计算求 式中， 为中间变量； = (r，I2，E，f)为t时刻r 解系统的有效增殖因子 ( )时，一般采用K静 处、能量为E、运动方向为 的单位立体角元内 态方法。该方法收敛迅速，在系统接近临界时(临 的中子角通量；S=S(r，，E，f)为初始源分布； 界系统在数学上为存在非负非零解的稳态无源问 (，．，E)、 (，E)和 (，．，E)分别为裂变、散射 题)结果准确，但在系统偏离临界较大时 (不再 和总的宏观截面；F(，．；Et，t2’ E，)为转移概率； 是稳态问题 )，往往会得出错误的结果。而动态 n(r，E)为每次裂变产生瞬发中子的平均数；v为 方法直接模拟中子的真实输运过程，不引入任何 中子速率。 近似，无论在临界附近还是偏离临界较远都可以 动态方法对时间相关 的中子输运问题求解 得到正确的结果 ；但是，动态方法消耗的计算资 时，从任意初始源分布开始 ，按照时间顺序，对 源相对要多 l【】。随着计算机技术的发展 ，动态方 所有源中子及次级中子发生的所有物理过程进行 法将会在临界计算中发挥更为重要的作用。 跟踪模拟，直至系统内中子通量分布演化为基模 本文从时间相关的中子输运方程和 fr的定 本征分布。这时即可对感兴趣的物理量 如【 义出发对动态方法进行了介绍，计算了 GodivaI 非伴随瞬发中子平均寿命 (r一 )和瞬发中子衰 反应堆进行基准实验时的 ff，并将计算结果与K 减常数 ( )等】进行统计求解。 静态方法的计算结果及实验结果进行了比较。为 对于给定系统，从中子平衡角度来看，有效 了说明K静态方法和动态方法在计算偏离临界较 增殖因子 ff定义如下：系统内中子产生率除以 远系统时的差别，分别