[fluka-discuss]: Gas production using new FLUKA functionality in 2011.2x.8

From: Paola Sala <paola.sala_at_mi.infn.it>
Date: Tue, 21 Jan 2020 09:29:37 +0100

From: Kriti Gupta <kriti.gupta_at_demcon.com>
To: fluka discuss <fluka-discuss_at_fluka.org>
Subject: Gas production using new FLUKA functionality in 2011.2x.8
Date: Mon, 20 Jan 2020 14:28:46 +0000

Dear FLUKA team and users,

I recently asked a question about estimating light ions production/depositi=
on (H,D,He etc.) using USRBIN+RES-NUCL+AUXSCORE introduced in FLUKA 2011.2x=
.8 and comparison of the results with the traditional way of estimating iso=
topes (RESNUCL card). I found a huge discrepancy in the numbers obtained fr=
om the two ways. Please find attached the input and output files for the sa=
me. It would be of great help if someone could point out where this discrep=
ancy is coming from.

Thanking you,

Best regards,
Kriti Gupta





-filename="I19.inp"
TITLE

* Set the defaults for precision simulations
DEFAULTS
PRECISIO
*EMF-BIAS 0.001 LPBEMF
* Define the beam characteristics
BEAM -0.014 0.0 0.0 0.4 0.0 -1.PROTON
* Define the beam position
BEAMPOS 0.0 0.0 -5.
IONTRANS -2
PHOTONUC 1. HYDROGEN _at_LASTMAT
PHYSICS 3.
EVAPORAT
PHYSICS 1.
COALESCE
GEOBEGIN
COMBNAME
    0 0
* Black body
SPH blkbody 0.0 0.0 0.0 100000.
* Void sphere
SPH void 0.0 0.0 0.0 10000.
* Cylindrical target
RPP degrader -1. 1. -1. 1. 0.0 0.095
* Cylindrical target
RPP target3 -1. 1. -1. 1. 5.095 5.105
END
* Black hole
BLKBODY 5 +blkbody -void
* Void around
VOID 5 +void -degrader -target3
* Degrader
DEGRADER 5 +degrader
* Target
TARGET3 5 +target3
END
GEOEND
MATERIAL 6. 1.88
Graphite
LOW-MAT Graphite 6. -3. 296. CARBON
MAT-PROP 20. Graphite
DPA-ENER
* Molybdenum-100
MATERIAL 42. 10.22 100.MOLY100
LOW-MAT MOLY100 42. -2. 296.
MOLYBDEN
* Molybdenum-98
MATERIAL 42. 10.22 98.MOLY98
LOW-MAT MOLY98 42. -2. 296.
MOLYBDEN
* Molybdenum-97
MATERIAL 42. 10.22 97.MOLY97
LOW-MAT MOLY97 42. -2. 296.
MOLYBDEN
* Molybdenum-96
MATERIAL 42. 10.22 96.MOLY96
LOW-MAT MOLY96 42. -2. 296.
MOLYBDEN
* Molybdenum-95
MATERIAL 42. 10.22 95.MOLY95
LOW-MAT MOLY95 42. -2. 296.
MOLYBDEN
* Molybdenum-94
MATERIAL 42. 10.22 94.MOLY94
LOW-MAT MOLY94 42. -2. 296.
MOLYBDEN
* Molybdenum-92
MATERIAL 42. 10.22 92.MOLY92
LOW-MAT MOLY92 42. -2. 296.
MOLYBDEN
* Material data - Molybdenum composition
MATERIAL 10.2 MOLY
MAT-PROP 65. MOLY
DPA-ENER
COMPOUND -0.95 MOLY100 -0.0404 MOLY98 -0.0026 MOLY97MOLY
COMPOUND -0.0023 MOLY96 -0.002 MOLY95 -0.001 MOLY94MOLY
COMPOUND -0.0017 MOLY92 MOLY
* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7..
ASSIGNMA BLCKHOLE BLKBODY
ASSIGNMA VACUUM VOID
ASSIGNMA Graphite DEGRADER
ASSIGNMA MOLY TARGET3
RESNUCLE 3. -25. 10. TARGET3 0.04resnuct
RESNUCLE 3. -26. 10. DEGRADER 0.38resnucd
USRBIN 10. RES-NUCL -41. 1. 1. 5.105Hdept
USRBIN -1. -1. 5.095 1. 1. 1. &
AUXSCORE USRBIN PROTON Hdept
USRBIN 10. RES-NUCL -42. 1. 1. 0.095Hdepd
USRBIN -1. -1. 0.0 1. 1. 1. &
AUXSCORE USRBIN PROTON Hdepd
USRBIN 10. RES-NUCL -42. 1. 1. 0.095Hedepd
USRBIN -1. -1. 0.0 1. 1. 1. &
AUXSCORE USRBIN 4-HELIUM Hedepd
USRBIN 10. RES-NUCL -42. 1. 1. 0.095Ddepd
USRBIN -1. -1. 0.0 1. 1. 1. &
AUXSCORE USRBIN DEUTERON Ddepd
* Set the random number seed
RANDOMIZ 1.
* Set the number of primary histories to be simulated in the run
*START 300000000.
* Set the number of primary histories to be simulated in the run
START 300000.
STOP


