The use of the LOW-PWSX card in FLUKA (FLAIR)

[32119110060 (jiefeng zhong)]

Versions

Please provide the used software versions.

FLUKA:4-5.0
Flair:3.4-3
Operating system:

Description

Dear Professor,
I am conducting low-energy neutron simulation with polyethylene configured as the moderator. My question pertains to the use of the LOW-PWSX card in FLUKA (FLAIR) – I was unable to load the Point-wise Materials successfully. As for my settings:

 after running the corresponding inp file, h_Polyethy was not loaded in the generated out file, with h_water being called by default. 

Nevertheless,

h_Polyethy can be loaded successfully with my alternative settings.

I am uncertain whether this method is valid and if it can accurately reflect the moderation effect of polyethylene.

Input files

The one with problem when loading the neutron Point Wise materials

TITLE

GLOBAL          1000         0         0         0         1       0.0
GEOBEGIN                                                              COMBNAME
0    0
SPH s1         0.0 0.0 0.0 1000
SPH s2         0.0 0.0 0.0 800
XYP c1d        0.0
XYP c1u        151.5
YZP c1l        0.0
YZP c1r        136.0
XZP c1b        0.0
XZP c1f        140.2
XYP c2d        4.0
XYP c2u        147.5
YZP c2l        4.0
YZP c2r        132.0
XZP c2b        4.0
XZP c2f        136.2
RCC tube1      136.0 70.1 75.7 -38.4 0.0 0.0 4.0
RCC tube2      136.0 70.1 75.7 -38.4 0.0 0.0 10.0
RCC tube3      97.6 70.1 75.7 -17.6 0.0 0.0 4.0
RCC tube4      69.1 140.2 75.7 0.0 -39.0 0.0 4.0
RCC tube5      69.1 140.2 75.7 0.0 -39.0 0.0 10.0
RCC tube6      69.1 101.2 75.7 0.0 -26.2 0.0 4.0
XYP c3d        39.0
XYP c3u        112.5
YZP c3l        43.6
YZP c3r        97.6
XZP c3b        37.2
XZP c3f        101.2
XYP c4d        46.0
XYP c4u        105.5
YZP c4l        49.0
YZP c4r        91.1
XZP c4b        43.7
XZP c4f        94.7
XYP c5d        56.7
XYP c5u        94.7
YZP c5l        60.0
YZP c5r        80.0
XZP c5b        65.0
XZP c5f        75.0
XYP c6d        70.0
XYP c6u        80.0
YZP c6l        156.0
YZP c6r        166.0
XZP c6b        65.0
XZP c6f        75.0
XYP c7d        71.0
XYP c7u        81.0
YZP c7l        64.0
YZP c7r        74.0
XZP c7b        160.0
XZP c7f        170.0
YZP mod1_l     80.2
YZP mod2_l     90.2
YZP mod3_l     100.2
YZP mod4_l     110.2
YZP mod5_l     120.2
YZP mod1_r     80.2
YZP mod2_r     90.2
YZP mod3_r     100.2
YZP mod4_r     110.2
YZP mod5_r     120.2
TRC cone1      136 70.1 75.7 -70.9 0.0 0.0 4.0 0.0
RCC det1       80.0 70.1 75.7 0.1 0.0 0.0 4.0
RCC det2       90.0 70.1 75.7 0.1 0.0 0.0 4.0
RCC det3       100.0 70.1 75.7 0.1 0.0 0.0 4.0
RCC det4       110.0 70.1 75.7 0.1 0.0 0.0 4.0
RCC det5       120.0 70.1 75.7 0.1 0.0 0.0 4.0
RCC det6       130.0 70.1 75.7 0.1 0.0 0.0 4.0
RCC det7       136.0 70.1 75.7 -0.1 0.0 0.0 4.0
RCC det8       156.0 70.1 75.7 1.0 0.0 0.0 10.0
RCC inside     136.0 70.1 75.7 -56 0.0 0.0 4.0
RCC inside1    136.0 70.1 75.7 -56 0.0 0.0 4
END
REG0         5 +s1 -s2
REG1         5 (+c2u -c2d +c2r -c2l +c2f -c2b) -(+c3u -c3d +c3r -c3l +c3f -c3b) -tube2 -tube5
REG2         5 (+c1u -c1d +c1r -c1l +c1f -c1b) -(+c2u -c2d +c2r -c2l +c2f -c2b) -tube2 -tube5
REG3         5 +s2 -(+c1u -c1d +c1r -c1l +c1f -c1b) -det8
REG4         5 +tube2 -tube1
#if 0
REG5         5 +tube1 -det3 -det4 -det5 -det6 -det7 -(+mod3_r +cone1 -mod3_l) -(+mod4_r +cone1 -mod4_l) -(+mod5_r
+cone1 -mod5_l)
#endif
REG6         5 (+c4u -c4d +c4r -c5l +c4f -c4b) -(c5u -c5d +c5r -c5l +c5f -c5b) -tube6 -tube3
REG11        5 +c5u -c5d +c5r -c5l +c5f -c5b
REG7         5 (+c3u -c3d +c3r -c3l +c3f -c3b) -(+c4u -c4d +c4r -c4l +c4f -c4b) -tube3 -tube6
#if 0
REG8         5 +tube3 -det1 -det2 -(+mod1_r +cone1 -mod1_l) -(+mod2_r +cone1 -mod2_l)
#endif
REG9         5 +tube5 -tube4
REG10        5 +tube4
REG13        5 +tube6
REG15        5 +det1
REG16        5 +det2
REG17        5 +det3
REG18        5 +det4
REG19        5 +det5
REG20        5 +det6
REG21        5 +det7
REG22        5 +det8
REG24        5 (+c4u -c4d -c4l +c5l +c4f -c4b)
REG12        5 +mod1_r +cone1 -mod1_l
REG14        5 +mod2_r +cone1 -mod2_l
REG23        5 +mod3_r +cone1 -mod3_l
REG25        5 +mod4_r +cone1 -mod4_l
REG26        5 +mod5_r +cone1 -mod5_l
REG27        5 +inside1-inside-det1-det2-det3-det4-det5-det6-det7
REG28        5 inside-det1-det2-det3-det4-det5-det6-det7
END
GEOEND
MATERIAL                            0.94                              Polyethy
COMPOUND         4.0  HYDROGEN       2.0    CARBON                    Polyethy
MATERIAL                           1.105                              D2O
COMPOUND         1.0      O-16       2.0       H-2                    D2O
MATERIAL           1                 1.0                             2H-2
MATERIAL           8                 1.0                            16O-16
MATERIAL                           1.157                              BOLENE
MATERIAL           5                2.37                              BORON
COMPOUND       -4.62    CARBON    -15.38     BORON     -80.0  POLYETHYBOLENE
MATERIAL                            7.85                              Q235b
MATERIAL          26                7.86                            56IRON
MATERIAL          14                2.57                            28SILICON
MATERIAL          25                7.21                            55MANGANES
MATERIAL          16                1.96                            32SULFUR
MATERIAL          15               1.823                            31PHOSPHO
COMPOUND      -97.96      IRON      -0.2    CARBON     -0.35   SILICONQ235b
COMPOUND        -1.4  MAGNESIU    -0.045   PHOSPHO    -0.045    SULFURQ235b
MATERIAL           6                2.26                              Graphite
ASSIGNMA    BLCKHOLE      REG0
ASSIGNMA      BOLENE      REG1
ASSIGNMA       Q235b      REG2
ASSIGNMA         AIR      REG3
ASSIGNMA    Polyethy      REG4
*ASSIGNMA         AIR      REG5
ASSIGNMA    Polyethy      REG6
 ASSIGNMA       Q235b      REG7
