Energy distribution and angle distribution

Scoring
1

Dear Developer
I want to use the Fluka program to calculate the energy and angle distribution of protons bombarding lithium targets. The following code gives the code I calculated using the mcnp program. How do I use the fluka program to complete the calculation of F21, F31, and F41 cards ?

c cell card
100 1 -0.534 1 -2 -3
200 0 -5 #100
300 0 5

c surface card
1 px 0
2 px 0.03
3 cx 8
4 px -10
5 so 100
30 kx 0 0.00765426624555 1
31 kx 0 0.0717967697 1
32 kx 0 0.217442832 1
33 kx 0 0.49029 1
34 kx 0 1 1
35 kx 0 2.03960673 1
36 kx 0 4.59891 1
37 kx 0 13.9282 1
38 kx 0 130.64609564 1
39 kx 0 130.64609564 -1
40 kx 0 13.9282 -1
41 kx 0 4.59891 -1
42 kx 0 2.03960673 -1
43 kx 0 1 -1
44 kx 0 0.49029 -1
45 kx 0 0.217442832 -1
46 kx 0 0.0717967697 -1
47 kx 0 0.00765426624555 -1
48 kx 0 0 -1

c date
mode h n
imp:h 1 1 0
imp:n 1 1 0
cut:h j 1e-3
sdef sur=4 pos=-10 0 0 rad=d1 vec=1 0 0 dir=1 erg=2.5 par=h
si1 0 5
m1 3007 0.925
3006 0.075
F11:n 5
F21:n 5
FT21 frv 1 0 0
*c21 0 5 15 25 35 45 55 65 75 85 95 105 115 125 135 145 155 165 175 T
F31:n 5
FS31 -30 -31 -32 -33 -34 -35 -36 -37 -38 39 40 41 42 43 44 45 46 47 48 T
F41:n 5
E41 0 999i 5
NPS 1E8

2

Dear Kaiwen,

in FLUKA you score surface current with the USRBDX card.

For F21, in FLUKA you can set up angular binning directly on the card, but in solid angles not in degrees, relative to the surface normal vector.
For F31, you will have to create separate scoring for all segments.
For F41, you can set the energy binning on the card itself.

Cheers,
David

1 Like
3

Dear David,
Thank you very much for your reply, F31 and F41 card calculation, I understand what you mean, F21 card implementation I still do not understand, the following is my Fluka program settings, I use 30 BIN successfully calculated the energy distribution, but I use 40 BIN calculation of angular distribution when there are some problems, only one data is not 0, other data is 0.
input

output
output887Ă—687 67.1 KB

li.flair (5.1 KB)
li.inp (1.8 KB)

4

Dear,

Kaiwen,

you only get result in only bin because you are scoring on the surface of a sphere, so the normal and the particle direction is almost parallel.

You need to score the on the surface between the TARGET and VOID, to get the angular distribution. You may need to separate the VOID into a “upstream” and “downstream” section, and use separate scorings to distinguish between the neutrons going forward and the backscattered ones.

Cheers,
David

1 Like
5

Thank you very much for your answer, probably because I just started using Fluka, so I don 't understand what you mean. If you can, can you help me to change the input card, so that I can understand, thank you very much.

6

Dear Kaiwen,

please, find the updated Flair project file here: li_dh.flair (3.7 KB)

Cheers,
David

7

Dear David
I understand what you mean, thank you very much.

image
image762Ă—620 64.5 KB

I also want to ask a question, the result file has two parts of the results, part is 2pi, part is 90 °, what do they mean ?

8

Dear Kaiwen,

one is showing the results with solid angles, the other now is the same, but the solid angles have been converted to degrees.

Cheers,
David

9

Is the calculation result of this solid angle corresponding to the calculation result of the F31 card in mcnp ? Is the calculation result of the angle corresponding to the F21 card in the MCNP ?

10

Dear Kaiwen,

I don’t looked at the geometry segments you are using for F31, and it is not possible to reproduce the linear angular binning in degrees with the FLUKA scoring, like in F21.
FLUKA can calculate the angular dependency only with linear binning along the solid angle.

Cheers,
David

1 Like
11

Dear David
I understand what you mean, thank you very much.