RMC在进行临界计算时,材料卡顺序发生变化,导致keff结果不同
如 92234.81c -1.32178481E-03
92235.81c -1.74035000E-01
92238.81c -7.05833089E-01
8016.81c -1.18527358E-01
8017.81c -4.51478478E-05
8018.81c -2.37620252E-04
改成 ,计算结果发生变化,请问这是什么原因导致的
8016.81c -1.18527358E-01
8017.81c -4.51478478E-05
8018.81c -2.37620252E-04
92234.81c -1.32178481E-03
92235.81c -1.74035000E-01
92238.81c -7.05833089E-01
请给出计算条件和keff及其标准偏差的值
计算条件是指哪些?
并行核数、计算粒子数、非活跃代和活跃代数等
如果输入卡可以提供的,建议附上输入卡
进程数:8;
线程数:10;
粒子数:10000;
总代数:500;
非活跃代数:250;
材料卡为:
92234.81c -1.32178481E-03
92235.81c -1.74035000E-01
92238.81c -7.05833089E-01
8016.81c -1.18527358E-01
8017.81c -4.51478478E-05
8018.81c -2.37620252E-04
所计算的keff为1.028815,标准偏差为0.000343;
材料卡为:
8016.81c -1.18527358E-01
8017.81c -4.51478478E-05
8018.81c -2.37620252E-04
92234.81c -1.32178481E-03
92235.81c -1.74035000E-01
92238.81c -7.05833089E-01
所计算的keff为1.029718,标准偏差为0.000348;
两者计算结果差别不大,但材料卡书写顺序的不同不应该对keff的计算结果产生影响;
输入卡为本人目前正在研究的内容,不方便给出,希望理解。以上是我使用RMC程序的问题反馈,并非求助。
问题的规模是?单棒,组件,还是全堆?
由于RMC并没有对材料顺序进行排序,所以实际计算碰撞核过程中,是按照用户输入顺序判断的,因此计算结果不一致是正常的,而你这个算例的偏差接近3sigma,怀疑可能有粒子数偏少的原因。
对于一个算例计算,建议每代粒子数相对多一些,而不是代数多一些