Xe135m microscopic cross section calculation

[motalab]

Hi Jaakko

Now we are using Serpent 2.1.31 and also use

set poi OPT [ VR XE135M ] to get the poison cross-section. We set the parameter 1 to treat Xe135 and Xe135m separately. Now in the output file, it is given microscopic and macroscopic cross-section and yield of Xe135 and Xe135m separately. But problem is that the microscopic and macroscopic cross-section of Xe135m is always given 0.000 for every burnup step.



INF_I135_YIELD (idx, [1: 4]) = [ 6.78544958E-02 1.4E-06 6.45136062E-02 3.3E-07 ];

INF_XE135_YIELD (idx, [1: 4]) = [ 6.02531656E-04 5.0E-05 2.07678679E-03 8.9E-06 ];

INF_XE135M_YIELD (idx, [1: 4]) = [ 1.35110647E-03 5.4E-05 4.93470490E-03 9.1E-06 ];

INF_I135_MICRO_ABS (idx, [1: 4]) = [ 1.32281258E+00 0.00019 2.66299691E+01 6.3E-05 ];

INF_XE135_MICRO_ABS (idx, [1: 4]) = [ 1.84533621E+02 0.00027 1.10822614E+06 7.2E-05 ];

INF_XE135M_MICRO_ABS (idx, [1: 4]) = [ 0.00000000E+00 0.0E+00 0.00000000E+00 0.0E+00 ];

INF_XE135_MACRO_ABS (idx, [1: 4]) = [ 1.76331260E-08 0.00027 1.05474452E-04 7.2E-05 ];

INF_XE135M_MACRO_ABS (idx, [1: 4]) = [ 0.00000000E+00 0.0E+00 0.00000000E+00 0.0E+00 ];


What is the main cause of this problem?

Regards

Motalab

[Jaakko Leppänen]

Try adding if-branch:

Code: Select all

  else if ((long)RDB[loc1 + NUCLIDE_ZAI] == 541351)
	    WDB[mat + MATERIAL_PTR_XE135M_ISO] = (double)loc0;

after branch:

Code: Select all

  else if ((long)RDB[loc1 + NUCLIDE_ZAI] == 541350)
	    WDB[mat + MATERIAL_PTR_XE135_ISO] = (double)loc0;

in putcompositions.c (around line 179).