[fluka-discuss]: Re: How to calculate activation cross sections with FLUKA

From: Francesco Cerutti <Francesco.Cerutti_at_cern.ch>
Date: Sat, 2 Nov 2013 12:11:02 +0100

Dear Oyeon,

please send your questions to the dedicated mailing list here in cc
(which you have to have registered) and not privately.

Your guess is not correct: to calculate specific reaction channels, FLUKA
does not use precomputed cross sections but its own models (PEANUT in case
of single hadron - like neutron - induced reactions). However, low energy
(<20MeV) neutrons represent an exception: in this case residue production
- and neutron transport - is calculated according to library data, that
you can access through the printing flag of the LOW-NEUT card.
Note that the use of FLUKA is not intended to extract microscopic cross
sections, rather to calculate - with respect to your example -
radionuclide yields/activity for a specific geometry.

Kind regards

FrancEsco

**************************************************
Francesco Cerutti
CERN-EN/STI
CH-1211 Geneva 23
Switzerland
tel. ++41 22 7678962
fax ++41 22 7668854

On Fri, 1 Nov 2013, Oyeon Kum wrote:

>
> Dear Francisco,
>
> I am going to bother you again. 
> I would appreciate it very much if you could let me know how to calculate
> activation cross sections with FLUKA, for example, 27Al(n,p)27Mg or
> 27Al(n,alpha)24Na, etc.
>
> I guess that FLUKA may use these activation cross sections inside.
> Then, we can catch it out in one way or the other. I looked for the FLUKA
> manual but I could not find.
>
> Thank you very much again in advance.
>
> Best regards,
> Oyeon
>
>
> 2013/10/29 Francesco Cerutti <Francesco.Cerutti_at_cern.ch>
>
> Dear Oyeon,
>
> if you want to get the isotope production per beam particle, use
> the RESNUCLEi card *with no association to any cooling time*
> (i.e. do no put your RESNUCLEi detector name in any DCYSCORE
> card). Concerning the normalization to the region volume (just
> used as a scaling factor), if you do not specify the latter in
> the RESNUCLEi card the resulting values actually will not be
> expressed per cmc (despite the unit label in the file, assuming
> in this case a default normalizing volume of 1 cmc) but will
> refer to the whole region (nuclei/p in that region).
> To process RESNUCLEi results, use usrsuw.
>
> (Bq is the unit for the activity of radioisotopes *at a given
> time*, which is the product of the radioisotope population at
> that time - depending on build-up [from direct production and
> parent decay] and decay - by its decay constant).
>
> Best regards
>
> Francesco
>
> **************************************************
> Francesco Cerutti
> CERN-EN/STI
> CH-1211 Geneva 23
> Switzerland
> tel. ++41 22 7678962
> fax  ++41 22 7668854
>
> On Tue, 29 Oct 2013, Oyeon Kum wrote:
>
> Dear Dr. Cerutti, 
> My name is Oyeon Kum, Ph.D. and is calculating
> radiation protection problems
> with FLUKA 2011 code at the Korean Institute of
> Radiological and Medical
> Sciences in Korea.
>
> I want to get Isotopic abundance (nuclei/p unit) of
> the produced nuclei as
> shown in L. Sarchiapone 's LNL Annual Report
> (attached here) about SPES
> project.
> I ran FLUKA code and supportive programs such as
> usrsuw and usrsuwev and
> tried to get the Isotopic abundance of the produced
> nuclei shown in Fig.2
> and Table 1 in his report, but I could not find it.
> Rather, I found the same
> values with different unit, nuclei/cmc/pr and
> Bq/cmc, from usrsuw and
> usrsuwev, respectively. The two results are also
> attached here. The file
> name, "output1_sum.lis" is from usrsuwev and
> "irr1day_61_sum.lis" from
> usrsuw.
>
> I would very much appreciate it if you could let me
> know how to get the
> isotopic abundance data from FLUKA code simulation
> and the meaning of the
> two different units with the same values, if
> possible. Thank you very much
> again in advance for your kindness.
>
> Sincerely,
> Oyeon Kum
> e-mail: okum1211_at_gmail.com; okum_at_kirams.re.kr
>
>
>
>
Received on Sat Nov 02 2013 - 13:07:17 CET