filename: I19_25_sum.lis
**** I19_25.rnc ****


  Total primaries run: 3000000

    Total weight of the primaries run: 3000000.00


  Detector n: 1 resnuct
     (Region 4 Volume: 3.99999991E-02 cmc,
      distr. type : 3 ,
      Z_max: 10 , N-Z_max: 22 , N-Z_min: -4 )

     Tot. response (n/cmc/pr) 4.26583327E-02 +/- 1.39916456 %
     ( --> Nuclei/pr 1.70633325E-03 +/- 1.39916456 % )

  **** Isotope Yield as a function of Mass Number **** ****
(nuclei / cmc / pr) ****


  A_min: 1 - A_max: 42

  A: 1 4.2658333E-02 +/- 1.3991646E+00 %

  **** Isotope Yield as a function of Atomic Number **** ****
(nuclei / cmc / pr) ****


  Z_min: 1 - Z_max: 10

  Z: 1 4.2658333E-02 +/- 1.3991646E+00 %

  **** Residual nuclei distribution **** **** (nuclei / cmc / pr)
  ****



  A \ Z 1 2 3 4 5 6
  7 8 9 10
   1 4.27E-02 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00
0.00E+00 0.00E+00 0.00E+00 0.00E+00
     +/- 1.4 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/-
0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 %


:
 filename : I19_26_sum.lis
 **** I19_26.rnc ****


  Total primaries run: 3000000

    Total weight of the primaries run: 3000000.00


  Detector n: 1 resnucd
     (Region 3 Volume: 0.379999995 cmc,
      distr. type : 3 ,
      Z_max: 10 , N-Z_max: 6 , N-Z_min: -4 )

     Tot. response (n/cmc/pr) 2.82517541E-02 +/- 1.09933329 %
     ( --> Nuclei/pr 1.07356664E-02 +/- 1.09933329 % )

  **** Isotope Yield as a function of Mass Number **** ****
(nuclei / cmc / pr) ****


  A_min: 1 - A_max: 26

  A: 13 9.8245614E-05 +/- 1.1739930E+01 %
  A: 12 9.2956140E-03 +/- 1.0847383E+00 %
  A: 10 4.2982456E-05 +/- 1.1958092E+01 %
  A: 6 1.2280702E-05 +/- 2.4281045E+01 %
  A: 4 9.1228071E-03 +/- 1.8669147E+00 %
  A: 2 1.2280702E-05 +/- 3.4006802E+01 %
  A: 1 9.6675440E-03 +/- 1.0937303E+00 %

  **** Isotope Yield as a function of Atomic Number **** ****
(nuclei / cmc / pr) ****


  Z_min: 1 - Z_max: 10

  Z: 7 6.1403510E-06 +/- 3.0491068E+01 %
  Z: 6 9.3877192E-03 +/- 1.1333162E+00 %
  Z: 5 4.2982456E-05 +/- 1.1958092E+01 %
  Z: 3 1.2280702E-05 +/- 2.4281045E+01 %
  Z: 2 9.1228071E-03 +/- 1.8669147E+00 %
  Z: 1 9.6798247E-03 +/- 1.0875541E+00 %

  **** Residual nuclei distribution **** **** (nuclei / cmc / pr)
  ****



  A \ Z 1 2 3 4 5 6
  7 8 9 10
  13 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 9.21E-05
