Re: AW: [fluka-discuss]: scoring neutron fluence & thermal neutrons

From: Mikhail Polkovnikov <pmk_at_ihep.ru>
Date: Tue, 16 May 2017 19:26:22 +0300

On 15.05.2017 14:04, Tom Jonas Janzen wrote:
> Dear Francesco,
>
> first of all, i want to thank you for your help with my problem. I already
> could fix the previous error with your given tipps on my input file.
>
> I used the 16O Oxygen variant because my topic of the Master Thesis is a
> follow up concept of another srtudents Master Thesis. My supervising
> Professor gave me his Input files so i could better understand the function
> of the LOW-MAT and MATERIAL cards. Therefore, i used his pattern as help.
> The student before me used this 16O Oxygen pattern for his simulation with
> PMMA so i thought it would also be good to use it too.
>
> The problem is, that currently an other error occured during the first test
> run after my adjustments ...
> Right now the error file says the following:
>
> *** (n,p) proton production activated for Xsec mat. # 1 ***
> **** Low energy neutron xsec not found for some media 2
> 5 ****
> HYDROGEN
> CARBON
> OXYGEN
>
> I can see, that the error written above must occure from the MATERIAL or
> LOW-MAT cards, but i don't know how to fix the error...
> Do you have any idea what is the problem ?
>
> I have written down the now changed INPUT file below.
> Thank you very much for your help and your time effort.
>
> Best regards,
> Tom Janzen
>
> Changed INPUT file:
> 04_Testsimulation_PMMA
> * Set the defaults for precision simulations
> DEFAULTS
> NEUTRONS
> * Define the beam characteristics
> BEAM -0.003 10000.0
> 1.0NEUTRON
> BEAMPOS 0.0 0.0 -5.0
> GEOBEGIN
> COMBNAME
> 3 0
> RCC blkbody 0.0 0.0 -100.0 0.0 0.0 200. 150.
> RCC void 0.0 0.0 -75.0 0.0 0.0 150.0 100.0
> RCC air 0.0 0.0 -10.0 0.0 0.0 20.0 15.0
> RCC PMMA 0.0 0.0 -5.0 0.0 0.0 1.0 10.0
> RCC PMMA2 0.0 0.0 -4.0 0.0 0.0 1.0 10.0
> RCC PMMA3 0.0 0.0 -3.0 0.0 0.0 1.0 10.0
> RCC PMMA4 0.0 0.0 -2.0 0.0 0.0 1.0 10.0
> RCC PMMA5 0.0 0.0 -1.0 0.0 0.0 1.0 10.0
> RCC PMMA6 0.0 0.0 0.0 0.0 0.0 1.0 10.0
> RCC PMMA7 0.0 0.0 1.0 0.0 0.0 1.0 10.0
> RCC PMMA8 0.0 0.0 2.0 0.0 0.0 1.0 10.0
> RCC PMMA9 0.0 0.0 3.0 0.0 0.0 1.0 10.0
> RCC PMMA10 0.0 0.0 4.0 0.0 0.0 1.0 10.0
> END
> * Black hole
> BLKBODY 5 +blkbody -void
> * Void around
> VOID 5 +void-air
> AIR 5 +air -PMMA -PMMA2 -PMMA3 -PMMA4 -PMMA5 -PMMA6 -PMMA7-PMMA8
> -PMMA9 -PMMA10
> PMMA 5 +PMMA
> PMMA2 5 +PMMA2
> PMMA3 5 +PMMA3
> PMMA4 5 +PMMA4
> PMMA5 5 +PMMA5
> PMMA6 5 +PMMA6
> PMMA7 5 +PMMA7
> PMMA8 5 +PMMA8
> PMMA9 5 +PMMA9
> PMMA10 5 +PMMA10
> END
> 9424777.964712388.9813823.0076314.159265314.159265314.159265314.159265
> 314.159265314.159265314.159265314.159265314.159265314.159265
> GEOEND
> * Polymethyl Methacrylate PMMA Photoresist
> * H
> * Chemical H H-C-H
> * Formula | |
> * Polymethyl ----- C - C -----
> * Methacrylate | | H
> * H C - O - C-H
> * C H O || H
> * 5 8 2 O
> *
> * atomic density calculated:
> * H-1
> * H-2
> * C
> * O-16
> MATERIAL 1.19 PMMA
> COMPOUND 0.5332533 HYDROGEN 0.3333333 CARBON 0.1333333 OXYGENPMMA
> MATERIAL 1. 8.99E-05
> 0.0HYDROGEN
> LOW-MAT HYDROGEN 1. 11. 296.
> HYDROGEN
> MATERIAL 6. 2.26 0.0CARBON
>
> LOW-MAT CARBON 6. -2. 293. CARBON
>
> MATERIAL 8. 0.001333 0.0OXYGEN
>
> LOW-MAT OXYGEN 8. 16. 293. OXYGEN
>
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7..
> ASSIGNMA BLCKHOLE BLKBODY
> ASSIGNMA VACUUM VOID
> ASSIGNMA AIR AIR
> ASSIGNMA PMMA PMMA
> ASSIGNMA PMMA PMMA2
> ASSIGNMA PMMA PMMA3
> ASSIGNMA PMMA PMMA4
> ASSIGNMA PMMA PMMA5
> ASSIGNMA PMMA PMMA6
> ASSIGNMA PMMA PMMA7
> ASSIGNMA PMMA PMMA8
> ASSIGNMA PMMA PMMA9
> ASSIGNMA PMMA PMMA10
> USRTRACK 1. NEUTRON -50. PMMA314.159265
> 10.detect1
> USRTRACK 1E-11 3E-12 &
> USRTRACK 1. NEUTRON -50. PMMA314.159265
> 10.detect2
> USRTRACK 3E-11 1E-11 &
> USRTRACK 1. NEUTRON -50. PMMA314.159265
> 10.detect3
> USRTRACK 1E-10 3E-11 &
> USRTRACK 1. NEUTRON -50. PMMA314.159265
> 10.detect4
> USRTRACK 3E-10 1E-10 &
> USRTRACK 1. NEUTRON -50. PMMA314.159265
> 10.detect5
> USRTRACK 1E-9 3E-10 &
> USRTRACK 1. NEUTRON -50. PMMA314.159265
> 10.detect6
> USRTRACK 3E-9 1E-9 &
> USRTRACK 1. NEUTRON -50. PMMA314.159265
> 10.detect7
> USRTRACK 1E-8 3E-9 &
> USRTRACK 1. NEUTRON -50. PMMA314.159265
> 10.detect8
> USRTRACK 3E-8 1E-8 &
> USRTRACK 1. NEUTRON -50. PMMA314.159265
> 10.detect9
> USRTRACK 1E-7 3E-8 &
> USRTRACK 1. NEUTRON -50. PMMA314.159265
> 10.detect10
> USRTRACK 3E-7 1E-7 &
> USRTRACK 1. NEUTRON -50. PMMA314.159265
> 10.detect11
> USRTRACK 1E-6 3E-7 &
> USRTRACK 1. NEUTRON -50. PMMA314.159265
> 10.detect12
> USRTRACK 3E-6 1E-6 &
> USRTRACK 1. NEUTRON -50. PMMA314.159265
> 10.detect13
> USRTRACK 1E-5 3E-6 &
> USRTRACK 1. NEUTRON -50. PMMA314.159265
> 10.detect14
> USRTRACK 3E-5 1E-5 &
> USRTRACK 1. NEUTRON -50. PMMA314.159265
> 10.detect15
> USRTRACK 1E-4 3E-5 &
> USRTRACK 1. NEUTRON -50. PMMA314.159265
> 10.detect16
> USRTRACK 3E-4 1E-4 &
> USRTRACK 1. NEUTRON -50. PMMA314.159265
> 10.detect17
> USRTRACK 1E-3 3E-4 &
> USRTRACK 1. NEUTRON -50. PMMA314.159265
> 10.detect18
> USRTRACK 3E-3 1E-3 &
> USRTRACK 1. NEUTRON -50. PMMA314.159265
> 10.detect19
> USRTRACK 1E-2 3E-3 &
> USRTRACK 1. NEUTRON -50. PMMA314.159265
> 10.detect20
> USRTRACK 1E-2 1E-11 &
> * Set the random number seed
> RANDOMIZ 1.0
> * Set the number of primary histories to be simulated in the run
> START 1000000.
> STOP
>
>
>
> -----Ursprüngliche Nachricht-----
> Von: owner-fluka-discuss_at_mi.infn.it [mailto:owner-fluka-discuss_at_mi.infn.it]
> Im Auftrag von Francesco Cerutti
> Gesendet: Montag, 15. Mai 2017 10:03
> An: Tom Jonas Janzen (tom.jonas.janzen_at_hs-furtwangen.de)
> <tom.jonas.janzen_at_hs-furtwangen.de>
> Cc: FLUKA discussion <fluka-discuss_at_fluka.org>
> Betreff: Re: [fluka-discuss]: scoring neutron fluence & thermal neutrons
>
>
>> Then, have you any reason to select 16O only instead of natural oxygen
>> composition? It's what you are doing with the 16. in the respective
>> MATERIAL card.
> well, since yours is a pure low energy neutron problem and the respective
> oxygen cross section set refers indeed to 16O, your monoisotopic material
> definition embeds in fact some coherence. In practice, the isotopic
> composition input choice is irrelevant, since low energy neutron transport
> will be perfomed with respect to 16O only.
>
> Regards
> Francesco
>
> **************************************************
> Francesco Cerutti
> CERN-EN/STI
> CH-1211 Geneva 23
> Switzerland
> tel. ++41 22 7678962
> fax ++41 22 7668854
>
> On Sun, 14 May 2017, Francesco Cerutti wrote:
>
>> Dear Tom,
>>
>> the problem is due to the presence of the word "compound" in your
>> LOW-MAT card. It's a character string totally unexpected and meaningless.
>> By the way, why are you asking for the obsolete 72-group neutron
>> library, which had a single thermal group and acutally is no longer
>> distributed? Just remove your LOW-NEUT card (this way you will
>> automatically use the current 260-group library).
>> Remove also the material index (WHAT(4)) from all MATERIAL cards, it's
>> useless and potentially dangerous.
>> Remove the AUXSCORE card that is pointless (in USRTRACK you are
>> scoring neutron fluence, so further filtering on NEUTRON makes no sense).
>> Then, have you any reason to select 16O only instead of natural oxygen
>> composition? It's what you are doing with the 16. in the respective
>> MATERIAL card.
>>
>> Best wishes
>>
>> Francesco
>>
>> **************************************************
>> Francesco Cerutti
>> CERN-EN/STI
>> CH-1211 Geneva 23
>> Switzerland
>> tel. ++41 22 7678962
>> fax ++41 22 7668854
>>
>> On Sun, 14 May 2017, Tom Jonas Janzen
>> (tom.jonas.janzen_at_hs-furtwangen.de)
>> wrote:
>>
>>> Hello the Fluka community!
>>>
>>> My name is Tom Janzen and i'm currently working with FLUKA for my
>>> Master Thesis at the Hochschule Furtwangen.
>>>
>>> My work experience with Fluka is scaled on an intermediate leval but
>>> i have a problem where i need some help of the community.
>>> For my master Thesis im need to score the Neutron Fluence of a
>>> isotropic beam in ten bodys of PMMA with the USRTRCK Card. Also i
>>> need to track the neutron fluence in the energy spectrum from 3E-12 GeV
> to 1E-2 GeV.
>>> The aim of this analysis should be results of the USRTRCK, which i
>>> can transform in a diagram, where i can visualize the Neutron
>>> Fluence in every PMMA Body at the Energy Range of 3E-12 GeV to 1E-2
>>> GeV. This means, that i am trying to track thermal neutrons and
>>> visualize the moderation of neutrons in the PMMA (moderator
>>> material).
>>>
>>> Right know i have a Input file which is looking fine to me. I am
>>> working with the Cards which should track the low energy neutrons
>>> and the thermal Neutrons as well. Also i took care of the Material
>>> cards and implemented in my Input file the LOW-MAT card in the way
>>> the manual refers to. But when i am trying to run the simulation, i
>>> get an Error with the following text:
>>> *** The 5th field compound of the following input card ***
>>> LOW-MAT HYDROGEN 1. 31. 293. compound
>>> HYDROG-1
>>> *** does not contain a valid formatted fortran real number ***
>>> *** when one is expected ***
>>> I tryed to fix it by my own beut i have no clue...
>>>
>>> I Hope the FLUKA users and experts can help me with this little problem.
>>> Thank you very much for your help and your time effort. The used
>>> Input file is written down below.
>>>
>>> Best regards,
>>> Tom Janzen
>>>
>>> My Input File:
>>> 04_Testsimulation_PMMA
>>> * Set the defaults for precision simulations DEFAULTS NEUTRONS
>>> * Define the beam characteristics
>>> BEAM -0.003 10000.0 1.0NEUTRON BEAMPOS
>>> 0.0 0.0 -5.0 GEOBEGIN COMBNAME
>>> 3 0
>>> RCC blkbody 0.0 0.0 -100.0 0.0 0.0 200. 150.
>>> RCC void 0.0 0.0 -75.0 0.0 0.0 150.0 100.0 RCC air 0.0
>>> 0.0 -10.0 0.0 0.0 20.0 15.0 RCC PMMA 0.0 0.0 -5.0 0.0 0.0 1.0
>>> 10.0 RCC PMMA2 0.0 0.0 -4.0 0.0 0.0 1.0 10.0 RCC PMMA3
>>> 0.0 0.0 -3.0 0.0 0.0 1.0 10.0 RCC PMMA4 0.0 0.0 -2.0 0.0 0.0
>>> 1.0 10.0 RCC PMMA5 0.0 0.0 -1.0 0.0 0.0 1.0 10.0 RCC PMMA6
>>> 0.0 0.0 0.0 0.0 0.0 1.0 10.0 RCC PMMA7 0.0 0.0 1.0 0.0 0.0 1.0
>>> 10.0 RCC PMMA8 0.0 0.0 2.0 0.0 0.0 1.0 10.0 RCC PMMA9 0.0
>>> 0.0 3.0 0.0 0.0 1.0 10.0 RCC PMMA10 0.0 0.0 4.0 0.0 0.0 1.0 10.0
>>> END
>>> * Black hole
>>> BLKBODY 5 +blkbody -void
>>> * Void around
>>> VOID 5 +void-air
>>> AIR 5 +air -PMMA -PMMA2 -PMMA3 -PMMA4 -PMMA5 -PMMA6
>>> -PMMA7-PMMA8
>>> -PMMA9 -PMMA10
>>> PMMA 5 +PMMA
>>> PMMA2 5 +PMMA2
>>> PMMA3 5 +PMMA3
>>> PMMA4 5 +PMMA4
>>> PMMA5 5 +PMMA5
>>> PMMA6 5 +PMMA6
>>> PMMA7 5 +PMMA7
>>> PMMA8 5 +PMMA8
>>> PMMA9 5 +PMMA9
>>> PMMA10 5 +PMMA10
>>> END
>>>
>>> 9424777.964712388.9813823.0076314.159265314.159265314.159265314.15926
>>> 5
>>> 314.159265314.159265314.159265314.159265314.159265314.159265
>>> GEOEND
>>> * Polymethyl Methacrylate PMMA Photoresist
>>> * H
>>> * Chemical H H-C-H
>>> * Formula | |
>>> * Polymethyl ----- C - C -----
>>> * Methacrylate | | H
>>> * H C - O - C-H
>>> * C H O || H
>>> * 5 8 2 O
>>> *
>>> * atomic density calculated:
>>> * H-1
>>> * H-2
>>> * C
>>> * O-16
>>> MATERIAL 1.19
>>> PMMA COMPOUND 0.5332533 HYDROGEN 0.3333333 CARBON 0.1333333
>>> OXYGENPMMA MATERIAL 1. 8.99E-05 3.
>>> 0.0HYDROGEN
>>> LOW-MAT HYDROGEN 1. 31. 293. compound
>>> HYDROG-1
>>> MATERIAL 6. 2.26 6. 0.0CARBON
>>> LOW-MAT CARBON 6. -2. 293. compound CARBON
>>> MATERIAL 8. 0.001333 8.
>>> 16.OXYGEN
>>> LOW-MAT OXYGEN 8. 16. 293. OXYGEN
> LOW-NEUT
>>> 72. 22. 0.0196 1. 1.
>>> *LOW-BIAS 20. 73. 0.95 PMMA PMMA 1.
>>> *SCORE NEUTRON THNEU-EQ
>>> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7..
>>> ASSIGNMA BLCKHOLE BLKBODY
>>> ASSIGNMA VACUUM VOID
>>> ASSIGNMA AIR AIR
>>> ASSIGNMA PMMA PMMA
>>> ASSIGNMA PMMA PMMA2
>>> ASSIGNMA PMMA PMMA3
>>> ASSIGNMA PMMA PMMA4
>>> ASSIGNMA PMMA PMMA5
>>> ASSIGNMA PMMA PMMA6
>>> ASSIGNMA PMMA PMMA7
>>> ASSIGNMA PMMA PMMA8
>>> ASSIGNMA PMMA PMMA9
>>> ASSIGNMA PMMA PMMA10
>>> AUXSCORE USRTRACK NEUTRON detect1 detect20
>>> AMB74
>>> USRTRACK 1. NEUTRON -50. PMMA314.159265
>>> 10.detect1
>>> USRTRACK 1E-11 3E-12
>>> & USRTRACK 1. NEUTRON -50. PMMA314.159265
>>> 10.detect2
>>> USRTRACK 3E-11 1E-11
>>> & USRTRACK 1. NEUTRON -50. PMMA314.159265
>>> 10.detect3
>>> USRTRACK 1E-10 3E-11
>>> & USRTRACK 1. NEUTRON -50. PMMA314.159265
>>> 10.detect4
>>> USRTRACK 3E-10 1E-10
>>> & USRTRACK 1. NEUTRON -50. PMMA314.159265
>>> 10.detect5
>>> USRTRACK 1E-9 3E-10
>>> & USRTRACK 1. NEUTRON -50. PMMA314.159265
>>> 10.detect6
>>> USRTRACK 3E-9 1E-9
>>> & USRTRACK 1. NEUTRON -50. PMMA314.159265
>>> 10.detect7
>>> USRTRACK 1E-8 3E-9
>>> & USRTRACK 1. NEUTRON -50. PMMA314.159265
>>> 10.detect8
>>> USRTRACK 3E-8 1E-8
>>> & USRTRACK 1. NEUTRON -50. PMMA314.159265
>>> 10.detect9
>>> USRTRACK 1E-7 3E-8
>>> & USRTRACK 1. NEUTRON -50. PMMA314.159265
>>> 10.detect10
>>> USRTRACK 3E-7 1E-7
>>> & USRTRACK 1. NEUTRON -50. PMMA314.159265
>>> 10.detect11
>>> USRTRACK 1E-6 3E-7
>>> & USRTRACK 1. NEUTRON -50. PMMA314.159265
>>> 10.detect12
>>> USRTRACK 3E-6 1E-6
>>> & USRTRACK 1. NEUTRON -50. PMMA314.159265
>>> 10.detect13
>>> USRTRACK 1E-5 3E-6
>>> & USRTRACK 1. NEUTRON -50. PMMA314.159265
>>> 10.detect14
>>> USRTRACK 3E-5 1E-5
>>> & USRTRACK 1. NEUTRON -50. PMMA314.159265
>>> 10.detect15
>>> USRTRACK 1E-4 3E-5
>>> & USRTRACK 1. NEUTRON -50. PMMA314.159265
>>> 10.detect16
>>> USRTRACK 3E-4 1E-4
>>> & USRTRACK 1. NEUTRON -50. PMMA314.159265
>>> 10.detect17
>>> USRTRACK 1E-3 3E-4
>>> & USRTRACK 1. NEUTRON -50. PMMA314.159265
>>> 10.detect18
>>> USRTRACK 3E-3 1E-3
>>> & USRTRACK 1. NEUTRON -50. PMMA314.159265
>>> 10.detect19
>>> USRTRACK 1E-2 3E-3
>>> & USRTRACK 1. NEUTRON -50. PMMA314.159265
>>> 10.detect20
>>> USRTRACK 1E-2 1E-11
>>> &
>>> * Neutron fluence scoring
>>> AUXSCORE USRTRACK NEUTRON detect1 detect20
>>> AMB74
>>> USRBIN 10. NEUTRON -70. 100. 100.
>>> 100.Neu-Flu USRBIN -100. -100. -100. 200.
>>> 200. 240. &
>>> * Set the random number seed
>>> RANDOMIZ 1.0
>>> * Set the number of primary histories to be simulated in the run
>>> START 1000000.
>>> STOP
>>>
>>>
> __________________________________________________________________________
> You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id¬c_info

Dear Tom,

Errors are in LOW-MAT parameters for CARBON and OXYGEN materials.
WHAT(4) must be 296 instead of 293.


Best regards,

Mikhail

__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Tue May 16 2017 - 19:49:54 CEST

This archive was generated by hypermail 2.3.0 : Tue May 16 2017 - 19:50:09 CEST