Skip to content

Conversation

@GuySten
Copy link
Contributor

@GuySten GuySten commented Jan 8, 2026

Description

This PR add an air dose response option to openmc.data.dose_coefficients

This will enable calculation of absorbed dose in air in subsequent PRs.

Checklist

  • I have performed a self-review of my own code
  • I have run clang-format (version 15) on any C++ source files (if applicable)
  • I have followed the style guidelines for Python source files (if applicable)
  • I have made corresponding changes to the documentation (if applicable)
  • I have added tests that prove my fix is effective or that my feature works (if applicable)

@vitmog
Copy link

vitmog commented Jan 8, 2026

@GuySten, as @marco-de-pietri noted in #3700, the mass energy-absorption coefficient $\mu_{en} / \rho$ defines absorbed dose, but not kerma. Also, these coefficients in cm2/g must be multiplied by the particle energy $E$ in MeV with the conversion coefficient 1.602E-10 Gy/(MeV/g) (or 160.2 pGy/(MeV/g)) to obtain absorbed dose in Gy to use them as a response function similarly to the fluence-to-dose conversion factors in pSv-cm2.

I really not sure that this response function deserves to be in the library, because there is not any known to me general problem of calculation of the absorbed dose specifically in air. The used in real world quantities are ICRP effective doses and ambient equivalents H*(10), H*(3) and H*(0.07) (two the last -- rarely), also they could be required equivalent doses ANSI-77 and ANSI-91 for comparisons with legacy measurements in the US and organ equivalent doses from ICRP for medical physics purposes.

EDIT.
The units conversion is implemented right, so my comment above is excessive.

@GuySten
Copy link
Contributor Author

GuySten commented Jan 8, 2026

I think I did convert the units correctly.

@vitmog, after thinking about it. Except maybe for comparison with FISPACT, I cannot think about a case where absorbed dose in air is relevant.
So unless @marco-de-pietri think otherwise I will close this PR.

@GuySten GuySten closed this Jan 8, 2026
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants