Re: [fluka-discuss]: thermal neutron induced radioactivity in natural tungsten

From: Stuart Ansell <stuart.ansell_at_maxiv.lu.se>
Date: Mon, 25 Jun 2018 14:31:19 +0000

Dear Saurabuh,


   Not sure that it will fix anything but you have a typo error in your source1.f file (I think). You write:


        ESAMPLE=XYZ*(ENEPRO(I)-ENEPRO(I-1))*(ENEPRO(I)+ENEPRO(I-1))
        ESAMPLE=SQRT(ESAMPLE+ENEPRO(I-1)**2)-ENEPRO(I-1)
        ESAMPLE=ESAMPLE/(ENEPRO(I)-ENEPRO(I-1))


attempting to get a linear energy based on XYZ being between 0.0 and 1.0.


and I believe that you intended that middle line should be

ESAMPLE=SQRT(ESAMPLE**2+ENEPRO(I-1)**2)-ENEPRO(I-1)


It leaves rather large chunks of the energy spectrum unsampled as it currently stands.



    Regards,


         Stuart Ansell






________________________________
From: owner-fluka-discuss_at_mi.infn.it <owner-fluka-discuss_at_mi.infn.it> on behalf of Saurabh Mukherjee <sxm3816_at_gmail.com>
Sent: Sunday, June 24, 2018 7:43:29 PM
To: fluka-discuss
Subject: [fluka-discuss]: thermal neutron induced radioactivity in natural tungsten

Dear experts,
I am trying to simulate residual activity due to thermal neutron incident on natural tungsten. The input file is attached, along with the source and spectrum file. The isotopes generated by FLUKA are - W-187, W-185 and Re-187.
 But when I compare it to online calculator (https://www.ncnr.nist.gov/resources/activation/) (the inputs to the online calculator are shown in attachment 2 and the results in attachment3), the results are different. I see that in online calculator, the W-187, W-185 isotopes are present but there are some extra isotopes present as well-W-181, W-188 and Ta-183. These are not present in FLUKA results. Is there some problem with my input file? Can some one please take a look?
sincerely
saurabh


__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Mon Jun 25 2018 - 17:57:43 CEST

This archive was generated by hypermail 2.3.0 : Mon Jun 25 2018 - 17:57:45 CEST