Dear fluka experts,
I have set the Emin、Emax、andEbins in the USRBDX card, but the output is not what I set. There are two kind of energy boundaries in the sum.lis file, and there are 310 energy intervals in the tab.lis file while I set 100 in the USRBDX card. Why is this happening?
50+50.inp (6.4 KB)
sum.lis&tab.lis.zip (135.1 KB)
This can be explained with FLUKA’s special treatment for low energy neutrons (E < 20 MeV). In this energy range FLUKA uses a group wise treatment of the neutrons based on evaluated data.
To be able to score low energy neutrons correctly the same group energy bins has to be used, which FLUKA does automatically for you.
The number and with of bins you get are calculated as the following:
First the bins from the USRBDX card are calculated, the this list is cut in half at 20 MeV, and the lower energy part is replaced with the relevant part of the neutron energy groups.
I want to sample the output neutron spectrum，but the result is not the same as I expected, there is a big difference.I’m afraid it’s related to the neutron energy gap.
I’ll upload my file. If you have time, please check it for me.Thank you in advance.
input&output.zip (61.4 KB)
The drop at 20 MeV you see is due to that your source routine interprets the 3rd column of the spectrum.dat file as probability of the of the bin, while the FLUKA output if in the units of [1/GeV/primary].
If I plot the flu_43.tab.lis, with the value option: Y*DX (Multiplying the value of the bin, with the width of the bin), I get the unit of [1/primary], and this plot matches with the spectrum specified in your source data file.
Should I replace the 3rd column of the spectrcum.dat with the “cumul. flux (part/cmq/pr)” data in the output_22_sum.lis file? If so, what bothers me is that I didn’t find the same output format as the. tab.lis file with the unit of part/cmq/pr. Is there any way for me to get output in this format？
Your understanding is very correct and my doubts and problems have been solved. Thank you again.