.. _section_ifp:
反复裂变几率法计算
================================
RMC反复裂变法可以在临界和燃耗计算模式下计算动态参数和几何微扰系数,其中几何微扰
系数仅用于连续能量截面,动态参数除此之外还可应用于多群截面。
反复裂变几率法模块输入卡
-----------------------------
.. code-block:: none
Adjoint
blocksize=
fisbranch=
OUTPUTINTERVAL=
kinetic
AdjointFluxMesh [Energy=] [Scope=] [Bound=]
GeoAdjointTally surf= type= parameters= constraint= universe= cellmat= cutoffcos=
其中,
- **Adjoint**\ 为反复裂变几率法模块关键词;
- **blocksize**\ 指定多少活跃代为一个块,默认为10;
- **fisbranch**\ 指定每次裂变能产生几个裂变中子,默认为1,数值越大,反复裂变
几率法初始分配内存越大;
- **OUTPUTINTERVAL**\ 指定每多少代输出伴随计算结果,默认所有活跃代结束才输出;
- **kinetic**\ 计算动态参数关键词;
- **AdjointFluxMesh**\ 为反复裂变几率法求解临界伴随通量输入卡的关键词。该选项卡不能和kinetic、GeoAdjointTally以及敏感性和不确定度分析功能一起使用。
- **Energy**\ 选项卡指定伴随通量的能量分界点。
- **Scope**\ 选项卡指定网格在x,y,z方向的数量。特别地,参数为“-1”表示该方向上
只有一层无限大网格 (注意:在Universe重复几何中的Scope选项卡当中,参数为1表示
该方向上只有一层无限大网格)。
- **Bound**\ 选项卡指定网格在x,y,z方向的边界范围,形如“Bound = x_min
x_max y_min y_max z_min z_max”。若某方向只有一层网格,
\ **Bound**\ 选项卡中对应的参数没有实际意义。
- **GeoAdjointTally**\ 为几何微扰输入卡的关键词。
- **surf**\ 计算几何微扰系数对应面在Surf卡中编号;
- **type**\ 几何形变类型。1:平移,2:均匀扩张,3:单向扩张,4:旋转;
- **parameters**\.描述几何形变的参数。平移需要三个参数,即指向平移方向的归一化向量;均匀扩张不需要参数;单向扩张需要三个参数,即指向扩张方向的归一化向量;旋转需要7个参数,分别是旋转轴的归一化方向向量,旋转不动点的坐标,和旋转方向(符合右手定则为正1,不符合为负1)。
- **constraint**\ 约束栅元。仅对约束栅元内的曲面进行扰动。输入格式通cell输入卡相同。
- **universe**\ 约束栅元所在的空间。
- **cellmat**\ 用于限制面的扰动部分,正数表示栅元,负数表示材料。均为程序内部编号。
- **cutoffcos**\ 穿面余弦的截断值。
动态参数输出文件
---------------------
输出结果给出以下动态参数结果的平均值和标准差:有效中子代时间,单位为微秒;rossi-alpha,
单位为1/微秒;有效缓发中子份额,单位为1。
.. code-block:: none
------------------ Current Cycle Number = 50 ------------------
------------------ kinetic parameters, gen. time, rossi-alpha, beta-eff ------------------
Ave std. dev.
gen. time (us) 2.39916E+01 1.23582E+00
rossi-alpha (us-1) -3.32090E-04 1.19101E-04
beta-eff 8.47770E-03 3.27803E-03
几何微扰系数输出文件
-------------------------
输出结果给出几何微扰系数结果的平均值和相对标准偏差:输运项;散射和裂变项。
总结果为两者之和,如实例中几何微扰系数为0.013019。
.. code-block:: none
------------------ Current Cycle Number = 300 ------------------
------------------ geometrical perturbation parameters ------------------
Ave Re
transport: 2.74153E-04 5.18707E-01
scatter + fission: 1.27445E-02 8.43237E-01