*ASSIGNMA         AIR      REG8
ASSIGNMA    Polyethy      REG9
ASSIGNMA         AIR     REG10
ASSIGNMA         AIR     REG13
ASSIGNMA         AIR     REG15
ASSIGNMA         AIR     REG11
ASSIGNMA         AIR     REG22
ASSIGNMA         AIR     REG16
ASSIGNMA         AIR     REG17
ASSIGNMA         AIR     REG18
ASSIGNMA         AIR     REG19
ASSIGNMA         AIR     REG20
ASSIGNMA         AIR     REG21
ASSIGNMA    Polyethy     REG24
ASSIGNMA        IRON     REG12
ASSIGNMA        IRON     REG14
ASSIGNMA        IRON     REG23
ASSIGNMA        IRON     REG25
ASSIGNMA        IRON     REG26
ASSIGNMA    Graphite     REG27
ASSIGNMA         AIR     REG28
BEAM          -0.014             10000.0       0.2         0      -1.0NEUTRON
BEAMPOS         65.1      70.1      75.7
SOURCE
LOW-PWXS
LOW-PWXS            h_polyethy            Polyethy  Polyethy
USERWEIG                               1
*USRBIN            10   NEUTRON       -21       160       160       160ntry1
*USRBIN             0       0.0       0.0       160       160       160 &
*USRBIN            10   NEUTRON       -28       157        80        85ntry8
*USRBIN           156        60        65         1        20        20 &
*USRBIN            10   NEUTRON       -29        79       161        86ntry9
*USRBIN            59       160        66        20         1        20 &
USRBDX           101   NEUTRON       -30     REG11     REG15       1.0ntry10
USRBDX         0.014       0.0     100.0                               &
*USRBDX           101   NEUTRON       -31      REG7      REG8     690.8ntry11
*USRBDX         0.014       0.0     100.0                               &
*USRBDX           101   NEUTRON       -32      REG4      REG5     690.8ntry10
*USRBDX         0.014       0.0     100.0                               &
*USRBDX           101   NEUTRON       -33      REG8      REG6     690.8ntry10
*USRBDX         0.014       0.0     100.0                               &
*USRBDX           101   NEUTRON       -34      REG8      REG7     690.8ntry10
*USRBDX         0.014       0.0     100.0                               &
*USRBDX           101   NEUTRON       -35      REG5      REG4     690.8ntry10
*USRBDX         0.014       0.0     100.0                               &
RANDOMIZ           1         2
START           1000
STOP

