FLUKA comparison with PHITS

Hi Volodymyr

Thank you for your help. I’m modelling a lightly more complex problem now with more complex source and geometry. I’m also getting discrepancy between PHITS and FLUKA and I dont suppose you have any suggestions as to what could be causing the difference?

screenshot13

snre_point_output.out (33.0 KB)
snre_spallation.inp (37.5 KB)
snre_gcr_spherical_source_21_tab.lis (87.8 KB)
snre_gcr_spherical_source_21_sum.lis (149.9 KB)
snre_gcr_spherical_source.flair (20.9 KB)

Best wishes,

Emma

Hi Emma,

First of all, let’s try to optimize your source definition.
In your case, I would suggest you to utilize source_newgen.f routine. Usage is described here:

Examples are from: 2023-fluka-beginner-course

There you can define shape, directions and file-based energy spectra.
I am attaching a modified example from this course.
histogram.txt (239 Bytes)
source_newgen.f (19.2 KB)
source_routine.flair (2.9 KB)
In addition, I included a uniform spherical shell generator into source_library.inc, so you need to replace it in /path/to/fluka/include
source_library.inc (42.5 KB)


image

Let me know if it works for you and then I will continue with further comparison.

Best wishes,
Volodymyr

Hi Volodymyr

I ran the model, and these are the plots I got


which look similar to yours.

Best wishes,

Emma

Hi Emma,

Then you can similarly define this source in your Fluka model.

I have obtained similar results with your geometry input and my sample histogram of protons.

Sample distribution 1-100 GeV protons:
image

Comparison of FLUKA (without PWXS/LOW-MAT cards, they activated by default) and PHITS (default physics in both codes otherwise).
image

One would need to investigate small differences in 1-100 MeV energy region with the PHITS team to understand the cause.

Best wishes,
Volodymyr