Neutron activation and decay products with time

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

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.

Also, this discussion might also interest you

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

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×1010 will give you the activity in Ci.