(no subject)

From: Sebastien WURTH (wurth@ipno.in2p3.fr)
Date: Wed Mar 05 2008 - 09:21:30 CET

  • Next message: Sebastien WURTH: "(no subject)"

    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=
    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=
    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.

    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=
    > 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=
    > >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=
    > >recommended to calculate dose equivalent?
    > >
    > >Best Regards
    > >Hamide
    > >
    > >
    > >

  • Next message: Sebastien WURTH: "(no subject)"

    This archive was generated by hypermail 2.1.6 : Wed Mar 05 2008 - 12:08:28 CET