How to rotating wide beam correctly?

Source Definition
1

Dear FLUKA experts,

In the x-z coordinate system, the angle between the direction of proton emission and positive z-axis is 80°. The radius of the source is 0.6 cm. The projection in the x-y coordinates is going to be an ellipse. However, the shape of the energy deposit calculated by USRBIN is a circle, which should be a ellipse. The input file and picture are below:

2021-10-08 17-18-32屏幕截图
2021-10-08 17-18-32屏幕截图2560×1440 221 KB

NI64-10.inp (4.0 KB)

So, what’s the problem about this phenomenon?

Thanks and best regards,
Xiaohe Wang

2

Your plot is in the x-y plane and reflects the annular beam shape you asked for in the BEAM card.

3

Dear Xiaohe,

Slide 15-17 of the following lecture may explain the issue you are facing:

Cheers,
David

4

Thank you very much .
According to the pdf, BEAMAXES is added. However, there is a mistake in FLAIR showed in the picture.

1
12560×1440 366 KB

So how to solve this problem?

5

Dear Xiaohe,

you need to input the numbers with more decimal places, as shown on the bottom of Slide 17.

Cheers,
David

6

Dear Cerutti:
Thank you very much.
But, in x-y plane, the shape of beam should be ellipse, not cycle.
I think, the source setting parameters have some mistake. May be BEAMAXES should be added. But, there are some problem. Please look at the picture in the next reply.

Thanks and best regards,
Xiaohe Wang

7

Why? You defined an annular beam in the x-y plane and this is what you get. Then, you have defined the beam direction at 80 deg in the x-z plane, therefore the circular annulus moves (not much, due to the low energy) along it, but the x-y shape remains circular as it should (as discussed on slide 15).
If you meant instead to get a circular beam shape in the plane perpendicular to the beam direction (that is not xy), such as to reflect into an ellipse on xy you expect, you should indeed make proper use of BEAMAXES as indicated by David.

8

Dear horvathd,
Yes, you are right.
But, in the inputfile, the length of what(5) and what(6) is 10 bytes. If decimal numbers less than 10 bytes, there is also a mistake.

Thanks and best regards,
Xiaohe Wang

9

Dear Horvathd,

The BEAMAXES card is added. But the results is still not right. Look at the next picture, when used the BEAMAXES, the radius of the beam become less than 0.6 cm. So, what’s the problem about this phenomenon?

4
42560×1440 186 KB

NI64-10.inp (4.1 KB)

Thanks and best regards,
Xiaohe Wang

10

Dear Xiaohe,

the issue is that you set the beam direction along the Y axis on the BEAMPOS card. So, the particles will propagate sideways from the starting surface.

When you are using BEAMAXES and you want to keep the beam direction perpendicular to the beam’s shape, they you should not change the beam direction on the BEAMPOS card.

Cheers,
David

13

Thank you very much.
Theo problem is solved.

Best regards,
Xiaohe Wang

14

Thank you very much.

Theo problem is solved.

Best regards,
Xiaohe Wang

15