Dear Fluka expert, I am currently conducting a simulation on decay heat. The content is to calculate the instantaneous decay heat of liquid lead bismuth target after a period of time when the proton beam stops bombarding. I plan to use the DCYTIMES card to set the cooling time, connect the DCYSCORE card to the USRBIN card, select the energy option in the USRBIN card, and use the RADDECAY card to enable decay. I’m not sure if the combination of these cards can achieve my simulation goals? Please ask all experts to clarify my doubts.
Hi @cuiyaoxuan
Yes, your idea is generally correct. Let me just make sure that I understand what you are trying to simulate. You have your target which is irradiated by a proton beam and then you would like to score the decay heat after a given cooling time, is that correct? If so, on top of the cards that you mentioned, you also have to set up the irradiation profile via IRRPROFI card.
The RADDECAY card should be set to Activation mode. Then it is also important to add two additional PHYSICS cards to enable COALESENCE and EVAPORATION.
I recommend you to look into two lectures from FLUKA courses, one about RP calculations given at the beginner’s course:
and the other about residual radiation exposure from the FLUKA RP topical course:
When you’re scoring ENERGY in a mesh-based USRBIN with an assigned DCYSCORE, keep in mind that the result will be expressed in GeV/cm^3/s (per second!). If you choose to go with a REGION scoring, the results will be given in GeV/s, so to get Watts for example you just need to convert GeV to Joules.
Hope this helps,
Cheers,
Jerzy
Thank you very much for your answer. It has been very helpful in guiding my work