Re: How to simulate URANIUM 235?

<alpine.LFD.2.02.1301101855421.17771_at_pceet073.cern.ch> <50EFB68A.1090701_at_gmail.com> <50EFCCA4.2090506_at_lnl.infn.it>
In-Reply-To: <50EFCCA4.2090506_at_lnl.infn.it>
Content-Type: multipart/mixed;
  boundary="------------030606040504000901040707"
X-Greylist: Sender IP whitelisted, not delayed by milter-greylist-4.2.7 (smtp1.mi.infn.it [192.84.138.69]); Fri, 11 Jan 2013
11:45:28 +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.
--------------030606040504000901040707
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
Content-Transfer-Encoding: 7bit

Thanks it worked! :)

but it doesn't work for RESNUCLEi card I still get errors: "empty data
set" why? :(

I also tried to implement second method ie. to bombard uranium target
with neutrons.
1) Is it possible to get Correct Equivalent dose not only from GAMMA
but other fission fragments (if alpha is not possible then at least
electromagnetic part) of URANIUM235 with the second method?
(Francesco's second note)
2) Generally how one can calculate other fission fragments contribution
of U235 into dose with FLUKA?
3)Can you list please what kind of similar limitations have fluka,
especially with dose calculations? (in order to be careful)

Thanks!

On 11.01.2013 12:26, sarchia wrote:
> Thank you Francesco for your answer, I missed many points!
>
> Dear Vahan, the way I suggested and that you are using in your input
> file only takes into account the gamma emission during uranium decay.
> You see zeros in the USRBIN file and empty RESNUCLEi file because in
> the DCYSCORE card WHAT(1) must be = -1.0 (check the manual).
>
> Lucia
>
>
>
> Il 11/01/2013 07:51, vahan petrosyan ha scritto:
>> Thanks for your Attentions!
>> I tried to implement Lucia's mentioned first method but as usual
>> happens to me I have problems:) I GET ALL ZEROES when plotting USRBIN
>> plotes and EMPTY data set for RESNUCLEI plot.
>>
>> I attached my input file, If you can please take a look.
>>
>> Thanks!
>>
>>
>>
>>
>>
>> On 10.01.2013 23:04, Francesco Cerutti wrote:
>>>
>>> Hi,
>>>
>>> be careful:
>>>
>>> i. for alpha emitters, as already pointed out
>>> (http://www.fluka.org/web_archive/earchive/new-fluka-discuss/4853.html),
>>> the present FLUKA version does NOT generate at all decay alphas, but
>>> lets the isotope decay and transport only electromagnetic particles
>>> emitted along the decay chain.
>>>
>>> ii. no low energy (<20MeV) induced fission is explicitly simulated,
>>> since the multigroup treatment applies. This means that fission
>>> fragments are just separately sampled for residual nuclei scoring
>>> purposes, but they are not at all transported and the interaction of
>>> a given energy group neutron with a given material isotope yields,
>>> in terms of energy deposition, a fixed kerma value averaged over all
>>> the open reaction channels
>>>
>>> Happy New Year
>>>
>>> Francesco
>>>
>>> **************************************************
>>> Francesco Cerutti
>>> CERN-EN/STI
>>> CH-1211 Geneva 23
>>> Switzerland
>>> tel. ++41 22 7678962
>>> fax ++41 22 7668854
>>>
>>> On Thu, 10 Jan 2013, sarchia wrote:
>>>
>>>> Il 10/01/2013 13:47, vahan petrosyan ha scritto:
>>>>> Hi everyone! :)
>>>>>
>>>>> I have may be simple but for me important questions.
>>>>> How one can simulate URANIUM 235 decay?
>>>>>
>>>>> If I have 1gramm Uranium 235 what Dose it will produce. Is it
>>>>> possible
>>>>> to do by fluka?
>>>>>
>>>>> Or I should take fission fragments as primaries and do for each them
>>>>> separate simulation?
>>>>>
>>>>> Is it possible to define two separate beams in fluka in the same
>>>>> input?
>>>>> Thanks!
>>>>>
>>>>>
>>>>
>>>> Hi,
>>>>
>>>> URANIUM 235 is an alpha emitter and FLUKA can simulate its decay using
>>>> the BEAM card with ISOTOPE as SDUM; with HI-PROPE card you specify
>>>> Z and
>>>> A of the element. Then you have to activate the decay mechanism
>>>> through
>>>> the RADDECAY card (semi-analogue mode), and add a DCYSCORE card to
>>>> associate the decay and the scoring. Results will be given normalized
>>>> per one disintegration.
>>>>
>>>> U-235 does not go spontaneous fission, so in case you want to simulate
>>>> the fission process, in my opinion, you should define a neutron beam
>>>> impinging on a Uranium target.
>>>>
>>>> In general, you can use separate beams in one simulation using an
>>>> external source routine edited according to your needs.
>>>> I hope I have answered your questions.
>>>> Cheers,
>>>>
>>>> Lucia
>>>>
>>>>
>>>>
>>
>

