Unit for dose calculation

[yieldz]

Hi all,
I am trying to do some simulation to get the dose rate in Sv/s on a 3d mesh. I tried to use openmc.data.dose_coefficients and openmc.EnergyFunctionFilter, but cannot figure out if I do the normalization correctly.
When searching on the forum I come accross these two posts:
[Tally Units, Material and Cell Meshes, Lethargy]
[Openmc.data.dose_coefficients]

I get confused with the discrepancy in reasonings of these discussions: Is the unit of dose coefficient always [pSv-cm2] or it becomes [pSv-cm3/source_particle] when using openmc.EnergyFunctionFilter? Do I need to perform a flux tally on the same 3d grid as the dose calculation for the normalization?

[cjwiguna]

Hello @yieldz,

Have you read:

dose calculation

I noticed this old post and wanted to add an example as I think dose for neutrons and photons can now be calculated using inbuilt dose coefficients from this code addition by @paulromano https://github.com/openmc-dev/openmc/pull/1558 I think this is nearly there but I appear to be getting strange numbers so I shall check the model a bit more energy_bins_n, dose_coeffs_n = openmc.data.dose_coefficients(particle=“neutron”, geometry=“ISO”) energy_bins_p, dose_coeffs_p = openmc.data.dose_coeffici…

You probably need to find the flux and normalize the flux first to get the correct dose.

[yieldz]

Hi @cjwiguna
Thanks for your reply, this helps.
Though I still have doubt on how the uncertainties are to be propagated, as the scores are probably correlated?