Re: Can't Calculate Residual Dose Rate, It always ZERO, Need Help!

From: Vahan Petrosyan <vpetrosyan_at_asls.candle.am>
Date: Thu, 29 Nov 2012 12:29:56 +0400

<0BC70B5D93E054469872FFD0FE07220EA7B2B176_at_CERNXCHG12.cern.ch> <50B5A857.4050900_at_gmail.com>
In-Reply-To: <50B5A857.4050900_at_gmail.com>
Content-Type: multipart/alternative;
  boundary="------------020701040801060402030403"
X-Antivirus: avast! (VPS 121128-1, 29.11.2012), Outbound message
X-Antivirus-Status: Clean
X-Greylist: Sender IP whitelisted, not delayed by milter-greylist-4.2.7 (smtp1.mi.infn.it [192.84.138.69]); Thu, 29 Nov 2012
09:30:16 +0100 (CET)
X-Virus-Scanned: clamav-milter 0.97.6 at smtp1.mi.infn.it
X-Virus-Status: Clean
Sender: owner-fluka-discuss_at_mi.infn.it

This is a multi-part message in MIME format.
--------------020701040801060402030403
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
Content-Transfer-Encoding: 7bit

many thanks!

1) I want to know If I calculate residual dose rate is USERWEIGHT card
required or it is optional card?\

For calculating Residual Dose Rate by fluka there is Stefan Roesler's
two-step method which gives more accurate results.
2) did it already implement in fluka or it should be implemented a part.?

>
> On 27.11.2012 12:00, Vasilis Vlachoudis wrote:
>> Hi Vahan,
>>
>> The results of eq-dose on the prompt part of the calculation are in
>> pSv/primary and on the decay part pSv/sec
>>
>> USERWEIGht allows the user to modify the scoring quantity before the
>> scoring takes place by enabling the call to fluscw and to comscw user
>> routines. e.g. apply the efficiency of your detector, make USRBIN
>> fluence plots with cuts in energy etc..
>>
>> concerning PHOTONUC your assumption is correct 20MeV is too low,
>> however I would leave 1.0.
>>
>> Regards
>> Vasilis
>>
>> ------------------------------------------------------------------------
>> *From:* owner-fluka-discuss_at_mi.infn.it
>> [owner-fluka-discuss_at_mi.infn.it] on behalf of Vahan Petrosyan
>> [vpetrosyan_at_asls.candle.am]
>> *Sent:* 20 November 2012 06:40
>> *To:* fluka-discuss_at_fluka.org
>> *Subject:* Can't Calculate Residual Dose Rate, It always ZERO, Need Help!
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>> many thanks! I was able to use it properly using your advice.
>>
>> one more question about residual dose.
>>
>> If I have for example 100pC charge per second and set it in IRRPROFI
>> (6.5*10^8) are this results for that charge already in units pSv/h?
>> I mean what is a unit of the residual (USERBIN) Eq. DOSE given by fluka?
>> is it pSv/hour ?
>>
>> what role does play card USERWEIG? when should I use it and with what
>> options?
>>
>> I want to ask question about PHOTONUC card option WHAT(1).
>> If I set What(1) 1, or 2 results are the same and I think that it is
>> because of my electron beam initial energy is about 0.02Gev ie. there is
>> no so high energy particles in order to activate it.
>>
>> My question is If I set WHAT(1) to 1/On, and I have low energy beam will
>> results be wrong? I mean activating photonuclear interactions ONLY (as
>> it said in fluka input options E > 0.7 Gev ) at higher energy ranges
>> calculations will miss something for 0.02Gev electron beam which is
>> mush less than 0.7 Gev.
>>
>>
>>
>> On 18.11.2012 21:40, Stefan Roesler wrote:
>> > Dear Vahan
>> >
>> > With an electron beam you should not kill the electromagnetic cascade for
>> > the prompt radiation... Thus, remove RADDECAY/What(5)=99999.
>> >
>> > Otherwise, note that your steel SS316LN has index 30 which is beyond
>> > lead.
>> > Thus, if you are interested in photonuclear interactions in steel please
>> > use instead of LEAD the generic variable @LASTMAT
>> >
>> > PHOTONUC 1010. VACUUM @LASTMAT
>> >
>> > Are you sure about the setting of What(1), i.e., excluding high-energy
>> > and
>> > quasi deuteron processes? I suggest to use PHOTONUC/What(1)=1
>> >
>> > Best regards
>> > Stefan
>> >
>> >
>> >
>> > On Fri, 16 Nov 2012, Vahan Petrosyan wrote:
>> >
>> >> Hi I tried to calculate residual dose rate according to this tutorial
>> >> but have problems.
>> >>http://www.fluka.org/content/course/NEA/lectures/Activation.pdf
>> >>
>> >> I switched on all necessary cards (RADDECCAY, IRRPROFI, DCYTIMES,
>> >> DCYSCORE ) and linked USRBIN, RESNUCLEI scoring cards with appropriate
>> >> dacay times, and FLUKA runs without errors!
>> >>
>> >> The problem that I get all ZEROS in reporting plots of USERBIN,
>> >> RESNUCLEI (in my case type Eq. Dose), and even more for any type of
>> >> scored quantity not only Eq. Dose
>> >>
>> >> I also attached my input file if You have time please take a look.
>> >>
>> >> Thanks!
>> >>
>> >>
>> >
>> >
>>
>>
>>
>>
>