--------------030606040504000901040707
Content-Type: application/x-fluka;
  name="uran-modified.inp"
Content-Transfer-Encoding: base64
Content-Disposition: attachment;
  filename="uran-modified.inp"

VElUTEUKCiogU2V0IHRoZSBkZWZhdWx0cyBmb3IgcHJlY2lzaW9uIHNpbXVsYXRpb25zCkRF
RkFVTFRTICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICAgICAgICAgICAgICBQUkVDSVNJTwoqIERlZmluZSB0aGUgYmVhbSBjaGFyYWN0ZXJpc3Rp
Y3MKQkVBTSAgICAgICAgICAgIC0wLjIgICAgIC0wLjAxICAgICAtMC4wMSAgICAgLTAuMDEg
ICAgIC0wLjAxICAgICAgIDEuMElTT1RPUEUKSEktUFJPUEUgICAgICAgICA5Mi4gICAgICAy
MzUuCiogRGVmaW5lIHRoZSBiZWFtIHBvc2l0aW9uCipCRUFNUE9TICAgICAgICAgIDAuMCAg
ICAgICAwLjAgICAgICAgMC4wICAgICAgIDAuMCAgICAgICAwLjAgICAgICAgICAgTkVHQVRJ
VkUKKiBEZWZpbmUgdGhlIGJlYW0gcG9zaXRpb24KKkJFQU1QT1MgICAgICAgICAgMC4xICAg
ICAgICAxLiAgICAgICAtMS4gICAgICAgMC4wICAgICAgIDAuMCAgICAgICAgICBGTE9PRApH
RU9CRUdJTiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICAgICAgICAgICAgICAgQ09NQk5BTUUKICAgIDAgICAgMCAgICAgICAgICAKKiBCbGFjayBi
b2R5ClNQSCBibGtib2R5ICAgIDAuMCAwLjAgMC4wIDEwMDAwMC4wCiogVm9pZCBzcGhlcmUK
U1BIIHZvaWQgICAgICAgMC4wIDAuMCAwLjAgMTAwMDAuMAoqIEN5bGluZHJpY2FsIHRhcmdl
dApSUFAgdGFyZ2V0ICAgICAtMC41IDAuNSAtMC41IDAuNSAtMC4xIDAuMQpFTkQKKiBCbGFj
ayBob2xlCkJMS0JPRFkgICAgICA1ICtibGtib2R5IC12b2lkCiogVm9pZCBhcm91bmQKVk9J
RCAgICAgICAgIDUgK3ZvaWQgLXRhcmdldAoqIFRhcmdldApUQVJHRVQgICAgICAgNSArdGFy
Z2V0CkVORApHRU9FTkQKTUFURVJJQUwgICAgICAgICA5Mi4gICAgMjM1LjA0ICAgICAxOS4w
NSAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIFVSQU5JVU0KTUFURVJJQUwgICAgICAg
ICAxNi4gICAgMzIuMDY2ICAgICAgMi4wNyAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
IFNVTEZVUgpNQVRFUklBTCAgICAgICAgIDE3LiAgIDM1LjQ1MjcgIDAuMDAzMjE0ICAgICAg
ICAgICAgICAgICAgICAgICAgICAgICAgQ0hMT1JJTkUKTUFURVJJQUwgICAgICAgICAxOS4g
ICAzOS4wOTgzICAgICAwLjg2MiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIFBPVEFT
U0lVCiogU2tpbiwgaHVtYW4gKFcmVyB0eXBlIDIpCiogZnJvbTogV29vZHdhcmQgJiBXaGl0
ZSwgQnJpdC4gSi4gUmFkaW9sb2d5LCB2b2wuIDU5LCAxMjE2ICgxOTg2KS4KTUFURVJJQUwg
ICAgICAgICAgICAgICAgICAgICAgICAgICAgMS4wOSAgICAgICAgICAgICAgICAgICAgICAg
ICAgICAgIFNraW5fSHVtCkNPTVBPVU5EICAgICAgIC0xMC4wICBIWURST0dFTiAgICAgLTIw
LjQgICAgQ0FSQk9OICAgICAgLTQuMiAgTklUUk9HRU5Ta2luX0h1bQpDT01QT1VORCAgICAg
ICAtNjQuNSAgICBPWFlHRU4gICAgICAtMC4yICAgIFNVTEZVUiAgICAgIC0wLjMgIENITE9S
SU5FU2tpbl9IdW0KQ09NUE9VTkQgICAgICAgIC0wLjIgICAgU09ESVVNICAgICAgLTAuMSAg
UE9UQVNTSVUgICAgICAgICAgICAgICAgICAgIFNraW5fSHVtCiogLi4rLi4uLjEuLi4uKy4u
Li4yLi4uLisuLi4uMy4uLi4rLi4uLjQuLi4uKy4uLi41Li4uLisuLi4uNi4uLi4rLi4uLjcu
LgpBU1NJR05NQSAgICBCTENLSE9MRSAgIEJMS0JPRFkKQVNTSUdOTUEgICAgICBWQUNVVU0g
ICAgICBWT0lECkFTU0lHTk1BICAgICAgVkFDVVVNICAgIFRBUkdFVApQSE9UT05VQyAgICAg
ICAgICAxLiAgICAgICAgICAgICAgICAgICAgICAgIFZBQ1VVTSAgQExBU1RNQVQKKiBSZXF1
ZXN0cyBzaW11bGF0aW9uIG9mIHJhZGlvYWN0aXZlIGRlY2F5cyBhbmQgc2V0cyB0aGUgY29y
cmVzcG9uZGluZwoqIGJpYXNpbmcgYW5kIHRyYW5zcG9ydCBjb25kaXRpb25zClJBRERFQ0FZ
ICAgICAgICAgIDIuICAgICAgICAgICAgICAgICAwLjAKKiBEZWZpbmUgYW4gaXJyYWRpYXRp
b24gcHJvZmlsZSBmb3IgcmFkaW9hY3RpdmUgZGVjYXkgY2FsY3VsYXRpb25zCipJUlJQUk9G
SSAgICAgICAzNjAwLiAgMTAwMDAwMC4gICAgIDM2MDAuICAgICAgIDAuMCAgICAgMzYwMC4g
IDEwMDAwMDAuCiogRGVmaW5lcyBkZWNheSB0aW1lcyBmb3IgcmFkaW9hY3RpdmUgcHJvZHVj
dCBzY29yaW5nCipEQ1lUSU1FUyAgICAgICAzNjAwLiAgICAgICAwLjAKKiBBc3NvY2lhdGVz
IHNlbGVjdGVkIHNjb3JpbmcgZGV0ZWN0b3JzIHdpdGggdXNlci1kZWZpbmVkIGRlY2F5IHRp
bWVzCkRDWVNDT1JFICAgICAgICAgLTEuICAgICAgICAgICAgICAgICAgICAgICAgIFJkb3Nl
ICAgICAgICAgICAgICAgICAgICBVU1JCSU4KVVNSQklOICAgICAgICAgICAxMC4gICBET1NF
LUVRICAgICAgLTUwLiAgICAgICAgMi4gICAgICAgIDIuICAgICAgICAyLlJkb3NlClVTUkJJ
TiAgICAgICAgICAgLTIuICAgICAgIC0yLiAgICAgICAtMi4gICAgICAxMDAuICAgICAgMTAw
LiAgICAgIDEwMC4mClVTUkJJTiAgICAgICAgICAgMTAuICAgRE9TRS1FUSAgICAgIC01MC4g
ICAgICAgIDIuICAgICAgICAyLiAgICAgICAgMi5Eb3NlClVTUkJJTiAgICAgICAgICAgLTIu
ICAgICAgIC0yLiAgICAgICAtMi4gICAgICAxMDAuICAgICAgMTAwLiAgICAgIDEwMC4mClVT
UkJJTiAgICAgICAgICAgMTAuICAgIFBIT1RPTiAgICAgIC01MC4gICAgICAgIDIuICAgICAg
ICAyLiAgICAgICAgMi5waApVU1JCSU4gICAgICAgICAgIC0yLiAgICAgICAtMi4gICAgICAg
LTIuICAgICAgMTAwLiAgICAgIDEwMC4gICAgICAxMDAuJgoqIEFzc29jaWF0ZXMgc2VsZWN0
ZWQgc2NvcmluZyBkZXRlY3RvcnMgd2l0aCB1c2VyLWRlZmluZWQgZGVjYXkgdGltZXMKRENZ
U0NPUkUgICAgICAgICAtMS4gICAgICAgICAgICAgICAgICAgICAgICAgIHJlc24gICAgICAg
ICAgICAgICAgICAgIFJFU05VQ0xFClJFU05VQ0xFICAgICAgICAgIDMuICAgICAgLTU0LiAg
ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICByZXNuCiogU2V0IHRoZSBy
YW5kb20gbnVtYmVyIHNlZWQKUkFORE9NSVogICAgICAgICAgMS4KKiBTZXQgdGhlIG51bWJl
ciBvZiBwcmltYXJ5IGhpc3RvcmllcyB0byBiZSBzaW11bGF0ZWQgaW4gdGhlIHJ1bgpTVEFS
VCAgICAgICAgIDUwMDAwLgpTVE9QCg==
--------------030606040504000901040707--
Received on Fri Jan 11 2013 - 12:52:17 CET