Neutron activation and decay products with time

[birsena]

Hi,

I am new. I have a 14 MeV neutrons and would like to activate an element (Al). I would like to have maybe 10 ns shot and after that I would like to tally the decay products with time.

Is there any examples out there I can use or start with?

Thank you for your help

Birsen

[wahidluthfi]

Hi Birsen, welcome back to the community.
Have you tried some examples from the fusion energy neutronic workshop? I think the example on transmutation isotope buildup might fit your cases.

neutronics-workshop/tasks/task_10_activation_transmutation_depletion at main ·...

Also, this discussion might also interest you

Shutdown dose rate calculations

Updated the shut down dose rate example on the fusion-energy neutronics_workshop with help from Paul If anyone is looking for a shut down dose rate simulation feedback on this one is most welcome

[birsena]

Thank you Wahid,

That was helpful.

May I have another question. In the “1-example_transmutation_isotope_buildup.ipynb” the results is obtained for number of atoms. Is there a way to get the result as “Ci” directly instead ?

results = openmc.deplete.ResultsList.from_hdf5(“depletion_results.h5”)

times, number_of_Ag110_atoms = results.get_atoms(my_material, ‘Ag110’)

Thank you

[paulromano]

Yes, you can use the Results.get_activity() method, which by default will give you the activity in Bq. Dividing the activity in Bq by 3.7×10¹⁰ will give you the activity in Ci.