Dose and dose rate estimation, some general questions

From: Andrey Smolyakov <Andrei.Smolyakov_at_itep.ru>
Date: Wed, 8 Sep 2010 13:13:15 +0400

Dear Fluka experts, I ask for your help and some clarification on how to
estimate dose and dose rate in Fluka.

The scope of my problem:

I need to calculate peak radiation level after one proton beam shot. In
order to do this I define a radiation profile by IRRPROFI, request
radiative decay simulation by RADDECAY, define decay time for scoring by
DCYTIMES and associate my estimator with this decay time by DCYSCORE.
This part of the code looks this way:

IRRPROFI 1.0E-6 1.0E+18
RADDECAY 1.0
DCYTIMES 0.0E+0
DCYSCORE 1.0 1.0
USRBIN

Here I define a beam with duration of 10mks and total number of
particles 1E12. And I request scoring right after the beam has passed.

And now I have some questions as to what data I will obtain in this way:

1. In my regular simulation WITHOUT the four lines mentioned above my
USRBIN estimator will produce data on dose deposited in it's binning by
all primary and secondary particles (for which the transport has been
requested, of course). I suppose it will not include dose deposited by
decay products. Is this true?

2. Now WITH these four lines, will the estimator produce data for dose
rate deposited by primary/secondary particles AND decay products (A)? Or
will it be dose rate due to ONLY decay products (B)?

3. Since the dose rate is defined as a dose over some period of time,
what (how long) period of time is taken in FLUKA to estimate dose rate?
Or if is not strictly defined, what are the guidelines to asses it?

4. In case of 1 and 2B are true, how to calculate total dose and dose
rate?

5. Extra question. If I want to score dose not in GeV but in Sv then in
case of difficult geometry it is better to use USERWEIG card and
corresponding user routine. According to the manual, I need to prepare
FLUSCW routine, compile and link to the Fluka library. But how can I
link multiple user routines, say, if apart from FLUSCW I need to define
my own MAGFLD routine? What is the correct syntax for command to do it?

Thanks in advance,

Andrey Smolyakov
Engineer
Institute for Theoretical and Experimental Physics
B.Cheremushkinskaya 25, Moscow, Russia
 
Received on Wed Sep 08 2010 - 15:14:44 CEST

This archive was generated by hypermail 2.2.0 : Wed Sep 08 2010 - 15:14:50 CEST