The one load successfully

TITLE

GLOBAL          1000         0         0         0         1       0.0
GEOBEGIN                                                              COMBNAME
0    0
SPH s1         0.0 0.0 0.0 1000
SPH s2         0.0 0.0 0.0 800
XYP c1d        0.0
XYP c1u        151.5
YZP c1l        0.0
YZP c1r        136.0
XZP c1b        0.0
XZP c1f        140.2
XYP c2d        4.0
XYP c2u        147.5
YZP c2l        4.0
YZP c2r        132.0
XZP c2b        4.0
XZP c2f        136.2
RCC tube1      136.0 70.1 75.7 -38.4 0.0 0.0 4.0
RCC tube2      136.0 70.1 75.7 -38.4 0.0 0.0 10.0
RCC tube3      97.6 70.1 75.7 -17.6 0.0 0.0 4.0
RCC tube4      69.1 140.2 75.7 0.0 -39.0 0.0 4.0
RCC tube5      69.1 140.2 75.7 0.0 -39.0 0.0 10.0
RCC tube6      69.1 101.2 75.7 0.0 -26.2 0.0 4.0
XYP c3d        39.0
XYP c3u        112.5
YZP c3l        43.6
YZP c3r        97.6
XZP c3b        37.2
XZP c3f        101.2
XYP c4d        46.0
XYP c4u        105.5
YZP c4l        49.0
YZP c4r        91.1
XZP c4b        43.7
XZP c4f        94.7
XYP c5d        56.7
XYP c5u        94.7
YZP c5l        60.0
YZP c5r        80.0
XZP c5b        65.0
XZP c5f        75.0
XYP c6d        70.0
XYP c6u        80.0
YZP c6l        156.0
YZP c6r        166.0
XZP c6b        65.0
XZP c6f        75.0
XYP c7d        71.0
XYP c7u        81.0
YZP c7l        64.0
YZP c7r        74.0
XZP c7b        160.0
XZP c7f        170.0
YZP mod1_l     80.2
YZP mod2_l     90.2
YZP mod3_l     100.2
YZP mod4_l     110.2
YZP mod5_l     120.2
YZP mod1_r     80.2
YZP mod2_r     90.2
YZP mod3_r     100.2
YZP mod4_r     110.2
YZP mod5_r     120.2
TRC cone1      136 70.1 75.7 -70.9 0.0 0.0 4.0 0.0
RCC det1       80.0 70.1 75.7 0.1 0.0 0.0 4.0
RCC det2       90.0 70.1 75.7 0.1 0.0 0.0 4.0
RCC det3       100.0 70.1 75.7 0.1 0.0 0.0 4.0
RCC det4       110.0 70.1 75.7 0.1 0.0 0.0 4.0
RCC det5       120.0 70.1 75.7 0.1 0.0 0.0 4.0
RCC det6       130.0 70.1 75.7 0.1 0.0 0.0 4.0
RCC det7       136.0 70.1 75.7 -0.1 0.0 0.0 4.0
RCC det8       156.0 70.1 75.7 1.0 0.0 0.0 10.0
RCC inside     136.0 70.1 75.7 -56 0.0 0.0 4.0
RCC inside1    136.0 70.1 75.7 -56 0.0 0.0 4
END
REG0         5 +s1 -s2
REG1         5 (+c2u -c2d +c2r -c2l +c2f -c2b) -(+c3u -c3d +c3r -c3l +c3f -c3b) -tube2 -tube5
REG2         5 (+c1u -c1d +c1r -c1l +c1f -c1b) -(+c2u -c2d +c2r -c2l +c2f -c2b) -tube2 -tube5
REG3         5 +s2 -(+c1u -c1d +c1r -c1l +c1f -c1b) -det8
REG4         5 +tube2 -tube1
#if 0
REG5         5 +tube1 -det3 -det4 -det5 -det6 -det7 -(+mod3_r +cone1 -mod3_l) -(+mod4_r +cone1 -mod4_l) -(+mod5_r
+cone1 -mod5_l)
#endif
REG6         5 (+c4u -c4d +c4r -c5l +c4f -c4b) -(c5u -c5d +c5r -c5l +c5f -c5b) -tube6 -tube3
REG11        5 +c5u -c5d +c5r -c5l +c5f -c5b
REG7         5 (+c3u -c3d +c3r -c3l +c3f -c3b) -(+c4u -c4d +c4r -c4l +c4f -c4b) -tube3 -tube6
#if 0
REG8         5 +tube3 -det1 -det2 -(+mod1_r +cone1 -mod1_l) -(+mod2_r +cone1 -mod2_l)
#endif
REG9         5 +tube5 -tube4
REG10        5 +tube4
REG13        5 +tube6
REG15        5 +det1
REG16        5 +det2
REG17        5 +det3
REG18        5 +det4
REG19        5 +det5
REG20        5 +det6
REG21        5 +det7
REG22        5 +det8
REG24        5 (+c4u -c4d -c4l +c5l +c4f -c4b)
REG12        5 +mod1_r +cone1 -mod1_l
REG14        5 +mod2_r +cone1 -mod2_l
REG23        5 +mod3_r +cone1 -mod3_l
REG25        5 +mod4_r +cone1 -mod4_l
REG26        5 +mod5_r +cone1 -mod5_l
REG27        5 +inside1-inside-det1-det2-det3-det4-det5-det6-det7
REG28        5 inside-det1-det2-det3-det4-det5-det6-det7
END
GEOEND
MATERIAL                            0.94                              Polyethy
COMPOUND         4.0  HYDROGEN       2.0    CARBON                    Polyethy
MATERIAL                           1.105                              D2O
COMPOUND         1.0      O-16       2.0       H-2                    D2O
MATERIAL           1                 1.0                             2H-2
MATERIAL           8                 1.0                            16O-16
MATERIAL                           1.157                              BOLENE
MATERIAL           5                2.37                              BORON
COMPOUND       -4.62    CARBON    -15.38     BORON     -80.0  POLYETHYBOLENE
MATERIAL                            7.85                              Q235b
MATERIAL          26                7.86                            56IRON
MATERIAL          14                2.57                            28SILICON
MATERIAL          25                7.21                            55MANGANES
MATERIAL          16                1.96                            32SULFUR
MATERIAL          15               1.823                            31PHOSPHO
COMPOUND      -97.96      IRON      -0.2    CARBON     -0.35   SILICONQ235b
COMPOUND        -1.4  MAGNESIU    -0.045   PHOSPHO    -0.045    SULFURQ235b
MATERIAL           6                2.26                              Graphite
ASSIGNMA    BLCKHOLE      REG0
ASSIGNMA      BOLENE      REG1
ASSIGNMA       Q235b      REG2
ASSIGNMA         AIR      REG3
ASSIGNMA    Polyethy      REG4
*ASSIGNMA         AIR      REG5
ASSIGNMA    Polyethy      REG6
ASSIGNMA       Q235b      REG7
*ASSIGNMA         AIR      REG8
ASSIGNMA    Polyethy      REG9
ASSIGNMA         AIR     REG10
ASSIGNMA         AIR     REG13
ASSIGNMA         AIR     REG15
ASSIGNMA         AIR     REG11
ASSIGNMA         AIR     REG22
ASSIGNMA         AIR     REG16
ASSIGNMA         AIR     REG17
ASSIGNMA         AIR     REG18
ASSIGNMA         AIR     REG19
ASSIGNMA         AIR     REG20
ASSIGNMA         AIR     REG21
ASSIGNMA    Polyethy     REG24
ASSIGNMA        IRON     REG12
ASSIGNMA        IRON     REG14
ASSIGNMA        IRON     REG23
ASSIGNMA        IRON     REG25
ASSIGNMA        IRON     REG26
ASSIGNMA    Graphite     REG27
ASSIGNMA         AIR     REG28
BEAM          -0.014             10000.0       0.2         0      -1.0NEUTRON
BEAMPOS         65.1      70.1      75.7
SOURCE
LOW-PWXS
LOW-PWXS      10010.h_polyethy            HYDROGEN  HYDROGEN          JEFF-3.3
USERWEIG                               1
*USRBIN            10   NEUTRON       -21       160       160       160ntry1
*USRBIN             0       0.0       0.0       160       160       160 &
*USRBIN            10   NEUTRON       -28       157        80        85ntry8
*USRBIN           156        60        65         1        20        20 &
*USRBIN            10   NEUTRON       -29        79       161        86ntry9
*USRBIN            59       160        66        20         1        20 &
USRBDX           101   NEUTRON       -30     REG11     REG15       1.0ntry10
USRBDX         0.014       0.0     100.0                               &
*USRBDX           101   NEUTRON       -31      REG7      REG8     690.8ntry11
*USRBDX         0.014       0.0     100.0                               &
*USRBDX           101   NEUTRON       -32      REG4      REG5     690.8ntry10
*USRBDX         0.014       0.0     100.0                               &
*USRBDX           101   NEUTRON       -33      REG8      REG6     690.8ntry10
*USRBDX         0.014       0.0     100.0                               &
*USRBDX           101   NEUTRON       -34      REG8      REG7     690.8ntry10
*USRBDX         0.014       0.0     100.0                               &
*USRBDX           101   NEUTRON       -35      REG5      REG4     690.8ntry10
*USRBDX         0.014       0.0     100.0                               &
RANDOMIZ           1         2
START           1000
STOP

[kltaylor (Kate Louise Taylor)]

Dear Jiefeng,

Thank you for your question.

In the first case you present the reason it is not showing h_polyethly and instead hydrogen is due to how the LOW-PWXS card works. The thermal scattering acts on an isotope basis and not on a compound basis and so the input card is searching for an isotope and cannot find one so is falling back on the default.

>= 0 - 7.22.44. LOW-PWXS — FLUKA Manual

And so your second card is quite correct as you are directing FLUKA to use to hydrogen in the h_polyethy data file.

If you want to add different hydrogen materials with different thermal scattering treatments you can use different FLUKA material definitions and different names.

I hope that helps,

Katie