--------------020701040801060402030403
Content-Type: text/html; charset=ISO-8859-1
Content-Transfer-Encoding: 7bit

<html>
   <head>
     <meta content="text/html; charset=ISO-8859-1"
       http-equiv="Content-Type">
   </head>
   <body bgcolor="#FFFFFF" text="#000000">
     many thanks!<br>
     <br>
     1) I want to know If I calculate residual dose rate is USERWEIGHT
     card required or it is optional card?\<br>
     <br>
     For calculating Residual Dose Rate by fluka there is Stefan
     Roesler's two-step method which gives more accurate results.<br>
     2) did it already implement in fluka or&nbsp; it should be implemented a
     part.? <br>
     <pre>
</pre>
     <br>
     <blockquote cite="mid:50B5A857.4050900_at_gmail.com" type="cite">
       <div class="moz-cite-prefix"> <br>
         On 27.11.2012 12:00, Vasilis Vlachoudis wrote:<br>
       </div>
       <blockquote
         cite="mid:0BC70B5D93E054469872FFD0FE07220EA7B2B176_at_CERNXCHG12.cern.ch"
         type="cite">
         <meta http-equiv="Content-Type" content="text/html;
           charset=ISO-8859-1">
         <style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>
         <div style="direction: ltr;font-family: Tahoma;color:
           #000000;font-size: 10pt;">Hi Vahan,<br>
           <br>
           The results of eq-dose on the prompt part of the calculation
           are in pSv/primary and on the decay part pSv/sec<br>
           <br>
           USERWEIGht allows the user to modify the scoring quantity
           before the scoring takes place by enabling the call to fluscw
           and to comscw user routines. e.g. apply the efficiency of your
           detector, make USRBIN fluence plots with cuts in energy etc..<br>
           <div><br>
             concerning PHOTONUC your assumption is correct 20MeV is too
             low, however I would leave 1.0.<br>
             <br>
             Regards<br>
             <div style="font-family:Tahoma; font-size:13px"><font
                 size="1"><span style="font-family:Helvetica"></span><span
                   style="font-family:Helvetica">Vasilis </span><span
                   style="font-family:Helvetica"><br>
                   <br>
                 </span></font></div>
           </div>
           <div style="font-family: Times New Roman; color: #000000;
             font-size: 16px">
             <hr tabindex="-1">
             <div style="direction: ltr;" id="divRpF715060"><font
                 color="#000000" face="Tahoma" size="2"><b>From:</b> <a
                   moz-do-not-send="true"
                   class="moz-txt-link-abbreviated"
                   href="mailto:owner-fluka-discuss_at_mi.infn.it">owner-fluka-discuss_at_mi.infn.it</a>
                 [<a moz-do-not-send="true"
                   class="moz-txt-link-abbreviated"
                   href="mailto:owner-fluka-discuss_at_mi.infn.it">owner-fluka-discuss_at_mi.infn.it</a>]
                 on behalf of Vahan Petrosyan [<a moz-do-not-send="true"
                   class="moz-txt-link-abbreviated"
                   href="mailto:vpetrosyan_at_asls.candle.am">vpetrosyan_at_asls.candle.am</a>]<br>
                 <b>Sent:</b> 20 November 2012 06:40<br>
                 <b>To:</b> <a moz-do-not-send="true"
                   class="moz-txt-link-abbreviated"
                   href="mailto:fluka-discuss_at_fluka.org">fluka-discuss_at_fluka.org</a><br>
                 <b>Subject:</b> Can't Calculate Residual Dose Rate, It
                 always ZERO, Need Help!<br>
               </font><br>
             </div>
             <div>
               <div class="moz-forward-container">
                 <table class="moz-email-headers-table" border="0"
                   cellpadding="0" cellspacing="0" height="61" width="98">
                   <tbody>
                     <tr>
                       <th align="RIGHT" nowrap="nowrap"
                         valign="BASELINE"><br>
                       </th>
                       <td><br>
                       </td>
                     </tr>
                     <tr>
                       <th align="RIGHT" nowrap="nowrap"
                         valign="BASELINE"><br>
                       </th>
                       <td><br>
                       </td>
                     </tr>
                     <tr>
                       <th align="RIGHT" nowrap="nowrap"
                         valign="BASELINE"><br>
                       </th>
                       <td><br>
                       </td>
                     </tr>
                     <tr>
                       <th align="RIGHT" nowrap="nowrap"
                         valign="BASELINE"><br>
                       </th>
                       <td><br>
                       </td>
                     </tr>
                   </tbody>
                 </table>
                 <br>
                 <br>
                 <pre>many thanks! I was able to use it properly using your advice.

