openmc icon indicating copy to clipboard operation
openmc copied to clipboard

Adding material decay gamma dose

Open shimwell opened this issue 2 months ago • 5 comments

Description

This PR adds a simplified gamma dose for a material. This allows users to find the dose at a distance from a gamma emitting material.

In this implementation I sample the gammas but I would be interested knowing if there is a way to do this without sampling the gammas.

The user can provide a distance from the material to get dose in pSv or the distance is assumed to be at the surface of a sphere of the material.

@tokasamwin and @samha-pf I think this feature will make comparing activated materials easier.

Fixes # (issue)

Checklist

  • [x] I have performed a self-review of my own code
  • [x] I have followed the style guidelines for Python source files (if applicable)
  • [x] I have made corresponding changes to the documentation (if applicable)
  • [x] I have added tests that prove my fix is effective or that my feature works (if applicable)

shimwell avatar Nov 20 '25 15:11 shimwell

You can avoid sampling in the special (but very common) case of a discrete gamma distribution.

GuySten avatar Nov 20 '25 17:11 GuySten

re: I would be interested knowing if there is a way to do this without sampling the gammas.

Can't we use Gamma-Ray Dose Constants to do that? (been thinking about it myself, and I think they'd be super useful for rapid evals) https://www.nrc.gov/docs/ML1121/ML11210B521.pdf page 10 http://www.iem-inc.com/information/tools/gamma-ray-dose-constants https://www.osti.gov/servlets/purl/1607209 <--- these guys use ICRP116 and AP https://www.osti.gov/biblio/6246345

This way it depends solely on the nuclide inventories.

yrrepy avatar Nov 21 '25 02:11 yrrepy

re: I would be interested knowing if there is a way to do this without sampling the gammas.

Can't we use Gamma-Ray Dose Constants to do that? (been thinking about it myself, and I think they'd be super useful for rapid evals) https://www.nrc.gov/docs/ML1121/ML11210B521.pdf page 10 http://www.iem-inc.com/information/tools/gamma-ray-dose-constants https://www.osti.gov/servlets/purl/1607209 <--- these guys use ICRP116 and AP https://www.osti.gov/biblio/6246345

This way it depends solely on the nuclide inventories.

This is super interesting, shall we put this in as a separate alternative PR. Could also be nice to compare answers for each option. I do really like your proposed idea, I guess it would be super quick to calculate. It looks like these Gamma-Ray Dose Constants use a distance of 1m.

shimwell avatar Nov 21 '25 08:11 shimwell

You can avoid sampling in the special (but very common) case of a discrete gamma distribution.

Thanks for this suggestion, I have had a go at this.

shimwell avatar Nov 21 '25 10:11 shimwell

I've create an 'issue' feature request for this I think you can change the distance and mass, just need data on air attenuation.

re: I would be interested knowing if there is a way to do this without sampling the gammas. Can't we use Gamma-Ray Dose Constants to do that? (been thinking about it myself, and I think they'd be super useful for rapid evals) https://www.nrc.gov/docs/ML1121/ML11210B521.pdf page 10 http://www.iem-inc.com/information/tools/gamma-ray-dose-constants https://www.osti.gov/servlets/purl/1607209 <--- these guys use ICRP116 and AP https://www.osti.gov/biblio/6246345 This way it depends solely on the nuclide inventories.

This is super interesting, shall we put this in as a separate alternative PR. Could also be nice to compare answers for each option. I do really like your proposed idea, I guess it would be super quick to calculate. It looks like these Gamma-Ray Dose Constants use a distance of 1m.

yrrepy avatar Nov 21 '25 18:11 yrrepy