Incorporating energy into score of scattering reaction rate

User Support
1

Hi everyone,

OpenMC score scattering rate as shown by equation 1 in the following image but i want result obtained by equation 2. I simply want to incorporate energy into scattering reaction rate (or equation 1), please let me know whether openMC has provision of it or not.

Auto Generated Inline Image 1.png
Auto Generated Inline Image 1.png994×518 535 KB

Best regards
Amna Zulfiqar

2

Hi Amna,

The score you’ve written there does not exist in OpenMC. However, it would be fairly straightforward to create a modified version of the code that changes the “scatter” score to tally what you want.

Best,
Paul

Auto Generated Inline Image 1.png
Auto Generated Inline Image 1.png994×518 535 KB

3

Hi Paul,

thank you very much Paul for always replying me on time, please tell me how to score required tally, i am a little bit confused about it.

Regards,
Amna Zulfiqar

4

If you make the following changes in tally.F90 then tallying the ‘scatter’ score should give you the quantity you desire:

diff --git a/src/tally.F90 b/src/tally.F90
index f7d8ab5…9d511ea 100644
— a/src/tally.F90
+++ b/src/tally.F90
@@ -105,6 +105,12 @@ contains
real(8) :: score ! analog tally score
real(8) :: E ! particle energy

  • if (t % estimator == ESTIMATOR_TRACKLENGTH) then
  • E = p % E
  • else
  • E = p % last_E
  • end if