From: Sebastien WURTH (wurth@ipno.in2p3.fr)
Date: Wed Mar 05 2008 - 09:21:30 CET
040910
X-Accept-Language: fr, en
MIME-Version: 1.0
To: hamide jalali <jalali.hb@gmail.com>
CC: fluka-discuss@fluka.org
Subject: Re: deq99c.f & dose equivalent
References: <200803040821.m248L7As029152@smtp2.mi.infn.it>=09 <47CD0CCB.409=
0509@ipno.in2p3.fr> <64c3dcd00803041240v4bc7ffbcmf89f4bec5898ce3e@mail.gmai=
l.com>
In-Reply-To: <64c3dcd00803041240v4bc7ffbcmf89f4bec5898ce3e@mail.gmail.com>
Content-Type: text/plain; charset=3DISO-8859-1; format=3Dflowed
Content-Transfer-Encoding: 8bit
Sender: owner-fluka-discuss@smtp2.mi.infn.it
Hello Hamide,
I reply to the all discussion group, I hope you don't mind, some people=20
could just have better answers than mine.
I'm not sure I really understood your question. In shielding problems=20
you want to calculate dose equivalent for alpha particles ?
The dose should be huge yes, but on very short distances !.
Then, if you want to calculate dose from secondaries induced by=20
reactions of these alpha, you have the case we discussed yesterday.
For heavy ions, I am no expert.
But, in general, you can always use FLUSCW routine without deq99c, and=20
convert your fluence into dose with your own conversion factors,
I did it before I started to use Stefan Roesler's routine.
For instance for photons (monoenergetic example here), you can find this=20
relation : D' =3D F'*E*=B5en/rh=F4, where :
D' is the absorbed dose (=3D equivalent dose for photons) rate in Gy per=20
second, E the energy in J and =B5en/rh=F4 the absorption coefficient in m=
=B2/kg.
F' is the fluence rate in particles per m=B2 per second : F' =3D A* Nu /=20
4*pi*d=B2 where A is the activity in Bq, Nu the branching ratio of the=20
energy considered, d the distance considered in meters. This is given by=20
FLUKA, in fact, FLUKA gives F'/A in particules/cm=B2/primary.
Of course, some work on the units must be done to use MeV or GeV for=20
energy, cm=B2/g for absorption coefficient, etc.
You calculate your factor to be applied to fluence during run time, you=20
modify FLUSCW routine as explained in the manual (user routines=20
section), you compile, link with your executable and run it. You could=20
use such relations for alpha or heavy ions...
Hope it helps.
Best regards.
Sebastien.
hamide jalali wrote :
> Dear Sebastien
> =20
> I am very thankful to you for your reply. Now I can be sure that my=20
> calculation is correct.
> I use deq99c.f for the same purpose as you, I mean shielding problem.=20
> In a part of my problem I should calculate dose equivalent when the=20
> source is alpha particles and also heavy ions. In this case FLUKA=20
> don't support my calculations because it hasn't any routine for this=20
> kind of source so my calculation will remain half-done.
> Have you ever faced this problem? What is your suggestion when source=20
> is heavy nuclei?
> =20
> Any help would be appreciated in advance.
>
> Kind Regards
> Hamide
> =20
> =20
> On Tue, Mar 4, 2008 at 12:48 AM, Sebastien WURTH <wurth@ipno.in2p3.fr=20
> <mailto:wurth@ipno.in2p3.fr>> wrote:
>
> Hello Hamide,
>
> The code to include in your detector's SDUM is explained at the
> beginning of the deq routine (maybe AMB74 here).
> Using particle code 201 to score the dose equivalent for all particle=
s
> included in the routine is the right thing to do, I used it very
> often.
> You can calculate separately and add up, if everything goes OK, it
> should be the same.
> Generalised particle is not an option in my opinion. We are talking
> about conversion from fluence to dose equivalent here (using FLUSCW
> routine), not from energy deposited to dose (for this one, COMSCW
> routine is the suitable one).
> You can use the three kinds of detectors to calculate dose, I usually
> use USRBIN to plot spatial distributions for shielding purposes, it
> depends of what you need.
>
> Hope it helps.
> Regards.
> Sebastien.
>
>
> hamide jalali wrote :
>
> >Dear FLUKA users
> >
> >I'd like to get ambient dose equivalent in ICRU sphere by using
> >deq99c.f while the source is Proton with a spectrum of energy
> which is
> >defined by source.f.
> >I want to know how can get total dose equivalent from all secondarie=
s
> >which are included deq99c.f ? Is it enough to set the particle type
> >to be scored=3D 201 (all particle) or I should calculate dose
> equivalent
> >for each secondary separately, then add them up? How about
> generalised
> >particle 208 (energy)?
> >And the other question, among USRBIN,USRBDX and USRTRACK which one i=
s
> >recommended to calculate dose equivalent?
> >
> >Best Regards
> >Hamide
> >
> >
> >
>
>
This archive was generated by hypermail 2.1.6 : Wed Mar 05 2008 - 12:08:28 CET