one more question about residual dose.

If I have for example 100pC charge per second and set it in IRRPROFI
(6.5*10^8) are this results for that charge already in units pSv/h?
I mean what is a unit of the residual (USERBIN) Eq. DOSE given by fluka?
is it pSv/hour ?

     what role does play card USERWEIG? when should I use it and with what
options?

I want to ask question about PHOTONUC card option WHAT(1).
If I set What(1) 1, or 2 results are the same and I think that it is
because of my electron beam initial energy is about 0.02Gev ie. there is
no so high energy particles in order to activate it.

My question is If I set WHAT(1) to 1/On, and I have low energy beam will
results be wrong? I mean activating photonuclear interactions ONLY (as
it said in fluka input options E &gt; 0.7 Gev ) at higher energy ranges
calculations will miss something for 0.02Gev electron beam which is
mush less than 0.7 Gev.

On 18.11.2012 21:40, Stefan Roesler wrote:
&gt; Dear Vahan
&gt;
&gt; With an electron beam you should not kill the electromagnetic cascade for
&gt; the prompt radiation... Thus, remove RADDECAY/What(5)=99999.
&gt;
&gt; Otherwise, note that your steel SS316LN has index 30 which is beyond
&gt; lead.
&gt; Thus, if you are interested in photonuclear interactions in steel please
&gt; use instead of LEAD the generic variable @LASTMAT
&gt;
&gt; PHOTONUC 1010. VACUUM @LASTMAT
&gt;
&gt; Are you sure about the setting of What(1), i.e., excluding high-energy
&gt; and
&gt; quasi deuteron processes? I suggest to use PHOTONUC/What(1)=1
&gt;
&gt; Best regards
&gt; Stefan
&gt;
&gt;
&gt;
&gt; On Fri, 16 Nov 2012, Vahan Petrosyan wrote:
&gt;
&gt;&gt; Hi I tried to calculate residual dose rate according to this tutorial
&gt;&gt; but have problems.
&gt;&gt; <a moz-do-not-send="true" class="moz-txt-link-freetext"
href="http://www.fluka.org/content/course/NEA/lectures/Activation.pdf"
target="_blank">http://www.fluka.org/content/course/NEA/lectures/Activation.pdf</a>
&gt;&gt;
&gt;&gt; I switched on all necessary cards (RADDECCAY, IRRPROFI, DCYTIMES,
&gt;&gt; DCYSCORE ) and linked USRBIN, RESNUCLEI scoring cards with appropriate
&gt;&gt; dacay times, and FLUKA runs without errors!
&gt;&gt;
&gt;&gt; The problem that I get all ZEROS in reporting plots of USERBIN,
&gt;&gt; RESNUCLEI (in my case type Eq. Dose), and even more for any type of
&gt;&gt; scored quantity not only Eq. Dose
&gt;&gt;
&gt;&gt; I also attached my input file if You have time please take a look.
&gt;&gt;
&gt;&gt; Thanks!
&gt;&gt;
&gt;&gt;
&gt;
&gt;

</pre>
                 <br>
                 <br>
               </div>
               <br>
             </div>
           </div>
         </div>
       </blockquote>
       <br>
     </blockquote>
     <br>
   </body>
</html>

--------------020701040801060402030403--
Received on Thu Nov 29 2012 - 10:29:41 CET

This archive was generated by hypermail 2.2.0 : Thu Nov 29 2012 - 10:30:12 CET