6.14E-06 0.00E+00 0.00E+00 0.00E+00
      +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/-12.2 %
+/-30.5 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
  12 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 9.30E-03
0.00E+00 0.00E+00 0.00E+00 0.00E+00
      +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 1.1 %
+/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
  10 0.00E+00 0.00E+00 0.00E+00 0.00E+00 4.30E-05 0.00E+00
0.00E+00 0.00E+00 0.00E+00 0.00E+00
      +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/-12.0 % +/- 0.0 %
+/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
   6 0.00E+00 0.00E+00 1.23E-05 0.00E+00 0.00E+00 0.00E+00
0.00E+00 0.00E+00 0.00E+00 0.00E+00
      +/- 0.0 % +/- 0.0 % +/-24.3 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
+/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
   4 0.00E+00 9.12E-03 0.00E+00 0.00E+00 0.00E+00 0.00E+00
0.00E+00 0.00E+00 0.00E+00 0.00E+00
      +/- 0.0 % +/- 1.9 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
+/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
   2 1.23E-05 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00
0.00E+00 0.00E+00 0.00E+00 0.00E+00
      +/-34.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
+/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
   1 9.67E-03 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00
0.00E+00 0.00E+00 0.00E+00 0.00E+00
      +/- 1.1 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 %
+/- 0.0 % +/- 0.0 % +/- 0.0 % +/- 0.0 %


 filename="I19_42.bnn.lis"; size=1896;

1
   Cartesian binning n. 1 "Hdepd " , generalized particle n. 251
      X coordinate: from -1.0000E+00 to 1.0000E+00 cm, 1 bins (
2.0000E+00 cm wide)
      Y coordinate: from -1.0000E+00 to 1.0000E+00 cm, 1 bins (
2.0000E+00 cm wide)
      Z coordinate: from 0.0000E+00 to 9.5000E-02 cm, 1 bins (
9.5000E-02 cm wide)
      Data follow in a matrix A(ix,iy,iz), format (1(5x,1p,10(1x,e11.4)))

      accurate deposition along the tracks requested
      this is a track-length binning
       9.5404E-07

      Percentage errors follow in a matrix A(ix,iy,iz), format
(1(5x,1p,10(1x,e11.4)))

       1.0963E+00
1
   Cartesian binning n. 2 "Hedepd " , generalized particle n. 251
      X coordinate: from -1.0000E+00 to 1.0000E+00 cm, 1 bins (
2.0000E+00 cm wide)
      Y coordinate: from -1.0000E+00 to 1.0000E+00 cm, 1 bins (
2.0000E+00 cm wide)
      Z coordinate: from 0.0000E+00 to 9.5000E-02 cm, 1 bins (
9.5000E-02 cm wide)
      Data follow in a matrix A(ix,iy,iz), format (1(5x,1p,10(1x,e11.4)))

      accurate deposition along the tracks requested
      this is a track-length binning
       8.9434E-07

      Percentage errors follow in a matrix A(ix,iy,iz), format
(1(5x,1p,10(1x,e11.4)))

       1.8249E+00

1
   Cartesian binning n. 3 "Ddepd " , generalized particle n. 251
      X coordinate: from -1.0000E+00 to 1.0000E+00 cm, 1 bins (
2.0000E+00 cm wide)
      Y coordinate: from -1.0000E+00 to 1.0000E+00 cm, 1 bins (
2.0000E+00 cm wide)
      Z coordinate: from 0.0000E+00 to 9.5000E-02 cm, 1 bins (
9.5000E-02 cm wide)
      Data follow in a matrix A(ix,iy,iz), format (1(5x,1p,10(1x,e11.4)))

      accurate deposition along the tracks requested
      this is a track-length binning
       1.2211E-09

      Percentage errors follow in a matrix A(ix,iy,iz), format
(1(5x,1p,10(1x,e11.4)))

       3.4017E+01


__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Tue Jan 21 2020 - 11:03:57 CET

This archive was generated by hypermail 2.3.0 : Tue Jan 21 2020 - 11:04:04 CET