有效增殖因子的动态模拟计算 calculation of effective multiplication factor with dynamic monte carlo method.pdfVIP

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