search

openmc-dev / openmc

OpenMC Monte Carlo Code
https://docs.openmc.org
Other
765 stars 493 forks source link

Heating tallies don't work with windowed multipole data #1697

Open paulromano opened 4 years ago

paulromano commented 4 years ago

Right now heating tallies rely on tabulated data in the MT=301 or MT=901 cross sections. When multipole is turned on, a check for the temperature index (which is -1 when cross sections were calculated using the WMP data) causes the heating values to be returned as zero. It's not immediately clear what the right solution is for how to calculate heating data when WMP data is in use. In the interim, we should probably give a warning/error if a user tries to tally heating with WMP on.

rockfool commented 4 years ago

I met with similar issue as weell. Take an example for B-10 (n, alpha) reaction rate at BOC of VERA 2E testing case, b10-alpha This will cause significantly increased B-11 production in the depletion. This is because the absorption is considered as gamma in current wmp implementation. No (n, alpha) cross section can be explicitely expressed by poles and residues.

paulromano commented 4 years ago

Thanks for pointing that out @rockfool. I'm going to create a separate issue for that since it is a distinctly different problem from that of the heating tallies.